Changeset 3996

Timestamp:

Jan 12, 2026, 9:49:50 AM (16 hours ago)

Author:

jbclement

Message:

Mars PCM:
Deletion of unused legacy code in "initracer.F90".
JBC

Location:

trunk/LMDZ.MARS

Files:

• changelog.txt (modified) (1 diff)
• libf/phymars/dyn1d/init_testphys1d_mod.F90 (modified) (1 diff)
• libf/phymars/initracer.F90 (modified) (3 diffs)
• libf/phymars/physiq_mod.F (modified) (1 diff)
• libf/phymars/planete_h.F90 (modified) (1 diff)
• libf/phymars/surfdat_h.F90 (modified) (1 diff)

trunk/LMDZ.MARS/changelog.txt

@@ -5006 +5006 @@
 == 21/11/2025 == JBC
 Simplification of the orbital parameters initialization (resolve ticket #109): removing redundant actions in subroutine 'iniorbit' in regard of what is done in "tabfi.F" and moving the subroutine 'lsp2solp' from "tabfi.F" to "planete_h.F90" + some cleanings.
+
+== 12/01/2026 == JBC
+Deletion of unused legacy code in "initracer.F90".

trunk/LMDZ.MARS/libf/phymars/dyn1d/init_testphys1d_mod.F90

@@ -431 +431 @@
 call ini_fillgeom(1,latitude,longitude,(/1.0/))
 call conf_phys(1,llm,nq)
-call initracer(ngrid,nq,qsurf)
+call initracer(ngrid,nq)
 
 ! In 1D model physics are called every time step

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

@@ -5 +5 @@
 contains
 
-SUBROUTINE initracer(ngrid,nq,qsurf)
+SUBROUTINE initracer(ngrid,nq)
 
 use tracer_mod
@@ -40 +40 @@
 integer, intent(in) :: ngrid ! number of atmospheric columns
 integer, intent(in) :: nq    ! number of tracers
-real, intent(out) :: qsurf(ngrid,nq) ! tracer on surface (e.g.  kg.m-2)
 
 integer :: iq,ig,count
@@ -548 +547 @@
 endif
 
-! if water cycle but iceparty=.false., there will nevertheless be
-! water ice at the surface (iceparty is not used anymore, but this
-! part is still relevant, as we want to stay compatible with the
-! older versions).
-if (water.and.(igcm_h2o_ice == 0)) then
-    igcm_h2o_ice=igcm_h2o_vap ! so that qsurf(i_h2o_ice) is identified
-                              ! even though there is no q(i_h2o_ice)
-else
-    ! surface ice qsurf(i_h2o_ice) was loaded twice by phyetat0,
-    ! as qsurf(i_h2o_vap) & as qsurf(i_h2o_ice), so to be clean:
-    if (igcm_h2o_vap /= 0) then
-        qsurf(1:ngrid,igcm_h2o_vap)=0
-    endif
-endif
-
 ! Additional test required for HDO
 ! We need to compute some things for H2O before HDO

trunk/LMDZ.MARS/libf/phymars/physiq_mod.F

@@ -711 +711 @@
 c        initialize tracers
 c        ~~~~~~~~~~~~~~~~~~
-         call initracer(ngrid,nq,qsurf)
+         call initracer(ngrid,nq)
 
 c        Initialize albedo and orbital calculation

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

@@ -17 +17 @@
 REAL,SAVE :: lsperi    ! Solar longitude of the perihelion, angle in rad
 REAL,SAVE :: e_elips   ! Orbit eccentricity
-REAL,SAVE :: p_elips   ! Ellipse semi-latus rectum
+REAL,SAVE :: p_elips   ! Ellipse parameter (semi-latus rectum)
 !$OMP THREADPRIVATE(lsperi,e_elips,p_elips)
 REAL,PARAMETER :: unitastr=149.597927 ! Astronomical unit AU, in Mkm

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

@@ -62 +62 @@
 
   implicit none
-  integer,intent(in) :: ngrid ! number of atmospheric columns
-  integer,intent(in) :: nq ! number of tracers
-  integer,intent(in) :: nslope ! number of sub-grid scale slope
+  integer,intent(in) :: ngrid  ! number of atmospheric columns
+  integer,intent(in) :: nq     ! number of tracers
+  integer,intent(in) :: nslope ! number of sub-grid scale slopes
 
     allocate(albedodat(ngrid))