2014年气象统计实习报告.pdf

实 习 报 告 书 课程名称：气象统计方法课程实践 姓名： 学号： 班级： 级气科 班 * 实习一 求 500hPa 高度场气候场、距平场和均方差场 实习时间：第 9 周周三 1、2 节 1. 资料介绍 有一 500hPa 高度场资料，文件名 h500.dat，范围：60~150E，0~40N. 时段： 1982.1~1985.12 共 48 个月。水平分辨率： 2.5*2.5 ，格点 数： 37*17 。 2．要求 编 fortran 程序，求 500hPa 高度场的 （1）气候场； （2）距平场； （3）均方差场。 并能用 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 do i=1,nx xd(i,j,t,1)=0 do it=1,4 xd(i,j,t,it)=var(i,j,t,it)-at(i,j,t) enddo enddo enddo enddo ! 求均方差场 sx do t=1,12 do j=1,ny do i=1,nx sx(i,j,t)=0 do it=1,4 sx(i,j,t)=sx(i,j,t)+(var(i,j,t,it)-at(i,j,t))**2 enddo sx(i,j,t)=sqrt(sx(i,j,t)/4) enddo enddo enddo ! 写入气候场 open(10,file=d:\study\form\shixione\at.grd,form=unformatted,access=dir ect,recl=nx*ny) irec=0 do t=1,12 irec=irec+1 write (10,rec=irec)((at(i,j,t),i=1,nx),j=1,ny) enddo close (10) ! 写入距平场 open(11,file=d:\study\form\shixione\xd.grd,form=unformatted,access=dir ect,recl=nx*ny) irec=0 do it=1,4 do t=1,12 irec=irec+1 write (11,rec=irec)((xd(i,