Fwd: [ncl-talk] Report bug of function "mixhum_ptrh"

From: Mary Haley <haley_at_nyahnyahspammersnyahnyah>
Date: Wed Jun 25 2014 - 12:39:53 MDT

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 Wed Jun 25 06:39:54 2014

This archive was generated by hypermail 2.1.8 : Wed Jul 23 2014 - 15:33:46 MDT