如何从wrfout数据中提取海拔差异较大地区的云量

cici破晓碎寒

大家好！

我想获取地形复杂区wrfout数据的云量结果(青藏高原，四川盆地一带)，使用NCL中的“wrf_user_getvar(cfrac)”函数可以获得低、中、高云。但是这个函数对低、中、高云的判断是通过高度判断的，即低于2km是低云等等，但是这个高度放在青藏高原一带，并没有加上地形海拔高度，得出的结果肯定是不准确的。与ERA5再分析结果相比，低、中、高三种云的空间分布不合理，这里我只展示低云的分布(图1)。我认为地形海拔高度应该加到对“高度”判断中，但是NCL这个函数的代码中似乎很难加(图2)。然后我又尝试用“气压”的标准来判断云量，但是和ERA5结果比较，这个分布也不对不上(图3)。请问大家如何从wrfout数据中获取地形起伏比较大的区域的云量?

感谢帮助！

Candymaker

Cloud fraction
This diagnostic is a special one that allows variable to either just be "cfrac", in which a series of default values are assumed, or it can be an array of six ordered string values that represent the following:

variable(0) - diagnostic name - "cfrac" (no other values possible)
variable(1) - vertical type - "height_agl", "height_msl", "pressure", "pres"
Default: "height_agl"
variable(2) - low threshold - the low cloud threshold value
Default: "300." for height_agl
variable(3) - mid threshold - The middle cloud threshold value
Default: "2000." for height_agl
variable(4) - high threshold - the high cloud threshold value
Default: "6000." for height_agl
variable(5) - missing value - the missing/fill value to use for non-cloud areas
Default: "9.9692099683868690E36"
这里的height_agl,就是离地高度吧uu，还需要加上地形吗

cici破晓碎寒

您好，我看了下，好像确实是离地高度，agl是离海平面高度减去地形高度，那这个结果不大对的上，可能就是模式的模拟和ERA5还是有差异的。谢谢！

yyyxxx666

您好，我想问一下wrf_user_getvar怎么把高度转成气压的啊？

Jelly不会ncl

请问楼主解决了吗？

cici破晓碎寒

yyyxxx666 发表于 2023-8-11 12:11
您好，我想问一下wrf_user_getvar怎么把高度转成气压的啊？

图片所示的这个函数有可以用气压表示的选项，默认的pressure分类是下面这个代码所示，和ERA5的中低高云气压分类标准差不多，这个代码是在NCL的"$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"里找到的，但是需要NCL版本在6.6.0以上才可以用这个函数        dim_vars = dimsizes(varin)
        do idims = 1,dim_vars-1
          if ( idims .eq. 1 ) then
            if ( varin(idims) .eq. "pressure" .or. varin(idims) .eq. "pres") then
              vert_type = varin(idims)
              low_thresh = 97000.
              mid_thresh = 80000.
              high_thresh = 45000.
              vert_inc_w_height = 0
            else
              vert_type = varin(idims)
            end if
          end if
          if ( idims .eq. 2 ) then
            low_thresh = stringtodouble(varin(2))
          end if
          if (idims .eq. 3) then
            mid_thresh = stringtodouble(varin(3))
          end if
          if (idims .eq. 4) then
            high_thresh = stringtodouble(varin(4))
          end if
你用6.6.0以上的NCL，然后直接写下面这个代码就可以用气压分类的方法区分低中高云了

cfracvar=(/"cfrac","pressure"/)
f1=addfile("/data/data_s12_1/zqfan/xscai/wrfout_four_groups/con/wrfout_d01_2020-06-25_00:00:00.nc", "r")
cfrtemp = wrf_user_getvar(f1, cfracvar, 0)

cici破晓碎寒

Jelly不会ncl 发表于 2023-8-22 00:19
请问楼主解决了吗？

没有完全解决，后来觉得可能ERA5给的低中高云的分布可能也并不准确，不好用来比较，只能大概看看，要是真的做评估，还是要看卫星资料的云，然后处理模式数据，可能要用其他更精确的区分低中高云的方法

cici破晓碎寒

Candymaker 发表于 2023-5-30 11:30
Cloud fraction
This diagnostic is a special one that allows variable to either just be "cfrac", in  ...

对，我看了下代码，这个应该确实是离地高度了

yyyxxx666

cici破晓碎寒 发表于 2023-8-22 16:46
图片所示的这个函数有可以用气压表示的选项，默认的pressure分类是下面这个代码所示，和ERA5的中低高云气 ...

谢谢楼主，马上试试

Jelly不会ncl

cici破晓碎寒 发表于 2023-8-22 16:49
没有完全解决，后来觉得可能ERA5给的低中高云的分布可能也并不准确，不好用来比较，只能大概看看，要是真 ...

我之前有看过用UPP处理云量，但是我还没尝试过