[ncl-talk] Help : calculation of Total precipitable water, overestimated values

[Dhirendra Kumar]

I am trying to calculate total precipitable water from 6 hourly data
containing specific humidity (kg/kg) at different pressure levels (hPa). I
am able to write the output successfully but I think the  values are
overestimated 10 times of the usual value.

I checked my code but could not figure out the mistake.

I will appreciate if people here can help me finding the error. The NCL
version for me is 6.4. The code written is following:

begin
f=addfile("ps_qas.nc","r")
sphum=f->qas
lat=f->lat
lon=f->lon
time=f->time
dp=f->plev*100                                                            ;
conversion into Pa
print(dp)
ps=f->ps(time|:,lat|:,lon|:)
re_sphum=sphum(time|:,lat|:,lon|:,plev|:)
pw = prcwater_dp (re_sphum, dp)
pw at long_name = "total column precipitable water"
pw at units     = "kg/m2"
printVarSummary(pw)
pw at coordinates = "lat lon"
copy_VarCoords(ps,pw)
out=addfile("pw.nc","w")
out->pw=pw
out->lat=lat
out->lon= lon
out->time= time
end


Thanks in anticipation
Dhirendra

-- 
========================================
Dhirendra Kumar
PhD Student
School of Environmental Sciences, JNU
New Delhi (India)
Contact: +91 9910778043
Alternate email: dhirendra.cub at outlook.com
=========================================
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170607/45c58b4e/attachment.html