<div dir="ltr">Hi:<div><br></div><div>From what you wrote, that would be precipitable water.</div><div><br></div><div>You can zero out the levels you don't want to keep, and then calculate pw.</div><div><br></div><div><a href="https://www.ncl.ucar.edu/Document/Functions/WRF_arw/wrf_user_getvar.shtml">https://www.ncl.ucar.edu/Document/Functions/WRF_arw/wrf_user_getvar.shtml</a><br></div><div><br></div><div>The difficulty is zeroing it out and then writing it back to the wrf output.  </div><div><br></div><div>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.</div><div><br></div><div>However, I have calculated precipitable water before and maybe able to find it for you.</div><div><br></div><div>I will try to take a look.</div><div><br></div><div>Barry</div><div><br></div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Mon, Dec 5, 2022 at 3:46 PM Barry Lynn <<a href="mailto:barry.h.lynn@gmail.com">barry.h.lynn@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">Hello Anie:<div><br></div><div>Could you please describe, in a bit more detail, what you are trying to do.</div><div><br></div><div>Barry</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Mon, Dec 5, 2022 at 3:12 PM Anie K Lal via ncl-talk <<a href="mailto:ncl-talk@mailman.ucar.edu" target="_blank">ncl-talk@mailman.ucar.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">I was able to find another method for vertical integration using <br>
"vibeta."<br>
<br>
------------------------------------------------------------------------------------------<br>
   filename = <br>
"/media/hp/Expansion/WRF/2019_simulation/EXP02/<a href="http://wrfout_d02.nc" rel="noreferrer" target="_blank">wrfout_d02.nc</a>"<br>
   f1        = addfile(filename,"r")<br>
<br>
   q1 = f1->QCLOUD(74,0:8,:,:)*1000<br>
   psfc  =  wrf_user_getvar(f1,"PSFC",74)*0.01                    ; <br>
surface pressure, turn to hPa (; unit of PSFC in WRF is Pa)<br>
   plev  =  wrf_user_getvar(f1,"pressure",74)                     ; air <br>
pressure  hPa (time,level,lat,lon)<br>
   plev1 =  plev(0:8,:,:)<br>
<br>
   printVarSummary(plev1)<br>
<br>
   copy_VarCoords(plev1, q1)<br>
<br>
   printVarSummary(q1)<br>
<br>
   ptop = 913.9985<br>
   pbot = 999.4952<br>
<br>
   q2      =   <br>
vibeta(plev1(south_north|:,west_east|:,bottom_top|:),q1(south_north|:,west_east|:,bottom_top|:),2,psfc,pbot,ptop)<br>
<br>
   printVarSummary(q2)<br>
------------------------------------------------------------------------------------------------<br>
<br>
I am not sure whether this is the right way to find the vertically <br>
integrated value of mixing ratio. Is this method correct?<br>
<br>
Kindly advise.<br>
<br>
Thank you<br>
Anie<br>
<br>
<br>
<br>
<br>
<br>
-------- Original Message --------<br>
Subject: Vertical Integration of cloud water mixing ratio from WRF out <br>
file<br>
Date: 05.12.2022 15:22<br>
 From: Anie K Lal <<a href="mailto:Anieklal@cas.iitd.ac.in" target="_blank">Anieklal@cas.iitd.ac.in</a>><br>
To: <a href="mailto:ncl-talk@mailman.ucar.edu" target="_blank">ncl-talk@mailman.ucar.edu</a><br>
<br>
Hi all,<br>
<br>
I want to calculate the liquid water path by vertically integrating<br>
cloud water mixing ratio for the first 9 levels from surface.<br>
<br>
Cloud water mixing ratio= QCLOUD(Time, bottom_top, south_north, <br>
west_east)<br>
<br>
This is the script I tried:<br>
------------------------------------------------------------------------------<br>
   filename = "/WRF/2019/EXP02/<a href="http://wrfout_d02.nc" rel="noreferrer" target="_blank">wrfout_d02.nc</a>"<br>
   f1        = addfile(filename,"r")<br>
<br>
   q1 = f1->QCLOUD(74,0:8,:,:)<br>
   P  = f1->P(74,0:8,:,:)<br>
   T1  = f1->T(74,0:8,:,:)<br>
   Pb = f1->PB(74,0:8,:,:)<br>
   P  = (P+Pb)   ;Total pressure in Pa<br>
<br>
   Ptop  =min(P)<br>
   Psurf =max(P)<br>
<br>
   dp   = dpres_plevel_Wrap(P, Psurf, Ptop, 0)  ; dp(time,plev,lat,lon)<br>
<br>
   printVarSummary(dp)<br>
--------------------------------------------------------------------------------<br>
<br>
I am getting the error:<br>
"fatal:Number of dimensions in parameter (0) of (dpres_plevel_Wrap) is<br>
(3), (1) dimensions were expected"<br>
<br>
I understand that 'P' should have only one dimension, but how to solve<br>
this in my case where P is three dimensional.<br>
<br>
Looking forward to any kind of help.<br>
<br>
Thank you.<br>
_______________________________________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@mailman.ucar.edu" target="_blank">ncl-talk@mailman.ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="https://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">https://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr"><div dir="ltr"><div><div dir="ltr"><br></div><div dir="ltr">Barry H. Lynn, Ph.D<div><div>Senior <font color="#000000">Scientist, Lecturer,</font></div><div><div><font color="#000000">The Institute of Earth Sciences, <br>The Hebrew University of Jerusalem, <br>Givat Ram, Jerusalem 91904, Israel <br></font></div><font color="#000000">Tel: 972 547 231 170<br>Fax: (972)-25662581</font></div></div><div><span style="color:rgb(136,136,136)"><br></span></div><div>Weather It Is, LTD<br>Weather and Climate Focus<br><a href="http://weather-it-is.com" target="_blank">https://weather-it-is.com</a><br>Jerusalem, Israel<br>Local: 02 930 9525<br>Cell: 054 7 231 170<br>Int-IS: x972 2 930 9525<br><br></div></div></div></div></div>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><br></div><div dir="ltr">Barry H. Lynn, Ph.D<div><div>Senior <font color="#000000">Scientist, Lecturer,</font></div><div><div><font color="#000000">The Institute of Earth Sciences, <br>The Hebrew University of Jerusalem, <br>Givat Ram, Jerusalem 91904, Israel <br></font></div><font color="#000000">Tel: 972 547 231 170<br>Fax: (972)-25662581</font></div></div><div><span style="color:rgb(136,136,136)"><br></span></div><div>Weather It Is, LTD<br>Weather and Climate Focus<br><a href="http://weather-it-is.com" target="_blank">https://weather-it-is.com</a><br>Jerusalem, Israel<br>Local: 02 930 9525<br>Cell: 054 7 231 170<br>Int-IS: x972 2 930 9525<br><br></div></div></div></div></div>