Fotran 气象笔记

注：2012-2016年抄，已不用fotran，停止更新

Fortran 求春季高度距平场

program Height_spring

!----------------------------------------------------------

integer :: nx=144,ny=73,mon=12,yr=68

integer i,j,iz,it,t  !it:yr循环变量，iz：mon循环变量

real,allocatable :: hgt850(:,:,:,:),hgt500(:,:,:,:),hgt200(:,:,:,:)

real,allocatable :: spring_850ave(:,:),spring_500ave(:,:),spring_200ave(:,:)

real,allocatable :: spring_850ano(:,:,:),spring_500ano(:,:,:),spring_200ano(:,:,:)

allocate (hgt850(nx,ny,mon,yr),spring_850ave(nx,ny),spring_850ano(nx,ny,64))

allocate (hgt500(nx,ny,mon,yr),spring_500ave(nx,ny),spring_500ano(nx,ny,64))

allocate (hgt200(nx,ny,mon,yr),spring_200ave(nx,ny),spring_200ano(nx,ny,64))

!----------------------------------------------------------

open(1,file='hgt850.grd',form='binary')

do it=1,yr

do iz=1,mon

do j=1,ny

read(1)(hgt850(i,j,iz,it),i=1,nx)

enddo;enddo;enddo

!print*,hgt850

!pause

 

open(2,file='hgt500.grd',form='binary')

do it=1,yr

do iz=1,mon

do j=1,ny

read(2)(hgt500(i,j,iz,it),i=1,nx)

enddo;enddo;enddo

 

open(3,file='hgt200.grd',form='binary')

do it=1,yr

do iz=1,mon

do j=1,ny

read(3)(hgt200(i,j,iz,it),i=1,nx)

enddo;enddo;enddo

 

!----------------------------------------------------------

!                        850hPa

!----------------------------------------------------------

!---求平均场-----------------------------------------------

!----------------------------------------------------------

do it=4,67 !1951-2014

do iz=3,5  !春

do j=1,ny

do i=1,nx

if (hgt850(i,j,iz,it)/=-9.96921e+36) then

spring_850ave(i,j)=spring_850ave(i,j)+hgt850(i,j,iz,it)/(64*3)

else 

spring_850ave(i,j)=-9.96921e+36

end if

enddo;enddo;enddo;enddo

!---求距平场-----------------------------------------------

!----------------------------------------------------------

do it=4,67 !1951-2014

do iz=3,5

t=it-3

do j=1,ny

do i=1,nx

if (hgt850(i,j,iz,it)/=-9.96921e+36) then

spring_850ano(i,j,t)=spring_850ano(i,j,t)+hgt850(i,j,iz,it)/3

else 

spring_850ano(i,j,t)=-9.96921e+36

end if

enddo;enddo;enddo;enddo

 

do t=1,64 !1951-2014

do j=1,ny

do i=1,nx

if (spring_850ano(i,j,t)/=-9.96921e+36 .and. spring_850ave(i,j)/=-9.96921e+36) then

spring_850ano(i,j,t)=spring_850ano(i,j,t)-spring_850ave(i,j)

else

spring_850ano(i,j,t)=-9.96921e+36

end if

enddo;enddo;enddo

 

!----------------------------------------------------------

 

 

open(11,file='D:\biye\output\850\hgt_spring_850ave.grd',form='binary')

do j=1,ny

write(11)(spring_850ave(i,j),i=1,nx)

enddo

 

open(12,file='D:\biye\output\850\hgt_spring_850ano.grd',form='binary')

do it=1,64

do j=1,ny

write(12)(spring_850ano(i,j,it),i=1,nx)

enddo;enddo

 

!open(13,file='D:\biye\output\850\