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

硕硕

攒金币。。。。

黄叁3

river 发表于 2017-12-20 07:48
建议你看看原始资料是不是正确的

好的谢谢回复，我的资料应该是没问题的，我在找找原因

叶绾棠

谢谢楼主分享，好好学习一下

爱喝柠檬汁

请问下面这几个语句是啥意思呢？
'set x 1'
'set y 1'

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

爱喝柠檬汁

还有个问题
'd coll2gr(1,-u)'
里的-u是啥呀？
计算的不是omega??

爱喝柠檬汁

下面是我的gs和图，不知道对不对，请@river 帮我看看呗
我用的是欧洲中心分辨率0.5*0.5一天四次的再分析资料
'reinit'
'sdfopen E:/era.nc'
'set grads off'
'set zlog on'
'set csmooth on'
'set t 109'
'set x 1'
'set y 1'
'set lev 1000 100'
lon1 = 110.0
lon2 = 130.0
lat1 = 60.0
lat2 = 40.0
lon = lon1
'collect 1 free'
while (lon <= lon2)
  lat = lat1 + (lat2-lat1)*(lon-lon1) / (lon2-lon1)
  'collect 1 gr2stn(w,'lon','lat')'
  lon = lon + 1
endwhile

'set map 15 1 9'
'set font 0'

'set x 14 16'
'set xaxis 'lon1' 'lon2
'set clab on'
'set gxout shaded'
'set cmax 0'
'd coll2gr(1,-w)'
'set gxout contour'
*'set cint .5'
'set cthick 11'
'd coll2gr(1,-w)'
'draw title omega along (60N,110W) to (40N,130W)'
'printim E:/wpou.png white'
;

爱喝柠檬汁

有几个问题，不懂
'set x 1'
'set y 1'有啥意义，为什么要这么设置，
**********
'set x 14 16'又是啥意思
************
原来的gs里
'd coll2gr(1,-u)'
这里为啥是u呢，算的不是垂直速度吗，我给改成w了，不知道对不对
@river

爱喝柠檬汁

运行上面那个gs是有警告的，如图
翻译了一下：
x轴xaxis标签覆盖集
标签可能不反映正确的扩展维度或数据
******
啥意思呢，咋改呢呢啊啊啊啊啊！！！

爱喝柠檬汁

@river 理理我呀！

爱喝柠檬汁

{:eb303:}{:eb303:}