任意两点间斜剖面做法总结【利用GRADS自带函数】

伯牙琴弦

感谢楼主！

酒国

很给力，谢谢楼主

酒国

谢谢楼主分享，学习一下

贪吃大王怕长胖

river 发表于 2014-10-16 07:15
那你就是只了解了grads的一半的功能啊，grads可以处理各种数据并画图，只要这些数据符合grads要求的格 ...

首先谢谢版主的耐心解答，小女子受益匪浅。可能是我表达的不清楚，我的意思是比如：grads可以根据grib这种格点数据画出等压面吗？请版主答疑解惑

river

贪吃大王怕长胖 发表于 2014-10-20 21:39
首先谢谢版主的耐心解答，小女子受益匪浅。可能是我表达的不清楚，我的意思是比如：grads可以根据grib这 ...

你为什么不试一下呢，我只能说是可以的

leig90

收藏一下。 经典

不骗你哦

lz能帮我看看么 。。我用的fnl1*1的资料，想做假相当位温
gs如下
'reinit'


'open d:\ffnnll\1.ctl'
'open d:\ffnnll\orog.ctl'
'set grads off'
'set zlog on'
'set csmooth on'
'set x 1'
'set y 1'
'set lev 1000 100'

'set t 31'
#'define xxx=hdivg(UGRDprs,VGRDprs)*100000'


lon1=102.0
lon2=106.0
lat1=32.0
lat2=28.0
lon=lon1
'collect 1 free'
'collect 2 free'
while(lon<=lon2)
  lat=lat1+(lat2-lat1)*(lon-lon1)/(lon2-lon1)


#计算假相当位温的计算
'define t0=tmpprs'
'define rh=rhprs'
'define prs=lev'
'define es=(6.1078*exp(17.2693882*(t0-273.16)/(t0-35.86)))'
'define qq=rh*(0.62197*es/(prs-0.378*es))/100.'
'define e=prs*qq/(0.62197+qq)+1e-10'
'define tlcl=55.0+2840.0/(3.5*log(t0)-log(e)-4.805)'
'define theta=t0*pow((1000/prs),(0.2854*(1.0-0.28*qq)))'
'define eqt=theta*exp(((3376./tlcl)-2.54)*qq*(1.0+0.81*qq))'


  'collect 1 gr2stn(eqt,'lon','lat')'
  'collect 2 gr2stn(orog.2(t=1),'lon','lat')'
  lon=lon+1
endwhile
'set map 15 1 9'
*'set parea 1 10.2 1.5 8'
'set xlopts 1 4 0.18'
'set ylopts 1 4 0.18'
'set font 0'

'set x 14 16'
'set xaxis 'lon1' 'lon2
'set clab on'

'set gxout shaded'
'set csmooth on'
'set clab off'
'set clevs 0'
'set ccols 0 1 0'
'd coll2gr(2,-u)'

'set gxout contour'
*'set cint 5'
'set clab forced'
'set cthick 11'
'd coll2gr(1,-u)'


'print'
'draw title omega along (55N,95W) to (15N,90W)'
'printim d:\ffnnll\section.gmf'
做出来的图问题有点大啊 就是要做运算的话这个改如何编gs

走路睡觉的猫WD

很厉害！赞一个！

river

不骗你哦 发表于 2014-10-21 20:59
lz能帮我看看么 。。我用的fnl1*1的资料，想做假相当位温
gs如下
'reinit'

先做运算，这个还用想啊。你要插值到一个斜剖面上，首先正常的点上得有假相当位温的值

不骗你哦

river 发表于 2014-10-21 21:21
先做运算，这个还用想啊。你要插值到一个斜剖面上，首先正常的点上得有假相当位温的值

嗯 我已经解决了 先算好在坐 为什么同样的方法咋涡度和散度上就行不通了呢？