[MITgcm-support] "blanklist" parameter in MITgcm

[Matthew Mazloff]

This may be helpful too

https://github.com/MITgcm/MITgcm/pull/382

Matt


> On Nov 8, 2022, at 7:54 AM, Daniel Goldberg <dan.goldberg at ed.ac.uk> wrote:
> 
> Thanks J-M. I did not know i needed to include the exch2 package. 
> 
> best
> dan
> 
> On Tue, Nov 8, 2022 at 3:46 PM Jean-Michel Campin <jmc at mit.edu> wrote:
> This email was sent to you by someone outside the University.
> You should only click on links or attachments if you are certain that the email is genuine and the content is safe.
> 
> Hi Dan,
> 
> To use blank-tiles, you need:
> 1) to compile pkg/exch2 (will be used if compiled, no run-time switch)
> 2) provide a list of blank-tile in "data.exch2".
>  the example I am aware of arei, in verification/:
>    adjustment.cs-32x32x1/input/data.exch2.mpi  (cs-grid, 4 blank-tiles)
>    global_ocean.90x40x15/input/data.exch2.mpi  (lat-lon grid, 1 blank-tile)
>  these 2 are intented to be copied to "data.exch2"
>  And note that you will have to set "dimsFacets" to the size of the full domain
>  (including the blank-tiles) like in the later example: dimsFacets = 90, 40,
> 3) to find the list of empty tiles to specify, the simplest is to run 1 time-step
>  with an executable compiled without blank-tile (but same sNx,sNy) and to surch for
>   "Empty tile: #" in all STDOUT ; those numbers will be the ones to list in "blankList".
> 
> Cheers,
> Jean-Michel
> 
> On Tue, Nov 08, 2022 at 03:15:33PM +0000, Daniel Goldberg wrote:
> > Hello All
> >
> > Can I ask about how to use the blankList parameter in data.exch2?
> >
> > Can it be used for cartesian as well as cubed sphere grids? If so, Im unsure of how exactly to use it without examples shown.
> >
> > the docs on the CS grid say:
> >
> >
> >   *   The tiles omitted are specified in the file blanklist.txt by their tile number in the topology, separated by a newline.
> >   *   The total number of tiles in the topology minus those listed in blanklist.txt must equal nSx*nPx.  (which i assume extends to the y direction -- but possibly not.)
> >
> > So, if i wanted to have a grid where nSx=nSy=1, and the grid structure
> >
> > nPx = 2
> > nPy = 4
> >
> > then i believe with a cartesian grid, the tile numbering would be
> >
> > 3 7
> > 2 6
> > 1 5
> > 0 4
> >
> > where left = west, up = north, and i wanted to make tiles 2 and 3 blank. Would i set
> >
> > blanklist = 2,3
> >
> > and then change nPy = 3? And would i run with e.g. "mpirun -n 6" ? And, importantly, could input files still be defined as arrays of dimensions (sNx*2, sNy*4), i.e. of the same size as if i were not blanking tiles?
> >
> > Many thanks!!!
> >
> > Dan
> >
> >
> >
> > --
> > Daniel Goldberg, PhD (he/him)
> > Reader in Glaciology
> > School of Geosciences, University of Edinburgh
> > Geography Building, Drummond Street, Edinburgh EH8 9XP
> >
> >
> > em: dan.goldberg at ed.ac.uk<mailto:dan.goldberg at ed.ac.uk>
> > web: https://www.geos.ed.ac.uk/homes/dgoldber
> > The University of Edinburgh is a charitable body, registered in Scotland, with registration number SC005336. Is e buidheann carthannais a th??? ann an Oilthigh Dhùn Èideann, clàraichte an Alba, àireamh clàraidh SC005336.
> 
> > _______________________________________________
> > MITgcm-support mailing list
> > MITgcm-support at mitgcm.org
> > http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support
> 
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> https://urldefense.com/v3/__http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support__;!!Mih3wA!E4-fwGJEfQtsA4B5eOhcwtEH3nR90Dx0uAt0TddrTNRxQj8zfGaI39hFxdEFckoJr82lqv3l8NMrWQ9YeUrWvr6fWg$