<div dir="ltr"><div><div>vibeta was created for rawindsonde data that had many missing values. <br></div>I would suggest *not* using it.<br>===============================================================================<br><div><div>The file you sent to me (offline) is a classic example of something ncl-talk constantly advises <br><b>***</b>ALWAYS, LOOK AT YOUR DATA<b>***</b><br>----<br><br>The file you sent has some confounding information. For example, I am not sure how one interprets the 'lev' vertical coordinate<br><br><br>[SNIP]<br> float lev(lev) ;<br> lev:units = "hPa" ;<br> lev:long_name = "height" ;<br> lev:standard_name = "height" ;<br> lev:axis = "Z" ;<br> lev:positive = "down" ;<br></div><div>[SNIP]<br><br></div><div>Clearly, the following are *not* pressure levels with units "hPa". They are 'height' values. Likely with units "m"<br><br></div></div> lev = <br> 250, 750, 1250, 1750, 2250, 2750, 3250, 3750, 4250, 4750, 5250, 5750, <br> 6250, 6750, 7250, 7750, 8250, 8750, 9250, 9750, 10250, 10750, 11250, <br> 11750, 12250 ;<br><br><div><div>===<br></div><div>Further, the <b>file</b>'s attributes include:<br></div><div><br>:Info2 = "(Height from the ground; in meter)" ;<br><br>===<br></div><div>To me, this means 'zsfc=0' See attached script. <br>It produces the following. PLEASE carefully look at the dimensions of the following:<br><br>Variable: x <br>Type: double<br>Total Size: 19660800 bytes<br> 2457600 values<br>Number of Dimensions: 4<br>Dimensions and sizes: [time | 12] x [lev | 25] x [lat | 64] x [lon | 128]<br>Coordinates: <br> time: [ 1..12]<br> lev: [250..12250]<br> lat: [-87.8638..87.8638]<br> lon: [ 0..357.1875]<br>Number Of Attributes: 2<br> units : (K/sec)<br> standard_name : Atmospheric heating (K/sec)<br>(0) Atmospheric heating (K/sec): min=0 max=0.0001740113802952692<br>(0) ---<br><br>Variable: dz<br>Type: float<br>Total Size: 9830400 bytes<br> 2457600 values<br>Number of Dimensions: 4<br>Dimensions and sizes: [12] x [25] x [64] x [128]<br>Coordinates: <br>Number Of Attributes: 1<br> _FillValue : 9.96921e+36<br>(0) ---<br><br>Variable: xint<br>Type: double<br>Total Size: 786432 bytes<br> 98304 values<br>Number of Dimensions: 3<br>Dimensions and sizes: [time | 12] x [lat | 64] x [lon | 128] <== <b>no</b> vertical dimension<br>Coordinates: <br> time: [ 1..12]<br> lat: [-87.8638..87.8638]<br> lon: [ 0..357.1875]<br>Number Of Attributes: 3<br> units : (K-m/sec)<br> long_name : Vertical Integral: (Absp3D)*dz<br> _FillValue : 9.969209968386869e+36<br>(0) Vertical Integral: (Absp3D)*dz: min=0 max=0.2500214350220631<br></div><div><br></div></div>===============================================================================<br>Note: 6.4.0 will have<br><br><a href="http://www.ncl.ucar.edu/Document/Functions/Contributed/wgt_vertical_n.shtml" target="_blank">http://www.ncl.ucar.edu/<wbr>Document/Functions/<wbr>Contributed/wgt_vertical_n.<wbr>shtml</a><br><br><br></div>Good Luck<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Feb 7, 2017 at 9:55 PM, Anahita Amiri Farahani <span dir="ltr"><<a href="mailto:aamir003@ucr.edu" target="_blank">aamir003@ucr.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hello all,<div><br></div><div>I have found this function in ncl to calculate vertical integration of atmospheric heating:</div><div><br></div><div><a href="http://www.ncl.ucar.edu/Document/Functions/Built-in/vibeta.shtml" target="_blank">http://www.ncl.ucar.edu/<wbr>Document/Functions/Built-in/<wbr>vibeta.shtml</a><br></div><div><br></div><div>are there any other functions?</div><div><br></div><div>Thank you,</div><div>Ana</div></div>
<br>______________________________<wbr>_________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/<wbr>mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>