[ncl-talk] wrf_cape_2d errors
Kai Zhang
kzkaizhang at gmail.com
Mon Sep 12 20:19:09 MDT 2016
Dear Bill,
Thanks so much for your response!
I've check in several ways that, all of the units are correct. And
non missing value, non strange value. The problem still exists.
Here are the result for printVarSummary on each of my input variables:
Variable: p
Type: float
Total Size: 595200 bytes
148800 values
Number of Dimensions: 4
Dimensions and sizes: [hour | 8] x [lev | 24] x [lat | 31] x [lon | 25]
Coordinates:
hour: [0..21]
lev: [10000..100000]
lat: [-48.31578947368421..8.526315789473685]
lon: [ 270.. 330]
Number Of Attributes: 1
units : Pa
Variable: T
Type: float
Total Size: 595200 bytes
148800 values
Number of Dimensions: 4
Dimensions and sizes: [hour | 8] x [lev | 24] x [lat | 31] x [lon | 25]
Coordinates:
hour: [0..21]
lev: [10000..100000]
lat: [-48.31578947368421..8.526315789473685]
lon: [ 270.. 330]
Number Of Attributes: 5
day : 1
year : 1
units : K
long_name : Temperature
_FillValue : 9.96921e+36
Variable: Q
Type: float
Total Size: 595200 bytes
148800 values
Number of Dimensions: 4
Dimensions and sizes: [hour | 8] x [lev | 24] x [lat | 31] x [lon | 25]
Coordinates:
hour: [0..21]
lev: [10000..100000]
lat: [-48.31578947368421..8.526315789473685]
lon: [ 270.. 330]
Number Of Attributes: 5
day : 1
year : 1
units : kg/kg
long_name : Specific humidity
_FillValue : 9.96921e+36
Variable: Z
Type: float
Total Size: 595200 bytes
148800 values
Number of Dimensions: 4
Dimensions and sizes: [hour | 8] x [lev | 24] x [lat | 31] x [lon | 25]
Coordinates:
hour: [0..21]
lev: [10000..100000]
lat: [-48.31578947368421..8.526315789473685]
lon: [ 270.. 330]
Number Of Attributes: 5
day : 1
year : 1
units : m
long_name : Geopotential Height (above sea level)
_FillValue : 9.96921e+36
Variable: zs
Type: float
Total Size: 24800 bytes
6200 values
Number of Dimensions: 3
Dimensions and sizes: [8] x [31] x [25]
Coordinates:
Number Of Attributes: 2
units : m
_FillValue : 1e+15
Variable: ps
Type: float
Total Size: 24800 bytes
6200 values
Number of Dimensions: 3
Dimensions and sizes: [8] x [31] x [25]
Coordinates:
Number Of Attributes: 1
units : Pa
capecalc3d: Outside of lookup table bounds. prs,thte= 925.000000000000
856.034700407579
Thanks again!
Best,
Kai
On Mon, Sep 12, 2016 at 10:41 AM, Bill Ladwig <ladwig at ucar.edu> wrote:
> Hi Kai,
>
> wrf_cape_2d requires the pressure argument (first argument) to be in units
> of Pa and the surface pressure argument (6th argument) to be in units of
> hPa, so double check that. Temperature needs to be in Kelvin. Height is
> in meters. Water vapor mixing ratio is in units of kg/kg. The error
> message you are getting is indicating that theta-e is being calculated to
> 856.03 K, which is obviously nonsensical, so one of your input arrays is in
> the wrong units, contains missing values, or has a bogus number in it.
> This algorithm relies on lookup tables, which is why things can go wrong in
> a hurry if something is not right. If you're still having problems, please
> reply with the result for printVarSummary on each of your input variables.
>
> Hope this helps,
>
> Bill
>
> On Fri, Sep 9, 2016 at 8:34 PM, Kai Zhang <kzkaizhang at gmail.com> wrote:
>
>> Dear NCL support,
>>
>> I am currently using function wrf_cape_2d to calculate CIN and CAPE. But
>> I cannot stop getting errors:
>>
>> "capecalc3d: Outside of lookup table bounds. prs,thte= 925.000000000000
>> 856.034700407579"
>>
>> I am sure that I got units for each variable correctly.
>>
>> It is weird that, if I change "p (unit: Pa)" to "p (unit: hPa)", I can
>> get it run smoothly. But the CIN and CAPE calculated are obviously wrong.
>>
>> Do you know why?
>>
>>
>> Thanks a lot for your help!
>>
>> Best,
>>
>> Kai
>>
>> _______________________________________________
>> ncl-talk mailing list
>> ncl-talk at ucar.edu
>> List instructions, subscriber options, unsubscribe:
>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>
>>
>
--
-------------------------------------------------------------------------------
Kai Zhang
The University of Texas at Austin, TX (2012 - )
Graduate Research Assistant,
Jackson School of Geosciences
BS, Peking University, China (2008 - 2012)
Dept. Atmospheric and Oceanic Science, School of Physics
Cellphone: +1 512 520 6902
https://sites.google.com/site/kaizhanghomepage/
Email: kzkaizhang at utexas.edu <zhangkai at utexas.edu>
kzkaizhang at gmail.com
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20160912/fe17f941/attachment.html
More information about the ncl-talk
mailing list