温度.涡度平流计算程序,fortran
REAL,PARAMETER:: Omega=7.292e-5,R=6371e+3,PI=3.1415926,Delta=1.00
INTEGER,PARAMETER :: nx=111,ny=61,nz=7,nt=4
INTEGER WtTopLevel(NX,NY,NT)
REAl lat(ny),lon(nx),sigmadeltD
CHARACTER timefile(4)*8 ,levelfile(5)*1,P(NZ)*2
real,allocatable :: temper(:,:,:,:),q(:,:,:,:),u(:,:,:,:),v(:,:,:,:),height(:,:,:,:),qu(:,:,:,:),qv(:,:,:,:),adq(:,:,:,:),adqv(:,:,:,:),vorg(:,:,:,:),voro(:,:,:,:),div(:,:,:,:),temp(:,:,:,:),vorop(:,:,:,:)
allocate (temper(NX,NY,NZ,NT),q(NX,NY,nz,NT),u(NX,NY,NZ,NT),v(NX,NY,NZ,NT),height(NX,NY,NZ,NT),qu(NX,NY,NZ,NT),qv(NX,NY,NZ,NT),adq(NX,NY,NZ,NT),adqv(NX,NY,NZ,NT),vorg(NX,NY,NZ,NT),voro(NX,NY,NZ,NT),div(NX,NY,NZ,NT),temp(NX,NY,NZ,NT),vorgp(NX,NY,NZ,NT) )
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
DATA P/'00','92','85','70','50','40','30'/
DATA timefile/'00041300','00041312','00041400', '00041412'/
do i=1,nx
lon(i)=50.0+(i-1)*Delta
lon(i)=lon(i)*pi/180.
enddo
do i=1,ny
lat(i)=10.0+(i-1)*Delta
lat(i)=lat(i)*pi/180.
enddo
1