Changeset 4170 for trunk/LMDZ.COMMON/libf/evolution/planet.F90

Timestamp:

Apr 3, 2026, 4:34:51 PM (10 days ago)

Author:

jbclement

Message:

PEM:

• Deletion of 'flux_ssice' ('zqdsdif_tot') from the "startfi.nc" as it is not needed.
• Using the yearly average flux for the sublimating subsurface ice instead of the value at last timestep.
• Keeping a clear separation between the subsurface ice flux and the surface ice tendency.
• Making sure that subsurface ice depth is well given to the PCM at the end of the PEM (ice table VS layering).

JBC

File:

• trunk/LMDZ.COMMON/libf/evolution/planet.F90 (modified) (5 diffs)

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

@@ -38 +38 @@
 
 ! Ice-related:
-real(dp),    dimension(:,:), allocatable :: h2o_ice                    ! H2O ice [kg.m-2]
-real(dp),    dimension(:,:), allocatable :: co2_ice                    ! CO2 ice [kg.m-2]
+real(dp),    dimension(:,:), allocatable :: h2o_ice                    ! H2O ice [kg/m2]
+real(dp),    dimension(:,:), allocatable :: co2_ice                    ! CO2 ice [kg/m2]
 real(dp)                                 :: h2oice_sublim_coverage_ini ! Initial surface area of sublimating H2O ice [m2]
 real(dp)                                 :: co2ice_sublim_coverage_ini ! Initial surface area of sublimating CO2 ice [m2]
@@ -73 +73 @@
 real(dp), dimension(:,:),   allocatable :: icetable_depth_old ! Old depth of the ice table [m]
 real(dp), dimension(:),     allocatable :: delta_icetable     ! Total mass of the H2O exchanged with the ice table [kg]
+real(dp), dimension(:,:),   allocatable :: flux_ssice_avg     ! Average of total flux exchanged with subsurface ice [kg/m2/y]
 
 ! Tracer-related:
@@ -153 +154 @@
 allocate(q_h2o_ts(ngrid,nday))
 allocate(q_co2_ts(ngrid,nday))
+allocate(flux_ssice_avg(ngrid,nslope))
 
 ps_avg(:) = 0._dp
@@ -163 +165 @@
 q_h2o_ts(:,:) = 0._dp
 q_co2_ts(:,:) = 0._dp
+flux_ssice_avg(:,:) = 0._dp
 
 END SUBROUTINE allocate_xios_state
@@ -422 +425 @@
 if (allocated(d_h2oice)) deallocate(d_h2oice)
 if (allocated(layerings_map)) deallocate(layerings_map)
+if (allocated(flux_ssice_avg)) deallocate(flux_ssice_avg)
 
 END SUBROUTINE end_planet