资源描述：
《数值天气预报实习》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
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）只做正平
3、滑和做正逆平滑的对比试验；2）地转风子程序中不同差分格式的数值试验；3）是否做边界平滑和内点平滑的数值试验；4）是否做时间平滑的数值试验。四、实习程序（自己编制的子程序）地转风初值的子程序subroutinecgw(ua,va,za,rm,f,d,m,n)dimensionua(m,n),va(m,n),za(m,n),rm(m,n),f(m,n)doi=1,mua(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*(
4、za(i,n)-za(i,n-1))doj=2,n-1ua(i,j)=-rm(i,j)/f(i,j)*9.8/2/d*(za(i,j+1)-za(i,j-1))enddoenddodoj=1,nva(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))doi=2,m-1va(i,j)=rm(i,j)/f(i,j)*9.8/2/d*(za(i+1,j)-za(i-1,j))enddoend
5、doreturnend区域内5点平滑(正逆平滑)的子程序，l=1为只执行正平滑，l=2为执行正逆平滑.subroutinessip(a,w,s,m,n,l)dimensiona(m,n),w(m,n)do10i=2,m-1do10j=2,n-110w(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)thendo20i=2,m-1do20j=2,n-120w(i,
6、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)endifreturnendctl文件：DSETH:shuzhishixiduizhaoh.grdTITLE500hpaheightUNDEF-9990.0XDEF20LINEAR853.5YDEF16LINEAR32.52.5ZDEF1linear11TDEF2LINEAR00Z29apr197324hrVARS1h09950
7、0ENDVARS绘图gs文件'reinit''openh:shuzhishixiduizhaoHUV.CTL''openh:shuzhishixishiyan1HUV.CTL''openh:shuzhishixishiyan2houHUV.CTL''openh:shuzhishixishiyan2qianHUV.CTL''openh:shuzhishixishiyan2wuHUV.CTL''openh:shuzhishixishiyan3bsHUV.CTL''openh:
8、shuzhishixishiyan3bnsHUV.CTL''openh:shuzhishixishiyan3sHUV.CTL''openh:shuzhishixishiyan4HUV.CTL''setccolor1''setcstyle1''setcthick6''dh.1(t=1)''printimH:shuzhishixid
