2014年气象统计实习报告

1、 960 实习报告 课程名称：气象统计方法课程实践 姓名: 学号: 班级： 级气科班 * 实习一 求 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做出图形，实习报告中气候场、距平场、均方差场任意 给出两张图 ，图注要清楚，即要注明是哪个

2、时间的图形，并做简单分 析。 注：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:studyformshixioneh500.dat) do iy=1,4 do m=1

3、,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,

4、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:studyformshixioneat.grd,form=unformatted,access=dir ect,recl=nx

5、*ny) irec=0 do t=1,12 irec=irec+1 write (10,rec=irec)(at(i,j,t),i=1,nx),j=1,ny) en ddo close (10) !写入距平场 open( 11,file=d:studyformshixi on exd.grd,form=un formatted,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,j,t,it),i=1,nx),j=1,ny) en ddo en ddo close

6、(11) !写入均方差场 open( 12,file=d:studyformshixi on esx.grd,form=un formatted,access=dir ect,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) en ddo close (12) 1000 2000 3000 format format format 4000 format (2i7) (37F6.2) (37f8.1) (37f7.2) end program sx1 运行结果: r )a

7、t at.gmf J atgrd ；|h500 H500.DAT IhSOO j) h 500 write JhSOOdi -JC9 Fort SojrCe SK0 E SSTPX 2011/9/16 11;02 CTL袁件 1 KB 201 l/3Zie 11:02 GKD文柠 4充KB G stdridard 2014/12/15 17:055目 5 71 KB 编写CTL文件 Evf.ctI dset d:studyformshixieightevf.grd title Coads SSTA E undef -999.0 xdef 18 linear 120 10 ydef 12 lin

8、ear -27.5 5 zdef 1 linear 1000 1 tdef 2 lin ear 1jan 1948 1month vars 1 S 0 99 Coads SST ano maly interp erated using En dvars Tcf.ctl dset d:studyformshixieighttcf.grd title Coads SSTA T undef -999.0 xdef 1 linear 120 10 ydef 1 linear -27.5 5 zdef 1 linear 1000 1 tdef 516 linear 1jan1948 1month var

9、s 2 a 0 99 time coefficient 1 b 0 99 time coefficient 2 endvars 根据 TCL 文件编写相应的 GS 文件 Evf.gs reinit enable print d:studyformshixieightevf.gmf open d:studyformshixieightevf.ctl day=1 while(day=2) set t day d s draw title dayevf print c day=day+1 endwhile disable print Tcfa.gs reinit enable print d:stu

10、dyformshixieighttcf.gmf open d:studyformshixieighttcf.ctl set t 1 516 d a draw title tcfa print disable print Tcfb.gs reinit en able print d:studyformshixieighttcfb.gmf open d:studyformshixieighttcf.ctl set t 1 516 d b draw title tcfb print disable print 效果图 空间场 2evf V： -n nfl. 7DN 丁 J 1* 曲: -P 口3 丁

11、 Y. I_ 心，” 7 由图中可以看出: 1949-1951、1953、1955-1956、1960-1963、1968、1971-1972、1974、1976、1979、1981、 而在 1952、1954、 1985、1989这些年份，赤道太平洋中部至东部的海温都是较正常水平偏低, 1958、1964、1969、1973、1977-1978、1980、1982-1983、1987、1991 这些年份，赤道太平 洋中部至东部的海温都是偏高的，尤其在 1982-1983年，据记载，这段时间出现了厄尔尼诺 现象。 空间场 1 evf J - 耳ii I 1 JQf 1仕 阳A 时间系数 1令爭& 1寸吕0 1?峙 5 19R i*tQ 1 $吕5 1990 run- y-ifl 5H eg 1曲i. ? -.士 1 1令耳 4:、 itpw tcfb 从图中可以看出: 1