[ncl-talk] Vertical integration in NCL

[Dennis Shea]

vibeta was created for rawindsonde data that had many missing values.
I would suggest *not* using it.
===============================================================================
The file you sent to me (offline) is a classic example of something
ncl-talk constantly advises
*****ALWAYS, LOOK AT YOUR DATA*****
----

The file you sent has some confounding information. For example, I am not
sure how one interprets the 'lev' vertical coordinate


[SNIP]
        float lev(lev) ;
                lev:units = "hPa" ;
                lev:long_name = "height" ;
                lev:standard_name = "height" ;
                lev:axis = "Z" ;
                lev:positive = "down" ;
[SNIP]

Clearly, the following are *not* pressure levels with units "hPa". They are
'height' values. Likely with units "m"

 lev =
    250, 750, 1250, 1750, 2250, 2750, 3250, 3750, 4250, 4750, 5250, 5750,
    6250, 6750, 7250, 7750, 8250, 8750, 9250, 9750, 10250, 10750, 11250,
    11750, 12250 ;

===
Further, the *file*'s attributes include:

:Info2 = "(Height from the ground; in meter)" ;

===
To me, this means 'zsfc=0' See attached script.
It produces the following. PLEASE carefully look at the dimensions of the
following:

Variable: x
Type: double
Total Size: 19660800 bytes
            2457600 values
Number of Dimensions: 4
Dimensions and sizes:    [time | 12] x [lev | 25] x [lat | 64] x [lon | 128]
Coordinates:
            time: [ 1..12]
            lev: [250..12250]
            lat: [-87.8638..87.8638]
            lon: [ 0..357.1875]
Number Of Attributes: 2
  units :    (K/sec)
  standard_name :    Atmospheric heating (K/sec)
(0)    Atmospheric heating (K/sec): min=0   max=0.0001740113802952692
(0)    ---

Variable: dz
Type: float
Total Size: 9830400 bytes
            2457600 values
Number of Dimensions: 4
Dimensions and sizes:    [12] x [25] x [64] x [128]
Coordinates:
Number Of Attributes: 1
  _FillValue :    9.96921e+36
(0)    ---

Variable: xint
Type: double
Total Size: 786432 bytes
            98304 values
Number of Dimensions: 3
Dimensions and sizes:    [time | 12] x [lat | 64] x [lon | 128]  <== *no*
vertical dimension
Coordinates:
            time: [ 1..12]
            lat: [-87.8638..87.8638]
            lon: [ 0..357.1875]
Number Of Attributes: 3
  units :    (K-m/sec)
  long_name :    Vertical Integral: (Absp3D)*dz
  _FillValue :    9.969209968386869e+36
(0)    Vertical Integral: (Absp3D)*dz: min=0   max=0.2500214350220631

===============================================================================
Note: 6.4.0 will have

http://www.ncl.ucar.edu/Document/Functions/Contributed/wgt_vertical_n.shtml


Good Luck

On Tue, Feb 7, 2017 at 9:55 PM, Anahita Amiri Farahani <aamir003 at ucr.edu>
wrote:

> Hello all,
>
> I have found this function in ncl to calculate vertical integration of
> atmospheric heating:
>
> http://www.ncl.ucar.edu/Document/Functions/Built-in/vibeta.shtml
>
> are there any other functions?
>
> Thank you,
> Ana
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at ucar.edu
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170212/98f73bf1/attachment.html 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: vertical_int.ncl
Type: text/x-ncl
Size: 1140 bytes
Desc: not available
Url : http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170212/98f73bf1/attachment.bin