[MITgcm-support] internal-wave tutorial: changing grid sizes

[Jeremy Miller]

Dear MITgcm forum,

Hope all are well. I have a couple of questions about the internal wave
tutorial.

*1.* I am trying to change Nr, the number of cells in the z-direction. When
changing it from 20 (the default) foe example to 22, at run-time I get the
error:

Fortran runtime error: Cannot match namelist object name -0.0440182

Error termination. Backtrace:
#0  0x7f0d178d1d3a
#1  0x7f0d178d2849

etc.

I changed Nr in input/data in the relevant lines:

Line 8: Tref =  0.0484200, ..... -0.0484200 (20 values),  -> Tref =
0.0486517,..., -0.0486517, (22 values)
Line 14: sRef= 20*35. ->  sRef= 22*35.,
Line 70: delZ=22*10., -> delZ=22*9.,

and I changed the relevant lines in code/SIZE.h and code/SIZE.h_mpi:

Line 55:      &           Nr  =  20) ->      &           Nr  =  22)

*2. *Changing Nx (the number of slices in the x direction) from the default
60 to, say 100, in data and SIZE.h as above, the code runs, but when trying
to read the temperature data back in from the file
'T.0000000100.002.001.meta', then resize the array in python with the
snippet (with Nz=20, Nx=100)

T = mds.rdmds('T',Nt);
T_vals = T.reshape(Nz,Nx)

I get the error

ValueError: cannot reshape array of size 1200 into shape (20,100)

i.e. the results data has not been reshaped according to the new Nx=100.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20210419/367f7c70/attachment.html>