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

钱磊037

楼主，请教一下：为什么我用fnl资料画不出来剖面图啊，grads总是提示：Collect Error:Z or T must be the only
varying dimension ?
PS：我并未做时间的剖面啊。。。

钱磊037

附上gs文件，拜托各位大神帮忙看看a
'reinit'
'open d:\data\1000-200hPa0811-17q.ctl'
'enable print d:\data\35N113E-42N123E11-17q.gmf'
i=1
while(i<=28)
'set vpage 0 8 0 8'
'set parea 1 6 1 8'
'set grads off'
'set zlog on'
'set csmooth on'

'set x 1'
'set y 1'
'set lev 1000 200'
'set t 'i''

'set lon 80 150'
'set lat 0 70'
'define prs=lev'
'define tc=tmpprs-273.16'
'define es=(6.1078*exp((17.27*tc)/(tc+237.3)))'
'define e=rhprs/100*es'
'define q=(0.62197*e/(prs-e*0.378)*1000)'
'define qtu=ugrdprs*q/9.8'
'define qtv=vgrdprs*q/9.8'
'define div=(hdivg(qtu,qtv)*1e6)'
'define qt=mag(qtu,qtv)'

lon1=113.0
lon2=123.0
lat1=35.0
lat2=42.0
lon=lon1
'collect 1 free'
while(lon<=lon2)
  lat=lat1+(lat2-lat1)*(lon-lon1)/(lon2-lon1)
  'collect 1 gr2stn(qt,'lon','lat')'
  lon=lon+1
endwhile

*'set lat 35 42'
'set lon 113 123'
'set z 1 19'
'set gxout shaded'
'set clevs 5 10 15 20 25 30 35'
'set ccols 0 4 5 3 7 8 2 6'
'd coll2gr(1,-u)'
'q time'
xx=subwrd(result,3)
'draw title 35N113E-42N123E QT VCS Field at 'xx''
'print'
'c'
i=i+1
endwhile
'disable print'
;

river

钱磊037 发表于 2014-3-10 21:32

                               
登录/注册后可看大图

附上gs文件，拜托各位大神帮忙看看a
'reinit'
'open d:\data\1000-200hPa0811-17q.ctl'

计算的先后顺序有问题，要先计算涡度、散度之类的量，然后再做剖面

wtf2ftw

学习了，标记下

liusx2007

请问river 楼主：是否有grib 1（FNL 1X1 资料）下使用的地形文件orog.ct'，可否能提供使用？

river

liusx2007 发表于 2014-3-15 10:19

                               
登录/注册后可看大图

请问river 楼主：是否有grib 1（FNL 1X1 资料）下使用的地形文件orog.ct'，可否能提供使用？

没有，自己把2.5*2.5的插值一下不就可以了吗

ml11178

river大大您好: 我使用wrf的out做成grads可以讀的資料之後

試用第一個程序畫了看看，但是一直出現:

Warning:  X axis labels overridden by SET XAXIS.
   Labels may not reflect correct scaling for dimensions or data.
Warning:  X axis labels overridden by SET XAXIS.
   Labels may not reflect correct scaling for dimensions or data.


實在再不懂為什麼... 資料我有內插到氣壓層，ctl檔也有寫對,我有先display別的變數也是沒問題的，

但是套用第一個程序就是一直出現以上錯誤，懇請大大幫我解答一下  ORZ

river

ml11178 发表于 2014-3-16 20:06

                               
登录/注册后可看大图

river大大您好: 我使用wrf的out做成grads可以讀的資料之後

試用第一個程序畫了看看，但是一直出現:

其实之前就解释过了的，你仔细看这个警告里面提到的 SET XAXIS 命令，这个是强制改变坐标轴标值的命令，所以肯定会有警告啊。就像fortran里面也时常有一些警告，但是你也可以忽略它继续运行这个程序。其实这个是一个道理，重要的是你强制改变的东西真的能和实际的对上就可以了

ml11178

恩~我畫出來的圖，只有外框出來，並沒有contour在框裡面，不知道是哪裡出了問題，訊息也只有給那幾行

river

ml11178 发表于 2014-3-16 21:30

                               
登录/注册后可看大图

恩~我畫出來的圖，只有外框出來，並沒有contour在框裡面，不知道是哪裡出了問題，訊息也只有給那幾行

如果只有那几行警告的话，那就说明其他的都运行通过了，怎么会没有contour呢。你还是把你用的资料的ctl 和你的gs都贴上来看一下吧