- 积分
- 1025
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2017-4-12
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
本帖最后由 zlh 于 2019-8-15 21:14 编辑
参考论坛里的帖子画的温度平流的垂直剖面图,画出来值咋这么小呢,是正常的吗?虽然是冬季,但看其他文献冬季的也没有这么低呀gs脚本:
'reinit'
'sdfopen E:\dalunwen\20110225\background\air.2011.nc'
'sdfopen E:\dalunwen\20110225\background\uwnd.2011.nc'
'sdfopen E:\dalunwen\20110225\background\vwnd.2011.nc'
'set grads off'
'set grid off'
'set xlpos 0.0 b'
'set ylpos 0.0 l'
'set csmooth on'
'set xlopts 1 4 0.15'
'set ylopts 1 4 0.15'
'set font 0'
'set t 220'
'set lev 1000 100'
'set lat 10 50'
'set lon -115 -65'
'set xlint 5'
'set ylint 5'
'define pi=3.1415726'
'define r=6.356e6'
'define vor=hcurl(uwnd.2,vwnd.3)'
'define dTx=cdiff(air,x)'
'define dTy=cdiff(air,y)'
'define dx=cdiff(lon,x)*cos(lat*pi/180)*pi*r/180'
'define dy=cdiff(lat,y)*pi*r/180'
'define advT=-(uwnd.2*dtx/dx+vwnd.3*dty/dy)'
'set gxout contour'
'set cint 1'
'set cthick 6'
'set clopts 1 0.1 0.15'
'set lat 33'
'set lon -110 -70'
'set xlpos -10'
'set ylpos -10'
'd advT*10000'
'run axis.gs -type b -position o -lable on -start -110 -end -70 -interval 5 -asuffix on -suffix `3. -size 0.1 -lsize 0.16 lthick 0.8'
'run axis.gs -type l -position o -lable on -start 1000 -end 100 -interval 100 -asuffix on -size 0.1 -lsize 0.16 lthick 0.8 -hoffset -0.1'
|
-
|