气象统计方法实习BD.doc
实习一：气候场、距平场、均方差场编程如下:parameter(ii=37,jj=17,mon=12,year=4) real var(ii,jj,mon,year),ave(ii,jj,mon),jp(ii,jj,mon,year) real s(ii,jj,mon) integer i,j,iy,m open(5,file=d:ex1h500.dat) open(6,file=d:ex1ave.grd,form=binary)open(7,file=d:ex1jp.grd,form=binary)open(8,file=d:ex1s.grd,form=binary)open(12,file=d:ex1outall.grd,form=binaryopen(9,file=d:ex1ave.txt)open(10,file=d:ex1jp.txt) open(11,file=d:ex1s.txt) ！读数据 DO iy=1,4 do m=1,12!ccc read h500 read(5,1000) read(5,2000) (var(i,j,m,iy),i=1,ii),j=1,jj) enddoenddo！计算气候场 do j=1,jj do i=1,ii do m=1,12 ave(i,j,m)=var(i,j,m,1)+var(i,j,m,2)+var(i,j,m,3)+var(i,j,m,4) ave(i,j,m)=ave(i,j,m)/4.0 enddo enddo enddo！计算距平场 do iy=1,4 do m=1,12 do j=1,jj do i=1,ii jp(i,j,m,iy)=var(i,j,m,iy)-ave(i,j,m) enddo enddo enddo enddo！计算均方差场 do j=1,jj do i=1,ii do m=1,12 s(i,j,m)=jp(i,j,m,1)*jp(i,j,m,1)+jp(i,j,m,2)*jp(i,j,m,2)+jp(i,j /,m,3)*jp(i,j,m,3)+jp(i,j,m,4)*jp(i,j,m,4) s(i,j,m)=s(i,j,m)/4.0 s(i,j,m)=sqrt(s(i,j,m) enddo enddo enddo do iy=1,4 do m=1,12 write(6)(ave(i,j,m),i=1,ii),j=1,jj) write(7)(jp(i,j,m,iy),i=1,ii),j=1,jj) write(8)(s(i,j,m),i=1,ii),j=1,jj) write(9,2000)(ave(i,j,m),i=1,ii),j=1,jj) write(10,2000)(jp(i,j,m,iy),i=1,ii),j=1,jj) write(11,2000)(s(i,j,m),i=1,ii),j=1,jj) write(12)(ave(i,j,m),i=1,ii),j=1,jj) write(12)(jp(i,j,m,iy),i=1,ii),j=1,jj) write(12)(s(i,j,m),i=1,ii),j=1,jj) enddoenddo1000 format(2i7)2000 format(37f8.1) close(5) close(6) close(7) close(8) close(9) close(10) close(11) close(12) end给ave配的ctl文件：dset d:ex1ave.grdundef -9.99E+33title NCEP/NCAR REANALYSIS PROJECTxdef 37 linear 60.000 2.500ydef 17 linear 0.000 2.500zdef 1 levels 500tdef 12 linear JAN1982 12movars 1ave 1 99 H500endvars给ave配的gs文件：reinitopen d:ex1ave.ctlenable print d:ex1ave.gmfmon=1while(mon<=12)set t mond avedraw title qihouchang of mon printcmon=mon+1endwhiledisable print;气候场图：一月份高度的气候场呈现南高北低的状态，陆地上的高度场比较稀疏，而在西太平洋上高度场比较密集。八月份高度的气候场呈现东高西低的状态，在我国东北部以北以及印度东北部出现低压中心，而在赤道西太平洋地区出现高压中心。35N以北高度分布很密集，而35N以南比较稀疏。给jp配的ctl文件：dset d:ex1jp.grdundef -9.99E+33title NCEP/NCAR REANALYSIS PROJECTxdef 37 linear 60.000 2.500ydef 17 linear 0.000 2.500zdef 1 levels 500tdef 48 linear JAN1982 1movars 1jp 1 99 H500endvars给jp配的gs文件：reinitopen d:ex1jp.ctlenable print d:ex1jp.gmfyear=1982while(year<=1985)mon=1while(mon<=12)set t mond jpdraw title jupingchang of year.monprintcmon=mon+1endwhileyear=year+1endwhiledisable print;距平场图：1983年6月距平场在日本地区出现低压中心，在我国南部出现高压中心,在亚洲西北部也有高压中心。赤道至25N间以及25N-40N，60E-100E间基本都是正距平，而在25N-40N，100E-150E间基本都是负距平。1984年7月距平场在亚洲大陆西部、日本地区、赤道西太平洋地区形成低压中心，太平洋西北部形成高压中心。给s配的ctl文件：dset d:ex1s.grdundef -9.99E+33title NCEP/NCAR REANALYSIS PROJECTxdef 37 linear 60.000 2.500ydef 17 linear 0.000 2.500zdef 1 levels 500tdef 12 linear JAN1982 12movars 1s 1 99 H500endvars给s配的gs文件：reinitopen d:ex1s.ctlenable print d:ex1s.gmfmon=1while(mon<=12)set t mond sdraw title junfangchachang of mon printcmon=mon+1endwhiledisable print;均方差场图：一月份高度的均方差场整体呈现南小北大的状态。说明低纬地区高度的波动幅度比较小，而中高纬地区高度的波动比较大。八月份高度的均方差场在亚洲大陆西部有极大值，在30N处包括赤道-30N 、60E-85E这些区域高度的波动幅度比较小，30N以南以北地区高度的波动幅度较大。实习二：相关系数Fortran程序如下：program ex2integer,parameter:n=20,p=10integer i,j,t1,t2,t3real a(n),b(n),jpa(n),jpb(n),zxfc1(p),zxgxs1(p),zxfc2(p),zxgxs2(p),lhxfc(p),lhxgxs(p)real:s1=0.0,s2=0.0,sum1=0.0,sum2=0.0,sum3=0.0,ave1,ave2,r,fc1,fc2data a/3.40,3.30,3.20,2.90,3.40,2.80,3.60,3.00,2.80,3.00,3.10,3.00,2.90,2.70,3.50,3.20,3.10,2.80,2.90,2.90/data b/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/data jpa/0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0/data jpb/0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0/data zxfc1/0,0,0,0,0,0,0,0,0,0/data zxgxs1/0,0,0,0,0,0,0,0,0,0/data zxfc2/0,0,0,0,0,0,0,0,0,0/data zxgxs2/0,0,0,0,0,0,0,0,0,0/data lhxfc/0,0,0,0,0,0,0,0,0,0/data lhxgxs/0,0,0,0,0,0,0,0,0,0/!求均值do i=1,ns1=s1+a(i)enddoave1=s1/ndo i=1,ns2=s2+b(i)enddoave2=s2/n!求距平do i=1,njpa(i)=a(i)-ave1enddodo i=1,njpb(i)=b(i)-ave2enddo!求相关系数do i=1,nsum1=sum1+jpa(i)*jpa(i)sum2=sum2+jpb(i)*jpb(i)sum3=sum3+jpa(i)*jpb(i)enddor=sum3/(sqrt(sum1*sum2)print*,中国1970-1989年年平均和冬季平均气温的相关系数为r=,r!求方差fc1=sum1/nfc2=sum2/n!求自
