用欧拉后差格式积分1小时 时间平滑是否算完 时 6小是否算完12小时 用中央差格式积分11小时 空间平滑 是否算完24小时
五、实习程序
运行出图：边界五点平滑的程序
c space smoothing for internal points 区域内5点平滑(正逆平滑) c 请同学编写区域内5点平滑(正逆平滑)的子程序！！！应用书中（4.126）式
打印 c 注：此程序必须设计成开关形式，保证既可选做正逆平滑，又可选做正平滑 l=1为只c 执行正平滑，l=2为执行正逆平滑.
subroutine ssip(a,w,s,m,n, l)
if(l==2) then do i=2,m-1
do j=2,n-1 enddo
w(i,j)=a(i,j)-s/4*(a(i+1,j)+a(i,j+1)+a(i-1,j)+a(i,j-1)-4*a(i,j))
enddo do i=2,m-1
do j=2,n-1
dimension a(m,n),w(m,n) do i=2,m-1
do j=2,n-1 enddo
w(i,j)=a(i,j)+s/4*(a(i+1,j)+a(i,j+1)+a(i-1,j)+a(i,j-1)-4*a(i,j))
enddo do i=2,m-1
do j=2,n-1
a(i,j)=w(i,j)
enddo enddo
a(i,j)=w(i,j)
enddo enddo endif return end
六、作图结果
原始场：
预报场：
2、在已编写地转风初值的子程序后，删去老师的数据文件ua.dat和va.dat后. 对原程序进行修改：
首先，将原程序ua.dat和va.dat删去。 然后，写入地砖风初始值编程程序。
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=2,m-1 do j=2,n-1
ua(i,j)=0-rm(i,j)*9.8/f(i,j)*((za(i,j+1)-za(i,j-1))/2/d) va(i,j)=rm(i,j)*9.8/f(i,j)*((za(i+1,j)-za(i-1,j))/2/d) enddo enddo
ua(m,n)=0-rm(m,n)*9.8/f(m,n)*((za(m,n)-za(m,n-1))/d) va(m,n)=rm(m,n)*9.8/f(m,n)*((za(m,n)-za(m-1,n))/d) ua(1,n)=0-rm(1,n)*9.8/f(1,n)*((za(1,n)-za(1,n-1))/d) va(1,n)=rm(1,n)*9.8/f(1,n)*((za(2,n)-za(1,n))/d) ua(m,1)=0-rm(m,1)*9.8/f(m,1)*((za(m,2)-za(m,1))/d) va(m,1)=rm(m,1)*9.8/f(m,1)*((za(m,1)-za(m-1,1))/d) ua(1,1)=0-rm(1,1)*9.8/f(1,1)*((za(1,2)-za(1,1))/d) va(1,1)=rm(1,1)*9.8/f(1,1)*((za(2,1)-za(1,1))/d) do i=2,m-1
ua(i,1)=0-rm(i,1)*9.8/f(i,1)*((za(i,2)-za(i,1))/d) va(i,1)=rm(i,1)*9.8/f(i,1)*((za(i+1,1)-za(i-1,1))/2/d)
ua(i,n)=0-rm(i,n)*9.8/f(i,n)*((za(i,n)-za(i,n-1))/d) va(i,n)=rm(i,n)*9.8/f(i,n)*((za(i+1,n)-za(i-1,n))/2/d) enddo do j=2,n-1
ua(1,j)=0-rm(1,j)*9.8/f(1,j)*((za(1,j+1)-za(1,j-1))/2/d) va(1,j)=rm(1,j)*9.8/f(1,j)*((za(2,j)-za(1,j))/d)
ua(m,j)=0-rm(m,j)*9.8/f(m,j)*((za(m,j+1)-za(m,j-1))/2/d) va(m,j)=rm(m,j)*9.8/f(m,j)*((za(m,j)-za(m-1,j))/d) enddo return end 最后作图： 原始场
预报场：
可见：两种作图方式作出的原始场、预报场都是一样的。 变高场：
七、结果分析
根据以上原始场和预报场，可以看出东北冷涡向东北方向移动，原东北有一高压脊570线，在预报时刻已经减弱，并在其东部生成更强的575线高压脊，我国东北地区有正变高，且出现一高压脊，，从585等高线大致可知西太平洋副热带高压可能向西北移动，我国华北地区形成575闭合等高线，从变高场可以看出俄罗斯远东地区出现强烈正变高，贝加尔湖西部正变高，东部出现负变高，槽线可能从西到东方向移动，进而影响我国东北地区。