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

Jeremy Miller jeremysharonmiller at gmail.com
Mon Apr 19 03:08:02 EDT 2021

Dear MITgcm forum,

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

*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

Fortran runtime error: Cannot match namelist object name -0.0440182

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


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>

More information about the MITgcm-support mailing list