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

yuyuxuan

您好我在做假相当位温剖面时无报错，但最后出图是这样子的，可以帮忙看看嘛，谢谢
'reinit'
'sdfopen D:/KY/ECMWF/201503.nc'
'set t 20'
'set grads off'
'set grid off'
'set zlog on'
'set csmooth on'
'set xlopts 1 6 0.14'
'set ylopts 1 6 0.14'
'set x 1'
'set y 1'
'set lev 1000 100'

'define tmp=t'
'define rh=r'
'define prs=lev'
'define es=(6.11*exp(17.27*(tmp-273.16)/(tmp-35.86)))'
'define q=rh*(0.62197*es/(prs-es))/100.'
'define e=prs*q/(0.62197+q)+1e-10'
'define tlcl=55.0+2840.0/(3.5*log(tmp)-log(e)-4.805)'
'define theta=tmp*pow((1000/prs),(0.2854*(1.0-0.28*q)))'
'define eqt=theta*exp(((3376./tlcl)-2.54)*q*(1.0+0.81*q))'

lon1=116
lon2=125
lat1=37
lat2=42
lon=lon1
'collect 1 free'
while(lon<=lon2)
lat=lat1+(lat2-lat1)*(lon-lon1)/(lon2-lon1)
'collect 1 gr2stn(eqt,'lon','lat')'
lon=lon+1
endwhile

'set lon 116 125'
*'set xaxis 116 125 2'
'set clab on'
'set gxout contour'
'd coll2gr(1,-u)'
'draw title fnl,'m_time
'gxprint D:/KY/tankong/20150307-11/T/jiaxiangdangweiwen/'t'.png white'
;

river

yuyuxuan 发表于 2020-4-19 17:05
您好我在做假相当位温剖面时无报错，但最后出图是这样子的，可以帮忙看看嘛，谢谢
'reinit'
'sdfopen D:/ ...

你一个大暴雨怎么会出现这种问题呢。你只计算了一个点的假相当位温，你还能指望它出来什么图？

yuyuxuan

river 发表于 2020-4-28 22:27
你一个大暴雨怎么会出现这种问题呢。你只计算了一个点的假相当位温，你还能指望它出来什么图？

握草   突然惊醒，谢了！

river

yuyuxuan 发表于 2020-5-1 08:27
握草   突然惊醒，谢了！

准备过五一，有些“轻敌”了吧，哈哈

ZHUTOU

学习中，学习楼主

李小毛123

'reinit'
'open F:\fnl\fnl072000.ctl'
'set grads off'
'set csmooth on'
'set zlog on'
'set x 1'
'set y 1'
'set lev 1000 100'
lon1 =113.8
lon2 =118.2
lat1 =36.5
lat2 = 41.8
lon = lon1
'collect 1 free'
while (lon <= lon2)
  lat = lat1 + (lat2-lat1)*(lon-lon1) / (lon2-lon1)
  'collect 1 gr2stn(UGRDprs,'lon','lat')'
  lon = lon + 1
endwhile

*set map 15 1 9'
'set xlopts 1 6 0.14'
'set ylopts 1 6 0.14'
'set font 0'

*'set xaxis 'lon1' 'lon2
'set xlabs 36.5N,113.8E|41.8N,118.2E'
'set clab on'
'set gxout shaded'
'set cmax 0'
'd coll2gr(1,-u)'
'set gxout contour'
'set cthick 11'
'd coll2gr(1,-u)'

'enable print F:\fnl\AB.gmf white'
'print white'
'disable print'
'c'
'reinit'

大神，这个是我的gs文件，运行之后报错，不知道什么原因，垂直坐标是气压

puck66

一直没搞明白遮地形的原理，不知道论坛哪里是否有介绍？另外就是，我用gtopo30_10min.nc的地面高程数据，是否能像orog的这样子用作地形遮挡？

puck66

下载用用，谢谢楼主分享

river

李小毛123 发表于 2020-8-8 20:41
'reinit'
'open F:\fnl\fnl072000.ctl'
'set grads off'

你这个问题应该是因为在出图之前少设置一个经度范围，你不设置他默认是最前面的set x 1 和 set y 1，只有一个点，所以函数就报错了。你需要在出图之前重新设置，比如 set lon 113.8 118.2，应该就可以了

river

puck66 发表于 2020-9-13 13:28
一直没搞明白遮地形的原理，不知道论坛哪里是否有介绍？另外就是，我用gtopo30_10min.nc的地面高程数据，是 ...

遮地形的原理很简单啊，就是有地形的数据和变量的数据，画图的时候用地形的图把变量的图盖住。至于你说的那个地形资料，我不了解它的格式，不清楚可不可以用