实?习?报?告?书
课程名称：气象统计方法课程实践
姓名： 学号：
班级： 级气科 班
\\\\*实习一 求500hPa高度场气候场、距平场和均方差场 实习时间：第9周周三1、2节 1. 资料介绍
有一500hPa高度场资料，文件名h500.dat，范围：60~150E，0~40N.
时段：1982.1~1985.12共48个月。水平分辨率：2.5*2.5，格点数：2．要求
编fortran程序，求500hPa高度场的 （1）气候场； （2）距平场； （3）均方差场。
37*17。 并能用Grads做出图形,实习报告中气候场、距平场、均方差场任意给出两张图，图注要清楚，即要注明是哪个时间的图形，并做简单分析。
注：h500.For给出了如何用fortran读取ASCII码资料h500.dat.
? FORTRAN program sx1 implicit none integer nx,ny,mo,yr
parameter(nx=37,ny=17,mo=12,yr=4) real var(nx,ny,mo,yr)
real at(nx,ny,mo),xd(nx,ny,mo,yr),sx(nx,ny,mo) integer i,j,m,t,it,iy,irec
open(5,file='d:\\study\\form\\shixione\\h500.dat') do iy=1,4 do m=1,12 read(5,1000)
read(5,3000)((var(i,j,m,iy),i=1,nx),j=1,ny) enddo
enddo close(5) !计算气候场at do t=1,12 do j=1,ny do i=1,nx at(i,j,t)=0 do it=1,4
at(i,j,t)=at(i,j,t)+var(i,j,t,it) enddo
at(i,j,t)=at(i,j,t)/4 enddo enddo enddo !求距平场xd do t=1,12 do j=1,ny