正压原始方程模式实习报告.docVIP

正压原始方程模式实习报告

一、实习目得与要求

1、实习目得:通过正压原始方程模式得实习,加深理解该章所学得基本内容；掌握当今制作数值天气预报得一般方法与主要步骤；并在数值计算、编制程序与上机操作等方面得到进一步训练.

2、实习要求:以１９73年４月29日08时（北京时）我国东北、华北地区5０0百帕等压面位势高度场及地转风场作为初值，采用固定得水平侧边界条件，应用正压原始方程二次守恒平流格式得模式，制作未来24小时有限区域5０0百帕位势高度场与风场得预报，并写一份实习报告.

二、模式得主要计算框图

三、实习任务

1、编写两个子程序：五点平滑子程序,地转风初值子程序

2、子程序编写完成并保证正确得前提下,做4个数值试验：

①只做正平滑与做正逆平滑得对比试验;

②地转风子程序中不同差分格式得数值试验;

③就是否做边界平滑得数值试验；

④就是否做时间平滑得数值试验。

四、实习程序

?1、五点平滑得子程序：?

ｓｕbroutｉnｅssｉｐ(ａ,w，s，m,n,ｋ,l)定义五点平滑子程序

dimｅnsiona(m，n）,ｗ(m，n）定义数组

ｉf(l==1）thｅn

doi=2，m—1

doj=2，n-１

ｗ(ｉ,j)=a（ｉ，ｊ)+s*(a(i—１,j)+ａ(ｉ+１,j)+a(i,ｊ-1)+ａ(ｉ，ｊ＋1)—４*a（ｉ，j）)／4、0执行正平滑

eｎddo

eｎｄdo

doi=2,m－1

ｄoj＝2,ｎ-1

ａ(i,j)=w（i，j）

ｅnddo

Ｅnddo

rｅturn

Eｌse

doｉ=2,m—1

doj=2，n—1

w(i，j)=a(i,ｊ）+ｓ*（a(i-１，j)+ａ（ｉ+1,j)+a(i,ｊ—１)+ａ(i,j+1）-4*ａ(i,ｊ）)/4、０执行正平滑

eｎddo

enｄdo

doi=2,ｍ—1

dｏj=2，n—１

a（ｉ,j）＝w(i,j）

eｎddｏ

enｄdｏ

ｄoｉ=2，m—１

dｏｊ=2,n—1??

w（ｉ,j）=a（i,ｊ)+（-s)*(a（i－1，j）＋a(i+１,j)+a（ｉ，j-1)+a(i，j+1)-４＊a（i,j))/4、０执行正逆平滑

ｅｎddo

enｄdo

doi=2，m－1

ｄoj=2,n-１

a(i,j)=w(i,j)

ｅnｄdo

enｄdo

ｅｎdif

reｔurn

ｅｎd结束子程序

2、地转风初值得子程序：?

suｂroutｉne?cgw(ua,vａ，zａ，ｒm,f,d，m，n）定义地转风子程序

implicit?ｎｏnｅ?

inteｇｅr?ｍ,n，p，q，i,j?定义整型变量m，n，ｐ，q,i,ｊ

ｒeａｌ?ｕa(m，n)，va（m，ｎ)，ｚa(m,n),ｒm（ｍ，n),ｆ(m,n）,d?定义数组

p＝ｍ-1?

q=n—1?

do??i＝1，m?循环ｘ方向上格点数

ｕa（i,1）＝—rm(i,1)＊9、8*(ｚa(i,2)—zａ(ｉ,1))/(f(i,1)*d）????利用公式（4、1３4)x方向风速

ua(i,n）=—rm（i,n)＊9、8*（za(i,n）-zａ(i，n—１））/（f（i,n)*d)?

do??ｊ=2，q?循环y方向上格点数

ｕa(i,j)＝－rm(i,ｊ)*9、８＊（ｚa(i,j＋１）—ｚａ(i，ｊ－１)）／（２、0*f(i，j）＊ｄ)?利用公式(４、1３4)x方向风速

enｄ?dｏ

end?do?

dｏ??j＝1,n循环y方向上格点数

ｖa（1,j)=rm（１，j)*9、8*(zａ（2,ｊ）-za（１，j）)/（f(１，j)*ｄ）?利用公式(4、134)ｙ方向风速

ｖa(ｍ，j)=rm（m,j）*９、８＊(ｚa(m,j)—za(ｍ-１，j))/(f（m,ｊ）*d)?

do??i=２,ｐ?循环x方向上格点数

vａ(i，j)=rｍ(i,ｊ）＊9、8＊(ｚa(i+１,j）-zａ（i-1，j）)/（2、0＊f（ｉ，j)＊ｄ)??利用公式（4、13４）ｙ方向风速

end?do?

eｎd?do?

end?suｂroutine?cgw结束子程序

五、结果(图形)及其分析

１、30日预报场分析

初始场

预报场

2、只做正平滑与做正逆平滑得对比试验:

预报场：五点正逆平滑