<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=us-ascii">
<meta name="Generator" content="Microsoft Word 15 (filtered medium)">
<style><!--
/* Font Definitions */
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:"Calibri Light";
panose-1:2 15 3 2 2 2 4 3 2 4;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
margin-bottom:.0001pt;
font-size:11.0pt;
font-family:"Calibri",sans-serif;
mso-fareast-language:EN-US;}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:#0563C1;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:#954F72;
text-decoration:underline;}
p.MsoPlainText, li.MsoPlainText, div.MsoPlainText
{mso-style-priority:99;
mso-style-link:"Plain Text Char";
margin:0cm;
margin-bottom:.0001pt;
font-size:11.0pt;
font-family:"Calibri",sans-serif;
mso-fareast-language:EN-US;}
span.EmailStyle17
{mso-style-type:personal-compose;
font-family:"Calibri",sans-serif;
color:windowtext;}
span.PlainTextChar
{mso-style-name:"Plain Text Char";
mso-style-priority:99;
mso-style-link:"Plain Text";
font-family:"Calibri",sans-serif;}
.MsoChpDefault
{mso-style-type:export-only;
font-family:"Calibri",sans-serif;
mso-fareast-language:EN-US;}
@page WordSection1
{size:612.0pt 792.0pt;
margin:72.0pt 72.0pt 72.0pt 72.0pt;}
div.WordSection1
{page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
</head>
<body lang="EN-GB" link="#0563C1" vlink="#954F72">
<div class="WordSection1">
<p class="MsoPlainText">Dear MITgcm Support Team,<o:p></o:p></p>
<p class="MsoPlainText"><o:p> </o:p></p>
<p class="MsoPlainText">I am trying to create a cube-sphere (CS96) MITgcm Ocean model setup which I plan to force with climatology derived from an Atmosphere-only GCM in a way similar to ECCO (but without the adjoint package and CS96 rather than LLC90 grid).
I believe (well, hope) that I can create correctly the appropriate bathymetry and forcing binary files.
<o:p></o:p></p>
<p class="MsoPlainText"><o:p> </o:p></p>
<p class="MsoPlainText">My problem lies in the initialisation of the model, in particular regarding the grid setup. I have searched the manual and verification experiments for code but I am still unsure how to create the grid_cs96_face00[1-6].bin files. Fortunately,
I found a set within <a href="http://mit.ecco-group.org/datasets/ecco_for_las/version_4/grids/grids_input/cs_nostretch/cube_96/">
http://mit.ecco-group.org/datasets/ecco_for_las/version_4/grids/grids_input/cs_nostretch/cube_96/</a>, which appear to be the correct size (8bytes*(96+1)*(96+1)*16records = 1204352 bytes). I have renamed the files to grid_cs96.face001.bin etc. and set up my
input/data to see read them. Additionally, I can read these grid_cs96.face00[1-6].bin binaries using Matlab GCMfaces and can create a mygrid structure (which I have used to create the bathymetry.bin). I believe these binary grid files are what I need. I have
within input/data<o:p></o:p></p>
<p class="MsoPlainText"><o:p> </o:p></p>
<p class="MsoPlainText"># Gridding parameters<o:p></o:p></p>
<p class="MsoPlainText">&PARM04<o:p></o:p></p>
<p class="MsoPlainText">horizGridFile='grid_cs96',<o:p></o:p></p>
<p class="MsoPlainText">usingCurvilinearGrid=.TRUE.,<o:p></o:p></p>
<p class="MsoPlainText"><o:p> </o:p></p>
<p class="MsoPlainText">On running, the model fails with the following error<o:p></o:p></p>
<p class="MsoPlainText">forrtl: severe (36): attempt to access non-existent record, unit 9 file
<i>[my_full_path]/</i>build/grid_cs96.face001.bin <o:p></o:p></p>
<p class="MsoPlainText"><o:p> </o:p></p>
<p class="MsoPlainText">I have searched MITgcm-support and it seems the model is trying to access records that are not within the grid_cs96.face00[1-6].bin files – which suggests an incorrect SIZE.h file. I have tried numerous configurations of SIZE.h. A simple
case I have for SIZE.h which (I believe) sets up a CS topology with 24 tiles running on 6 processors.<o:p></o:p></p>
<p class="MsoPlainText"><o:p> </o:p></p>
<p class="MsoPlainText"> PARAMETER (<o:p></o:p></p>
<p class="MsoPlainText"> & sNx = 48,<o:p></o:p></p>
<p class="MsoPlainText"> & sNy = 48,<o:p></o:p></p>
<p class="MsoPlainText"> & OLx = 4,<o:p></o:p></p>
<p class="MsoPlainText"> & OLy = 4,<o:p></o:p></p>
<p class="MsoPlainText"> & nSx = 4,<o:p></o:p></p>
<p class="MsoPlainText"> & nSy = 1,<o:p></o:p></p>
<p class="MsoPlainText"> & nPx = 6,<o:p></o:p></p>
<p class="MsoPlainText"> & nPy = 1,<o:p></o:p></p>
<p class="MsoPlainText"> & Nx = sNx*nSx*nPx, <o:p></o:p></p>
<p class="MsoPlainText"> & Ny = sNy*nSy*nPy, <o:p></o:p></p>
<p class="MsoPlainText"> & Nr = 50 )<o:p></o:p></p>
<p class="MsoPlainText"><o:p> </o:p></p>
<p class="MsoPlainText">Looking at the cs32 verification experiments which I can run, the experiment uses grid_cs32 .face00[1-6].bin files that are 156816bytes in size (e.g. within global_oce_cs32/input_fields/) whereas I would expect them to be 8bytes*(32+1)*(32+1)*16records=139392
bytes. Hence I think that I am misunderstanding something. I believe that my compiler optfile is correct (the verification experiments work).<o:p></o:p></p>
<p class="MsoPlainText"><o:p> </o:p></p>
<p class="MsoNormal">Kind regards,<o:p></o:p></p>
<p class="MsoNormal">Stephen Hunter<o:p></o:p></p>
<p class="MsoNormal"><span style="font-size:10.5pt;font-family:"Calibri Light",sans-serif"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:10.5pt;mso-fareast-language:EN-GB">School of Earth & Environment<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:10.5pt;mso-fareast-language:EN-GB">University of Leeds<o:p></o:p></span></p>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
</body>
</html>