<html>
<head>
<meta http-equiv="content-type" content="text/html; charset=ISO-8859-1">
</head>
<body bgcolor="#FFFFFF" text="#000000">
<font face="Courier New, Courier, monospace">Hi.<br>
<br>
A few years ago I had a problem with GMT/GDAL interactions such
that the GMT netCDF files were reflected when read into Qgis (for
example). This was worked around by running it through
'gdal_translate -of GMT'. There was also a discussion about this
being a bug in something in GDAL (since Qgis was using GDAL if I
remember correctly) that was supposed to have been fixed in some
version (obviously don't remember the details). At any rate, it
stopped being a problem for me but now something else is causing a
problem in files that have the following workflow:<br>
<br>
--> mbgrid (from MB-System) --> gdal_translate -of GMT
--> ArcGIS netCDF reader<br>
<br>
Now I find that gdal_translate appears to be changing the header
in a way such that ArcGIS cannot read the netCDF. The first
'ncdump -h' below is what the header looks like before running
through gdal_translate and the second one is after gdal_translate.<br>
<br>
I'm sending this note to MB-System, GMT and GDAL groups since it
seems like an issue of import to all three. Any thoughts on this?<br>
<br>
J.<br>
</font><br>
<font face="Courier New, Courier, monospace"><font face="Courier
New, Courier, monospace">=======================================Output
from mbgrid before gdal_translate
==================================================</font><br>
ncdump -h MCBCP_BeachModel_UTM11N_Baseline.grdclip.nc<br>
netcdf MCBCP_BeachModel_UTM11N_Baseline.grdclip {<br>
dimensions:<br>
x = 12626 ;<br>
y = 12251 ;<br>
variables:<br>
double x(x) ;<br>
x:long_name = "Easting (meters)" ;<br>
x:actual_range = 444057.827834802, 469309.046567392 ;<br>
double y(y) ;<br>
y:long_name = "Northing (meters)" ;<br>
y:actual_range = 3670183.00468016, 3694684.53612513 ;<br>
float z(y, x) ;<br>
z:long_name = "Topography (m)" ;<br>
z:_FillValue = NaNf ;<br>
z:actual_range = -1995.541015625, 599.818359375 ;<br>
<br>
// global attributes:<br>
:Conventions = "COARDS/CF-1.0" ;<br>
:title = "Topography Grid" ;<br>
:history = "grdclip -V
/media/EXT4_Projects_01/SLR_Data/level1/MCBCP/Basemap/MCBCP_BeachModel_UTM11N_Baseline.grd
-G/media/EXT4_Projects_01/SLR_Data/level1/MCBCP/Basemap/MCBCP_BeachModel_UTM11N_Baseline.grdclip.nc
-Sa600/NaN -Sb-2000/NaN" ;<br>
:description = "\n",<br>
"\tProjection: UTM11N\n",<br>
"\tGrid created by mbgrid\n",<br>
"\tMB-system Version 5.3.1982\n",<br>
"\tRun by <hellyj> on <StonestepsLinux> at
<Thu Apr 25 19:44:41 2013>" ;<br>
:GMT_version = "4.5.9_r9909 [64-bit]" ;<br>
}<br>
=====================================Output from mbgrid after
gdal_translate
====================================================<br>
ncdump -h MCBCP_BeachModel_UTM11N_Baseline.grdclip.flipped.nc<br>
netcdf MCBCP_BeachModel_UTM11N_Baseline.grdclip.flipped {<br>
dimensions:<br>
side = 2 ;<br>
xysize = 154681126 ;<br>
variables:<br>
double x_range(side) ;<br>
x_range:units = "meters" ;<br>
double y_range(side) ;<br>
y_range:units = "meters" ;<br>
double z_range(side) ;<br>
z_range:units = "meters" ;<br>
double spacing(side) ;<br>
int dimension(side) ;<br>
float z(xysize) ;<br>
z:scale_factor = 1. ;<br>
z:add_offset = 0. ;<br>
z:node_offset = 1 ;<br>
<br>
// global attributes:<br>
:title = "" ;<br>
:source = "" ;<br>
}</font><br>
<br>
<pre class="moz-signature" cols="72">--
John Helly, University of California, San Diego / San Diego Supercomputer Center / Scripps Institution of Oceanography / 760 840 8660 mobile / stonesteps (Skype) / stonesteps7 (iChat) / <a class="moz-txt-link-freetext" href="http://www.sdsc.edu/~hellyj">http://www.sdsc.edu/~hellyj</a></pre>
</body>
</html>