帮忙看看这个CTL有没有错，怎么画不出来图啊？？

勇敢的悲伤

给出了全国160站降水资料，我把6、7、8月份的专为GRD格式。下面的是程序：
program main
implicit none
integer,parameter:: m=160,n=58,y=60
integer i,j,yr,a
integer::anum=0,bnum=0
integer rain(m,n,1),rain1(m,n,2),rain2(m,n,3),rain3(m,n),rjuping(m),rain4(m),zrain(m),prain(m)
real atemp(n)
open(21,file='f:\1\rain6.grd',form='binary')
open(31,file='f:\1\rain7.grd',form='binary')
open(23,file='f:\1\rain8.grd',form='binary')
read(21) ((rain(i,j,1),i=1,m),j=1,n)
read(31) ((rain1(i,j,2),i=1,m),j=1,n)
read(23) ((rain2(i,j,3),i=1,m),j=1,n)
open(11,file='f:\1\r1606.txt')
write(11,*)((rain(i,j,1),i=1,m),j=1,n)
open(12,file='f:\1\r1607.txt')
write(12,*)((rain1(i,j,2),i=1,m),j=1,n)
open(13,file='f:\1\r1608.txt')
write(13,*)((rain2(i,j,3),i=1,m),j=1,n)
end
然后我给rain6（就是六月份的降水）配了一个CTL文件
dset F:\1\rain6.grd
undef -9.99E+08  
title NCEP/NCAR REANALYSIS PROJECT
xdef  63 linear   73.000  1
ydef  41 linear   15.000  1
zdef    1 levels  1000
tdef 1 linear   JAN1951      1yr
vars 1
p  0 99 u
endvars
这么写对吗？可以画出来六月份各站点的降水图吗？

lotina

是站点资料  应该在CTL里边加一个dtype    station

勇敢的悲伤

lotina 发表于 2012-4-5 22:05

                               
登录/注册后可看大图

是站点资料  应该在CTL里边加一个dtype    station

可是加了之后，更不能用 了，根本就打不开GRD了

mofangbao

看看站点资料绘图详解的帖子吧：
http://bbs.06climate.com/forum.php?mod=viewthread&tid=4903
得按照GrADS中的格式来写，不能想当然的写了

mofangbao

看看站点资料绘图详解的帖子吧：
http://bbs.06climate.com/forum.php?mod=viewthread&tid=4903
得按照GrADS中的格式来写，不能想当然的写了

勇敢的悲伤

mofangbao 发表于 2012-4-5 22:33

                               
登录/注册后可看大图

看看站点资料绘图详解的帖子吧：
http://bbs.06climate.com/forum.php?mod=viewthread&tid=4903
得按照Gr ...

好的。谢啦~~