Changeset 3208 for trunk/LMDZ.MARS/libf

Timestamp:

Feb 8, 2024, 4:32:31 PM (11 months ago)

Author:

jbclement

Message:

Mars PCM:
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.
JBC

File:

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

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

@@ -84 +84 @@
 ! 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 +410 @@
 ! 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 +479 @@
     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