[ncl-talk] Strange problem related to linint1_Wrap

[Barry Lynn]

Dear All (& Dave):

I still have the strange problem where there is a lot of white space above
the mapped radar cross-section.

If, however, I reduce the scope of the data (i_z is reduced), then the top
of the contours get cut off.

So, I have a feeling that the gsn plotting program is assuming some kind of
non-linear height field, while my data is in height coordinates.

Is this correct?  If not, what might be the solution?  Note, I'm not using
any WRF functions and I converted the WRF data to equally spaced data
levels.

Thank you,

Barry

On Mon, Oct 30, 2023 at 6:12 PM Barry Lynn <barry.h.lynn at gmail.com> wrote:

> Hi Dave:
>
> I thought those quotes were "suspicious."
>
> The labeling is the same regardless of which variable is called for the
> plot, so that wasn't the problem.
>
> *Regardless*, your solution worked!  (delete (ff1&HEIGHT)).
>
> Thank you!
>
> On Mon, Oct 30, 2023 at 5:20 PM Dave Allured - NOAA Affiliate <
> dave.allured at noaa.gov> wrote:
>
>> Correction.  Do not use quotes here:  delete (ff1&HEIGHT)
>>
>>
>> On Mon, Oct 30, 2023 at 9:16 AM Dave Allured - NOAA Affiliate <
>> dave.allured at noaa.gov> wrote:
>>
>>> Barry, please notice that there is a HEIGHT coordinate variable attached
>>> to ff1 after interpolation, but there is no such coordinate variable
>>> attached to ff1_s before interpolation.  It looks like the attached HEIGHT
>>> coordinate variable is misleading the plot function in some way.  Try
>>> delete ("ff1&HEIGHT").
>>>
>>> I am not familiar with the cross section plotting function.  But it
>>> looks like in the first plot, you are setting some plot resources to set
>>> the range, spacing, and units of the height axis.  I suspect that these
>>> settings are in conflict with the attached HEIGHT coordinates which may
>>> have been generated by linint1_Wrap.  The units are different.  Instead of
>>> deleting the HEIGHT coordinate variable, you could try deleting or
>>> adjusting the conflicting resource settings.
>>>
>>> If you want to rescale the attached HEIGHT coordinates from meters to
>>> km, then simply "ff1&HEIGHT = ff1&HEIGHT / 1000".
>>>
>>>
>>> On Mon, Oct 30, 2023 at 8:08 AM Barry Lynn via ncl-talk <
>>> ncl-talk at mailman.ucar.edu> wrote:
>>>
>>>> Hello:
>>>>
>>>> I wanted to use linint1 to interpolate WRF data to height surfaces with
>>>> constant spacing between levels.  But, there are two mysteries.
>>>>
>>>> You can see the WRF data in a cross-section of radar in (a; the
>>>> variable ff1_s below).  When I interpolate the data (b, ff1) is compacted
>>>> towards the bottom of the cross-sectional plot.  Moreover, the height
>>>> information on the y-axis will not plot.
>>>>
>>>> I understand why the data is compacted, but not why I can't eliminate
>>>> the white space above.  If I reduce the range of WRF data sampled in the z
>>>> direction, I just have a truncated cross-section (it doesn't get any closer
>>>> to the top of the cross-section).
>>>>
>>>> I don't see why using ff1 (from the interpolate) would make the
>>>> labeling disappear.
>>>>
>>>> Variable: ff1_s
>>>>
>>>> Type: float
>>>>
>>>> Total Size: 33728 bytes
>>>>
>>>>             8432 values
>>>>
>>>> Number of Dimensions: 2
>>>>
>>>> Dimensions and sizes: [height | 31] x [lon | 272]
>>>>
>>>> Coordinates:
>>>>
>>>> Number Of Attributes: 1
>>>>
>>>>   _FillValue : 9.96921e+36
>>>>
>>>>
>>>> Variable: ff1
>>>>
>>>> Type: float
>>>>
>>>> Total Size: 33728 bytes
>>>>
>>>>             8432 values
>>>>
>>>> Number of Dimensions: 2
>>>>
>>>> Dimensions and sizes: [HEIGHT | 31] x [lon | 272]
>>>>
>>>> Coordinates:
>>>>
>>>>             HEIGHT: [500..15500]
>>>>
>>>> Number Of Attributes: 1
>>>>
>>>>   _FillValue : 9.96921e+36
>>>>
>>>> --
>>>> Barry H. Lynn, Ph.D
>>>> Senior Scientist, Lecturer,
>>>> The Institute of Earth Sciences,
>>>> The Hebrew University of Jerusalem,
>>>> Givat Ram, Jerusalem 91904, Israel
>>>> Tel: 972 547 231 170
>>>> Fax: (972)-25662581
>>>>
>>>> Weather It Is, LTD
>>>> Weather and Climate Focus
>>>> https://weather-it-is.com <http://weather-it-is.com>
>>>> Jerusalem, Israel
>>>> Local: 02 930 9525
>>>> Cell: 054 7 231 170
>>>> Int-IS: x972 2 930 9525
>>>>
>>>
>
> --
>
> Barry H. Lynn, Ph.D
> Senior Scientist, Lecturer,
> The Institute of Earth Sciences,
> The Hebrew University of Jerusalem,
> Givat Ram, Jerusalem 91904, Israel
> Tel: 972 547 231 170
> Fax: (972)-25662581
>
> Weather It Is, LTD
> Weather and Climate Focus
> https://weather-it-is.com <http://weather-it-is.com>
> Jerusalem, Israel
> Local: 02 930 9525
> Cell: 054 7 231 170
> Int-IS: x972 2 930 9525
>
>

-- 

Barry H. Lynn, Ph.D
Senior Scientist, Lecturer,
The Institute of Earth Sciences,
The Hebrew University of Jerusalem,
Givat Ram, Jerusalem 91904, Israel
Tel: 972 547 231 170
Fax: (972)-25662581

Weather It Is, LTD
Weather and Climate Focus
https://weather-it-is.com <http://weather-it-is.com>
Jerusalem, Israel
Local: 02 930 9525
Cell: 054 7 231 170
Int-IS: x972 2 930 9525
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20231030/0a9dada6/attachment.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: radar_vert_crossection_10UTC.pdf
Type: application/pdf
Size: 58313 bytes
Desc: not available
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20231030/0a9dada6/attachment.pdf>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: test2.ncl
Type: application/octet-stream
Size: 11759 bytes
Desc: not available
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20231030/0a9dada6/attachment.obj>