[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