Index: /LMDZ6/trunk/libf/phylmd/conf_phys_m.f90
===================================================================
--- /LMDZ6/trunk/libf/phylmd/conf_phys_m.f90	(revision 6178)
+++ /LMDZ6/trunk/libf/phylmd/conf_phys_m.f90	(revision 6179)
@@ -535,5 +535,7 @@
     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)
 
Index: /LMDZ6/trunk/libf/phylmd/limit_read_mod.f90
===================================================================
--- /LMDZ6/trunk/libf/phylmd/limit_read_mod.f90	(revision 6178)
+++ /LMDZ6/trunk/libf/phylmd/limit_read_mod.f90	(revision 6179)
@@ -259,5 +259,5 @@
 
        !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,5 +344,5 @@
              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,5 +358,5 @@
           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,4 +414,7 @@
 ! 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,5 +512,5 @@
           CALL Scatter(sst_glo,sst)
           !GG
-          IF (iflag_seaice==1) THEN
+          IF (iflag_seaice==2) THEN
              CALL Scatter(sih_glo,sih)
           END IF