[ncl-talk] missing values when finding depths of isotherms

[Denise Worthen - NOAA Affiliate]

Hi All,

I am using the example of iso_1.ncl to find the depth of a series of
isotherms. I am not getting errors, but I am not finding values where I
think there should be values.

I've attached a test version of my script and the output. If I look at two
nearby points, I get an isothermal depth at one point but not at a nearby
point with a similar T(z). I actually get values for isothermal depths over
large areas but they are "patchy". I don't understand why I'm not finding
depths in some places where I think I should be.

In the attached output for example, the T(z) at two nearby points look like
this:

index      Pt1    depth   Pt2
(12)    23.6228  125  22.7855
(13)    22.68  135  21.8488
(14)    21.717  145  20.8878
(15)    20.7226  155  19.8876
(16)    19.689  165  18.8448
(17)    18.6207  175  17.7681

The 20 deg isotherm should be at about 160m at Pt1 and 150m at Pt2. What I
get is a missing value at Pt1 and 153.6m at Pt2.

Any help appreciated,
Denise.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20190509/84ca02e1/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: depth_isotherms2.ncl
Type: application/octet-stream
Size: 931 bytes
Desc: not available
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20190509/84ca02e1/attachment.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: out
Type: application/octet-stream
Size: 2272 bytes
Desc: not available
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20190509/84ca02e1/attachment-0001.obj>