资源描述：
《马嘉理的数值预报》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。
1、南京信息工程大学正压原始方程模式实验（实习）报告实验（实习）名称正压原始方程模式实验（实习）日期2010.6.1得分指导教师邵海燕系大气科学专业大气科学年级08班次2姓名马嘉理学号20081301057一、实验目的利用五点正逆平滑共式和地转风初始值的公式来计算初始场和预报值并且利用所学的数值预报中有关的相关知识利用FORTRAN编写出计算地转风初值的程序和计算五点正逆平滑的程序并且利用GRADS来画图二、实习要求mbpm1.for为提供学生上机实习用的正压原始方程模式，程序中删去了计算地转风的子程序，也删去了五点平滑（正逆平滑）的子程序。同学们需要先编写五点平滑的子程序，在保证能够得到正确结果
2、的前提下（是否正确得作图给老师看），再去编写地转风初值的子程序。在未编写地转风初值的子程序前，地转风初值可由数据文件ua.dat和va.dat提供（见程序说明）。编写完地转风初值子程序后，需将我提供的数据文件ua.dat和va.dat删去，此时的地转风初值由自己的程序提供。huv.ctl和huv.gs分别为用于Grads做图的数据描述文件和批处理文件。三模式的基本假定，基本方程，技术方法，计算流程等的介绍（1）基本方程地图投影坐标下的正压原始方程模式组空间差分格式应用二次平流守恒格式——有限差分近似方程其中，，这就是正压原始方程模式的预报方程组理论分析和预报实践表明，由于观测的风场与高度场的不
3、平衡，，以及风场，高度场与模式之间的不协调，直接用观测的风场和高度场作为原始方程模式的初始值容易产生高频振荡，使数值积分变为不稳定。因此在应用原始方程模式作数值预报之前，必须对资料加以处理，即所谓的资料初始化，动力初始化和变分初始化。2）基本假设正压原始方程模式是在一定的基本假设下建立起来的、能够比较准确地描述大气运动状态的方程组，假定的条件有：1、大气为正压大气；2、大气为均匀不可压缩的流体；3、大气满足静力平衡；4、大气上界为自由面；5、不考虑摩擦和非绝热加热作用。3）技术方法1、模式的初始条件：理论分析和预报实践表明，由于观测的风场与高度场之间的不平衡，以及风场、高度场与模式之间的不协调
4、，直接用观测的风场和高度场作为原始方程模式的初始值容易产生高频振荡，使数值积分变得不稳定。为了解决这个问题就必须对资料进行初始化，所谓初始化就是指用一些已知的风压场平衡关系，或运用运动方程等求得诊断方程来处理初值，使风场同气压场平衡或近似平衡的方法。资料初始化的方法主要有：静力初始化、动力初始化和变分初始化。此处我们采用的是静力初始化，初始条件采用地转风初值2、边界条件：模式中采用的是固定边界条件，其中β表示预报区域的水平侧边界。3、时间积分方案（1）先采用欧拉—后差格式；（2）随后，采用三步法起步的时间中央差格式。另外，在数值积分的过程中，为了阻尼高频振荡，抑制计算解的增长，可以穿插进行时间
5、平滑，其对应公式为：式中S为时间平滑系数；同时，可以穿插进行空间平滑，来滤除短波扰动，抑制非线性计算不稳定。四计算流程输入初始位势高度场计算mi,j和fi,j计算初始风场u0i,j和v0i,j计算预报要素的变化倾向Ei,jGi,j和Hi,j用欧拉后差格式积分1小时用中央差格式积分11小时是否算完12小时是否算完6小时空间平滑是否算完24小时时间平滑打印五、实习程序运行出图：边界五点平滑的程序cspacesmoothingforinternalpoints区域内5点平滑(正逆平滑)c请同学编写区域内5点平滑(正逆平滑)的子程序！！！应用书中（4.126）式c注：此程序必须设计成开关形式，保证既可
6、选做正逆平滑，又可选做正平滑l=1为只c执行正平滑，l=2为执行正逆平滑.subroutinessip(a,w,s,m,n,l)dimensiona(m,n),w(m,n)doi=2,m-1doj=2,n-1w(i,j)=a(i,j)+s/4*(a(i+1,j)+a(i,j+1)+a(i-1,j)+a(i,j-1)-4*a(i,j))enddoenddodoi=2,m-1doj=2,n-1a(i,j)=w(i,j)enddoenddoif(l==2)thendoi=2,m-1doj=2,n-1w(i,j)=a(i,j)-s/4*(a(i+1,j)+a(i,j+1)+a(i-1,j)+a(i,j
7、-1)-4*a(i,j))enddoenddodoi=2,m-1doj=2,n-1a(i,j)=w(i,j)enddoenddoendifreturnend六、作图结果原始场：预报场：2、在已编写地转风初值的子程序后，删去老师的数据文件ua.dat和va.dat后.对原程序进行修改：首先，将原程序ua.dat和va.dat删去。然后，写入地砖风初始值编程程序。subroutinecgw(ua,va
