《数值天气预报实习报告》由会员分享，可在线阅读，更多相关《数值天气预报实习报告（10页珍藏版）》请在金锄头文库上搜索。
1、 课程名称:数值天气预报 姓 名： 班 级: 学 号： 正压原始方程实习报告1.实习目的与规定通过正压原始方程模式的实习，加深理解该章所学的基本内容；掌握当今制作数值天气预报的一般措施和重要环节;并在数值计算、编制程序和上机操作等方面得到进一步训练。以193年4月9日08时（北京时)国内东北、华北地区00hPa等压面位势高度场及地转风场作为初值,采用固定的水平侧边界条件,应用正压原始方程二次守恒平流格式的模式,制作将来24小时有限区域500hPa位势高度场和风场的预报,并写一份实习报告。2.模式的重要计算框图： 基本假定： a大气是均匀不可压流体 b大气处在静力平衡状态 .上边界为自由面 d.
2、无摩擦和绝热过程 .正压大气基本方程： 地转风公式: 技术措施:1.为了减小模式中重力惯性外波的波速,增长差分格式的计算稳定性，把持续方程中水平散度项的系数z改为(z- z0）,并取50pm。2.时间积分方案选择：先采用欧拉后差格式数值积分1小时，然后采用三步法起步的时间中央差格式数值积分11小时。3.边界采用点空间平滑,积分2h后采用5点平滑对个要素场进行内点空间滤波。 五点平滑公式为:4.为了阻尼高频振荡、克制计算解的增长,在积分6小时之后的两个持续时间层应用时间平滑公式对个要素进行时间平滑。5应用二次平流守恒格式的正压原始方程模式。3.实习任务：1）编写两个子程序：五点平滑子程序，地转风
3、初值子程序2)子程序编写完毕并保证对的的前提下,做4个数值实验： 只做正平滑和做正逆平滑的对比实验； 地转风子程序中不同差分格式的数值实验； 与否做边界平滑的数值实验; 与否做时间平滑的数值实验。注意:在原始方程模式中既包具有缓慢移动的大气长波，又有迅速移动的重力惯性波。因此,这种模式不仅可以模拟准地转演变过程,并且还可以模拟地转适应过程。虽然此模式比准地转模式能更好的描述实际大气中的物理过程,但是在这种模式的数值积分中也浮现了某些新的待解决的问题。原始方程模式中包具有快波解，为保证计算稳定,时间步长必须获得很短，这就大大增长了计算工作量。因此,选用既能保证计算稳定、又能节省计算机时的时间积分
4、方案，便成为一种很重要的问题。由于原始方程模式中时间步长必须获得很短,因此非线性计算不稳定的问题就显得尤为突出。因此,构造性能良好的空间差分格式也就成为一种很重要的问题。原始方程模式规定给出多种初始条件。除了要给出初始气压场(位势高度场）外,还必须给出初始风场等资料,如果初始风压场不协调,则在积分过程中就会产生虚假的重力惯性波。这种虚假的重力惯性波会迅速增长，把天气尺度的波动掩盖掉,使预报遭到破坏。因此，资料的初始化也是一种相称重要的问题。原始方程模式对边界条件也尤为敏感，规定给出合适的边界条件也很重要。 4.程序如下）地转风：surouti cw(ua,a,rm，f,m,n)dimena（m
5、,n),v（m,n)，za(m,)，r（m,n)，(m，n） d i=2，m-1do j=2,n-1 u（,j）=(1）*r(i,j）*9(a（i，j+1)-a(i,j-)/2/f(i,j） a(i,)m(i,j)*8*(z(i+1,j)-(i1,j)2/df（，j) eddo enddrtun2）5点平滑:subroutn sp（a，s，,n,k，）dimenin a（m，n),(m,n）(l=1)hend i=2,1do j=2,n-w（i，j)=a(i，j）*(a（i-1，j)+a（+1,j)+a(i，j1)+a(i,j+1）-4*a（,j)/.0enddonodoi=2,m-1 d 2
6、,-1 a(i，j）=w（i，j）endoendoretrnelse doi=2,m- d j=2,n- (,j）a(,j)+*(a(-1，j)(i,j)a(i,-1）a(i,j+1）-4a（i,j)）/40enddonddo 2,-1 do j,- (,j)=w(,j) endenddod i=2,m-1 do j=2,n1 w（i，)=a(i，j)+(-）(a(i-1,j）+a（i+,j)+a(i,-1)+(i,j+1）4*a(，j)/. edoenddodo i=,m- doj=2,-1 a（i，j)=w(i,j) edenddoendifrtu nd五、实验成果(图形）及其分析1、30
7、日预报场分析 图 初始场 图 预报场 2、只做正平滑和做正逆平滑的对比实验： 图 预报场:五点正逆平滑 图 4预报场：五点正平滑、地转风子程序中不同差分格式的数值实验: 图初始场：后差地转风子程序 图 初始场：前插地转风子程序 图预报场:前插地转风子程序 图8预报场：后插地转风子程序 图 10预报场： 无时间平滑 图 9预报场: 无边界平滑成果分析: 由图一和图二,图三和图四,我们可以清晰的看到：4月29日08时在鄂霍次克海的西岸，有一种大的低压中心,其中心强度达到5250，在贝加尔湖附近有一条大的低压槽;月30日08时,在预报场上,鄂霍次克海西岸的低压向东移动,达到鄂霍次克海的上空,中心强度
8、达到200，比2日加强,贝加尔湖的低压槽向南,向东移动。也许导致一次冷空气的南下过程。 我们拿30日的预报场与30日的实况场(课本上)做一下对比，可以看到预报场的低压中心偏北大概偏个纬距左右，有较大的偏差，模拟效果很不好；但是，贝加尔湖附近的低压槽预报场与实况场非常的近似，故而模拟效果较好。而我们需要做对中国的预报，影响系统重要是贝加尔湖地区的低压槽影响，因此对整个预报来说，还是有较好的参照价值的。5.小结 通过本次实习加强了对Fotran语言编程和GrD绘图的熟悉以及使用，是后来的工作生活中可以更加纯熟的使用它们,同步也加强了对正压原是方程的理解,以及五点平滑,与地转风的计算,加深了对数值天气预报这门课的理解与掌握。本次接触到的正压原始方程模式,算是个比较简朴的模式，由于在00hPa高度上,大气是正压的，准地转,准静力平衡的,因此大气环流的形势比较简朴、稳定,所得成果基本与现实差距不大。但是大气并没有假设的那样简朴,大气是斜压的，复杂的系统。在如今,模式正在发展中,仍然无法完全描述目前的大气变化趋势,大气内部仍然有诸多机制对我们而言是未知的。这需要我们对大气继续摸索，并发展模式。
