Changeset 3996 for trunk/LMDZ.MARS/libf/phymars/initracer.F90

Timestamp:

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

Author:

jbclement

Message:

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

File:

• trunk/LMDZ.MARS/libf/phymars/initracer.F90 (modified) (3 diffs)

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