Changeset 6179

Timestamp:

Apr 20, 2026, 9:52:01 AM (2 weeks ago)

Author:

evignon

Message:

petits changements d'initialisation des leads pour les nouvelles params

Location:

LMDZ6/trunk/libf/phylmd

Files:

• conf_phys_m.f90 (modified) (1 diff)
• limit_read_mod.f90 (modified) (5 diffs)

LMDZ6/trunk/libf/phylmd/conf_phys_m.f90

@@ -535 +535 @@
     CALL getin_p('amax_n', amax_n)
 
-    amax_s = 0.95
+    !amax_s = 0.95
+    ! by default, set amax_s (south) = amax_n (north) for consistency in LAM simulations
+    amax_s = amax_n
     CALL getin_p('amax_s', amax_s)
 

LMDZ6/trunk/libf/phylmd/limit_read_mod.f90

@@ -259 +259 @@
 
        !GG
-       IF (iflag_seaice==1) THEN
+       IF (iflag_seaice==2) THEN
              ALLOCATE(sih(klon), stat=ierr)
              IF (ierr /= 0) CALL abort_physic(modname, 'PB in allocating sih',1)
@@ -344 +344 @@
              IF (is_omp_master) CALL xios_recv_field("sst_limin",sst_mpi)
              !GG
-             IF (iflag_seaice == 1) THEN
+             IF (iflag_seaice == 2) THEN
                IF (is_omp_master) CALL xios_recv_field("sih_limin",sih_mpi)
              ENDIF
@@ -358 +358 @@
           CALL Scatter_omp(sst_mpi,sst)
           !GG
-          IF (iflag_seaice==1) THEN
+          IF (iflag_seaice==2) THEN
             CALL Scatter_omp(sih_mpi,sih)
           ENDIF
@@ -414 +414 @@
 ! GG
 ! Account for leads
+! pay attention, amax_n and amax_s stands for sea ice cover threshold values different
+! for the northern and southern hemispheres.
+! This is misleading for LAM simulations
              IF (iflag_seaice>0) THEN
                DO ii=1,klon_glo/2
@@ -509 +512 @@
           CALL Scatter(sst_glo,sst)
           !GG
-          IF (iflag_seaice==1) THEN
+          IF (iflag_seaice==2) THEN
              CALL Scatter(sih_glo,sih)
           END IF