Changeset 3130 for trunk/LMDZ.MARS/libf/phymars/phyredem.F90

Timestamp:

Nov 20, 2023, 1:25:31 PM (14 months ago)

Author:

jbclement

Message:

PEM:
The perennial co2 ice is now taken into account with co2 frost (qsurf) to compute the tendency and to make the update + Rework of how co2 frost is converted to perennial co2 ice at the end of the PEM run + Correction of the value of 'threshold_co2_frost2perennial' to correspond to 10 m + Perennial co2 ice is now handled outside 'paleoclimate' in "phyetat0_mod.F90" of the Mars PCM + Some cleanings.

/!\ Commit for the PEM management of co2 ice before a rework of ice management in the PEM!
JBC

File:

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

trunk/LMDZ.MARS/libf/phymars/phyredem.F90

@@ -168 +168 @@
                     albedo,emis,q2,qsurf,qsoil,&
                     tauscaling,totcloudfrac,wstar, &
-                    watercap,perenial_co2ice)
+                    watercap,perennial_co2ice)
   ! write time-dependent variable to restart file
   use iostart, only : open_restartphy, close_restartphy, & 
@@ -205 +205 @@
   real,intent(in) :: wstar(ngrid)
   real,intent(in) :: watercap(ngrid,nslope)
-  real,intent(in) :: perenial_co2ice(ngrid,nslope)
+  real,intent(in) :: perennial_co2ice(ngrid,nslope)
 
   integer :: iq
@@ -226 +226 @@
   call put_field("watercap","H2O ice cover",watercap,time)
  
- ! Perenial CO2 ice layer
-  if (paleoclimate) call put_field("perenial_co2ice","CO2 ice cover",perenial_co2ice,time)
+  ! Perennial CO2 ice layer
+  call put_field("perennial_co2ice","CO2 ice cover",perennial_co2ice,time)
 
   ! Surface temperature