[ncl-talk] Strange problem related to linint1_Wrap

Barry Lynn barry.h.lynn at gmail.com
Tue Oct 31 00:36:02 MDT 2023 

Hello:

Strangely, I was able to solve the problem with the extra white space this
way (by limiting the vertical extent of the data passed into the
gsm_contour_map program).

    *ff1_x = ff1(0:25,:)*

    plot  = gsn_csm_contour(wks, ff1, res )

    plots(i_file) = plot

    plot  = gsn_csm_contour(wks, ff1_x, res )

    plots(i_file+1) = plot


I say strangely because reducing the vertical extent of the data read into
the linear interpolation program just led to a cropping of the radar data.
That makes no sense, as the map window clearly shows white space above (a).

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

> 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
>
>

-- 

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/20231031/01481c0e/attachment.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: radar_vert_crossection_10UTC.pdf
Type: application/pdf
Size: 84958 bytes
Desc: not available
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20231031/01481c0e/attachment.pdf>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: test3.ncl
Type: application/octet-stream
Size: 11521 bytes
Desc: not available
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20231031/01481c0e/attachment.obj>

More information about the ncl-talk mailing list