enddo
enddo
!Calculating the Climatological Field
do it=1,nt
do i=1,nx
do j=1,ny
sum=0
do iz=1,nz
sum=sum+var(i,j,iz,it)
enddo
cl(i,j,it)=sum/4
enddo
enddo
enddo
! tput climate-file
冬季（此处以1月为代表）等高线分布整体平缓，表明高度场分布相对均匀，且北部接近极地位势高度低，赤道地区位势高度高，这与太阳直射点在1月在南半球，使北半球整体的辐射吸收随纬度增加而减小有关。北半球的气块受热随纬度递增而递减，因而膨胀率递减，故南方接近赤道地区的气体膨胀大，位势高，而北方近极地地区气体膨胀小，位势低。
do iz=1,nz
do it=1,nt
write(13)((jp(i,j,iz,it),i=1,nx),j=1,ny)
enddo
enddo
!Calculating the Mean-square Deviation
do it=1,nt
do i=1,nx
do j=1,ny
sum=0
do iz=1,nz
气象统计方法实习报告
实习一求500hPa高度场气候场、距平场和均方差场
1、资料介绍
有一500hPa高度场资料，文件名h500.dat，范围：60~150E，0~40N.
时段：1982.1~1985.12共48个月。水平分辨率：2.5*2.5，格点数：37*17。
2．要求
编fortran程序，求500hPa高度场的
夏季（此处以多年平均的8月气候场为代表）在太平洋上有副热带高压，副高是深厚的系统，所以可以一直延伸到500hPa高度。东亚地区存在明显的西风带长波槽，即东亚大槽，东亚大槽的移动和变化配合副热带高压和夏季风的影响，会使我国大部分地区天气造成强烈变化，如形成大范围降水或强对流天气。同时，孟加拉湾处存在一低涡。由于高原的阻挡作用，这一系统对我国的影响并不显著。图中副高脊线8月的多年平均位置位于30°N以北，且东部长波槽位于110°E附近，故由气候场分析，华北地区位于长波槽前，又槽线受到副高的阻挡作用，因而华北地区容易形成降水。
（1）气候场；
（2）距平场；
（3）均方差场。
并能用Grads做出图形,实习报告中气候场、距平场、均方差场任意给出两张图，图注要清楚，即要注明是哪个时间的图形，并做简单分析。
注：h500.For给出了如何用fortran读取ASCII码资料h500.dat.
3、实习结果
1）、FORTRAN源程序
program ex_grads
从图中还可以看出冬夏季的位势高度分布存在明显差异，这与太阳直射点的年纪变化密切相关。
②、气候场
2月
8月
结果分析：
气候场是多年数据中同时段的平均值序列，表征了区域内多年平均的位势高度变化。
从气候场图形可以看出多年平均的500hPa高度场中，冬季（此处以多年平均的2月气候场为代表）等高线较平直，大陆上等高线较稀疏，而海洋上等高线较密集，这表明大陆上空的位势高度变化率较海洋小。此外，冬季等位势高度线分布平直，还说明冬季的天气变化显著。
sum=sum+(jp(i,j,iz,it))**2
enddo
jfc(i,j,it)=sqrt(sum/4)
enddo
enddo
enddo
!Output mean-square deviation-file
open(14,file='g:\gradsdata\deviation.grd',form='binary')
open(10,file='g:\gradsdata\h500.dat')
do iz=1,nz
do it=1,nt
read(10,1000)
read(10,3000) ((var(i,j,iz,it),i=1,nx),j=1,ny)
enddo
enddo
1000 format(2i7)
2000 format(37f6.2)
do i=1,nx
do j=1,ny
jp(i,j,iz,it)=var(i,j,iz,it)-cl(i,j,it)
enddo
enddo
enddo
enddo
open(13,file='g:\gradsdata\anomaly.grd',form='binary')
!Output anomaly-file
夏季（此处以7月为代表），大洋上空出现副热带高压（588线位置），在东亚地区存在较为明显的位势高度槽，即东亚大槽。夏季在30°N以北的地区位于西风带中，从图中可看出明显的西风带长波特征。除东亚大槽外，在中亚地区也存在一长波槽，这些槽线发生长波调整时，会在部分地区发生较剧烈的天气变化。此外，从图中可以看出，1982年7月副热带高压脊线的平均位置位于25°N ,125°E附近，我国华东地区位于副高北侧西南气流控制下，西南气流为水汽输送的主要通道，该地区发生降水较为频繁。
open(12,file='g:\gradsdata\climate.grd',form='binary')
do it=1,nt
write(12)((cl(i,j,it),i=1,nx),j=1,ny)
enddo
!Calculating the Anomaly
do iz=1,nz
do it=1,nt
'enable print g:\gradsdata\***.gmf'（***为所求场名称）
i=1
while(i<=48(或12))
'set t 'i
'd h'
'print'
'c'
i=i+1
endwhile
'disable print'
;
（3）、实习结果
①、原始场
1982年1月
1982年7月
结果分析：
implicit none
integer,parameter::nx=37,ny=17,nz=4,nt=12
integer i,j,iz,it
real var(nx,ny,nz,nt),cl(nx,ny,nt),sum,jp(nx,ny,nz,nt),jfc(nx,ny,nt)
! Opening file
do it=1,nt
write(14)((jfc(i,j,it),i=1,nx),j=1,ny)
enddo
end
(2)、grads文件
'open g:\gradsdata\***.ctl' (***为所求场对应的ctl文件名)
'set lat 0 40'
'set lon 60 150'
'set lev 500'
3000 format(37f8.1)
4000 format(37f7.2)
close(10)
!Output
open(16,file='g:\gradsdata\h500.grd',form='binary')
do iz=1,nz
do it=1,nt
write(16)((var(i,j,iz,it),i=1,nx),j=1,ny)