正压原始方程模式实习报告课件.docVIP

正压原始方程模式实习报告 一、实习目的与要求 1.实习目的：通过正压原始方程模式的实习，加深理解该章所学的基本内容；掌 握当今制作数值天气预报的一般方法和主要步骤； 并在数值计算、 编制程序和上 机操作等方面得到进一步训练。 2.实习要求： 以 1973 年 4 月29 日 08 时(北京时)我国东北、 华北地区 500 百帕等 压面位势高度场及地转风场作为初值， 采用固定的水平侧边界条件， 应用正压原 始方程二次守恒平流格式的模式，制作未来 24 小时有限区域 500 百帕位势高度 场和风场的预报，并写一份实习报告。 二、模式的主要计算框图 三、实习任务 1. 编写两个子程序：五点平滑子程序，地转风初值子程序 2. 子程序编写完成并保证正确的前提下，做 4 个数值试验： ①只做正平滑和做正逆平滑的对比试验； ②地转风子程序中不同差分格式的数值试验； ③是否做边界平滑的数值试验； ④是否做时间平滑的数值试验。 四、实习程序 1、五点平滑的子程序： subroutine ssip(a,w,s,m,n,k,l) 定义五点平滑子程序 dimension a(m,n),w(m,n) 定义数组 if(l==1)then do i=2,m-1 do j=2,n-1 w(i,j)=a(i,j)+s*(a(i-1,j)+a(i+1,j)+a(i,j-1)+a(i,j+1)-4*a(i,j))/4.0 执行正平滑 enddo enddo do i=2,m-1 do j=2,n-1 a(i,j)=w(i,j) enddo Enddo return Else do i=2,m-1 do j=2,n-1 w(i,j)=a(i,j)+s*(a(i-1,j)+a(i+1,j)+a(i,j-1)+a(i,j+1)-4*a(i,j))/4.0 执行正平滑 enddo enddo do i=2,m-1 do j=2,n-1 a(i,j)=w(i,j) enddo enddo do i=2,m-1 do j=2,n-1 w(i,j)=a(i,j)+(-s)*(a(i-1,j)+a(i+1,j)+a(i,j-1)+a(i,j+1)-4*a(i,j))/4.0 执行正逆平滑 enddo enddo do i=2,m-1 do j=2,n-1 a(i,j)=w(i,j) enddo enddo endif return end 结束子程序 2.地转风初值的子程序： subroutine cgw(ua,va,za,rm,f,d,m,n) 定义地转风子程序 implicit none integer m,n,p,q,i,j 定义整型变量 m,n,p,q,i,j real ua(m,n),va(m,n),za(m,n),rm(m,n),f(m,n),d 定义数组 p=m-1 q=n-1 do i=1,m 循环 x 方向上格点数 ua(i,1)=-rm(i,1)*9.8*(za(i,2)-za(i,1))/(f(i,1)*d) 利用公式（4.134）x 方向风速 ua(i,n)=-rm(i,n)*9.8*(za(i,n)-za(i,n-1))/(f(i,n)*d) do j=2,q 循环 y方向上格点数 ua(i,j)=-rm(i,j)*9.8*(za(i,j+1)-za(i,j-1))/(2.0*f(i,j)*d) 利用公式（4.134）x 方向风速 end do end do do j=1,n 循环 y方向上格点数 va(1,j)=rm(1,j)*9.8*(za(2,j)-za(1,j))/(f(1,j)*d) 利用公式（4.134）y 方向风速 va(m,j)=rm(m,j)*9.8*(za(m,j)-za(m-1,j))/(f(m,j)*d) do i=2,p 循环 x方向上格点数 va(i,j)=rm(i,j)*9.8*(za(i+1,j)-za(i-1,j))/(2.0*f(i,j)*d) 利用公式（4.134）y 方向风速 end do end do end subroutine cgw 结束子程序 五、结果（图形）及其分析 1、30 日预报场分析 初始场 预报场 2. 只做正平滑和做正逆平滑的对比试验： 预报场：五点正逆平滑 预报场：五点正平滑 3. 地转风子程序中不同差分格式的数值试验： 初始场：后差地转风子程序 初始场：前差地转风子程序 预报场：前插地转风子程序 预报场：后插地转风子程序 预报场：无时间平滑 预报场：无边界平滑 结果分析： 由前四张图片，我们可以发现， 4 月 29 日 08 时在鄂霍次克海的西岸，有一个大 的低压中心，其中心强度达到 5250，在贝加尔湖附近有一条大的低压槽； 4 月 30 日 08 时，在预报场上，鄂霍次克海西岸的低压向