[MITgcm-support] "blanklist" parameter in MITgcm

[Daniel Goldberg]

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
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20221108/588159de/attachment-0001.html>