[ncl-talk] Vertical Integration of cloud water mixing ratio from WRF out file

[Anie K Lal]

Hi Barry

Thank you for the response. I had checked 'pw'. But it will not serve my 
need. pw has only three dimensions (time, lat, lon).

Could you kindly share the calculations for pw so that I can select the 
vertical levels also.

Anie.




-------- Original Message --------
Subject: Re: [ncl-talk] Vertical Integration of cloud water mixing ratio 
from WRF out file
Date: 05.12.2022 20:27
 From: Barry Lynn <barry.h.lynn at gmail.com>
To: Anie K Lal <Anieklal at cas.iitd.ac.in>

Hi:

 From what you wrote, that would be precipitable water.

You can zero out the levels you don't want to keep, and then calculate
pw.

https://www.ncl.ucar.edu/Document/Functions/WRF_arw/wrf_user_getvar.shtml

The difficulty is zeroing it out and then writing it back to the wrf
output.

The easiest thing to do would be to request the function for pw and
then modify it the way you want, and then call it.

However, I have calculated precipitable water before and maybe able to
find it for you.

I will try to take a look.

Barry

On Mon, Dec 5, 2022 at 3:46 PM Barry Lynn <barry.h.lynn at gmail.com>
wrote:

> 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 [1]"
>> 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 [1]"
>> 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 [2]
> Jerusalem, Israel
> Local: 02 930 9525
> Cell: 054 7 231 170
> Int-IS: x972 2 930 9525

-- 

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 [2]
Jerusalem, Israel
Local: 02 930 9525
Cell: 054 7 231 170
Int-IS: x972 2 930 9525

CAUTION: This email originated from outside of IIT Delhi. Do not click
links or open attachments unless you recognize the sender and know the
content is safe.


Links:
------
[1] http://wrfout_d02.nc
[2] http://weather-it-is.com