数值天气预报实习

1、数值天气预报实习题 目：正压原始方程模式实习报告院 系： xxxx 学院 .专业年级： xxxx 专业 .班 级： xxxx级x班 .姓 名： 步留情 .一、实习目的与要求 通过正压原始方程模式的实习，加深理解该章所学的基本内容；掌握当今制作数值天气预报的一般方法和主要步骤；并在数值计算、编制程序和上机操作等方面得到进一步训练。 以1973年4月29日08时(北京时)我国东北、华北地区500百帕等压面位势高度场及地转风场作为初值，采用固定的水平侧边界条件，应用正压原始方程二次守恒平流格式的模式，制作未来24小时有限区域500百帕位势高度场和风场的预报，并写一份实习报告。二、模式的主要计算框图输

2、入初始位势高度场Zi,j0计算初始风场ui,j和vi,j计算预报要素的变化倾向Ei,j、Gi,j和Hi,j用欧拉后差积分格式1h用中央差积分格式1h是否算完6h是否算完12h计算mi,j和fi,j时间平滑停机空间平滑是否算完24h打印否否是是否是三、实习任务编写2个子程序：五点平滑子程序、地转风初值子程序先编写五点平滑的子程序，在保证能够得到正确结果的前提下，再去编写地转风初值的子程序。所用公式：两个子程序编写完成并保证正确的前提下，还需要做4个数值试验：1）只做正平滑和做正逆平滑的对比试验；2）地转风子程序中不同差分格式的数值试验；3）是否做边界平滑和内点平滑的数值试验；4）是否做时间平滑的

3、数值试验。四、实习程序（自己编制的子程序）地转风初值的子程序subroutine cgw(ua,va,za,rm,f,d,m,n)dimension ua(m,n),va(m,n),za(m,n),rm(m,n),f(m,n)do i=1,m ua(i,1)=-rm(i,1)/f(i,1)*9.8/d*(za(i,2)-za(i,1) ua(i,n)=-rm(i,n)/f(i,n)*9.8/d*(za(i,n)-za(i,n-1) do j=2,n-1 ua(i,j)=-rm(i,j)/f(i,j)*9.8/2/d*(za(i,j+1)-za(i,j-1) end do end do do j

4、=1,n va(1,j)=rm(1,j)/f(1,j)*9.8/d*(za(2,j)-za(1,j) va(m,j)=rm(m,j)/f(m,j)*9.8/d*(za(m,j)-za(m-1,j) do i=2,m-1 va(i,j)=rm(i,j)/f(i,j)*9.8/2/d*(za(i+1,j)-za(i-1,j) end do end do returnend区域内5点平滑(正逆平滑)的子程序，l=1为只执行正平滑，l=2为执行正逆平滑.subroutine ssip(a,w,s,m,n,l)dimension a(m,n),w(m,n) do 10 i=2,m-1 do 10 j=2

5、,n-1 10 w(i,j)=a(i,j)+s*0.25*(a(i-1,j)+a(i+1,j)+a(i,j-1)+a(i,j+1)-4*a(i,j) a(2:m-1,2:n-1)=w(2:m-1,2:n-1) if(l=2) then do 20 i=2,m-1 do 20 j=2,n-1 20 w(i,j)=a(i,j)-s*0.25*(a(i-1,j)+a(i+1,j)+a(i,j-1)+a(i,j+1)-4*a(i,j) a(2:m-1,2:n-1)=w(2:m-1,2:n-1) end if return endctl文件：DSET H:shuzhishixiduizhaoh.grdT

6、ITLE 500hpa heightUNDEF -9990.0XDEF 20 LINEAR 85 3.5YDEF 16 LINEAR 32.5 2.5ZDEF 1 linear 1 1TDEF 2 LINEAR 00Z29apr1973 24hrVARS 1h 0 99 500ENDVARS绘图gs文件'reinit''open h:shuzhishixiduizhaoHUV.CTL''open h:shuzhishixishiyan1HUV.CTL''open h:shuzhishixishiyan2houHUV.CTL'

7、9;open h:shuzhishixishiyan2qianHUV.CTL''open h:shuzhishixishiyan2wuHUV.CTL''open h:shuzhishixishiyan3bsHUV.CTL''open h:shuzhishixishiyan3nbnsHUV.CTL''open h:shuzhishixishiyan3nsHUV.CTL''open h:shuzhishixishiyan4HUV.CTL''set ccolor 1''set cstyle

8、 1''set cthick 6''d h.1(t=1)''printim H:shuzhishixiduizhao1.png white''c''set ccolor 1''set cstyle 1''set cthick 6''d h.1(t=2)''printim H:shuzhishixiduizhao2.png white''c''set ccolor 1''set cstyle 1'&

9、#39;set cthick 6''d h.1(t=1)''set ccolor 1''set cstyle 6''set cthick 5''d h.1(t=2)''printim H:shuzhishixiduizhao1-2.png white''c'i=2while(i<10)'set gxout contour''set ccolor 1''set cstyle 1''set cthick 6'&

10、#39;d h.'i'(t=2)''printim H:shuzhishixi'i'.png white''set ccolor 1''set cstyle 6''set cthick 5''d h.1(t=2)''printim H:shuzhishixi'i'-duizhao.png white''c''set gxout shaded''d abs(h.'i'(t=2)-h.1(t=

11、2)''cbar 1.0 0''set gxout contour''set ccolor 1''set cstyle 1''set cthick 6''d h.'i'(t=2)''printim H:shuzhishixi'i'wucha-duizhao.png white''c'i=i+1endwhile'set fwrite h:shuzhishixiwucha.grd'i=2while(i<10)&

12、#39;set gxout fwrite''d aave(abs(h.'i'(t=2)-h.1(t=2),lon=85,lon=155,lat=32.5,lat=72.5)'i=i+1endwhile'disable fwrite''reinit'读取fwrite的f90程序program wuchagrdreal wucha(8)open(1,file='H:shuzhishixiwucha.grd',form='binary') read(1)wuchaclose(1)print*,wu

13、chaend五、结果（图形）及其分析 初始场图预报结果图与初始场之间的对比图（实线初始场，点划线为预报场）从图上可以看出，系统整体的在东移。29日鄂霍次克海西安有一个低压中心，强度在5250左右，贝加尔湖以北有低槽，贝加尔湖位于槽前。30日预报场，东边的低压中心向东移动大约5度，到达鄂霍次克海上空，并且中心强度加强，下降到5200左右，贝加尔湖槽向南、向东移动，贝加尔湖位于槽后。贝湖槽的移动可能会引起冷空气南下。实验图与预报场之间差值绝对值的分布图（阴影填色部分为差值大小（实验场减对照场），实线为实验预报图，长虚线为对照图）1） 只做正平滑和做正逆平滑的对比试验；只做正平滑由1）图可知，整体偏

14、差不大，在槽脊处的差值相对较大。最大差值绝对值达到50位势米。特点一：在气旋性曲率处相对正逆平滑的结果而言要大，反气旋性曲率处相对正逆平滑的结果而言要小。这说明正逆平滑对长波振幅的衰减比只做一次正平滑得对长波整幅的衰减要小。特点二：而且在图中可以看出，波长较长的波的振幅差别不是很大，波长在2倍格距（5个经度）到8倍格距（20个经度）之间的波的振幅差别就比较大了。上述两个特点符合响应函数的图像特征2）地转风子程序中不同差分格式的数值试验；1、内点全部用后差格式计算 2、内点全部用前差格式计算 3、次外层点用中央差格式计算，内层点用四阶精度的五点差分格式计算。由2）1、2图可以看出由后差格式与前差

15、格式得出的初始风场作出的预报场与对照场之间的差距还是比较明显的，后差在-40到+100之间，前差在-100到+40之间，而五点差分格式得出的初始风场作出的预报场与对照场之间的差距几乎为零。可以看出后差与前差的结果中，由于精度不高，截断误差较大，激发了许多快波，在边界处尤其明显。值得一提的是，后差、前差分别与对照场的差值，在空间上的分布相当一致，在数值上，在对应格点的绝对值大小几乎相同，符号相反。由四阶精度的五点差分方案作出的初值场得出的实验结果与二阶精度的中央差格式的得出的对照场结果可以说是完全一样，说明在风场初始化方面，二阶精度的差分格式已经足够了。3）是否做边界平滑和内点平滑的数值试验；1

16、、只做边界平滑，不做内点平滑2、不做内点平滑也不做边界平滑3、只做内点平滑，不做边界平滑从3）1、2、3的图可以看出，不做内点平滑的影响相对较小，不做边界平滑则影响非常大。内点平滑是为了滤除短波波动，抑制计算解增长，从3）1、的图可以看出，不做内点平滑，对整体环流形势影响较小，但差值阴影图中会有布满全图的相互之间间隔5五个经纬度的小块正负值区域，且正值区与负值区相间分布，相邻的正负值区域数值的绝对值很接近。越靠近途中央，正负值区数值的绝对值大小越大。从3）2、3可以看出，仅仅不做边界平滑，由固定边界条件引起的边界点与内点之间的要素梯度对结果的影响非常大，激发的快波振幅较大，数量较多，差值分布从

17、-120到+80都有。不做边界平滑也不做内点平滑得出的结果与仅仅不做边界平滑得出的结果相似，在差值大小向上，边界、内点平滑都不做的结果更大一些，从-150到+120.4）是否做时间平滑的数值试验。不做时间平滑时间平滑的目的是阻尼高频振荡，抑制计算解增长，从4）图可以看出，不做时间平滑，对整体形势影响较小，差值分布从-9到+15。各个实验的 差值/差值绝对值 以面积为权重的平均结果（单位：位势米）1）只做正平滑和做正逆平滑的对比试验； 0.7436748/5.1711942）地转风子程序中不同差分格式的数值试验； 1、内点全部用后插格式计算 5.635431/17.05854 2、内点全部用前插

18、格式计算 -5.526467/17.96309 3、次外层点用中央差格式计算，内层点用四阶精度的五点差分格式计算。 0.0000000E+00/0.0000000E+003）是否做边界平滑和内点平滑的数值试验； 1、只做边界平滑，不做内点平滑 -0.061916102/5.058856 2、不做内点平滑也不做边界平滑 6.993883/34.08777 3、只做内点平滑，不做边界平滑 5.259184/27.857114）是否做时间平滑的数值试验。 -1.323557/2.804939六、小结与讨论 （按整体差值大小从小到大描述）五点差分方案的结果几乎没有差距；不做时间平滑的影响在短时间内（因为积分时间跨度不是很长，只能说在短时间内影响很小）很小；只做正平滑，对长坡振幅有一定的衰减，影响较小；不做内点平滑，整体影响较小，其中差值以面积为权重的平均结果只有0.061916102，在数值上说明了