ncl计算假相当位温的函数pot_temp_equiv_tlcl

小懒C

• iounits(2)=0 input w are kg/kg
• iounits(2)=1 input w are g/kg
  

在函数的源代码中，iounits(2)=1表示w是kg/kg，0则是g/kg
;===== 以下为该部分代码 ==========                                                                                        ; Bolton (2000); Eq 16
  xcon = 1000.                          ; water vapor pressure
  if (iounits(2).eq.1) then
      rmix = mixr*xcon                  ; kg/kg ==> g/kg
  else
      rmix = mixr
  end if                                ;evp@units = "hPa"
  evp  = (p*pcon)*rmix/(622+rmix)

;======================================================
附上看代码的方式，ncl官网一些函数在写说明时会告诉各位调用了什么库，这个库便是函数所在位置，例如
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"  
function pot_temp_equiv_tlcl (
        p  : numeric,  
        t  : numeric,  
        tlcl : numeric,
        w  : numeric,  
        iounits [4] : integer           )     
   return_val [dimsizes(t)] :  float or double

pot_temp_equiv_tlcl位于$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl中

白衣风月倾天下

dei 官网写错了 实际应用的时候根据结果也应该是0表示g/kg

眼里有星辰

确实很坑，文档说明和下面的例子都不一致

versatile

楼主，这个源代码不是说iounits(2)=1的话，w是g/kg吗，因为语句的意思应该是把原先的kg/kg乘以1000换算成了g/kg，

kongyueyiyun

感谢楼主，我算的结果不对劲，排查了半天都不知道哪里有问题。变量的单位已经检查了无数次，没看到这个帖子的话不知道还要蒙蔽多久。。。。。

小懒C

versatile 发表于 2021-9-30 10:22
楼主，这个源代码不是说iounits(2)=1的话，w是g/kg吗，因为语句的意思应该是把原先的kg/kg乘以1000换算成了 ...

因为他0-1指示的是函数里虚参（mixr）的单位，最终计算ncl是规定必须转化为g/kg，所以你看到的是要kg/kg===>g/kg

我超胖

用了这个公式出现一堆nan是什么问题呢

                               
登录/注册后可看大图

                               
登录/注册后可看大图

                               
登录/注册后可看大图

来自安卓APP客户端

小懒C

我超胖 发表于 2022-6-1 12:03
用了这个公式出现一堆nan是什么问题呢

我用不会，检查一下数据

白衣风月倾天下

我超胖 发表于 2022-6-1 12:03
用了这个公式出现一堆nan是什么问题呢

我偶尔也会出现这个问题 如果使用不含抬升凝结高度的函数就不会 我检查过原数据 但是这个问题没法解决 建议找出nan的位置 然后对部分有问题的数据进行处理

shayingying

我就说算出来怎么严重低估了！！！我的天，我改一下再试试