数值天气预报实习.docVIP

标准文档 实用文案 数值天气预报实习 题 目：正压原始方程模式实习报告 院 系： xxxx 学院 . 专业年级： xxxx 专业 . 班 级： xxxx级x班 . 姓 名： 步留情 . 一、实习目的与要求 通过正压原始方程模式的实习，加深理解该章所学的基本内容；掌握当今制作数值天气预报的一般方法和主要步骤；并在数值计算、编制程序和上机操作等方面得到进一步训练。 以1973年4月29日08时(北京时)我国东北、华北地区500百帕等压面位势高度场及地转风场作为初值，采用固定的水平侧边界条件，应用正压原始方程二次守恒平流格式的模式，制作未来24小时有限区域500百帕位势高度场和风场的预报，并写一份实习报告。 二、模式的主要计算框图 输入初始位势高度场 输入初始位势高度场Z 计算初始风场u 计算预报要素的变化倾向E 用欧拉后差积分格式1h 用中央差积分格式1h 是否算完6h 是否算完12h 计算m 时间平滑 停机 空间平滑 是否算完24h 打印 否 否 是 是 否 是 三、实习任务 编写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) 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=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 return end 区域内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,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 end ctl文件： DSET H:\shuzhishixi\duizhao\h.grd TITLE 500hpa height UNDEF -9990.0 XDEF 20 LINEAR 85 3.5 YDEF 16 LINEAR 32.5 2.5 ZDEF 1 linear 1 1 TDEF 2 LINEAR 00Z29apr1973 24hr VARS 1 h 0 99 500 ENDVARS 绘图gs文件 reinit open h:\shuzhishixi\duizhao\HUV.CTL open h:\shuzhishixi\shiyan1\HUV.CTL open h:\shuzhishixi\shiyan2\hou\HUV.CTL open h:\shuzhishixi\