程序及程序运行出来的图.docVIP

以下是所编程序及程序运行出来的图 1、程序 program EP real a(160,60),ind(61),r(160),aa(61),bb(61),lat(160),lon(160),hh(144,73,12,61),eu(61),aveeu,aveh(144,73),avet(160),sx,sy(144,73),sz(160),sum,rxy(144,73),l(320),rxz(160),t(144,73) character*8 id(160) open(2,file=g:\shixi\3\data\t1601.txt) open(3,file=g:\shixi\3\ind.grd,form=binary) open(4,file=g:\shixi\3\data\lat_lon.txt) open(5,file=g:\shixi\3\r160.grd,form=binary) open(10,file=g:\shixi\3\data\hgt500.grd,form=binary) !ccccccccccccccc 读数据（指数、经纬度、站温度） do it=1,61 do k=1,12 do j=1,73 do i=1,144 read(10)hh(i,j,k,it) enddo;enddo;enddo;enddo write(*,*)ok1 !ccc a:160站气温(1951~2010年) ind：指数序列(1948~2008年) iyb=1951 read(2,*)((a(i,j),i=1,160),j=1,60) do i=1,160 do j=1,60 a(i,j)=a(i,j)/10 end do end do close(2) print*,ok2 read(4,*) (l(i),i=1,320) do i=1,160 lat(i)=l(2*i-1) lon(i)=l(2*i) end do print*,ok3 !ccccccccccccccc 编程求EU（欧亚）遥相关指数 k=1 j=59 do it=1,61 eu(it)=(-1)*hh(9,j,k,it)/4.0+(hh(31,j,k,it)/2.0)-(hh(59,53,k,it)/4.0) end do !cccccccccc将EU（欧亚）遥相关指数进行标准化 aveeu=0.0 do it=1,61 aveeu=aveeu+eu(it) end do aveeu=aveeu/61 sx=0.0 do it=1,61 sx=sx+(eu(it)-aveeu)**2 end do sx=sqrt(sx/61) do it=1,61 eu(it)=(eu(it)-aveeu)/sx end do !ccccccccccccc存储标准化的EU（欧亚）遥相关指数数据 do it=1,61 write(3) eu(it) enddo close(3) write(*,*)ok 4 !ccccccccccc将高度场进行标准化 k=1 do i=1,144 do j=1,73 do it=1,61 aveh(i,j)=aveh(i,j)+hh(i,j,k,it) end do aveh(i,j)=aveh(i,j)/61 end do end do do i=1,144 do j=1,73 do it=1,61 sy(i,j)=sy(i,j)+(hh(i,j,k,it)-aveh(i,j))**2 end do sy(i,j)=sqrt(sy(i,j)/61) end do end do do i=1,144 do j=1,73 do it=1,61 hh(i,j,k,it)=(hh(i,j,k,it)-aveh(i,j))/sy(i,j) end do end do end do !cccccccccccEU遥相关指数与同期环流场（hPa高度场或海平面气压场）的相关系数 k=1 do