- 积分
- 55609
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2015-1-21
- 最后登录
- 1970-1-1
|
发表于 2016-12-25 15:46:51
|
显示全部楼层
'reinit'
'open e:\nc\fnl_20160420_00_00.ctl'
'set lev 850'
*'set mpvals 60 180 20 80'
'set mpdset hires'
*'set mproj nps'
'set lon 100 150'
'set lat 10 50'
'set map 1 1 6'
i=1
while(i<9)
'set grid on 5 1'
'set grads off'
'define_colors'
'set t 'i
'set clevs -36 -32 -28 -24 -20 -16 -12 -8 -4 0 4 8 12 16 20 24 28 32'
'set ccols 49 48 47 46 45 44 43 42 41 0 21 22 23 24 25 26 27 28 29'
'set gxout shaded'
'set cthick 6'
'set cstyle 2'
'set cint 4'
'define t=TMPprs-273.15'
'd t'
'cbarn'
'set gxout contour'
'set ccolor 2'
'set cstyle 3'
'set cthick 6'
'set cint 4'
*'set clab forced'
'd t'
'set gxout barb'
'set ccolor 1'
'd skip(UGRDprs,2);skip(VGRDprs,2)'
'set gxout contour'
'set cthick 6'
'set cint 4'
'set ccolor 1'
*'set clab forced'
'd HGTprs/10'
'q time'
t=subwrd(result,3)
'draw title 850hPa_HGT_TEMP_WIND`b(barb)`n_'i'_'t''
'printim E:\nc\FNL_850hPa_HGT_WIND_TEMP_'i'_'t'.png white x1000 y800'
i=i+1
'c'
endwhile
'reinit'
|
|