Re: Report bug of function "mixhum_ptrh"

Hello,

You are correct. This will be fixed in the next release of NCL (6.2.1 or
6.3.0).

Thank you very much.

Regards
Dennis Shea


On Wed, Jun 25, 2014 at 12:39 PM, Mary Haley <haley@ucar.edu> wrote:

> Did you see this bug report?
>
> ---------- Forwarded message ----------
> From: huangyj@mail.iap.ac.cn <huangyj@mail.iap.ac.cn>
> Date: Mon, Jun 23, 2014 at 5:38 PM
> Subject: Report bug of function "mixhum_ptrh"
> To: ncl-talk <ncl-talk@ucar.edu>
>
>
> Hello,
>
> Today, I found a bug of funtion "mixhum_ptrh" (
> http://www.ncl.ucar.edu/Document/Functions/Built-in/mixhum_ptrh.shtml).
>
> On the website,
>
> p = 1000. ; hPa
> tk = 18. + 273.15 ; K
> rh = 46.5 ; %
> q = *mixhum_ptrh* (p, tk, rh, 2) ; q = 5.9e-3 (kg/kg)
> q = *mixhum_ptrh* (p, tk, rh,-2) ; q = 5.9 ( g/kg)
>
> Actually, *mixhum_ptrh* (p, tk, rh,-2) return 6.01846.
>
> I checked the source code of this funtion "mixhum_ptrh", and I found a
> bug:
>
>
> Line 46 of Fortran source code "mixhum_ptd.f", should be *IF
> (ABS(ISWIT).EQ.2) THEN.*
> So, when the value of *iswit* equal to *-2*, the function of
> "mixhum_ptrh" will return wrong values.
>
>
> Best Regards,
> Yongjie Huang
>
> ------------------------------
> HUANG Yongjie
> Key Laboratory of Cloud-Precipitation Physics and Severe Storm(LACS),
> Chinese Academy of Sciences
> Email Address: *huangyj@mail.iap.ac.cn* <huangyj@mail.iap.ac.cn> OR
> *huangynj@gmail.com* <huangynj@gmail.com>
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
>

Received on Thu Jun 26 12:24:48 2014