栖霞市长春湖地块可研报告.ppt

气候变化及其诊断实习报告
姓名:
学号:20091349026
班级:2009级气候一班
实习一
(1)实习目的:编fortran程序,求500hPa高度场的气候场、距平场、标准化场,标准差场,并能用Grads做出图形。
(2)程序:
program dl
implicit none
parameter N=37,M=17,nyr=4,nmon=12
real h(N,M,nmon,nyr),s(N,M,nmon),qh(N,M,nmon),bzc(N,M,nmon,nyr),jp(N,M,nmon,nyr)
integer ix,iy,it,itt
open(1,file='f:/sx/h500/',form='binary')
open(2,file='f:/sx/h500/')
open(3,file='f:/sx/h500/',form='binary')
open(4,file='f:/sx/h500/',form='binary')
open(5,file='f:/sx/h500/',form='binary')
open(6,file='f:/sx/h500/',form='binary')
do itt=1,nyr
do it=1,nmon

read(1) ((h(ix,iy,it,itt),ix=1,N),iy=1,M)
write(2,1000) ((h(ix,iy,itt,itt),ix=1,N),iy=1,M)
end do
end do
call qhc(h,qh)
call jpc(h,s,qh,jp)
call (h,s,qh,jp,bzc)
call bzhh(h,qh,s,jp,bzc)
1000 FORMAT()
end program
!气候场
subroutine qhc(h,qh)
implicit none
parameter N=37,M=17,nyr=4,nmon=12
real h(N,M,nmon,nyr),s(N,M,nmon),qh(N,M,nmon)
integer ix,iy,it,itt
s=0
do it=1,nmon
do ix=1,N
do iy=1,M

do itt=1,nyr
s(ix,iy,it)=s(ix,iy,it)+h(ix,iy,it,itt)
end do
qh(ix,iy,it)=s(ix,iy,it)/nyr
end do
end do
end do
do it=1,nmon
write(3) ((qh(ix,iy,it),ix=1,N),iy=1,M)
end do
1000 FORMAT()
end subroutine qhc
!求距平场
subroutine jpc(h,s,qh,jp)
implicit none
parameter N=37,M=17,nyr=4,nmon=12
real h(N,M,nmon,nyr),s(N,M,nmon),qh(N,M,nmon),jp(N,M,nmon,nyr)
integer ix,iy,it,itt
call qhc(h,qh)

do it=1,nmon
do iy=1,M
do ix=1,N
do itt=1,nyr
jp(ix,iy,it,itt)=h(ix,iy,it,itt)-qh(ix,iy,it)

end do
end do
end do
end do
do itt=1,nyr
do it=1,nmon
write(4) ((jp(ix,iy,itt,itt),ix=1,N),iy=1,M)
end do
end do
1000 FORMAT()
end subroutine jpc
!求标准差场
subroutine (h,s,qh,jp,bzc)
implicit none
parameter N=37,M=17,nyr=4,nmon=12
real h(N,M,nmon,nyr),s(N,M,nmon),qh(N,M,nmon),jp(N,M,nmon,nyr),bzc(N,M,nmon)
integer ix,iy,it,itt
call qhc(h,qh)
call jpc(h,s,qh,jp)
bzc=0

do it=1,