数值天气预报实习.docx

数值天气预报实习题目：正压原始方程模式实习报告院系：xxxx学院 .专业年级：xxxx专业 .班级：xxxx级x班 .姓名：步留情 .一、实习目的与要求通过正压原始方程模式的实习，加深理解该章所学的基本内容；掌握当今制作数值天气预报的一般方法和主要步骤；并在数值计算、编制程序和上机操作等方面得到进一步训练。以1973年4月29日08时(北京时)我国东北、华北地区500百帕等压面位势高度场及地转风场作为初值，采用固定的水平侧边界条件，应用正压原始方程二次守恒平流格式的模式，制作未来24小时有限区域500百帕位势高度场和风场的预报，并写一份实习报告。二、模式的主要计算框图三、实习任务编写2个子程序：五点平滑子程序、地转风初值子程序先编写五点平滑的子程序，在保证能够得到正确结果的前提下，再去编写地转风初值的子程序。所用公式：两个子程序编写完成并保证正确的前提下，还需要做4个数值试验：1）只做正平滑和做正逆平滑的对比试验；2）地转风子程序中不同差分格式的数值试验；3）是否做边界平滑和内点平滑的数值试验；4）是否做时间平滑的数值试验。四、实习程序（自己编制的子程序）地转风初值的子程序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)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*(za(i,n)-za(i,n-1))do j=2,n-1ua(i,j)=-rm(i,j)/f(i,j)*9.8/2/d*(za(i,j+1)-za(i,j-1))end do end dodo j=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))end do end doreturnend区域内5点平滑(正逆平滑)的子程序，l=1为只执行正平滑，l=2为执行正逆平滑.subroutinessip(a,w,s,m,n,l)dimension a(m,n),w(m,n)do 10 i=2,m-1do 10 j=2,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) thendo 20 i=2,m-1do 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 ifreturnendctl文件：DSET H:\shuzhishixi\duizhao\h.grdTITLE 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文件reinitopen h:\shuzhishixi\duizhao\HUV.CTLopen h:\shuzhishixi\shiyan1\HUV.CTLopen h:\shuzhishixi\shiyan2\hou\HUV.CTLopen h:\shuzhishixi\shiyan2\qian\HUV.CTLopen h:\shuzhishixi\shiyan2\wu\HUV.CTLopen h:\shuzhishixi\shiyan3\bs\HUV.CTLopen h:\shuzhishixi\shiyan3\nbns\HUV.CTLopen h:\shuzhishixi\shiyan3\ns\HUV.CTLopen h:\shuzhishixi\shiyan4\HUV.CTLsetccolor 1setcstyle 1setcthick 6d h.1(t=1)printim H:\shuzhishixi\duizhao1.png whitecsetccolor 1setcstyle 1setct