<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style type="text/css" style="display:none;"><!-- P {margin-top:0;margin-bottom:0;} --></style>
</head>
<body dir="ltr">
<div id="divtagdefaultwrapper" style="font-size:12pt;color:#000000;font-family:Calibri,Helvetica,sans-serif;" dir="ltr">
<p style="margin-top:0;margin-bottom:0">Hello,</p>
<p style="margin-top:0;margin-bottom:0"><br>
</p>
<p style="margin-top:0;margin-bottom:0">I am working on modelling water flow through a strait using MITgcm. I am new to this and I have some questions setting up the bathymetry file for upload to MITgcm. I would appreciate any help.</p>
<p style="margin-top:0;margin-bottom:0"><br>
</p>
<p style="margin-top:0;margin-bottom:0">To create my bathymetry file, I downloaded the bedrock in a geotiff format from NOAA. The coordinates were in latitude/longitude and I warped it to metres. Then I cropped the geotiff to remove all areas without data and
obtained a rectangle. This I converted to a Numpy array on python. My array is of shape: (points in x, points in y) and the elevation at each point is the value in the matrix. I converted the array to a bin using the following code on Python: <span style="font-size: 12pt;">
</span></p>
<p style="margin-top:0;margin-bottom:0"><span style="font-size: 12pt;">dt='>f8'</span></p>
<p>fid = open("bathymetry.bin, "wb")</p>
<p>bandarray.astype(dt).ravel(order='F').tofile(fid)</p>
<p>fid.close()</p>
<p><br>
</p>
<p>My question is, when I provide info in the SIZE.h file, are the following correct:</p>
<p></p>
<ol style="margin-bottom: 0px; margin-top: 0px;">
<li><span style="font-size: 12pt;"></span><span style="font-size: 12pt;"></span><span style="font-size: 12pt;">sNx*nPx should equal the number of X points on the bathymetry file.</span></li><li><span style="font-size: 12pt;"></span><span style="font-size: 12pt;">sNy*nPy should equal the number of Y points on the bathymetry file. </span></li><li><span style="font-size: 12pt;"></span><span style="font-size: 12pt;">The resolution should equal the spacing between each point in the X/Y and in metres? </span></li><li><span style="font-size: 12pt;"></span><span style="font-size: 12pt;">Then what is Nr? We are not providing a grid for the Z coordinates only values of elevation. </span></li><li><span style="font-size: 12pt;">Also, is my process above correct?</span></li></ol>
<br>
<p></p>
<p>Regards,</p>
<p>Keerthijan</p>
</div>
</body>
</html>