- 积分
- 88
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 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
|
|