【完美】解决FNL资料任意两点间斜剖面叠加地形的问题

songww

请问我要同时添加地形和变量填色，但是最后出来地形只有一条线，是什么原因呢？ctl如下：
'reinit'
'sdfopen e:\2020\ec2.nc'
'open e:\2020\xpm\orog0.25.ctl'
'enable print  e:\2020\xpm\vor+div-terrain.gmf'

'set dfile 2'
'set t 1'
'set x 1'
'set y 1'
'set z 1 17'

lon1=102.0
lon2=106.0
lat1=28.0
lat2=33.0

lon=lon1
'collect 3 free'
while(lon<=lon2)
lat=lat1+(lat2-lat1)*(lon-lon1)/(lon2-lon1)
'collect 3 gr2stn(orog,'lon','lat')'
lon=lon+1
endwhile


'set dfile 1'
'set grads off'
'set zlog on'
'set csmooth on'
'set x 1'
'set y 1'
'set lev 1000 100'
'set t 7'
lon1 =102
lon2 =106
lat1 = 28.0
lat2 = 33.0
lon = lon1
'collect 1 free'
'collect 2 free'
while (lon <= lon2)
  lat = lat1 + (lat2-lat1)*(lon-lon1) / (lon2-lon1)
  'collect 1 gr2stn(d*100000,'lon','lat')'
  'collect 2 gr2stn(vo*100000,'lon','lat')'
  lon = lon + 1
endwhile



'set grads off'
'set zlog on'
'set grid off'
'set lon 102 106'
*'set z 1 27'
*'set ylab off'
*'set xlab off'
'set csmooth on'


*'set parea 1.5 9.0 1.5 6.5'
*'set xlopts 1 6 0.14'
*'set ylopts 1 6 0.14'
*'set font 0'

'set x 14 16'
*'set xaxis 'lon1' 'lon2
'set xaxis 'lat1' 'lat2
'set clab on'

'set gxout shaded'
'set xlint 1'
*'set ccols 0 0 7 0'
*'set ccolor 1'
'd coll2gr(1,-u)'
'set gxout contour'
*'set cint 0.5'
'set cthick 6'
'set ccolor 1'
'd coll2gr(2,-u)'

'set lon 102 106'
'set z 1 27'
'set gxout shade2'
'set clevs 0'
'set ccols -1 1'
'd coll2gr(3,-u)'

*'cbarn 1 0'
'print'

'disable print'
'reinit'

songww

songww 发表于 2021-5-25 17:21
请问我要同时添加地形和变量填色，但是最后出来地形只有一条线，是什么原因呢？ctl如下：
'reinit'
'sdfo ...

已解决，grads版本问题，换2.1或者2.2就行

红河岸

lz的资料是ncep的吗

river

红河岸 发表于 2021-7-18 12:35
lz的资料是ncep的吗

用的是fnl 的grib2 格式的资料，你用NC的话道理是一样的。而且我帖子的第一句就说之前我发过用NCEP资料画任意剖面的帖子，你去看一下不就行了吗

YonghangX

想请教一下楼主，运行后出现下面的反馈，怎么解释。

Cannot use log scaling when coordinates <= 0
Linear scaling used

river

YonghangX 发表于 2021-7-27 14:09
想请教一下楼主，运行后出现下面的反馈，怎么解释。

Cannot use log scaling when coordinates

信息不全没办法回答。我只能说提示的是坐标有小于零的值，所以对数坐标轴报错了 set zlog on

红河岸

river 发表于 2021-7-18 15:55
用的是fnl 的grib2 格式的资料，你用NC的话道理是一样的。而且我帖子的第一句就说之前我发过用NCEP资料画 ...

不好意思我可能说的不够清楚，我的意思是您用的ctl里面配的是ncep单个时次的数据，还是您把所有的时次都写到一个文件里？

river

红河岸 发表于 2021-8-7 18:24
不好意思我可能说的不够清楚，我的意思是您用的ctl里面配的是ncep单个时次的数据，还是您把所有的时次都 ...

这个不影响吧，都可以。如果要用到很多时次的资料我会直接使用批量描述，不会把资料写成一个文件，那样太麻烦了。

逝者如斯

太牛了！！！

逝者如斯

，如果沿某一经度做剖面，是否也能使用该gs？