计算eta层对应垂直高度的方法

中尺度风速模拟

首先感谢各路大神的帮忙，这里根据自己的理解说两个方法，给初学者一些建议，也希望大家批评指正。
一、用ncl很简单：
      用这个命令： Z=wrf_user_getvar(a,"z",0)
                            print (Z)
        得出来的Z就是高度
二、用matlab
       1、通过matlab提取PH、PHB、HGT三个值
       2、根据公式gmp=(PH+PHB)/9.81-HGT计算，找出gmp为60时所对应的PH的值  
       3、找出这个PH值所对应的垂直层，这一层的风速值即为60米高度的风速

       4、（这一段出自天目大神的帖子）先按照eta值运行下去（在跑wrf时，可以先只算一个小时的数据，因为现在还属于试验阶段），得到模式结果后，你就可以画出每个eta值对应的height的等值线图了，这样你就可以看出每个eta对应的离地高度值。这样，得到的结果，可能是eta=0.993对应的位势高度是57m,  eta=0.9922对应的是66m,这样一来，你就知道eta值对应的60m的范围了（0.993-0.9922）。第一次试验结束，现在你就把0.993改为0.9928，大概差不多，下次你的结果就基本是60m了，当然你也可能第二次试验。由于eta坐标是等气压值坐标，所以每个eta面上的位势高度都不可能相等。也就是说等气压面并不是一个平面，所以无论你怎么调节都不可能全部等于60m。 所以在60m左右就可以了，一般会在59-61m。相当不大。

绿衣茶爽

有点乱呢，用NCL得到的z是什么高度？位势高度吗？

中尺度风速模拟

绿衣茶爽 发表于 2016-1-5 11:17
有点乱呢，用NCL得到的z是什么高度？位势高度吗？

嗯嗯 对 是位势高度

piaopiao1992

楼主可以帮我看一下吗http://bbs.06climate.com/forum.php?mod=viewthread&tid=48917

Sunshining

中尺度风速模拟 发表于 2016-1-5 19:47
嗯嗯 对 是位势高度

这个z单位是m吧，不应该是位势高度吧？今天也碰到了问题，纯粹讨论没有杠的意思

杜不苦

嗯嗯，学习了，回顾一下

2868559739

Sunshining 发表于 2020-4-3 10:45
这个z单位是m吧，不应该是位势高度吧？今天也碰到了问题，纯粹讨论没有杠的意思

我也觉得不是位势高度，但我感觉是上层与下层的中心高度