Changeset 3159 for trunk/LMDZ.COMMON/libf/evolution/pemetat0.F90

Timestamp:

Dec 14, 2023, 6:56:40 PM (12 months ago)

Author:

jbclement

Message:

PEM:

• A number of corrections related to r3149 for CO2 ice management:

  'co2_ice' was not transferred to 'perennial_co2ice' at the end of the PEM run;
  In the Mars PCM, 'perennial_co2ice' is now correctly handled regarding the frost in case of CO2 sublimation/condensation.
  

• Addition of (CO2 and H2O) ice metamorphism: if the frost is above a given value, then the excess frost is transferred to the perennial ice.
• Thresholds value for ice management can now be set in the "run_PEM.def".

JBC

File:

• trunk/LMDZ.COMMON/libf/evolution/pemetat0.F90 (modified) (2 diffs)

trunk/LMDZ.COMMON/libf/evolution/pemetat0.F90

@@ -82 +82 @@
 #ifdef CPP_STD
     logical, dimension(ngrid) :: watercaptag
-    watercaptag(:) = .false.
+    watercaptag = .false.
 #endif
 
@@ -152 +152 @@
                             TI_PEM(ig,iloop,islope) = TI_breccia
                         enddo
-                    else ! we keep the high ti values
+                    else ! we keep the high TI values
                         do iloop = index_breccia + 1,index_bedrock
                             TI_PEM(ig,iloop,islope) = TI_PEM(ig,index_breccia,islope)