<div dir="ltr"><div>There are no specific functions or examples to calculate vertical vapor flux divergence:</div><div> <br>
(1/g*[d(qw)/dz]) <br></div><div><br></div><div>from a WRF output.</div><div>=================<br></div><div><br></div><div>I do not remember much about WRF output. <br></div><div><br></div><div>Generically:<br></div><div><br></div><div>QW = Q*W<br></div><div><pre>dQWdz = <a href="https://www.ncl.ucar.edu/Document/Functions/Built-in/center_finite_diff_n.shtml"><b>center_finite_diff_n</b></a> (QW,'z',False,0,0)</pre></div><div>=======================<br></div><div><br></div><div>There are several examples for doing this type of activity on 'regular/ rectilinear' grids with pressure [also, hybrid levels]: <br></div><div><a href="https://www.ncl.ucar.edu/Applications/wind.shtml"><b>https://www.ncl.ucar.edu/Applications/wind.shtml</b></a></div><div><br></div><div>========================</div><div>Good Luck<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, Jul 13, 2023 at 9:35 PM Anie K Lal via ncl-talk <<a href="mailto:ncl-talk@mailman.ucar.edu">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">Hi all,<br>
<br>
I want to calculate the vertical vapor flux divergence <br>
(1/g*[d(qw)/dz])from a WRF output.<br>
<br>
I found the vapor flux as below:<br>
   shum        =   mixhum_ptrh(pres,temp,rh,2)  ;units: <br>
(hPA,K,%,kg/kg);calculate the specific humidity(2),(1)means mixing ratio<br>
   ;printVarSummary(shum)<br>
   uq          =   ua * shum                        ;vapor flux<br>
   vq          =   va * shum<br>
   wq          =   wa * shum<br>
<br>
Not sure how to find the divergence.<br>
I have gone through few examples in:<br>
<a href="https://www.ncl.ucar.edu/Applications/wind.shtml" rel="noreferrer" target="_blank">https://www.ncl.ucar.edu/Applications/wind.shtml</a><br>
<a href="https://www.ncl.ucar.edu/Applications/ESMF.shtml" rel="noreferrer" target="_blank">https://www.ncl.ucar.edu/Applications/ESMF.shtml</a><br>
<br>
All consider only u and v. How do I do it?<br>
<br>
Looking forward to any kind of help or advise.<br>
<br>
Thank you.<br>
Anie K Lal<br>
PhD Student<br>
IIT Delhi<br>
India<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>