温度.涡度平流计算程序，fortran

温度.涡度平流计算程序，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