<div dir="ltr">Hi Rick,<div><br></div><div>thanks for getting back to me. I do get a plot which however looks distorted (see attached). </div><div>Yes the xgrid and ygrid are the x-offset and y-offset <span style="color:rgb(71,71,71);font-family:'Open Sans','Segoe UI',Arial,Helvetica,sans-serif;font-size:14px;line-height:25.200000762939453px;background-color:rgb(247,248,250)">in meters of the projection grid centers (taken from </span><a href="http://nsidc.org/data/G02202#daily-projection-var" target="_blank">http://nsidc.org/data/<wbr>G02202#daily-projection-var</a>, under the tab 'user guide'). </div><div><br></div><div>For example reading in the latitude gives the following:</div><div><br></div><div><div>Variable: latitude</div><div>Type: double</div><div>Total Size: 1089536 bytes</div><div> 136192 values</div><div>Number of Dimensions: 2</div><div>Dimensions and sizes:<span class="m_6928723214790830998gmail-Apple-tab-span" style="white-space:pre-wrap"> </span>[ygrid | 448] x [xgrid | 304]</div><div>Coordinates:</div><div> ygrid: [5837500..-5337500]</div><div> xgrid: [-3837500..3737500]</div><div>Number Of Attributes: 5</div><div> standard_name :<span class="m_6928723214790830998gmail-Apple-tab-span" style="white-space:pre-wrap"> </span>latitude</div><div> long_name :<span class="m_6928723214790830998gmail-Apple-tab-span" style="white-space:pre-wrap"> </span>latitude</div><div> units :<span class="m_6928723214790830998gmail-Apple-tab-span" style="white-space:pre-wrap"> </span>degrees_north</div><div> valid_range :<span class="m_6928723214790830998gmail-Apple-tab-span" style="white-space:pre-wrap"> </span>( 0, 90 )</div><div> _FillValue :<span class="m_6928723214790830998gmail-Apple-tab-span" style="white-space:pre-wrap"> </span>-999</div></div><div><br></div><div>and equivalent for longitude. Printing the single values gives reasonable results:</div><div><div>(0)<span class="m_6928723214790830998gmail-Apple-tab-span" style="white-space:pre-wrap"> </span>31.9602770108575</div><div>(1)<span class="m_6928723214790830998gmail-Apple-tab-span" style="white-space:pre-wrap"> </span>32.1115587186697</div><div>(2)<span class="m_6928723214790830998gmail-Apple-tab-span" style="white-space:pre-wrap"> </span>32.2628753734599</div><div>(3)<span class="m_6928723214790830998gmail-Apple-tab-span" style="white-space:pre-wrap"> </span>32.4142244216501</div><div>(4)<span class="m_6928723214790830998gmail-Apple-tab-span" style="white-space:pre-wrap"> </span>32.5656032790591</div><div><br></div></div><div>However there remains a problem with the plotting:</div><div> "warning:ScalarFieldSetValues: 2d coordinate array sfXArray has an incorrect dimension size: defaulting sfXArray" comes when I set 'res@sfXArray' and 'res@sfYArray' </div><div> Without setting these resources:</div><div>"(0)<span class="m_6928723214790830998gmail-Apple-tab-span" style="white-space:pre-wrap"> </span>check_for_y_lat_coord: Warning: Data either does not contain a valid latitude coordinate array or doesn't contain one at all." </div><div>So definitely there is a problem with the coordinates. I don't which function or resource can handle this problem in NCL.</div><div>I attached also the printVarSummary of the sea ice concentration variable. </div><div><br></div><div>thanks in advance</div><div>Anne</div><div><br></div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On 23 October 2017 at 17:36, Rick Brownrigg <span dir="ltr"><<a href="mailto:brownrig@ucar.edu" target="_blank">brownrig@ucar.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div>Hi,<br><br></div>Are you having trouble georeferencing the plot, or not getting any sort of plot at all. From the printVarSummary, it looks like the x/y coordinates are in meters on the ground, relative to some polar stereographic cartographic projection. Is there any sort of other metadata describing the projection? If you are able to share the file, I can to take a look. Otherwise, what is the output of a print-statement on the file-variable that reads this grid?<br><br></div>Rick<br></div><div class="gmail_extra"><br><div class="gmail_quote"><div><div class="h5">On Mon, Oct 23, 2017 at 9:27 AM, Anne <span dir="ltr"><<a href="mailto:anne.seidenglanz@unive.it" target="_blank">anne.seidenglanz@unive.it</a>></span> wrote:<br></div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div class="h5"><div dir="ltr">Hello everybody,<div><br></div><div>I have a problem plotting data that have 2D latitude/longitude dimensions (in particular monthly sea ice concentration in the Arctic from the NSIDC, National Snow and Ice Data Center), even though I have used the </div><div><br></div><div><span style="color:rgb(80,0,80);font-family:Arial;font-size:12px">res@sfXArray = lon2d</span><br style="color:rgb(80,0,80);font-family:Arial;font-size:12px"><span style="color:rgb(80,0,80);font-family:Arial;font-size:12px">res@sfYArray = lat2d</span><br></div><div><span style="color:rgb(80,0,80);font-family:Arial;font-size:12px"><br></span></div><div><font color="#500050" face="Arial">resources as I have done earlier for another case (traingular ocean grid) which worked fine, however not so in this case.This is the information I got about the grid:</font></div><div><font color="#2b2e2f" face="Lucida Sans Unicode, Lucida Grande, Tahoma, Verdana, sans-serif">"The grid is in the NSIDC Sea Ice Polar Stereographic [North/South] projection. This is an irregular grid in that the latitude and longitude are 2-dimensional arrays, where the x and the y component are specified by xgrid and ygrid component (see below)"</font><font color="#500050" face="Arial"><br></font></div><div><font color="#2b2e2f" face="Lucida Sans Unicode, Lucida Grande, Tahoma, Verdana, sans-serif"><br></font></div><div><font color="#2b2e2f" face="Lucida Sans Unicode, Lucida Grande, Tahoma, Verdana, sans-serif">This is part of the the printVarSummary of the variable:</font></div><div><font color="#2b2e2f" face="Lucida Sans Unicode, Lucida Grande, Tahoma, Verdana, sans-serif"><br></font></div><div><div>Variable:<b> seaice_conc_monthly_cdr</b></div><div>Type: byte</div><div>Total Size: 46441472 bytes</div><div> 46441472 values</div><div>Number of Dimensions: 3</div><div>Dimensions and sizes:<span class="m_5438400607328969483m_3409168039474519372gmail-Apple-tab-span" style="white-space:pre-wrap"> </span>[time | 341] x <b>[ygrid | 448]</b> x <b>[xgrid | 304]</b></div><div>Coordinates:</div><div> time: [141195..151544]</div><div> <b> ygrid: [5837500..-5337500]</b></div><div><b> xgrid: [-3837500..3737500]</b></div><div>Number Of Attributes: 14</div></div><div>...</div><div><br></div><div><font size="1"><font color="#2b2e2f" face="Lucida Sans Unicode, Lucida Grande, Tahoma, Verdana, sans-serif"><br></font></font></div><div>I managed to convert time into a useful date using the cd_calender function, but was not successful with the latitude and longitude (neither with ascribing explicitly latitude and longitude to the 2nd and 3rd coordinate dimension, respectively).</div><div><br></div><div>Did I overlook something?</div><div><br></div><div>Thanks in advance,</div><div>Anne</div><div><br></div><div><br></div></div>
<br></div></div>______________________________<wbr>_________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailma<wbr>n/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>
</blockquote></div><br></div>