*a(i,j))/4.0 enddo enddo doi=2,m-1 doj=2,n-1 a(i,j)=w(i,j) enddo enddo doi=2,m-1 doj=2,n-1 w(i,j)=a(i,j)+(-s)*(a(i-1,j)+a(i+1,j)+a(i,j-1)+a(i,j+1)-4*a(i,j))/4.0 enddo enddo doi=2,m-1 doj=2,n-1 a(i,j)=w(i,j) enddo enddo endif return end根据得到的数据文件画出的高度场与风场?图一.原始高度场图二.预报高度场图三.原始风场图四.预报风场图五.预报高度场与风矢量场合成上面的图均是500hpa08时的要素场;结果分析:由图一和图二,图三和图四,我们可以清楚的看到:4月29日08时在鄂霍次克海的西岸,有一个大的低压中心,其中心强度达到5250,在贝加尔湖附近有一条大的低压槽;4月30日08时,在预报场上,鄂霍次克海西岸的低压向东移动,到达鄂霍次克海的上空,中心强度达到5200,比29日加强,贝加尔湖的低压槽向南,向东移动。可能造成一次冷空气的南下过程。我们拿30日的预报场与30日的实况场(书本上)做一下对比,可以看到预报场的低压中心偏北大概偏5个纬距左右,有较大的偏差,模拟效果很不好;但是,贝加尔湖附近的低压槽预报场与实况场非常的近似,故而模拟效果较好。而我们需要做对中国的预报,影响系统主要是贝加尔湖地区的低压槽影响,所以对整个预报来说,还是有很好的参考价值的。5.小结通过这次实习,加强了我对于整个正压原始模式的整个计算流程的认识,而且通过这次实习我的编程与分析能力也得到了很大的加强,与此同时,也加深了对数值预报理论的理解和认识。我也了解到数值预报需要处理的一些问题,以及我们需要面对的问题和如何解决这些问题.并且可以知道我们需要在哪些方面服出更多的努力.
全文预览
