数值天气预报实习 题 目：正压原始方程模式实习报告 院 系： xxxx 学院 . 专业年级： xxxx 专业 . 班 级： xxxx 级 x 班 . 姓 名： 步留情 . 一、实习目的与要求 通过正压原始方程模式的实习，加深理解该章所学的基本内容；掌握当今制作数值天气预报的一般方法和主要步骤；并在数值计算、编制程序和上机操作等方面得到进一步训练。 以 1973 年 4 月 29 日 08 时(北京时)我国东北、华北地区 500 百帕等压面位势高度场及地转风场作为初值，采用固定的水平侧边界条件，应用正压原始方程二次守恒平流格式的模式，制作未来 24 小时有限区域 500 百帕位势高度场和风场的预报，并写一份实习报告。 二、模式的主要计算框图 输入初始位势高度场ᵄᵅ,ᵅ0 计算初始风场ui,j和vi,j 计算预报要素的变化倾向Ei,j、Gi,j和Hi,j 用欧拉后差积分格式 1h 用中央差积分格式 1h 是否算完 6h 是否算完 12h 计算mi,j和ᵅᵅ,ᵅ 时间平滑 停机 空间平滑 是否算完 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 jijijijijijiyxjiFFFFFSFF,1,,11,,1,,,44~ 区域内5 点平滑(正逆平滑)的子程序，l=1 为只执行正平滑，l=2 为执行正逆平滑. subroutine ssip(a,w,s,m,n,l) dimension a(m,n),w(m,n)...
