Ocean Modeling Discussion

You need to flip Gx and Gy when you call pyroms.grid.Gridgen.
hgrd = pyroms.grid.Gridgen(lonp, latp, beta, (Gy, Gx), proj=maproms)
python use the C matrix indexing and ordering convention.

Etienne,
Can you please check that the parameters (theta_b, theta_s, ...) in your gridid file are consistant with the parameters used to generate the vertical grid?
Thanks,
Frederic

Hi lleoex,

What are the horizontal grid dimension?
Can you add print statment in script in order to return the shape of z_u?

Thanks,
Fred