fnl处理资料出错

linyu临宇

'reinit'
'open d:/shuju/201407/fnl_20140712_06_00.ctl'
'set gxout fwrite'
'set fwrite d:/shuju/201407wind.dat'
'set lat 10 60'
'set lon 90 160'
tt=1
while (tt<=40)
'set t 'tt
zz=1
while(zz<=26)
'set z'zz
'd UGRDprs'
zz=zz+1
endwhile
zz=1
while(zz<=26)
'set z'zz
'd VGRDprs'
zz=zz+1
end while
tt=tt+1
endwhile
zz=zz+1
endwhile
'disable fwrite'
'reinit'
;
这是我对应兰溪给的“fnl资料认识”写的gs文件，但是运行的时候总是出错，这是为什么呢？

传说中的谁

该空格的地方空格

river

不熟悉简写就不要用简写，这句'set z'zz都是错的
写全了'set z 'zz''

linyu临宇

river 发表于 2015-3-24 20:42
不熟悉简写就不要用简写，这句'set z'zz都是错的
写全了'set z 'zz''

'reinit'
'open d:/shuju/201407/fnl_20140712_06_00.ctl'
'set gxout fwrite'
'set fwrite d:/shuju/201407wind.dat'
'set lat 10 60'
'set lon 90 160'
tt=1
while (tt<=40)
'set t 'tt''
zz=1
while(zz<=26)
'set z 'zz''
'd UGRDprs'
zz=zz+1
endwhile
zz=1
while(zz<=26)
'set z 'zz''
'd VGRDprs'
zz=zz+1
endwhile
tt=tt+1
endwhile
'disable fwrite'
'reinit'
;
我这样改也不对，什么哪里还有错？

river

linyu临宇 发表于 2015-3-24 21:28
'reinit'
'open d:/shuju/201407/fnl_20140712_06_00.ctl'
'set gxout fwrite'

错误提示呢，看脚本好像没什么问题了

linyu临宇

river 发表于 2015-3-24 20:42
不熟悉简写就不要用简写，这句'set z'zz都是错的
写全了'set z 'zz''

zdef 26 levels
1000 975 950 925 900 850 800 750 700 650 600 550 500 450 400 350 300 250 200 150 100 70 50 30 20 10
vars 121
no4LFTXsfc  0 132,1,0  ** surface Best (4-layer) lifted index [K]
no5WAVAprs  0 230,100,500  ** 5-wave geopot. height anomaly [gpm]
no5WAVHprs  0 222,100,500  ** 5-wave geopotential height [gpm]
ABSVprs 26 41,100,0 ** Absolute vorticity
ACPCPsfc  0 63,1,0  ** surface Convective precipitation [kg/m^2]
ALBDOsfc  0 84,1,0  ** surface Albedo [%]
APCPsfc  0 61,1,0  ** surface Total precipitation [kg/m^2]
CAPEsfc  0 157,1,0  ** surface Convective Avail. Pot. Energy [J/kg]
CAPE180_0mb  0 157,116,46080 ** 180-0 mb above gnd Convective Avail. Pot.
这是fnl的ctl文件，我想画变量为APCPsfc这个变量，编写gs文件如下：
'reinit'
'open d:/shuju/201407/fnl_20140712_06_00.ctl'
'set lat -30 90'
'set lon 30 180'
'set lev 1000'
'enable print d:/shuju/tu/tp.gmf'
'set gxout shaded'
'set mpdset hires'
'set mpdset cnworld'
'set grid off'
'set grads on'
tt=1
while(tt<=40)
'set t 'tt
'd APCPsfc'
'q time'
x1=subwrd(result,3)
'draw title 'x1'tp'
'print'
'c'
tt=tt+1
endwhile
'disable print'
;
运行之后出来的图是d:/1.png

linyu临宇

river 发表于 2015-3-24 20:42
不熟悉简写就不要用简写，这句'set z'zz都是错的
写全了'set z 'zz''

zdef 26 levels
1000 975 950 925 900 850 800 750 700 650 600 550 500 450 400 350 300 250 200 150 100 70 50 30 20 10
vars 121
no4LFTXsfc  0 132,1,0  ** surface Best (4-layer) lifted index [K]
no5WAVAprs  0 230,100,500  ** 5-wave geopot. height anomaly [gpm]
no5WAVHprs  0 222,100,500  ** 5-wave geopotential height [gpm]
ABSVprs 26 41,100,0 ** Absolute vorticity
ACPCPsfc  0 63,1,0  ** surface Convective precipitation [kg/m^2]
ALBDOsfc  0 84,1,0  ** surface Albedo [%]
APCPsfc  0 61,1,0  ** surface Total precipitation [kg/m^2]
CAPEsfc  0 157,1,0  ** surface Convective Avail. Pot. Energy [J/kg]
CAPE180_0mb  0 157,116,46080 ** 180-0 mb above gnd Convective Avail. Pot.
这是fnl的ctl文件，我想画变量为APCPsfc这个变量，编写gs文件如下：
'reinit'
'open d:/shuju/201407/fnl_20140712_06_00.ctl'
'set lat -30 90'
'set lon 30 180'
'set lev 1000'
'enable print d:/shuju/tu/tp.gmf'
'set gxout shaded'
'set mpdset hires'
'set mpdset cnworld'
'set grid off'
'set grads on'
tt=1
while(tt<=40)
'set t 'tt
'd APCPsfc'
'q time'
x1=subwrd(result,3)
'draw title 'x1'tp'
'print'
'c'
tt=tt+1
endwhile
'disable print'
;
运行之后出来的图是d:/1.pngd:/1.png

river

没看见图啊

风之精灵

你把经纬度设置（set lon  lon1 lon2 ;set lat lat1 lat2）放到循环的最里面，再跑一遍你的gs。grads存储和读取数据是按照t>lev>vars>lat>lon的顺序，可能是这个原因造成的。当然，该空格的时候也要留出空格，以免系统识别出错。

linyu临宇

风之精灵 发表于 2015-3-30 08:57
你把经纬度设置（set lon  lon1 lon2 ;set lat lat1 lat2）放到循环的最里面，再跑一遍你的gs。grads存储和 ...

我现在想用这个fnl数据绘制台风路径，但是我不知道怎样用grads画，我的想法是先用fortran读取fnl数据，找到最低气压的点，，但是怎么用fortran读取fnl数据呢？是先用fwrite提取对应的气压么转为grd格式的么？还是？？？求解