Index: /trunk/LMDZ.MARS/changelog.txt
===================================================================
--- /trunk/LMDZ.MARS/changelog.txt	(revision 3207)
+++ /trunk/LMDZ.MARS/changelog.txt	(revision 3208)
@@ -4468,5 +4468,5 @@
 
 == 07/02/2024 == JBC
-- Addition of the possibility to use subslopes parametization in 1D. To do so, put 'nslope=x' in "run.def" where x (1, 3, 5 or 7) is the number of slopes you want to, or modify it in an appropriate "startfi.nc". Then, default subslopes definition and distribution are used by the model. The already existing case of using one slope whose inclination and orientation are defined in "run.def" is still possible.
+- Addition of the possibility to use subslopes parametrization in 1D. To do so, put 'nslope=x' in "run.def" where x (1, 3, 5 or 7) is the number of slopes you want to, or modify it in an appropriate "startfi.nc". Then, default subslopes definition and distribution are used by the model. The already existing case of using one slope whose inclination and orientation are defined in "run.def" is still possible.
 - Some cleanings throughout the code, in particular related to unused variables.
 
@@ -4476,2 +4476,5 @@
 == 08/02/2024 == JBC
 Modification of r3200 for the 1D model: the 'CO2cond_ps' coefficient (now defined in "co2condens_mod.F") is applied not only to update the surface pressure but also all the tracers and the winds from now on. It is done in "physiq_mod.F" just after the call to 'co2condens'.
+
+== 08/02/2024 == JBC
+Following r3203, initialization cases for subslopes parametrization in 1D have been settled: it can be done with 'nslope' in "run.def" or by using an appropriate "startfi.nc". Addition of error/warning messages for the different situations to be user-friendly.
Index: /trunk/LMDZ.MARS/libf/phymars/dyn1d/init_testphys1d_mod.F90
===================================================================
--- /trunk/LMDZ.MARS/libf/phymars/dyn1d/init_testphys1d_mod.F90	(revision 3207)
+++ /trunk/LMDZ.MARS/libf/phymars/dyn1d/init_testphys1d_mod.F90	(revision 3208)
@@ -84,5 +84,5 @@
 ! Local variables
 !=======================================================================
-integer                   :: ierr, iq, j, ilayer, isoil, nlevel, nsoil, flagthermo, flagh2o
+integer                   :: ierr, iq, j, ilayer, isoil, nlevel, nsoil, flagthermo, flagh2o, nslope_tmp
 integer                   :: dayn ! Final date
 real, dimension(nlayer)   :: zlay ! altitude estimee dans les couches (km)
@@ -410,10 +410,10 @@
 ! Initialize local slope parameters
 ! ---------------------------------
-nslope = 1 ! default: one slope
+nslope_tmp = nslope
+write(*,*)'Number of slopes?' ! Number of subgrid scale slopes
+call getin("nslope",nslope_tmp)
+write(*,*) " nslope = ", nslope_tmp
 if (.not. therestartfi) then
-    ! Number of subgrid scale slopes
-    write(*,*)'Number of slopes?'
-    call getin("nslope",nslope)
-    write(*,*) " nslope = ", nslope
+    nslope = nslope_tmp
     call end_comslope_h
     call ini_comslope_h(1,nslope)
@@ -479,9 +479,10 @@
     call end_paleoclimate_h
     call ini_paleoclimate_h(1,nslope)
-
-    ! Fraction of subslopes in mesh
-    !write(*,*)'What subslope distribution do you want to use?'
-    !call getin("slope_distribution",subslope_dist)
-    !write(*,*) " subslope_dist = ", subslope_dist(1,:)
+else
+    if (nslope /= nslope_tmp) then
+        error stop 'Number of slopes not compatible between "startfi.nc" and "run.def".'
+    else
+        write(*,*) 'Be careful: subslopes distribution and definition are taken from "startfi.nc" and not initialized by default!'
+    endif
 endif