grads如何画地面等压线图？

小河淌水

各位大神，想请教grads如何画地面等压线图？

独奏古筝到天明

FNLgrib2资料：
'reinit'
'open e:\nc\fnl_20160415_00_00.ctl'
*'set mpvals 60 180 20 90'
'set lon 80 160'
'set lat 10 60'
'set mpdset hires'
*'set mproj nps'
'set map 1 1 6'
i=1
while(i<27)
'set grid on 5 1'
'set grads off'
'set t 'i
'set lev 1000'
'set gxout barb'
'set cmin 6'
'set ccolor 5'

'd skip(ugrdprs(lev=1000),2,2);skip(vgrdprs(lev=1000),2,2);mag(ugrdprs(lev=1000),vgrdprs(lev=1000))'

'set gxout contour'
'set cthick 6'
'set cint 5'
'set clskip 2'
'set ccolor 1'
'q time'
t=subwrd(result,3)
'set clab forced'
'd PRMSLmsl/100'


'draw title Pressure_1000hPa_WIND`b(barb)`n_'i'_'t''
'printim E:\nc\Sea_Surface_Pressure_1000hPa_TEMP_'i'_'t'.png white x1000 y800'
i=i+1
'c'
endwhile

独奏古筝到天明

nc资料：
'reinit'
'sdfopen e:\grads\2014\mslp.2014.nc'
'sdfopen e:\grads\2014\uwnd.2014.nc'
'sdfopen e:\grads\2014\vwnd.2014.nc'
'sdfopen e:\grads\2014\air.2014.nc'

'set mpvals 60 180 20 90'
'set mpdset hires'
'set mproj nps'
'set map 1 1 6'
i=1300
while(i<1360)
'set grid on 5 1'
'set grads off'
'set t 'i
'set gxout barb'
'set lev 1000'
'set cthick 6'
'set cstyle 2'
'set cint 5'
'set cmin 5'
'set ccolor 0'
'd skip(uwnd.2(lev=1000),2,2);skip(vwnd.3(lev=1000),2,2);mag(uwnd.2(lev=1000),vwnd.3(lev=1000))'
'set gxout contour'
'set cthick 6'
'set cint 4'
'set clskip 2'
'q time'
t=subwrd(result,3)
'd mslp.1/100'
'set gxout contour'
'set cthick 6'
'set cstyle 2'
'set ccolor 2'
'set cint 4'
*'d air.4(lev=1000)-273.15'

'draw title MSLP_1000hPa_WIND`b(barb)`n_'i'_'t''
'printim E:\grads\surface\SFC-different\Sea_Surface_Pressure_1000hPa_TEMP_'i'_'t'.png white x1000 y800'
i=i+1
'c'
endwhile

小河淌水

感谢大神啊

好多鱼

独奏古筝到天明 发表于 2016-9-5 19:43
FNLgrib2资料：
'reinit'
'open e:\nc\fnl_20160415_00_00.ctl'

为什么在运行循环的时候出现问题了呢？

好多鱼

楼主，请问一下，你用FNL资料运行的时候有问题吗？我的出现了点错误。

好多鱼

独奏古筝到天明 发表于 2016-9-5 19:43
nc资料：
'reinit'
'sdfopen e:\grads\2014\mslp.2014.nc'

unable to locate endwhile statement for the while statement at line 12

嘟嘟的哒哒

大神抱抱大腿