新手再次求教，偏高年的500hpa高度距平场

勇敢的悲伤

我要画500hpa高度距平场。图画出来了，但是和原图不一样。不知道数据处理的问题还是什么，大家帮看看FORTRAN写的对不对。
program main
implicit none
integer,parameter:: iix=145,iiy=73,mt=696,m=12,y=58
integer i,j,t,yr,mo,kk
integer::anum=0,bnum=0
real atemp(mt),hgt(iix,iiy,mt),zhgt(iix,iiy),phgt(iix,iiy),zzhgt(iix,iiy),jphgt(iix,iiy),hgt1(iix,iiy,m,y),ghgt(iix,iiy),num(iix,iiy)
open(13,file='f:\1\hgt.grd',form='binary')
open(14,file='f:\1\temp1.grd',form='binary')
open(15,file='f:\1\ghgt.txt')
open(16,file='f:\1\ghgt.grd',form='binary')
do yr=1,y
read(14) (atemp(yr))
enddo
do t=1,mt
read(13) ((hgt(i,j,t),i=1,iix),j=1,iiy)
enddo
do yr=1,y
do mo=1,m
kk=(yr-1)*12+mo
do j=1,iiy
do i=1,iix
hgt1(i,j,mo,yr)=hgt(i,j,kk)
enddo
enddo
enddo
enddo
do yr=1,y
if(atemp(yr)>1.0)then
anum=anum+1
endif
enddo
do yr=1,y
if(atemp(yr)<1.0.and.atemp(yr)>-1.0)then
bnum=bnum+1
endif
enddo
do j=1,iiy
do i=1,iix
num(i,j)=0
do yr=1,y
do mo=6,8
if(hgt1(i,j,mo,yr)<=-9.9e+08)then
num(i,j)=num(i,j)+1
endif
enddo
enddo
enddo
enddo
do j=1,iiy
do i=1,iix
zhgt(i,j)=0
do mo=6,8
do yr=1,y
if(atemp(yr)>1.0.and.hgt1(i,j,mo,yr)>-9.9e+08)then
zhgt(i,j)=zhgt(i,j)+hgt1(i,j,mo,yr)
endif
enddo
enddo
enddo
enddo
do j=1,iiy
do i=1,iix
phgt(i,j)=zhgt(i,j)/(3*anum-num(i,j))
enddo
enddo
do j=1,iiy
do i=1,iix
zzhgt(i,j)=0
do mo=6,8
do yr=1,y
if(hgt1(i,j,mo,yr)>-9.9e+08)then
zzhgt(i,j)=hgt1(i,j,mo,yr)+zzhgt(i,j)
endif
enddo
enddo
enddo
enddo

do j=1,iiy
do i=1,iix
zzhgt(i,j)=zzhgt(i,j)/(3*y-num(i,j))
enddo
enddo
do j=1,iiy
do i=1,iix
jphgt(i,j)=phgt(i,j)-zzhgt(i,j)
enddo
enddo

write(15,*) ((jphgt(i,j),i=1,iix),j=1,iiy)
write(16) ((jphgt(i,j),i=1,iix),j=1,iiy)
end
令再附上CTL 和 GS文件
dset F:\1\ghgt.grd
undef -9.9e+08
xdef  145 linear   0.000  2.50   
ydef 73 linear   -90.000  2.50
zdef   1 levels   500
tdef 1 linear   JAN1951      1yr
vars   1
ghgt  0 99
endvars

'reinit'
'open f:\1\ghgt.ctl'
'enable print f:\1\ghgt.gmf'
'set xflip on'
'set lon 180 300'
'set lat -20 70'
'set lev 500'
'set gxout contour'
'set clevs 0 10 15 25'
'd ghgt'
'print'
'disable print'

陌小默的

你是怎么个不对法哦？

勇敢的悲伤

陌小默的 发表于 2012-4-7 14:23

                               
登录/注册后可看大图

你是怎么个不对法哦？

画出来的图和文献给的图不一样，差别很大。。。

陌小默的

勇敢的悲伤 发表于 2012-4-7 14:47

                               
登录/注册后可看大图

画出来的图和文献给的图不一样，差别很大。。。

同样的方法和数据？

勇敢的悲伤

陌小默的 发表于 2012-4-7 17:05

                               
登录/注册后可看大图

同样的方法和数据？

数据是一样的，就是比文献给的多了一年。方法文献没说，直接给的图。就是夏季温度偏高年合成的500hpa高度距平场。距平场不都应该是一样的算法吗？

cui111yy

谢谢{:5_275:}{:5_275:}

守望之月/;-D

能不能把文献的图和你画的图都放上来？

纯粹猜测，
1）夏季都是JJA吗？
2）距平场合成的话，都是用的1981-2010的气候态平均？
3）比文献多的那一年，这一年对整个合成的距平场有多大影响？