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

[Jim Means]

Dhirendra,

Since I don't have access to the netcdf file that you're using, perhaps 
I'm mistaken, but aren't the plev the pressure levels? I think the dp 
that's used in the function prcwater_dp are the differences between the 
levels (the layer thicknesses), if you're trying to do an integration to 
get the precipitable water.

I think if you used the levels instead of their thicknesses that would 
definitely give you a value that was too large.

Hope that helps.

Jim


On 6/7/2017 12:37 AM, Dhirendra Kumar wrote:
> 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 <http://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 <http://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 
> <mailto:dhirendra.cub at outlook.com>
> =========================================
>
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at ucar.edu
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk

-- 
Untitled Document <http://www.weatherextreme.com>

*James D. Means, Ph.D.*
Atmospheric and Climate Scientist
jim at weatherextreme.com <mailto:jim at weatherextreme.com>|vcard 
<http://www.weatherextreme.com/vcards/James%20MEANS.vcf>
930 Tahoe Blvd., Suite 802-560
Incline Village, Nevada 89451
www.weatherextreme.com <http://www.weatherextreme.com>

<https://www.facebook.com/wxextreme> <https://twitter.com/wxextreme> 
<http://www.linkedin.com/company/weatherextreme>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170607/db3dc67b/attachment.html 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: signaturelogo.png
Type: image/png
Size: 8813 bytes
Desc: not available
Url : http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170607/db3dc67b/attachment.png 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: facebookicon.png
Type: image/png
Size: 3261 bytes
Desc: not available
Url : http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170607/db3dc67b/attachment-0001.png 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: twittericon.png
Type: image/png
Size: 3524 bytes
Desc: not available
Url : http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170607/db3dc67b/attachment-0002.png 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: linkedin.png
Type: image/png
Size: 3508 bytes
Desc: not available
Url : http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170607/db3dc67b/attachment-0003.png