<div dir="ltr"><div>Looking at previous examples with RAP files, I found the following defined in resources as <br></div><div>res@mpLambertParallel1F = lon2d@Latin1<br> res@mpLambertParallel2F = lon2d@Latin2<br> res@mpLambertMeridianF = lon2d@Lov<br> res@mpCenterLonF = lon2d@Lo1<br> res@mpCenterLatF = lon2d@La1</div><div><br></div><div>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 ? <br></div><div><br></div><div>Thanks</div><div>Deba<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, Jun 29, 2023 at 2:27 PM Debasish Hazra <<a href="mailto:debasish.hazra5@gmail.com">debasish.hazra5@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div>Hi Rick,</div><div>Thanks for your reply. When I turned <br></div><div>res@tfDoNDCOverlay = False</div><div>got the following error</div><div><br></div><div> check_for_y_lat_coord: Warning: Data either does not contain<br>(0) a valid latitude coordinate array or doesn't contain one at all.<br>(0) A valid latitude coordinate array should have a 'units'<br>(0) attribute equal to one of the following values: <br>(0) 'degrees_north' 'degrees-north' 'degree_north' 'degrees north' 'degrees_N' 'Degrees_north' 'degree_N' 'degreeN' 'degreesN' 'deg north'<br>(0) check_for_lon_coord: Warning: Data either does not contain<br>(0) a valid longitude coordinate array or doesn't contain one at all.<br>(0) A valid longitude coordinate array should have a 'units'<br>(0) attribute equal to one of the following values: <br>(0) 'degrees_east' 'degrees-east' 'degree_east' 'degrees east' 'degrees_E' 'Degrees_east' 'degree_E' 'degreeE' 'degreesE' 'deg east'<br>warning:ContourPlotDraw: out of range coordinates encountered; standard AreaFill rendering method may be unreliable;<br> consider setting the resource trGridType to "TriangularMesh" if coordinates contain missing values</div><div><br></div><div>Debasish<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, Jun 29, 2023 at 2:13 PM Rick Brownrigg <<a href="mailto:brownrig@ucar.edu" target="_blank">brownrig@ucar.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div>Hi Debasish,</div><div><br></div><div>I'm a little confused by this line:</div><div><br></div><div> res@tfDoNDCOverlay = True</div><div><br></div><div>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.</div><div><br></div><div>Rick</div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, Jun 29, 2023 at 12:01 PM Debasish Hazra via ncl-talk <<a href="mailto:ncl-talk@mailman.ucar.edu" target="_blank">ncl-talk@mailman.ucar.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div>Hi,</div><div>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. <br></div><div><br></div><div> file global attributes:<br> dimensions:<br> ygrid_0 = 3268 <br> xgrid_0 = 5200</div><div><br></div><div> float gridrot_0 ( ygrid_0, xgrid_0 )<br> long_name : vector rotation angle<br> GridType : Polar Sterographic Projection (North or South)<br> units : radians<br> formula_u : Uearth = sin(rot)*Vgrid + cos(rot)*Ugrid<br> formula_v : Vearth = cos(rot)*Vgrid - sin(rot)*Ugrid<br> note1 : u and v components of vector quantities are resolved relative to grid<br> note2 : apply formulas to derive u and v components relative to earth<br> <br> float gridlat_0 ( ygrid_0, xgrid_0 )<br> corners : ( -4, -5.706548, 18.53903, 22.48842 )<br> long_name : latitude<br> grid_type : Polar Sterographic Projection (North or South)<br> units : degrees_north<br> Dy : 3.17<br> Dx : 3.17<br> Lov : 245 <br> Lo1 : 206.5<br> La1 : -4 <br> <br> float gridlon_0 ( ygrid_0, xgrid_0 )<br> corners : ( -153.5, -74.4254, -22.46207, 152.2681 )<br> long_name : longitude<br> grid_type : Polar Sterographic Projection (North or South)<br> units : degrees_east<br> Dy : 3.17<br> Dx : 3.17 <br> Lov : 245<br> Lo1 : 206.5<br> La1 : -4</div><div><br></div><div>Any help on this is appreciated. <br></div><div>Thanks</div><div>Deba<br></div></div>
_______________________________________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@mailman.ucar.edu" target="_blank">ncl-talk@mailman.ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="https://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">https://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
</blockquote></div>
</blockquote></div>
</blockquote></div>