实习一(2)

endwhile year=year+1 endwhile

'disable print' ;

距平场图：

1983年6月距平场在日本地区出现低压中心，在我国南部出现高压中心,在亚洲西北部也有高压中心。赤道至25°N间以及25°N-40°N，60°E-100°E间基本都是正距平，而在25°N-40°N，100°E-150°E间基本都是负距平。

1984年7月距平场在亚洲大陆西部、日本地区、赤道西太平洋地区形成低压中心，太平洋西北部形成高压中心。

给s配的ctl文件： dset ^d:\\ex1\\s.grd undef -9.99E+33

title NCEP/NCAR REANALYSIS PROJECT xdef 37 linear 60.000 2.500 ydef 17 linear 0.000 2.500 zdef 1 levels 500

tdef 12 linear JAN1982 12mo vars 1

s 1 99 H500 endvars

给s配的gs文件： 'reinit'

'open d:\\ex1\\s.ctl'

'enable print d:\\ex1\\s.gmf' mon=1

while(mon<=12) 'set t 'mon'' 'd s'

'draw title junfangchachang of 'mon' ' 'print' 'c'

mon=mon+1 endwhile

'disable print' ;

均方差场图：

一月份高度的均方差场整体呈现南小北大的状态。说明低纬地区高度的波动幅度比较小，而中高纬地区高度的波动比较大。

八月份高度的均方差场在亚洲大陆西部有极大值，在30°N处包括赤道-30°N 、60°E-85°E这些区域高度的波动幅度比较小，30°N以南以北地区高度的波动幅度较大。

实习二：相关系数 Fortran程序如下：

program ex2

integer,parameter::n=20,p=10 integer i,j,t1,t2,t3 real

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,fc2 data

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,n s1=s1+a(i) enddo ave1=s1/n do i=1,n s2=s2+b(i) enddo ave2=s2/n !求距平 do i=1,n

jpa(i)=a(i)-ave1 enddo do i=1,n

jpb(i)=b(i)-ave2 enddo

!求相关系数 do i=1,n

sum1=sum1+jpa(i)*jpa(i) sum2=sum2+jpb(i)*jpb(i) sum3=sum3+jpa(i)*jpb(i) enddo

r=sum3/(sqrt(sum1*sum2))

print*,'中国1970-1989年年平均和冬季平均气温的相关系数为r=',r !求方差 fc1=sum1/n fc2=sum2/n

!求自协方差 do j=1,p do i=1,n-j

zxfc1(j)=zxfc1(j)+jpa(i)*jpa(i+j) enddo

zxfc1(j)=zxfc1(j)/(n-j) enddo do j=1,p do i=1,n-j

zxfc2(j)=zxfc2(j)+jpb(i)*jpb(i+j) enddo

zxfc2(j)=zxfc2(j)/(n-j) enddo

!自相关系数 do i=1,p

zxgxs1(i)=zxfc1(i)/fc1 zxgxs2(i)=zxfc2(i)/fc2 enddo

!落后交叉协方差 do j=1,p do i=1,n-j

lhxfc(j)=lhxfc(j)+jpa(i)*jpb(i+j) enddo

lhxfc(j)=lhxfc(j)/(n-j) enddo

!落后相关系数 do i=1,p

lhxgxs(i)=lhxfc(i)/(sqrt(fc1*fc2)) enddo

print*,'年平均气温不同滞后时刻所对应的自相关系数为：' print *,((zxgxs1(i),','),i=1,P)

百度搜索“77cn”或“免费范文网”即可找到本站免费阅读全部范文。收藏本站方便下次阅读，免费范文网，提供经典小说综合文库实习一(2)在线全文阅读。