气象统计实习报告[共29页].doc

  实 习 报 告 书课程名称：气象统计方法课程实践姓名：学号： 班级： 级气科 班\\*实习一 求500hPa高度场气候场、距平场和均方差场 实习时间：第9周周三1、2节1. 资料介绍有一500hPa高度场资料，文件名h500.dat，范围：60~150E，0~40N.时段：1982.1~1985.12共48个月水平分辨率：2.5*2.5，格点数：37*172．要求编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:\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)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:\study\form\shixione\at.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:\study\form\shixione\xd.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:\study\form\shixione\sx.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文件气候场'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=1while(mon<=12)'set t 'mon'''d h''draw title 1982year'mon'month''print''c'mon=mon+1endwhile'disable print';距平场'reinit''enable print d:\study\form\shixione\sx.gmf''open d:\study\form\shixione\sx.ctl''set grid off''set grads off''set lat 0 40''set lon 60 150''set lev 500'year=1982while(year<=1984)mon=1while(mon<=12)'set t 'mon'''d h''draw title 500hPa 'year'year'mon'month anomaly''print''c'mon=mon+1endwhileyear=year+1endwhile'disable print';均方差'reinit''enable print d:\study\form\shixione\sx.gmf''open d:\study\form\shixione\sx.ctl''set grid off''set grads off''set lat 0 40''set lon 60 150''set lev 500'mon=1while(mon<=12)'set t 'mon'''d h''draw title 500hPa 1982year'mon'month Mean-square Deviation''print''c'mon=mon+1endwhile'disable print';n 运行结果：上面图中只展示了部分而未全部添加*实习二 计算给定数据资料的简单相关系数和自相关系数 实习时间：第10周周三1、2节根据下表中年平均气温和冬季平均气温的等级数据进行下列计算：1）计算两个气温之间的简单相关系数。

2）分别找出两个气温数据自相关系数绝对值最大的滞后时间长度滞后长度τ最大取10）要求：实习报告中附出简单相关系数或自相关系数程序答案：r=0.47 年平均气温在滞后长度j=7、冬季序列在j=4最大n FORTRAN计算相关系数rPROGRAM EXAMIMPLICIT NONEINTEGER,PARAMETER::N=20INTEGER i,j,k,ty,tw,tywREAL::avr_y=0,avr_w=0,sy=0,sw=0,rxy=0,max_y=0,max_w=0,max_yw=0REAL y(N),w(N)DATA y/3.4,3.3,3.2,2.9,3.4,2.8,3.6,3.0,2.8,3.0,3.1,3.0,2.9,2.7,3.5,3.2,3.1,2.8,2.9,2.9/DATA w/3.24,3.14,3.26,2.38,3.32,2.71,2.84,3.94,2.75,1.83,2.80,2.81,2.63,3.20,3.60,3.40,3.07,1.87,2.63,2.47/REAL syy(N),sww(N),r(N),rty(N),rtw(N),rtyw(N),rxy_ty(N),rxy_tw(N),rxy_tyw(N)!求两数组平均值 DO i=1,N avr_y=avr_y+y(i) avr_w=avr_w+w(i) END DO avr_y=avr_y/N avr_w=avr_w/N!简单相关系数DO j=1,N syy(j)=(y(j)-avr_y)**2 sy=sy+syy(j) sww(j)=(w(j)-avr_w)**2 sw=sw+sww(j)END DO sy=sqrt(sy/N) sw=sqrt(sw/N)DO j=1,N r(j)=((y(j)-avr_y)/sy)*((w(j)-avr_w)/sw) rxy=rxy+r(j)END DO rxy=rxy/NPRINT "(/'1970-1989年全年平均气温与冬季平均气温的简单相关系数rxy=',f5.2)",rxyk=0!自相关系数DO ty=1,N/2 DO i=1,N-ty rty(i)=((y(i)-avr_y)/sy)*((y(i+ty)-avr_y)/sy) rxy_ty(ty)=rxy_ty(ty)+rty(i) END DO rxy_ty(ty)=rxy_ty(ty)/(N-ty) rxy_ty(ty)=ABS(rxy_ty(ty)) IF(rxy_ty(ty)>max_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 运行成果：*实习四。