Changeset 2285 for trunk/LMDZ.MARS/libf/phymars

Timestamp:

Apr 14, 2020, 10:41:14 AM (5 years ago)

Author:

adelavois

Message:

MARS GCM:
Minor update of revision 2281 in order to run martian physics with callsoil option when no startfi.nc

File:

• trunk/LMDZ.MARS/libf/phymars/physiq_mod.F (modified) (5 diffs)

trunk/LMDZ.MARS/libf/phymars/physiq_mod.F

@@ -36 +36 @@
      &                      igcm_topdust_number
       use comsoil_h, only: inertiedat, ! soil thermal inertia
-     &                     tsoil, nsoilmx ! number of subsurface layers
+     &                     tsoil, nsoilmx,!number of subsurface layers
+     &                     mlayer,layer ! soil mid layer depths
       use geometry_mod, only: longitude, latitude, cell_area,
      &                        longitude_deg 
@@ -479 +480 @@
       REAL hsummit(ngrid)
 
+c when no startfi file is asked for init
+      real alpha,lay1 ! coefficients for building layers
+      integer iloop
+
 !      LOGICAL startphy_file
 
@@ -517 +522 @@
      &         mem_Mccn_co2,mem_Nccn_co2,
      &         mem_Mh2o_co2,watercap)
+
 
 #else
@@ -541 +547 @@
            ! starting without startfi.nc and with callsoil
            ! is not yet possible as soildepth default is not defined
-           if (callsoil) then
-              call abort_physic("physiq","callsoil option is not", 
-     &                     "yet available without startfi",1)
+           if (callsoil) then                     
+              ! default mlayer distribution, following a power law:
+              !  mlayer(k)=lay1*alpha**(k-1/2)
+              lay1=2.e-4
+              alpha=2
+              do iloop=0,nsoilmx-1
+                 mlayer(iloop)=lay1*(alpha**(iloop-0.5))
+              enddo
+              lay1=sqrt(mlayer(0)*mlayer(1))
+              alpha=mlayer(1)/mlayer(0)
+              do iloop=1,nsoilmx
+                 layer(iloop)=lay1*(alpha**(iloop-1))
+              enddo
            endif
            ! additionnal "academic" initialization of physics
@@ -558 +574 @@
         endif
 #endif
-
          if (pday.ne.day_ini) then
            write(*,*) "PHYSIQ: ERROR: bad synchronization between ",