[MITgcm-support] Internal-Wave: Whether I could enable MPI with using only one tile001.mitgrid (LLC grid)

[Jean-Michel Campin]

Hi Chen,

I think you are right, but might not be easy to fix for something
not many people uses.

I would like to suggest to compile pkg/exch2 (which means you will
use EXCH2) and try again with just one grid file: tile001.mitgrid
that contains the full grid (with all fields 61x2).
With default pkg/exch2 parameters, you don't need to provide a 
parameter file "data.exch2" and it should work just like with EXCH1.

Cheers,
Jean-Michel

On Sat, Jun 04, 2016 at 11:48:14AM -0500, ?????? wrote:
> Hi All,
> 
> I am interested in the experiment Internal_wave, and want to play with the
> input grid file. Here is what I have already tried and the problem I met.
> 
> step0: input grid info at /input/data, &parm4, delxvar.bin (Initially
> given)  ==> run serial succeed.
> 
> step1: input grid info at /input/data, &parm4, tile001.mitgrid(Generated)
> ==> run serial succeed.
> 
>      /code/SIZE.h
> 
>      &           sNx =  60,
>      &           sNy =   1,
>      &           OLx =   2,
>      &           OLy =   2,
>      &           nSx =   1,
>      &           nSy =   1,
>      &           nPx =   1,
>      &           nPy =   1,
>      &           Nx  = sNx*nSx*nPx,
>      &           Ny  = sNy*nSy*nPy,
>      &           Nr  =  20)
> 
> step2: input tile001~004.mitgrid(Generated) ==> run serial succeed.
> 
>      /code/SIZE.h
> 
>      &           sNx =   15,
>      &           sNy =    1,
>      &           OLx =    2,
>      &           OLy =    2,
>      &           nSx =    4,
>      &           nSy =    1,
>      &           nPx =    1,
>      &           nPy =    1,
>      &           Nx  = sNx*nSx*nPx,
>      &           Ny  = sNy*nSy*nPy,
>      &           Nr  =  20)
> 
> step3: input tile001~004.mitgrid(Generated) ==> run parallel succeed.
> 
>      /code/SIZE.h
> 
>      &           sNx =   15,
>      &           sNy =    1,
>      &           OLx =    2,
>      &           OLy =    2,
>      &           nSx =    1,
>      &           nSy =    1,
>      &           nPx =    4,
>      &           nPy =    1,
>      &           Nx  = sNx*nSx*nPx,
>      &           Ny  = sNy*nSy*nPy,
>      &           Nr  =  20)
> 
> =====================================
> Here comes the problem. Since eventually I will have a grid with
> (nSx,nSy,nPx,nPy)=(1,1,20,21), or (nSx,nSy,nPx,nPy)=(20,21,1,1). Generate
> 400 tiles grid seems crazy. So as in the experiment I would like to try
> things like :: enable SIZE.h of step2 (or step3) with only one input file
> tile001.mitgrid. Following is something I have already tried:
> 
> 1.  Try in Serial with only one tile001.mitgrid (each 16 fields in this
> grid file have size XC(61,2),YC(61,2), etc )
> /input/data/&parm01  :  add useSingleCpuIO=.TRUE.,  and/or
> useSingleCpuInput=.TRUE., with SIZE.h in step2:
> /code/SIZE.h
> 
>      &           sNx =   15,
>      &           sNy =    1,
>      &           OLx =    2,
>      &           OLy =    2,
>      &           nSx =    4,
>      &           nSy =    1,
>      &           nPx =    1,
>      &           nPy =    1,
>      &           Nx  = sNx*nSx*nPx,
>      &           Ny  = sNy*nSy*nPy,
>      &           Nr  =  20)
> 
> ===============> error: At line 696 of file mdsio_facef_read.f (unit = 9)
> Fortran runtime error: Cannot open file 'tile002.mitgrid': No such file or
> directory
> 
> 
> 
> 2.  Try in Parallel with only one tile001.mitgrid (each 16 fields in this
> grid file have size XC(61,2),YC(61,2), etc )
> /input/data/&parm01  :  add useSingleCpuIO=.TRUE.,  and/or
> useSingleCpuInput=.TRUE., with SIZE.h in step3:
> /code/SIZE.h
> 
>      &           sNx =   15,
>      &           sNy =    1,
>      &           OLx =    2,
>      &           OLy =    2,
>      &           nSx =    1,
>      &           nSy =    1,
>      &           nPx =    4,
>      &           nPy =    1,
>      &           Nx  = sNx*nSx*nPx,
>      &           Ny  = sNy*nSy*nPy,
>      &           Nr  =  20)
> 
> ===============> error: At line 696 of file mdsio_facef_read.f (unit = 9)
> Fortran runtime error: Cannot open file 'tile002.mitgrid': No such file or
> directory
> 
> 
> 3. add globalFiles=.TRUE., to &Parm01. Try both serial and parallel, gives
> the same error.
> 
> 
> 4. compared with the case on
> MIYgcm_contrib/atnguyen/llc_270/aste_270X450X180, of its Size.h and
> /input/data. Didn't notice there is anyone Flag which is related to grid
> input.
> 
> So Could anyone provide me any ideas about how to enable this?
> 
> Thank you in advance! I do appreciate it!
> 
> 
> Chen CHEN
> Graduate Student in ICES&EM
> POB 3.402J
> The University of Texas at Austin
> 
> Phone: +1 512-968-2126
> Email : chen636489 at utexas.edu

> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support