- 积分
- 631
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2013-4-26
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
楼主有个小的问题,就是通过wrf的结果提出了n天的逐小时的站点的边界层高高度的值,单位应该是m
放在了pblv的数组(一个n*24小时个数的一位数组)里面,现在我想画一个图,横坐标是时间,纵坐标是压力
这需要垂直方向上差值
lev = (/ 3000.,4000.,5000.,7000.,8500.,10000.,20000.,25000.,30000.,40000.,42500.,\
50000.,60000.,70000.,75000.,77500.,80000.,82500.,85000.,87000.,90000.,92500.,95000.,97500.,99000.,100000./)
用了pbls=int2p(lev,pblv,lev,linlog)
但是这句话是有问题的,因为它需要lev的个数(26)和pblv的个数(N*24个)是一致的?
所以是应该怎写这句话让它差值到压力方向上。
还有有没有大神画过要素随时间的变化图,纵坐标是压力的,温压湿风的我可以花,但是边界层高度的具体咋画啊?毕竟这个是个一维的,我希望在水平风场上叠加一条曲线,就是边界层高度的。
我之前写的代码是叫h_time_point.ncl
这里面是filname是已经通过wrf的结果差值以后的文件。
|
|