任意两点间斜剖面做法:coll2gr 異常

蕎蕎安

各位大大好:

看了river大的帖子(任意两点间斜剖面做法总结【利用GRADS自带函数】)，
試著做了，但是Open GrADS(版本2.0.2),
屏幕顯示coll2gr欄位1是空的(empty),

看了官網coll2gr(http://cola.gmu.edu/grads/gadoc/gadocindex.html)的
Usage Notes:The option"-u" will only work for data in pressure vertical coordinates. If your vertical coordinate is height, num should be the actual number of vertical levels in the collection of profiles. Note: This restriction has been removed as of version 2.1.0.“-u”選項僅適用於壓力垂直坐標中的數據。 如果您的垂直坐標是height，num應該是配置文件集合中垂直級別的實際數量。 注意：此限制已從2.1.0版本刪除。

所以把gs檔內的coll2gr(1.-u)改成coll2gr(19.-u)，
還是不能跑
不知道是我哪裡弄錯了,
麻煩各位大大幫忙了。


以下是gs檔
'reinit'

'open g:\FNL/20160413.ctl'
'set grads off'
'set zlog on'
'set x 1'
'set y 1'
'set lev 1000 100'
lon1 = 1220
lon2 = 124.0
lat1 = 25.0
lat2 = 20.0
lon = lon1
'collect 1 free'
while (lon <= lon2)
  lat = lat1 + (lat2-lat1)*(lon-lon1) / (lon2-lon1)
  'collect 1 gr2stn(20160413,'lon','lat')'
  lon = lon + 1
endwhile
'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 font 0'
'set x 14 16'
'set xaxis 'lon1' 'lon2
'set clab on'
'set gxout shaded'
'set cmax 0'
*'set ccols 0 0 7 0'
*'set clevs 0 .5 15'
*'set ccols 0 0 7 0'
'd coll2gr(19,-u)'
'set gxout contour'
'set cint .5'
'set cthick 11'
'd coll2gr(19,-u)'
'draw title omega along (55N,95W) to (15N,90W)'
'printim g:/picture/20160413.jpg x1280 y768'
*'c'

river

你还是自己去官网好好看看coll2gr这个函数的用法吧······

蕎蕎安

好的

李小毛123

蕎蕎安 发表于 2018-1-5 09:14
好的

请问解决了嘛，我也是这个问题