| 
 
	积分1078贡献 精华在线时间 小时注册时间2017-4-12最后登录1970-1-1 
 | 
 
| 
本帖最后由 zlh 于 2019-8-15 21:14 编辑
x
登录后查看更多精彩内容~您需要 登录 才可以下载或查看,没有帐号?立即注册 
  
 参考论坛里的帖子画的温度平流的垂直剖面图,画出来值咋这么小呢,是正常的吗?虽然是冬季,但看其他文献冬季的也没有这么低呀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'
 
 
 
 | 
 
  |