温度平流随时间的剖面图。

不骗你哦

做一个区域平均的温度平流随时间变化的高度剖面图，gs貌似有点毛病，帮看看把
这个温度平流的计算公式是我在这边帖子里看到的，不过貌似 经纬度定了求区域平均就不行了 我试过把经纬度改成区域计算就可以出图。帮帮解决下。。
'reinit'
'open d:/fff.ctl'
'set gxout contour'
'set lev 1000 500'
'set lon 180'
'set lat -90'
'set grid off'
'set grads off'
'set t 1 100'
'enable print d:\tpingliu.gmf'
'define pi=3.1415726/180.0'
'define vor=hcurl(ugrdprs,vgrdprs)'
'define r=6.356e6'
'define dtx=cdiff(tmpprs,x)'
'define dty=cdiff(tmpprs,y)'
'define dx=cdiff(lon,x)*cos(lat*pi/180)*pi*r/180'
'define dy=cdiff(lat,y)*pi*r'
'define advT=-(ugrdprs*dtx/dx+vgrdprs*dty/dy)'
'define n=aave(advT,lon=113,lon=120,lat=30,lat=40)'
'd n*1e3'
'disable print'
'reinit'

koala54

ls正解，求区域平均要设置成区域，而画时间-高度剖面要固定经纬度，建议把计算步骤和画图分开。
关于楼主的问题，吧里有帖子讲过，点我看看
我试了下，可以出图的。我觉得还有一个问题就是求温度平流时，pi定义时除以了180，后面就不用再除了。
下面是我的gs，给楼主参考一下，欢迎指正

'reinit'
'open d:\data\ht.ctl'
'enable print d:\wanwan\advT.gmf'
*--------------GrADS设置-----------------------
'set parea 1 10 1 8'
'set mpdset cnworld'
'set map 1 1 5'
'set xlopts 1 6 0.12'
'set ylopts 1 6 0.12'
'set grid off'
'set grads off'
'set csmooth on'
*-----------------计算----------------------------
'set lon 100 140'
'set lat 10 40'
'set lev 1000 200'
'set t 8 28'
*----------------温度平流的计算-------------------
'define pi=3.1415726'
'define r=6.356e6'
'define dTx=cdiff(TMPprs,x)'
'define dTy=cdiff(TMPprs,y)'
'define dx=cdiff(lon,x)*cos(lat*pi/180)*pi/180*r'
'define dy=cdiff(lat,y)*pi/180*r'
'define advT=-(ugrdprs*dTx/dx+vgrdprs*dTy/dy)'
'define n=aave(advT,lon=112,lon=126,lat=22,lat=35)'
*---------------画图-------------------------------
*画图前先把经纬度固定（好像任意值都行）
'set lon 110'
'set lat 20'
'set zlog on'
'd n*1e5'      
*量级为10的-5,单位是10-5K/s
'draw title avdT'
'print'
'disable print'
'close 1'
;

mofangbao

你把经纬度都定了还能叫区域平均么?

不骗你哦

koala54 发表于 2013-3-14 09:59

                               
登录/注册后可看大图

ls正解，求区域平均要设置成区域，而画时间-高度剖面要固定经纬度，建议把计算步骤和画图分开。
关于楼主的 ...

谢谢 就是像你所说 已经解决！

Mr.Gao

koala54 发表于 2013-3-14 09:59

                               
登录/注册后可看大图

ls正解，求区域平均要设置成区域，而画时间-高度剖面要固定经纬度，建议把计算步骤和画图分开。
关于楼主的 ...

请问为什么要除以dx和dy呢？温度平流不是速度乘以温度梯度吗

river

楼主就是进入了一个误区，咱们刚开始学习剖面图的时候基本不会有什么计算过程，所以直接固定经纬度什么的就行了。但是一旦涉及到计算过程，先固定经纬度计算过程就无法进行。所以要分开，先计算，再固定经纬度画剖面

1417920146

学习一下 很有帮助

koala54

Mr.Gao 发表于 2013-3-30 17:30

                               
登录/注册后可看大图

请问为什么要除以dx和dy呢？温度平流不是速度乘以温度梯度吗

你说的对，只要用温度的差分近似温度的梯度就行。我是这样理解的，首先这里的中央差函数cdiff(expr,dim)是针对格点计算的，具体就是用前一个格点的温度值减后一个给点的温度值，再除以两个值的格点间隔（也就是2）。而我们需要的中央差是要针对实际的格点距离而言，所以后面又定义了dx、dy，为所对应的2个格点纬向和经向的实际间距。于是dTx/dx=（温度差/格点距2）/（两格点间纬向实际距离/格点距2），这就是用中央差表示温度的纬向梯度。

Mr.Gao

koala54 发表于 2013-3-31 18:14

                               
登录/注册后可看大图

你说的对，只要用温度的差分近似温度的梯度就行。我是这样理解的，首先这里的中央差函数cdiff(expr,dim)是 ...

噢 ，这样，那是我对cdiff函数理解有错误，谢谢了

wc19900615

谢谢楼主分享~~~