<div dir="ltr"><div class="gmail_default" style="font-size:small">Michele,</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">All NCL questions should be directed to ncl-talk, and not ncarg-talk. This forum is for people who are using the older NCAR Graphics C and Fortran routines to do graphics.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">I see you've already posted to ncl-talk, so thanks for doing this.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">--Mary</div><div class="gmail_default" style="font-size:small"><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Jun 24, 2015 at 1:56 AM, Michele Petrini <span dir="ltr"><<a href="mailto:mpetrini139@yahoo.it" target="_blank">mpetrini139@yahoo.it</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div text="#000000" bgcolor="#FFFFFF">
Dear helping desk,<br>
I am a beginner with NCL and I have the following issue. I have a
netcdf file with some data obtained from an ice sheet model run. The
model transforms latitude-longitude coordinates into native
Cartesian distances coordinates via the Lambert equal Area
projection, using as a center of projection the North Pole (0,90).
Hence my data in the netcdf file are geographically dependent on
native Cartesian distances coordinates. I am trying now to plot a
contour over a map of the Northern Emisphere of one of the netcdf
file variables (let's say the ice thickness). First of all, I
imported succesfully the variables I am interested to, depending on
time, x and y, along with the latitude and longitude variables, on x
and y (where x and y are defined in native Cartesian distances from
(0,90). In order to go back from the native Cartesian distances to
latitude and longitude, I define a new array <i>without </i>attributes
(if the original variable is H, I just insert Hlatlon= (\ H \)).
Then I name the three dimensions of the new variable Hlatlon (time,
lat, lon), and I assign to the dimension number 1 and 2 (using the
NCL enumeration from 0) the coordinate variables lat and lon.
Finally, I plot the contour of ice thickness over the map. What I
get when I run the script is wrong (I'm sure as, for instance,
Greenland is ice free). Where I did wrong? Is there a better (faster
or more correct) way to plot my contour with such variables?<br>
Thank you for your help,<br>
<br>
Michele<br>
<pre cols="72">--
***
Michele Petrini
Ph.D. student in Earth Science and Fluid Mechanics
Università degli studi di Trieste,
Dipartimento di Matematica e Geoscienze
Palazzina C - via Weiss 1, 34128 Trieste, Italy
Email: <a href="mailto:mpetrini139@yahoo.it" target="_blank">mpetrini139@yahoo.it</a>
Skype: michele.petrins
Mobile: <a href="tel:%2B39%203398367372" value="+393398367372" target="_blank">+39 3398367372</a>
</pre>
</div>
<br>_______________________________________________<br>
ncarg-talk mailing list<br>
<a href="mailto:ncarg-talk@ucar.edu">ncarg-talk@ucar.edu</a><br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncarg-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncarg-talk</a><br>
<br></blockquote></div><br></div>