Changeset 4004 for trunk

Timestamp:

Jan 14, 2026, 8:29:52 PM (10 days ago)

Author:

aslmd

Message:

planete_h module: replace a couple of deprecated functions dtan dsqrt dsin by their modern equivalent tan sqrt sin

File:

• trunk/LMDZ.MARS/libf/phymars/planete_h.F90 (modified) (1 diff)

trunk/LMDZ.MARS/libf/phymars/planete_h.F90

@@ -42 +42 @@
   ! Compute orbit geometrical parameters
   e_elips = (aphelie - periheli)/(aphelie + periheli)
-  zx0 = -2.0*datan(dtan(0.5*lsp/degrad)*dsqrt((1. - e_elips)/(1. + e_elips)))
+  zx0 = -2.0*atan(tan(0.5*lsp/degrad)*sqrt((1. - e_elips)/(1. + e_elips)))
   if (zx0 <= 0.) zx0 = zx0 + 2.*pi
 
   ! Compute sol at perihelion
-  solp = year_day*(1. - (zx0 - e_elips*dsin(zx0))/(2.*pi))
+  solp = year_day*(1. - (zx0 - e_elips*sin(zx0))/(2.*pi))
 
   END SUBROUTINE lsp2solp