内容发布更新时间 : 2024/12/22 14:48:41星期一 下面是文章的全部内容请认真阅读。

实?习?报?告?书

课程名称：气象统计方法课程实践

姓名： 学号：

班级： 级气科 班

\\\\*实习一 求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