NCL计算整层可降水量

342901wzx

;-----------------气柱可降水量------------------------
      qtemp    = q(time|:,latitude|:,longitude|:,level|:)  ;Specific humidity,units = "kg kg**-1"
      plev          = lev*100
      plev@units    = "Pa"
      ptop          = min(plev)   
      dp1           = dpres_plevel(plev, psfc, ptop, 0)  ; dp(9) 计算等压层厚度 ;psfc:地面气压 psfc@units    = "Pa"
      dp1!0    = "tau"
      dp1!1    = "lev"
      dp1!2    = "lat"
      dp1!3    = "lon"
      dp       =  dp1(tau|:,lat|:,lon|:,lev|:)  ;调整数据结构（与比湿数据结构一致）
      delete(dp1)
      
      ;printVarSummary(dp)
      ;printVarSummary(psfc)
      ;printVarSummary(qtemp)
      copy_VarCoords(qtemp, dp)
      prcwater      = prcwater_dp(qtemp, dp)
      prcwater!0    = "tau"
      prcwater!1    = "lat"
      prcwater!2    = "lon"
      prcwater&tau  = tau
      prcwater&lat  = lat
      prcwater&lon  = lon
      prcwater@long_name  = "prcwater"
      prcwater@units      = "mm"
      
      ;printVarSummary(dp)
      ;printVarSummary(psfc)
      ;printVarSummary(qtemp)
      ;printMinMax(prcwater(2,:,:), 0)
      delete(qtemp)
            
      ;-----------------------------------------------

woshiwzm

能够把计算公式列出来就更好了

三岁最可爱

谢谢楼主分享，原来ncl有内置的函数，省得自己写了

mier

https://www.ncl.ucar.edu/Document/Functions/Built-in/prcwater_dp.shtml#:~:text=prcwater_dp%20Computes%20total%20column%20precipitable,water%20of%20the%20rightmost%20dimension.
NCL对应的可降水量pw的函数说明