资源描述：
《2014年气象统计实习报告.doc》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、实 习 报 告 书课程名称：气象统计方法课程实践姓名：学号：班级：级气科班\*实习一求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做出图形,实习报告中气候场、距平场、均方差场任意给出两张图，图注要清楚，即
2、要注明是哪个时间的图形，并做简单分析。注：h500.For给出了如何用fortran读取ASCII码资料h500.dat.nFORTRANprogramsx1implicitnoneintegernx,ny,mo,yrparameter(nx=37,ny=17,mo=12,yr=4)realvar(nx,ny,mo,yr)realat(nx,ny,mo),xd(nx,ny,mo,yr),sx(nx,ny,mo)integeri,j,m,t,it,iy,irecopen(5,file='d:studyformshixione
3、h500.dat')doiy=1,4dom=1,12read(5,1000)read(5,3000)((var(i,j,m,iy),i=1,nx),j=1,ny)enddoenddoclose(5)!计算气候场atdot=1,12doj=1,nydoi=1,nxat(i,j,t)=0doit=1,4at(i,j,t)=at(i,j,t)+var(i,j,t,it)enddoat(i,j,t)=at(i,j,t)/4enddoenddoenddo!求距平场xddot=1,12doj=1,nydoi=1,nxxd(i,j,t,1)=0
4、doit=1,4xd(i,j,t,it)=var(i,j,t,it)-at(i,j,t)enddoenddoenddoenddo!求均方差场sxdot=1,12doj=1,nydoi=1,nxsx(i,j,t)=0doit=1,4sx(i,j,t)=sx(i,j,t)+(var(i,j,t,it)-at(i,j,t))**2enddosx(i,j,t)=sqrt(sx(i,j,t)/4)enddoenddoenddo!写入气候场open(10,file='d:studyformshixioneat.grd',form='u
5、nformatted',access='direct',recl=nx*ny)irec=0dot=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=0doit=1,4dot=1,12irec=irec+1write(11,rec=ir
6、ec)((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=0dot=1,12irec=irec+1write(12,rec=irec)((sx(i,j,t),i=1,nx),j=1,ny)enddoclose(12)1000format(2i7)2000format(37F6.2)3
7、000format(37f8.1)4000format(37f7.2)endprogramsx1n运行结果：nGrads文件气候场'reinit''enableprintd:studyformshixioneat.gmf''opend:studyformshixioneat.ctl''setgridoff''setgradsoff''setlat040''setlon60150''setlev500'mon=1while(mon<=12)'sett'mon'''dh''drawtitle1982year'mon'm
8、onth''print''c'mon=mon+1endwhile'disableprint';距平场'reinit''enableprintd:studyformshixionesx.gmf''opend:studyformshixi
