[ncl-talk] A missing plot

Barry Lynn barry.h.lynn at gmail.com
Tue Oct 16 12:59:53 MDT 2018


HI:

Thank you, one can see the problem below.

Here is the print summary when I don't interpolate to the GEFS grid.  When
compared to when I pass through hilowres (below), you can see that the
lat_0, lon_0 are different.  It doesn't map in the first case because lon_0
is outside my mapping area. .

Variable: fgen_gefs

Type: float

Total Size: 15620 bytes

            3905 values

Number of Dimensions: 2

Dimensions and sizes: [lat_0 | 55] x [lon_0 | 71]

Coordinates:

            lat_0: [52..25]

            lon_0: [90..125]

Number Of Attributes: 3

  _FillValue : 0

  units : K/100 km/3 hr

  short_name : Petterssen Frontogenesis


Here is the printVarSummary when I interpolate first using hi2lowres


Variable: fgen_gefs

Type: float

Total Size: 15620 bytes

            3905 values

Number of Dimensions: 2

Dimensions and sizes: [lat_0 | 55] x [lon_0 | 71]

Coordinates:

            lat_0: [25..52]

            lon_0: [-90..-55]

Number Of Attributes: 3

  _FillValue : 0

  units : K/100 km/3 hr

  short_name : Petterssen Frontogenesis

On Tue, Oct 16, 2018 at 9:10 PM Mary Haley <haley at ucar.edu> wrote:

> Hi Barry,
>
> It helps to include the error message exactly as it appears, instead of
> paraphrasing it, because the wording is very important.
>
> I would do a printVarSummary on fgen_gefs right before you plot it, so you
> can see what exactly is attached to it.
>
> My guess is that you need to add a units to the lat/lon coordinate arrays
> of fgen_gefs. Something like this:
>
> fgen_gefs&lat at units = "degrees_north"
> fgen_gefs&lon at units = "degrees_east"
>
> Of course, without seeing your data, I have no idea if fgen_gefs even has
> coordinate arrays. This depends on what th_regrid looks like, so you should
> also do a printVarSummary on that.
>
> --Mary
>
>
>
>
>
>
> On Tue, Oct 16, 2018 at 10:15 AM, Barry Lynn <barry.h.lynn at gmail.com>
> wrote:
>
>> Hello:
>>
>> I am am trying to make some panel plots.
>>
>> Strangely, when I copy the coordinates to fgen_gefs the data has defined
>> coordinates, but the map disappears. WIthout the coordinates, it appears,
>> but I get a warning that the coordinates north south, etc, are not defined.
>>
>> copy_VarCoords(th_regrid,fgen_gefs)
>>
>> I use the same coordinates for the other maps which appear fine.
>>
>> The only difference between these plots is that I haven't interpolated
>> the data in the first map  to the GEFS grid because it is the GEFS data.
>>
>> I don't see any error messages.
>>
>> I am wondering what I should print out in order to find the error.
>>
>> Barry
>> P.S. By the way, if I interpolate the gefs grid to the gefs grid the map
>> appears.
>>
>>
>>
>>
>>
>> --
>> Barry H. Lynn, Ph.D
>> Senior Associate Scientist, Lecturer,
>> The Institute of the Earth Science,
>> The Hebrew University of Jerusalem,
>> Givat Ram, Jerusalem 91904, Israel
>> Tel: 972 547 231 170
>> Fax: (972)-25662581
>>
>> C.E.O, Weather It Is, LTD
>> Weather and Climate Focus
>> http://weather-it-is.com
>> Jerusalem, Israel
>> Local: 02 930 9525
>> Cell: 054 7 231 170
>> Int-IS: x972 2 930 9525
>>
>>
>> _______________________________________________
>> ncl-talk mailing list
>> ncl-talk at ucar.edu
>> List instructions, subscriber options, unsubscribe:
>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>
>>
>

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

C.E.O, Weather It Is, LTD
Weather and Climate Focus
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: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20181016/db8b3cc8/attachment.html>


More information about the ncl-talk mailing list