| 
 
	积分86贡献 精华在线时间 小时注册时间2014-8-26最后登录1970-1-1 
 | 
 
| 
x
登录后查看更多精彩内容~您需要 登录 才可以下载或查看,没有帐号?立即注册 
    grads里先将站点数据转为二进制数据,之后运行时出现了这个错误,不知道是什么原因,求指导啊
 这是ctl
 dset F:\VIS\VIS-PIC\VIS-1980-1989\vis-1980-1989.grd
 dtype station
 stnmap F:\VIS\VIS-PIC\VIS-1980-1989\AH.map
 undef -9999
 title rain data sample
 tdef 1 linear 01may1980 1mo
 vars 1
 day 0 99 rain data
 endvars
 
 这是gs
 'reinit'
 'open F:\VIS\VIS-PIC\VIS-1980-1989\AH.ctl'
 'open F:\VIS\VIS-PIC\VIS-1980-1989\first.ctl'
 'set mpdset ahs'
 'set lon 114.5 123';'set lat 27 35.5'
 'define a=oacres(mask.1,day.2,50)'
 'define a1=maskout(a,mask.1-1.0)'
 'define aa=smth9(a1)'
 'set grads off'
 'set grid off'
 'set xlopts 1 6 0.18';'set ylopts 1 6 0.18'
 'set xlint 2';'set ylint 2'
 
 
 'set rgb 20 255 244 197'
 'set rgb 21 255 227 101'
 'set rgb 22 255 198 42'
 'set rgb 23 255 167 42'
 'set rgb 24 255 144 42'
 'set rgb 25 255 100 42'
 'set rgb 26 255 74 0'
 'set rgb 27 227 32 23'
 'set rgb 28 194 23 15'
 
 'set rbcols 20 21 22 23 24 25 26 27 28 '
 
 'set clevs 15 20 25 30 35 40 45 50 '
 'set gxout shaded'
 'd aa'
 'cbar 1.0 1'
 'printim F:\VIS\VIS-PIC\VIS-1980-1989\AH-1980-1989.png white'
 
 这是转换数据的fortran
 Character*8 stid
 Real lat,lon,day,tim
 integer nlev,flag
 !变量声明结束
 !程序开始
 tim=0.0
 nlev=1
 flag=1
 PRINT*,111
 Open(1,file='F:\VIS\VIS-PIC\VIS-1980-1989\vis-1980-1989.txt',status='old')
 PRINT*,111
 open(2,file='F:\VIS\VIS-PIC\VIS-1980-1989\vis-1980-1989.grd',status='replace',form='binary')
 10 Read(1,*,end=100)stid,lat,lon,day
 print*,stid,lat,lon,day
 Write(2)stid,lat,lon,tim,nlev,flag,day
 !write(2)stid,lat,lon,day
 Goto 10
 100 continue
 Close(1)
 nlev=0
 Write(2)stid,lat,lon,tim,nlev,flag
 close(2)
 !程序结束
 End
 
 
 | 
 |