<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style type="text/css" style="display:none"><!-- P { margin-top: 0px; margin-bottom: 0px; }--></style>
</head>
<body dir="ltr" style="font-size:12pt;color:#000000;background-color:#FFFFFF;font-family:Calibri,Arial,Helvetica,sans-serif;">
<div>Dear Ncl-Team</div>
<div> I'm using NCL - version 6.3.0.</div>
<div> I have netcdf file representing the bathymetry of all greenLand</div>
<div> I would like to do a countor map of a specific region of the bed topography (below there is a file description)</div>
<div> but this file is mapped in polar_stereographic and the</div>
<div> the variable bed is define on Cartesian x and y coordinate and not in lon and lat.</div>
<div><br>
</div>
<div> I applied polar plot like this</div>
<div><br>
</div>
<div> res = True ; plot mods desired</div>
<div> res@gsnPolar = "NH" ; specify the hemisphere</div>
<div> plot = gsn_csm_contour_map_polar(wks,u,res) ; create the plot</div>
<div><br>
</div>
<div> but it is not working. It give me this error:</div>
<div><br>
</div>
<div><br>
</div>
<div>(0)<span class="Apple-tab-span" style="white-space:pre"> </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>(0)<span class="Apple-tab-span" style="white-space:pre"> </span>A valid latitude coordinate array should have a 'units' attribute equal to one of the following values: </div>
<div>(0)<span class="Apple-tab-span" style="white-space:pre"> </span> 'degrees_north' 'degrees-north' 'degree_north' 'degrees north' 'degrees_N' 'Degrees_north' 'degree_N' 'degreeN' 'degreesN' 'deg north'</div>
<div>(0)<span class="Apple-tab-span" style="white-space:pre"> </span>check_for_lon_coord: Warning: Data either does not contain a valid longitude coordinate array or doesn't contain one at all.</div>
<div>(0)<span class="Apple-tab-span" style="white-space:pre"> </span>A valid longitude coordinate array should have a 'units' attribute equal to one of the following values: </div>
<div>(0)<span class="Apple-tab-span" style="white-space:pre"> </span> 'degrees_east' 'degrees-east' 'degree_east' 'degrees east' 'degrees_E' 'Degrees_east' 'degree_E' 'degreeE' 'degreesE' 'deg east'</div>
<div><br>
</div>
<div> How can I then convert the the cartesian x,y coordinate in a regular lat,lon degree grid?<br>
</div>
<div><br>
</div>
<div> I would be grateful is somebody could assist me.</div>
<div> Thanks.</div>
<div><br>
</div>
<div>Francesco</div>
<div><br>
</div>
<div><br>
</div>
<div>/===============================</div>
<div> netcdf BedMachineGreenland-2017-02-01 {</div>
<div>dimensions:</div>
<div><span class="Apple-tab-span" style="white-space:pre"></span>x = 10018 ;</div>
<div><span class="Apple-tab-span" style="white-space:pre"></span>y = 17946 ;</div>
<div>variables:</div>
<div><span class="Apple-tab-span" style="white-space:pre"></span>char mapping ;</div>
<div><span class="Apple-tab-span" style="white-space:pre"></span>mapping:geoid = "eigen-6c4" ;</div>
<div><span class="Apple-tab-span" style="white-space:pre"></span>mapping:false_easting = 0. ;</div>
<div><span class="Apple-tab-span" style="white-space:pre"></span>mapping:false_northing = 0. ;</div>
<div><span class="Apple-tab-span" style="white-space:pre"></span>mapping:grid_mapping_name = "polar_stereographic" ;</div>
<div><span class="Apple-tab-span" style="white-space:pre"></span>mapping:latitude_of_projection_origin = 90. ;</div>
<div><span class="Apple-tab-span" style="white-space:pre"></span>mapping:standard_parallel = 70. ;</div>
<div><span class="Apple-tab-span" style="white-space:pre"></span>mapping:straight_vertical_longitude_from_pole = -45. ;</div>
<div><span class="Apple-tab-span" style="white-space:pre"></span>int x(x) ;</div>
<div><span class="Apple-tab-span" style="white-space:pre"></span>x:long_name = "Cartesian x-coordinate" ;</div>
<div><span class="Apple-tab-span" style="white-space:pre"></span>x:standard_name = "projection_x_coordinate" ;</div>
<div><span class="Apple-tab-span" style="white-space:pre"></span>x:units = "meter" ;</div>
<div><span class="Apple-tab-span" style="white-space:pre"></span>int y(y) ;</div>
<div><span class="Apple-tab-span" style="white-space:pre"></span>y:long_name = "Cartesian y-coordinate" ;</div>
<div><span class="Apple-tab-span" style="white-space:pre"></span>y:standard_name = "projection_y_coordinate" ;</div>
<div><span class="Apple-tab-span" style="white-space:pre"></span>y:units = "meter" ;</div>
<div><span class="Apple-tab-span" style="white-space:pre"></span>byte mask(y, x) ;</div>
<div><span class="Apple-tab-span" style="white-space:pre"></span>mask:long_name = "mask (0 = ocean, 1 = ice-free land, 2 = grounded ice, 3 = floating ice, 4 = non-Greenland land)" ;</div>
<div><span class="Apple-tab-span" style="white-space:pre"></span>mask:units = "none" ;</div>
<div><span class="Apple-tab-span" style="white-space:pre"></span>mask:grid_mapping = "mapping" ;</div>
<div><span class="Apple-tab-span" style="white-space:pre"></span>mask:valid_range = 0b, 4b ;</div>
<div><span class="Apple-tab-span" style="white-space:pre"></span>mask:flag_values = 0b, 1b, 2b, 3b, 4b ;</div>
<div><span class="Apple-tab-span" style="white-space:pre"></span>mask:flag_meanings = "ocean ice_free_land grounded_ice floating_ice non_greenland_land" ;</div>
<div><span class="Apple-tab-span" style="white-space:pre"></span>mask:source = "gimpmask v2.0 (http://bprc.osu.edu/GDG/icemask.php) combined with coastline from Jeremie Mouginot" ;</div>
<div><span class="Apple-tab-span" style="white-space:pre"></span>short bed(y, x) ;</div>
<div><span class="Apple-tab-span" style="white-space:pre"></span>bed:long_name = "bed topography" ;</div>
<div><span class="Apple-tab-span" style="white-space:pre"></span>bed:standard_name = "bedrock_altitude" ;</div>
<div><span class="Apple-tab-span" style="white-space:pre"></span>bed:units = "meters" ;</div>
<div><span class="Apple-tab-span" style="white-space:pre"></span>bed:grid_mapping = "mapping" ;</div>
<div><span class="Apple-tab-span" style="white-space:pre"></span>bed:source = "Mass conservation (Mathieu Morlighem)" ;</div>
<div><span class="Apple-tab-span" style="white-space:pre"></span>bed:_FillValue = -9999s ;</div>
<div><span class="Apple-tab-span" style="white-space:pre"></span>short geoid(y, x) ;</div>
<div><span class="Apple-tab-span" style="white-space:pre"></span>geoid:long_name = "EIGEN-EC4 Geoid - WGS84 Ellipsoid difference" ;</div>
<div><span class="Apple-tab-span" style="white-space:pre"></span>geoid:standard_name = "geoid_height_above_reference_ellipsoid" ;</div>
<div><span class="Apple-tab-span" style="white-space:pre"></span>geoid:units = "meters" ;</div>
<div><span class="Apple-tab-span" style="white-space:pre"></span>geoid:grid_mapping = "mapping" ;</div>
<div><span class="Apple-tab-span" style="white-space:pre"></span>geoid:geoid = "eigen-6c4 (Forste et al 2014)" ;</div>
<div><br>
</div>
<div>// global attributes:</div>
<div><span class="Apple-tab-span" style="white-space:pre"></span>:Conventions = "CF-1.4" ;</div>
<div><span class="Apple-tab-span" style="white-space:pre"></span>:Title = "BedMachine Greenland: dataset based on mass conservation" ;</div>
<div><span class="Apple-tab-span" style="white-space:pre"></span>:Author = "Mathieu Morlighem" ;</div>
<div><span class="Apple-tab-span" style="white-space:pre"></span>:version = "01-Feb-2017 (v3.2)" ;</div>
<div><span class="Apple-tab-span" style="white-space:pre"></span>:nx = 17946. ;</div>
<div><span class="Apple-tab-span" style="white-space:pre"></span>:ny = 10018. ;</div>
<div><span class="Apple-tab-span" style="white-space:pre"></span>:Projection = "Polar Stereographic North (70N, 45W)" ;</div>
<div><span class="Apple-tab-span" style="white-space:pre"></span>:proj4 = "+init=epsg:3413" ;</div>
<div><span class="Apple-tab-span" style="white-space:pre"></span>:xmin = -637925 ;</div>
<div><span class="Apple-tab-span" style="white-space:pre"></span>:ymax = -657675 ;</div>
<div><span class="Apple-tab-span" style="white-space:pre"></span>:spacing = 150 ;</div>
<div><span class="Apple-tab-span" style="white-space:pre"></span>:no_data = -9999. ;</div>
<div><span class="Apple-tab-span" style="white-space:pre"></span>:Data_citation = "M. Morlighem, E. Rignot, J. Mouginot, H. Seroussi and E. Larour, Deeply incised submarine glacial valleys beneath the Greenland Ice Sheet, Nat. Geosci., 7, 418-422, 2014, doi:10.1038/ngeo2167,
http://www.nature.com/ngeo/journal/vaop/ncurrent/full/ngeo2167.html" ;</div>
<div><span class="Apple-tab-span" style="white-space:pre"></span>:Notes = "Data processed at the Department of Earth System Science, University of California, Irvine" ;</div>
<div>}</div>
<div><br>
<br>
</div>
<p><br>
</p>
</body>
</html>