[ncl-talk] getting latitude writing on axis wrong

[Barry Lynn]

Hi:

You can try to estimate these values from the domain configuration.  Each
RAP forecast requires these values to set up the input files.  (e.g., 33
and 45, respectively).

Barry

On Fri, Jun 30, 2023 at 2:54 PM Debasish Hazra via ncl-talk <
ncl-talk at mailman.ucar.edu> wrote:

> Looking at previous examples with RAP files, I found the following defined
> in resources as
> res at mpLambertParallel1F    = lon2d at Latin1
>    res at mpLambertParallel2F    = lon2d at Latin2
>    res at mpLambertMeridianF     = lon2d at Lov
>    res at mpCenterLonF           = lon2d at Lo1
>    res at mpCenterLatF           = lon2d at La1
>
> There is no Latin1 or Latin2 info in the grib2 file I am dealing with, but
> has Lov, Lo1 and La1 . Any suggestion where to get that info ?
>
> Thanks
> Deba
>
> On Thu, Jun 29, 2023 at 2:27 PM Debasish Hazra <debasish.hazra5 at gmail.com>
> wrote:
>
>> Hi Rick,
>> Thanks for your reply. When I turned
>> res at tfDoNDCOverlay = False
>> got the following error
>>
>>  check_for_y_lat_coord: Warning: Data either does not contain
>> (0)     a valid latitude coordinate array or doesn't contain one at all.
>> (0)     A valid latitude coordinate array should have a 'units'
>> (0)     attribute equal to one of the following values:
>> (0)         'degrees_north' 'degrees-north' 'degree_north' 'degrees
>> north' 'degrees_N' 'Degrees_north' 'degree_N' 'degreeN' 'degreesN' 'deg
>> north'
>> (0)     check_for_lon_coord: Warning: Data either does not contain
>> (0)     a valid longitude coordinate array or doesn't contain one at all.
>> (0)     A valid longitude coordinate array should have a 'units'
>> (0)     attribute equal to one of the following values:
>> (0)         'degrees_east' 'degrees-east' 'degree_east' 'degrees east'
>> 'degrees_E' 'Degrees_east' 'degree_E' 'degreeE' 'degreesE' 'deg east'
>> warning:ContourPlotDraw: out of range coordinates encountered; standard
>> AreaFill rendering method may be unreliable;
>>  consider setting the resource trGridType to "TriangularMesh" if
>> coordinates contain missing values
>>
>> Debasish
>>
>> On Thu, Jun 29, 2023 at 2:13 PM Rick Brownrigg <brownrig at ucar.edu> wrote:
>>
>>> Hi Debasish,
>>>
>>> I'm a little confused by this line:
>>>
>>>     res at tfDoNDCOverlay = True
>>>
>>> Should that be False?  If that doesn't help, if you can send me the
>>> dataset, I'll try to figure out what is going on.
>>>
>>> Rick
>>>
>>>
>>> On Thu, Jun 29, 2023 at 12:01 PM Debasish Hazra via ncl-talk <
>>> ncl-talk at mailman.ucar.edu> wrote:
>>>
>>>> Hi,
>>>> I am trying to plot grib2 file with Lambert projection and getting
>>>> wrong latitude markers (0-25 N) on the y-axis. Attached is figure and code
>>>> and also highlighting attributes from the file.
>>>>
>>>>  file global attributes:
>>>>    dimensions:
>>>>       ygrid_0 = 3268
>>>>       xgrid_0 = 5200
>>>>
>>>>       float gridrot_0 ( ygrid_0, xgrid_0 )
>>>>          long_name :    vector rotation angle
>>>>          GridType :     Polar Sterographic Projection (North or South)
>>>>          units :        radians
>>>>          formula_u :    Uearth = sin(rot)*Vgrid + cos(rot)*Ugrid
>>>>          formula_v :    Vearth = cos(rot)*Vgrid - sin(rot)*Ugrid
>>>>          note1 :        u and v components of vector quantities are
>>>> resolved relative to grid
>>>>          note2 :        apply formulas to derive u and v components
>>>> relative to earth
>>>>
>>>>       float gridlat_0 ( ygrid_0, xgrid_0 )
>>>>          corners :      ( -4, -5.706548, 18.53903, 22.48842 )
>>>>          long_name :    latitude
>>>>          grid_type :    Polar Sterographic Projection (North or South)
>>>>          units :        degrees_north
>>>>          Dy :   3.17
>>>>          Dx :   3.17
>>>>          Lov :  245
>>>>          Lo1 :  206.5
>>>>          La1 :  -4
>>>>
>>>>       float gridlon_0 ( ygrid_0, xgrid_0 )
>>>>          corners :      ( -153.5, -74.4254, -22.46207, 152.2681 )
>>>>          long_name :    longitude
>>>>          grid_type :    Polar Sterographic Projection (North or South)
>>>>          units :        degrees_east
>>>>          Dy :   3.17
>>>>          Dx :   3.17
>>>>          Lov :  245
>>>>          Lo1 :  206.5
>>>>          La1 :  -4
>>>>
>>>> Any help on this is appreciated.
>>>> Thanks
>>>> Deba
>>>> _______________________________________________
>>>> ncl-talk mailing list
>>>> ncl-talk at mailman.ucar.edu
>>>> List instructions, subscriber options, unsubscribe:
>>>> https://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>>
>>> _______________________________________________
> ncl-talk mailing list
> ncl-talk at mailman.ucar.edu
> List instructions, subscriber options, unsubscribe:
> https://mailman.ucar.edu/mailman/listinfo/ncl-talk
>


-- 

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/20230630/f6f4471a/attachment.htm>