open TELEMAC-MASCARET The mathematically superior suite of solvers

Hello,
I want to run ocean model to figure out tidal currents.
Firstly, I make the domain which coordinate system is WGS84. The whole domain is
consisted of three(3) UTM zone.
My domain is quite big, so I want to consider 'tide generating force' and Coriolis force.
The commands can be run only in spherical coordinate system.
So, I set as follows.

'GEOGRAPHIC SYSTEM = 1 (WGS84 LONGITUDE/LATITUDE IN REAL DEGREES)'
'SPHERICAL COORDINATES = true'
'SPATIAL PROJECTION TYPE = 3(LATITUDE LONGITUDE)'
'CORIOLIS = true'
'TIDE GENERATING FORCE = true'

Is it right way to set up the model for running spherical coordinate?

And in reference to the telemac2d manual, basically, the spherical coordinate will be changed to Cartesian coordinate using 'LATITU.f (subroutine)'.

When using this subroutine file, 'LATITU.f', what coordinate system would change?


Thank you.

MJ Lee

Hello,

Yes, your settings seem to be right.

If you provide your mesh in real degrees (SPATIAL PROJECTION TYPE = 3), the code automatically projects your grid in a mercator projection for telemac (done in read_mesh_coord.f, and then becomes a SPATIAL PROJECTION TYPE = 2 in fact).
The latitu.f subroutine only computes arrays depending on the latitude for each grid point (sinus/cosinus).

I hope it helps.

Regards,
Laurent