实习报告书课程名称：气象统计方法课程实践姓名：学号： 班级： 级气科 班*实习一 求500hPa高度场气候场、距平场和均方差场 实习时间：第9周周三1、2节1. 资料介绍有一500hPa高度场资料，文件名h500.dat，范围：60150E，040N.时段：1982.11985.12共48个月。水平分辨率：2.5*2.5，格点数：37*17。2要求编fortran程序，求500hPa高度场的（1）气候场；（2）距平场；（3）均方差场。并能用Grads做出图形,实习报告中气候场、距平场、均方差场任意给出两张图，图注要清楚，即要注明是哪个时间的图形，并做简单分析。注：h500.For给出了如何用fortran读取ASCII码资料h500.dat.n FORTRANprogram sx1implicit noneinteger nx,ny,mo,yrparameter(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,irecopen(5,file=d:studyformshixioneh500.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)enddoenddo close(5)!计算气候场atdo 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 enddoenddo!求距平场xddo 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 enddoenddo!求均方差场sxdo 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 enddoenddo!写入气候场open(10,file=d:studyformshixioneat.grd,form=unformatted,access=direct,recl=nx*ny)irec=0do t=1,12irec=irec+1write(10,rec=irec)(at(i,j,t),i=1,nx),j=1,ny)enddoclose(10)!写入距平场open(11,file=d:studyformshixionexd.grd,form=unformatted,access=direct,recl=nx*ny)irec=0do it=1,4 do t=1,12irec=irec+1write(11,rec=irec)(xd(i,j,t,it),i=1,nx),j=1,ny) enddoenddoclose(11)！写入均方差场open(12,file=d:studyformshixionesx.grd,form=unformatted,access=direct,recl=nx*ny)irec=0do t=1,12irec=irec+1write(12,rec=irec)(sx(i,j,t),i=1,nx),j=1,ny)enddoclose(12)1000 format(2i7)2000 format(37F6.2)3000 format(37f8.1)4000 format(37f7.2)end program sx1n 运行结果：n Grads文件气候场reinitenable print d:studyformshixioneat.gmfopen d:studyformshixioneat.ctlset grid offset grads offset lat 0 40set lon 60 150set lev 500mon=1while(mon=12)set t mond hdraw title 1982yearmonmonthprintcmon=mon+1endwhiledisable print;距平场reinitenable print d:studyformshixionesx.gmfopen d:studyformshixionesx.ctlset grid offset grads offset lat 0 40set lon 60 150set lev 500year=1982while(year=1984)mon=1while(mon=12)set t mond hdraw title 500hPa yearyearmonmonth anomalyprintcmon=mon+1endwhileyear=year+1endwhiledisable print;均方差reinitenable print d:studyformshixionesx.gmfopen d:studyformshixionesx.ctlset grid offset grads offset lat 0 40set lon 60 150set lev 500mon=1while(monmax_y) THEN max_y=rxy_ty(ty) k=ty END IFEND DOPRINT (全年平均气温绝对值最大自相关系数rxy_ty=,f7.4,/,滞后时间长度 k=,I2),rxy_ty(k),kk=0DO tw=1,N/2 DO i=1,N-tw rtw(i)=(w(i)-avr_w)/sw)*(w(i+tw)-avr_w)/sw) rxy_tw(tw)=rxy_tw(tw)+rtw(i) END DO rxy_tw(tw)=rxy_tw(tw)/(N-tw) rxy_tw(tw)=ABS(rxy_tw(tw) IF(rxy_tw(tw)max_w) THEN max_w=rxy_tw(tw) k=tw END IFEND DOPRINT (冬季平均气温绝对值最大自相关系数rxy_tw=,f7.4,/,滞后时间长度 k=,I2),rxy_tw(k),kk=0ENDn 运行成果：*实习四