grads处理fnl资料绘制500hPa位势高度+850hPa风场和相对湿度

云树

'reinit'

'!grib2ctl d:/fnl/fnl/20000624/fnl_20000624_06_00.grib1 > d:/fnl/fnl/20000624/fnl_20000624_06_00.ctl'

'!gribmap -v -i d:/fnl/fnl/20000624/fnl_20000624_06_00.ctl > d:/fnl/b.txt'

'open d:/fnl/fnl/20000624/fnl_20000624_06_00.ctl'

'set lon 70 140'

'set lat 25 60'

'set lev 850'

'set gxout shaded'

'd RHclm'

'cbar'

'set gxout barb'

'd ugrd3658m;vgrd3658m'

'set gxout contour'

'set ccolor 0'

'set strmden 10'

'set lev 500'

'd hgtprs'

'q w2xy 114 37'

x1=subwrd(result,3)

y1=subwrd(result,6)

'set string 1 c 10 0'

'set strsiz 0.2'

'draw string 'x1' 'y1' s'

'printim d:/fnl/fnl/20000624/2000062406.png white'

这个gs是可以成功画出图的，但是改时次后画下一张偶尔程序会报错（再运行一次又没错了，想要知道为什么）

云树

应该发到脚本编辑区的，但是点了脚本编辑后老是显示为未选择发帖主题

278803532

，话说，你先把所有的时次的fnl文件用一个ctl描述，先整出来，再画图不就行了哇，估计是你描述文件的问题。。。。

云树

278803532 发表于 2020-5-6 08:55
，话说，你先把所有的时次的fnl文件用一个ctl描述，先整出来，再画图不就行了哇，估计是你描述文 ...

因为日期都不连续，不知道怎么才能批量处理fnl，所以应该怎么操作呢

278803532

多下几天啊，，文件很多的么