- 积分
- 57088
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2015-1-21
- 最后登录
- 1970-1-1
![未绑定新浪微博用户 新浪微博达人勋](source/plugin/sina_login/img/gray.png)
|
发表于 2016-7-19 08:00:18
|
显示全部楼层
'reinit'
'open e:\nc\fnl_20160415_00_00.ctl'
i=1
while(i<=27)
'set t 'i
'set lev 850'
'set lon 90 140'
'set lat 10 50'
'set grads off'
'set mpdset hires'
'set map 1 1 6'
'define_colors'
'set gxout shaded'
'set clevs -32 -28 -24 -20 -16 -12 -8 -4 0 4 8 12 16 20 24 28 32 36'
'set ccols 49 48 47 46 45 44 43 42 41 0 21 22 23 24 25 26 27 28 29'
'set lon 90 140'
'set lat 10 50'
'set t 'i
'define vor = hcurl(ugrdprs,vgrdprs)'
'd 1e5*vor'
'cbarn '
'set gxout contour'
'set cthick 4'
'set ccolor 8'
'd 1e5*vor'
'q time'
t=subwrd(result,3)
'draw title 850hPa_Vort_'i'_'t''
'printim E:\nc\850hPa_Vort_'i'_'t'.png white x1000 y800'
i=i+1
'c'
endwhile
这是我的代码,应该设置两遍经纬度,才能画散度和涡度,我之前遇到过这个问题。 |
|