爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 26150|回复: 8

[其他] 请问高度混合坐标怎样垂直插值到气压坐标下?

[复制链接]

新浪微博达人勋

发表于 2021-2-3 22:03:14 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
如题,想要把垂直方向是高度混合坐标的数据插值到气压坐标下,nc文件中关于高度的信息见图1使用vinth2p这个函数,系数a是每层高度lev,系数b文件中有单独变量,ps对应的是地表高度orog,这样计算出来画图看是不对的,程序部分见图2

求助一下论坛的朋友们这种该用哪个函数、怎么做呢?
非常感谢大家!!!!

图1 nc文件有关垂直坐标的信息

图1 nc文件有关垂直坐标的信息

图2 使用vinth2p函数计算的程序部分

图2 使用vinth2p函数计算的程序部分
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2023-10-2 22:24:10 | 显示全部楼层
杰西卡·芬 发表于 2021-3-27 17:16
具体的我也不是很懂原理是什么,你可以看看这个函数的解释

我看到函数解释要求参数hyam,hybm,ps的单位是帕或者百帕,请问高度坐标(单位m)也可以直接用吗?
密码修改失败请联系微信:mofangbao
回复 支持 0 反对 1

使用道具 举报

新浪微博达人勋

发表于 2021-3-19 10:10:48 | 显示全部楼层
请问楼主解决了吗?可以讨论一下?我插值出来以后1000hPa好多地方都是缺测值
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2021-3-20 09:49:39 | 显示全部楼层
杰西卡·芬 发表于 2021-3-19 10:10
请问楼主解决了吗?可以讨论一下?我插值出来以后1000hPa好多地方都是缺测值

我最后用了stdatmus_z2tdp这个函数将高度换算成了气压值,换出来看起来差别不大,也不知道对不对
然后用int2p_n_Wrap这个函数,插值到需要的气压高度层上
1000hPa很多缺测可能是这些地方本身地势高,1000hPa插在地下去了?我的猜测,我也不是很明白,可以再去问一下
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2021-3-22 09:16:04 | 显示全部楼层
我插值出来了,用Vnew = vinth2p(v,hyam,hybm,pnew,ps,2,P0mb,1,True)这个函数
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2021-3-25 14:18:30 | 显示全部楼层
杰西卡·芬 发表于 2021-3-22 09:16
我插值出来了,用Vnew = vinth2p(v,hyam,hybm,pnew,ps,2,P0mb,1,True)这个函数

请问您是用这个函数将高度混合坐标插值到气压坐标下吗,可以具体看一下您的这一段代码吗,我一直没尝试成功这个函数怎样对高度混合坐标进行插值
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2021-3-27 17:16:16 | 显示全部楼层
某某蒙 发表于 2021-3-25 14:18
请问您是用这个函数将高度混合坐标插值到气压坐标下吗,可以具体看一下您的这一段代码吗,我一直没尝试成 ...

具体的我也不是很懂原理是什么,你可以看看这个函数的解释
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2021-4-14 19:25:09 | 显示全部楼层
杰西卡·芬 发表于 2021-3-27 17:16
具体的我也不是很懂原理是什么,你可以看看这个函数的解释

好吧,谢谢您!
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2025-1-16 20:40:11 | 显示全部楼层
请问您解决了吗?我用这个插值位势高度500hpa,出来后竟然不能画5880线,最大值只有5622
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表