mixhum_ptd missing values

Dear all,

It seems that the function mixhum_ptd ignores a whole column of data if
there includes one missing value, instead of ignoring just the missing
value itself.
The following is a simple example:
;======================================
pp = (/500.,100000.,2000./)
tt = (/-32768,276.,-32768/)
tt@_FillValue = -32768
rr = mixhum_ptd(pp,tt,1)
print(pp)
print(tt)
print(rr)
;======================================
In this way, rr contains 3 missing values.
But if I comment off tt@_FillValue=-32768, rr is (/0,0.004718332,0/). Seems
more reasonable.

Any idea where I went wrong?
Thank you!

NCL version: 6.1.2
System: Linux grele 3.5.0-39-generic #60~precise1-Ubuntu SMP Wed Aug 14
15:38:41 UTC 2013 x86_64 x86_64 x86_64 GNU/Linux

Best,
Wenshan

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Thu Aug 22 22:30:13 2013