求助 关于绘制某测站温度平流随时间变化的垂直剖面图问题

qiuba

        本来要绘制某固定测站附近温度平流随时间变化的垂直剖面图，剖面图已绘制出，但总感觉不对劲，似乎每个地方均适用，这还有意义吗？请各位高人指点。
        具体脚本如下：
'reinit'
'sdfopen d:\grads\data\uwnd.2020.nc'
'sdfopen d:\grads\data\vwnd.2020.nc'
'sdfopen d:\grads\data\air.2020.nc'
*--------------GrADS设置-----------------------
'set grads off'
'set grid off'
'set csmooth on'
'set map 15 1 9'
'set parea 1 10.2 1.5 8'
'set xlopts 1 6 0.14'
'set ylopts 1 6 0.14'

*--------计算----------
'set lev 1000 300'
'set lon 70 150'
'set lat 10 60'
'set t  51 71'
*--------计算温度平流----------
'define pi=3.1415726'
'define r=6.356e6'
'define dtx = cdiff(air.3,x)'
'define dty = cdiff(air.3,y)'
'define dx=cdiff(lon,x)*cos(lat*pi/180)*pi/180*r'
'define dy=cdiff(lat,y)*pi/180*r'
'define advT=-(uwnd.1*dtx/dx+vwnd.2*dty/dy)'
'define n=aave(advT,lon=100,lon=130,lat=30,lat=50)'
*---------输出图形------------
'set lon 100'
'set lat 10'
*固定经纬度任意值均可以，图形似乎没变化，这合理吗？？？
'set zlog on'
*'set xlint 6'
'set ylevs 1000 925 850 700 600 500 400 300'
'set xlabs 20|02|08|14|20|02|08|14|20|02|08|14|20|02|08|14|20|02|08|14|20'
'set gxout shaded'
'set cint 1'
'set cthick 3'
'set clskip 1'
'd n*1e5'      
*量级为10的-5,单位是10-5K/s
'run cbarn.gs'
'set ylpos 0 r'
'set gxout contour'
'set cint 1'
'set cthick 3'
*'set clskip 1'
'd n*1e5'      
'set font 10'
'writeslhz 0.5 4 气压(hPa) fsj 2 1 1.3 0.75 90'
'writeslhz 1 8.1 温度平流 FSJ 2.5 1 1.3 0.6 0'
'writeslhz 0.5 0.9 1月13日 FSJ 2 1 1.3 0.6 0'
'writeslhz 2.4 0.9 1月14日 FSJ 2 1 1.3 0.6 0'
'writeslhz 4.3 0.9 1月15日 FSJ 2 1 1.3 0.6 0'
'writeslhz 6.1 0.9 1月16日 fsj 2 1 1.3 0.6 0'
'writeslhz 8.0 0.9 1月17日 fsj 2 1 1.3 0.6 0'
'writeslhz 9.8 0.9 1月18日 fsj 2 1 1.3 0.6 0'
'printim d:\grads/1/advT-poumian.png white'
;

lightmoon

'define n=aave(advT,lon=100,lon=130,lat=30,lat=50)'
‘’d n*1e5'   
就相当于你把这个区域平均值都赋值给你所圈定纬度里边的所有格点，那当然改经纬度图像不变啦。

qiuba

lightmoon 发表于 2020-3-16 12:42
'define n=aave(advT,lon=100,lon=130,lat=30,lat=50)'
‘’d n*1e5'   
就相当于你把这个区域平均值都赋 ...

        厉害！一针见血！去掉区域平均就ok了！非常感谢！向您致敬!        这是修改后的效果图！完全可以勾绘出测站上空温度平流的垂直变化情况，也具有预报意义！感谢您的指点，向您点赞！

杨豆

厉害了学到了谢谢

布拉风666

感谢

杨玲玲1014

你好，楼主，请问横坐标的时间怎么标出来？非常感谢

huyaoxing

你好，图画的很好看，线条均匀，你是用的哪个版本的grads呢？

来自苹果APP客户端