[ncl-talk] Vertical Integration of cloud water mixing ratio from WRF out file
Barry Lynn
barry.h.lynn at gmail.com
Mon Dec 5 06:46:14 MST 2022
Hello Anie:
Could you please describe, in a bit more detail, what you are trying to do.
Barry
On Mon, Dec 5, 2022 at 3:12 PM Anie K Lal via ncl-talk <
ncl-talk at mailman.ucar.edu> wrote:
> I was able to find another method for vertical integration using
> "vibeta."
>
>
> ------------------------------------------------------------------------------------------
> filename =
> "/media/hp/Expansion/WRF/2019_simulation/EXP02/wrfout_d02.nc"
> f1 = addfile(filename,"r")
>
> q1 = f1->QCLOUD(74,0:8,:,:)*1000
> psfc = wrf_user_getvar(f1,"PSFC",74)*0.01 ;
> surface pressure, turn to hPa (; unit of PSFC in WRF is Pa)
> plev = wrf_user_getvar(f1,"pressure",74) ; air
> pressure hPa (time,level,lat,lon)
> plev1 = plev(0:8,:,:)
>
> printVarSummary(plev1)
>
> copy_VarCoords(plev1, q1)
>
> printVarSummary(q1)
>
> ptop = 913.9985
> pbot = 999.4952
>
> q2 =
>
> vibeta(plev1(south_north|:,west_east|:,bottom_top|:),q1(south_north|:,west_east|:,bottom_top|:),2,psfc,pbot,ptop)
>
> printVarSummary(q2)
>
> ------------------------------------------------------------------------------------------------
>
> I am not sure whether this is the right way to find the vertically
> integrated value of mixing ratio. Is this method correct?
>
> Kindly advise.
>
> Thank you
> Anie
>
>
>
>
>
> -------- Original Message --------
> Subject: Vertical Integration of cloud water mixing ratio from WRF out
> file
> Date: 05.12.2022 15:22
> From: Anie K Lal <Anieklal at cas.iitd.ac.in>
> To: ncl-talk at mailman.ucar.edu
>
> Hi all,
>
> I want to calculate the liquid water path by vertically integrating
> cloud water mixing ratio for the first 9 levels from surface.
>
> Cloud water mixing ratio= QCLOUD(Time, bottom_top, south_north,
> west_east)
>
> This is the script I tried:
>
> ------------------------------------------------------------------------------
> filename = "/WRF/2019/EXP02/wrfout_d02.nc"
> f1 = addfile(filename,"r")
>
> q1 = f1->QCLOUD(74,0:8,:,:)
> P = f1->P(74,0:8,:,:)
> T1 = f1->T(74,0:8,:,:)
> Pb = f1->PB(74,0:8,:,:)
> P = (P+Pb) ;Total pressure in Pa
>
> Ptop =min(P)
> Psurf =max(P)
>
> dp = dpres_plevel_Wrap(P, Psurf, Ptop, 0) ; dp(time,plev,lat,lon)
>
> printVarSummary(dp)
>
> --------------------------------------------------------------------------------
>
> I am getting the error:
> "fatal:Number of dimensions in parameter (0) of (dpres_plevel_Wrap) is
> (3), (1) dimensions were expected"
>
> I understand that 'P' should have only one dimension, but how to solve
> this in my case where P is three dimensional.
>
> Looking forward to any kind of help.
>
> Thank you.
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at mailman.ucar.edu
> List instructions, subscriber options, unsubscribe:
> https://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
--
Barry H. Lynn, Ph.D
Senior Scientist, Lecturer,
The Institute of Earth Sciences,
The Hebrew University of Jerusalem,
Givat Ram, Jerusalem 91904, Israel
Tel: 972 547 231 170
Fax: (972)-25662581
Weather It Is, LTD
Weather and Climate Focus
https://weather-it-is.com <http://weather-it-is.com>
Jerusalem, Israel
Local: 02 930 9525
Cell: 054 7 231 170
Int-IS: x972 2 930 9525
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20221205/424963aa/attachment.htm>
More information about the ncl-talk
mailing list