<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>