[ncl-talk] Calculation of vertical vapor flux divergence from WRF out file

[Anie K Lal]

Hi all,

I want to calculate the vertical vapor flux divergence 
(1/g*[d(qw)/dz])from a WRF output.

I found the vapor flux as below:
   shum        =   mixhum_ptrh(pres,temp,rh,2)  ;units: 
(hPA,K,%,kg/kg);calculate the specific humidity(2),(1)means mixing ratio
   ;printVarSummary(shum)
   uq          =   ua * shum                        ;vapor flux
   vq          =   va * shum
   wq          =   wa * shum

Not sure how to find the divergence.
I have gone through few examples in:
https://www.ncl.ucar.edu/Applications/wind.shtml
https://www.ncl.ucar.edu/Applications/ESMF.shtml

All consider only u and v. How do I do it?

Looking forward to any kind of help or advise.

Thank you.
Anie K Lal
PhD Student
IIT Delhi
India