[MITgcm-support] to prepare horizGridFile for the curvilinear coordinates

[Yi-Chih Huang]

Martin,

    I got the following error messages "Transport endpoint is not
connected" in the beginning of execution.  I searched this error online,
but ended up nothing.  Could you indicate what this error about?

/bin/bash: /cvmfs/
soft.computecanada.ca/nix/var/nix/profiles/16.09/lmod/lmod/init/bash:
Transport endpoint is not connected
slurmstepd: error: execve():
/project/6001267/MITgcm/MITgcm_c66k/verification/rotating_tank2/runM/./mitgcmuv:
Transport endpoint is not connected
......................................................................................................................................

    My topic is about atmospheric convection so that I don't need any sea
ice or ocean packages.  I have successfully run the rotating_tank
example in spherical coordinates.  To extend to the polar region, I need to
apply the curvilinear coordinates.  At first, I thought if only I prepare
the horizGridFile file, I should be able to run this model in the
curvilinear coordinates.  If to reproduce the example in MITgcm_contrib is
the best way to learn the curvilinear coordinates, could you recommend an
atmospheric case?  I don't want finally I have to remove many sea ice and
ocean packages.

    Thanks much,

                    Yi-Chih

######################################################################################
Hi Yi-Chih,

now I am not quite sure what you are trying to do. Once you have specified
the curvilinear grid in one way (grid files) or the other (OLD_GRID_IO
method), you do not have to specify delX/delY in data. In fact, it does not
make sense to specify spherical coordinate deltaX/Y for a curvilinear grid,
does it? The best that can happen, is that these parameters are ignored,
and I think that is what you are seeing. Have look here:

http://wwwcvs.mitgcm.org/viewvc/MITgcm/MITgcm_contrib/arctic/

for examples of a curvilinear grid across the Arctic using the OLD_GRID_IO
method. The grid files contain longitude and latitude and there?s no
ambiguity about the grid.

Another option is to rotate a spherical grid so that the converging
longitudes are moved out of your region of interest. This has been done
before and there are parameters for this (from PARAMS.h:
C     rotateGrid      :: rotate grid coordinates to geographical coordinates
C                        according to Euler angles phiEuler, thetaEuler,
psiEuler
)

Only on a spherical grid (or cartesian grid), the parameters ygOrigin,
delx/dely have a meaning. f0 can only be used with a cartesian grid. It
does not make too much sense to have an f- or beta-plane in a non-cartesian
grid (as far a I know).

Martin
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20181004/6d40c71c/attachment.html>