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

? ? 实?习?报?告?书 课程名称：气象统计方法课程实习一 求500hPa高度场气候场、距平场和均方差场 实习时间：第周周、节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=direct,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=direct,recl=nx*ny) irec=0 do it=1,4 do t=1,12 irec=irec+1 write(11,rec=irec)((xd(i,j,t,it),i=1,nx),j=1,ny) enddo enddo close(11) ！写入均方差场 open(12,file=d:\study\form\shixione\sx.grd,form=unformatted,access=direct,recl=nx*ny) irec=0 do t=1,12 irec=irec+1 write(12,rec=irec)((sx(i,j,t),i=1,nx),j=1,ny) enddo close(12) 1000 format(2i7) 2000 format(37F6.2) 3000 format(37f8.1) 4000 format(37f7.2) end program sx1 运行结果： Grads文件 reinit enable print d:\study\form\shixione\at.gmf open d:\study\form\shixione\at.ctl set grid off set grads off set lat 0 40 set lon 60 150 set lev 500 mon=1 while(mon=12) set t mon d h draw title 1982yearmonmonth print c mon=mon+1 endwhile disable print ; 距平场 reinit enable print d:\study\form\shixione\sx.gmf open d:\study