- 积分
- 335
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2016-2-27
- 最后登录
- 1970-1-1
|
发表于 2021-10-13 22:03:33
|
显示全部楼层
回帖奖励 +5 金钱
得到skewt_func.ncl去改吧
在这个目录下$NCARG_ROOT/lib/ncarg/nclscripts/csm/skewt_func.ncl
大概157行左右的位置
; --- declare pressure values [hPa] and x coordinates of the endpoints of each
; --- isobar. these x,y values are computed from the equations in the
; --- transform functions listed at the beginning of this program.
; --- refer to a skew-t diagram for reference if necessary.
pres = (/1050.,1000.,850.,700.,500.,400.,300.,250.,200.,150.,100./)
xpl = (/-19.0,-19.0,-19.0,-19.0,-19.0,-19.0,-19.0,-19.0 \
,-19.0,-19.0,-19.0/)
xpr = (/27.1,27.1,27.1,27.1,22.83,18.6,18.6,18.6,18.6,18.6,18.6/)
npres= dimsizes (pres) ; pressures
if (dimsizes(xpl).ne.dimsizes(xpr) .or. dimsizes(xpl).ne.npres) then
print ("Error: dimsizes pres, xpl, xpr do not match")
exit
end if
这附近看看呢(虽然回答过时了,但可以给后面需要的人点方法)
没有实践,刚好最近用到t-logp来回答一波 |
|