任意两点的温度平流剖面gs

547168768

'reinit'
'open E:\hua\wrfout_d01.ctl'
'enable E:\hua\Tflow11.gmf'
'set mpdset  cnworld'
'set grads off'
'set grid off'
i=1
while(i<=13)
'set t 'i''
'set lon 110 118'
'set lat 35  43'
'set lev 1000 100'
'define dtx=cdiff(tk,x)'
'define dty=cdiff(tk,y)'
'define dx=cdiff(lon,x)*3.1416/180'
'define dy=cdiff(lat,y)*3.1416/180'
'define dit=-1*((U*dtx)/(cos(lat*3.1416/180)*dx)+V*dty/dy)/6.37e6'
'set zlog on'
'set x 1'
'set y 1'
'set z 1  19'
lon1 = 110.0
lon2 = 118.0
lat1 = 35.5
lat2 = 40.2
lon = lon1
'collect 1 free'
*'collect 2 free'
while (lon <= lon2)
lat = lat1 + (lat2-lat1)*(lon-lon1) / (lon2-lon1)
'collect 1 gr2stn(dit,'lon','lat')'
lon = lon + 0.054
endwhile
'set lon 110 118'
'set xlabs 35N,110E|36N,111E|37N,112E|38N,113E|39N,114E|40N,115E|41N,116E|42N,117E|43N,118E'
'set ylevs 1000 900 800 700 600 500 400 300 200 100'
'set gxout contour'
'set cterp on'
'set ylpos  0 l'
'set black -0.01 0.01'
'd coll2gr(1,-u)*1e4'
'q time'
    time=subwrd(result,3)
    'q dim'
   
    'draw title  Tflow in 'time' along (35.5N,110E) to (40.2N,118E)'
'draw string 9.87 0.46 unit:1e-4 K/s'
'print'
'c'
i=i+1
endwhile
'disable print'
;

547168768

欢迎大家提出宝贵意见！

/xin村儿/

贴个图..................

547168768

/xin村儿/ 发表于 2015-7-31 17:40
贴个图..................

不会贴图啊，请指导

/xin村儿/

你的代码左下方有个【编辑】，点击后，选择【图片】就会弹出图片上传对话框。

river

上个图就更完美了

hxyj

感谢楼主分享。

547168768

/xin村儿/ 发表于 2015-7-31 20:46
你的代码左下方有个【编辑】，点击后，选择【图片】就会弹出图片上传对话框。

图片地址是啥？我在自己电脑D盘存着

jack_ln

谢谢分享很实用

lifei910

贴个图，看下效果