- 积分
- 167
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2012-3-16
- 最后登录
- 1970-1-1
|
楼主 |
发表于 2012-10-27 11:19:37
|
显示全部楼层
'reinit'
'open e:\data\fnl\721.ctl'
'set grads off'
'set grid off'
'set mpdset hires'
'set mpdset cnworld cnriver'
'set map 15 1 2'
'set xlopts 1 4 0.2'
'set ylopts 1 4 0.2'
'set parea 1.5 10 1 7.5'
'set csmooth on'
nt=13
while(nt<=16)
'set t 'nt''
'set lev 850'
'set ylint 5'
'set xlint 5'
'set lat 24.999 45.001'
'set lon 99.999 130.001'
'define es=(6.112*exp((17.67*(TMPprs-273.15))/(TMPprs-29.65)))'
'define qs=(0.62197*es/(lev-0.378*es))'
'define q=qs*rhprs/100'
'define sqdiv=hdivg(q*ugrdprs/9.8,q*vgrdprs/9.8)*1e8'
'set gxout shaded'
*'run define_color.gs'
*'set ccols 9 14 4 11 5 13 3 10 12 8 2 27 68'
*'set clevs -1.6 -1.2 -0.8 -0.4 0 0.4 0.8 1.2 1.6'
'd sqdiv'
'cbarn 0.85 0 5.75 0.35'
这是画水汽通量散度的gs文件,请高手查看有错误的地方么?尤其是计算部分,后面的一些收尾程序没贴上来 |
|