open TELEMAC-MASCARET The mathematically superior suite of solvers

Hi,

In subroutine readgeo3 in the "TelemacV7\v7p0r1\sources\utils\bief\readgeo3.f ", i came across this formulation for
! LATITUDE-LONGITUDE COORDINATES (DEGREES)
! CHANGED INTO MERCATOR PROJECTION.

to me it seems that there is a factor of 2 missing in the lines 146, 147, 159 and 160.

these lines are:
146 LONGIRAD=LONGI0*PS4/90.D0
147 TAN2=TAN(LATI0*PS4/180.D0+PS4)
159 X(I)=R*(X(I)*PS4/90.D0-LONGIRAD)
160 TAN1=TAN(Y(I)*PS4/180.D0+PS4)

whereas i think they should be :

146 LONGIRAD=LONGI0*PS4*2/90.D0
147 TAN2=TAN(LATI0*PS4*2/180.D0+PS4)
159 X(I)=R*(X(I)*PS4*2/90.D0-LONGIRAD)
160 TAN1=TAN(Y(I)*PS4*2/180.D0+PS4)

could it be?
has anyone faced this issue before?

Bests,
Homayoon

Hello,

Yes, this is the bug we mentioned, that is now corrected in the new release 7.1 that is now available. The new release has also a better treatment of spherical equations proposed by HRW, so that they could mesh the world and propagate a wave on the sphere.

With best regards,

Jean-Michel Hervouet