Changeset 3328 for trunk/LMDZ.COMMON/libf/evolution

Timestamp:

May 15, 2024, 6:51:20 PM (6 months ago)

Author:

llange

Message:

PEM

• Remove a wrong line which modified the initial surface sublimating during

the main time loop

• Cleaning the previous cleaning commit (rage pas JB)

LL

Location:

trunk/LMDZ.COMMON/libf/evolution

Files:

• changelog.txt (modified) (1 diff)
• ice_table_mod.F90 (modified) (1 diff)
• pem.F90 (modified) (6 diffs)

trunk/LMDZ.COMMON/libf/evolution/changelog.txt

@@ -307 +307 @@
 - Correction of the computation for the change of sublimating ice necessary to the stopping criteria (the mistake was introduced in r3149) + simplification of the algorithm + renaming of variables with more explicit names;
 - Cleaning of "soil_thermalproperties_mod.F90".
+
+== 15/05/2024 == LL
+- Remove a wrong line which modified the  initial surface sublimating during
+the main time loop
+- Cleaning the previous cleaning commit.
+

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

@@ -178 +178 @@
 do ig = 1,ngrid
     do islope = 1,nslope
-        print*, ig, islope, ice_table_depth(ig,islope)
         call compute_Tice_pem(nsoil_PEM,tsoil(ig,:,islope),tsurf(ig,islope),ice_table_depth(ig,islope),Tice(ig,islope))
         rho(ig,islope) = -3.5353e-4*Tice(ig,islope)**2 + 0.0351* Tice(ig,islope) + 933.5030 ! Rottgers, 2012

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

@@ -223 +223 @@
 real                                  :: totmassh2o_adsorbded               ! Total mass of H2O that is exchanged because of adsorption / desoprtion over the planets [kg]
 logical                               :: bool_sublim                        ! logical to check if there is sublimation or not
+logical, dimension(:,:),  allocatable :: co2ice_has_disappeared             ! logical to check if a co2ice reservoir has already disappeared at a previous timestep
 real, dimension(:,:),     allocatable :: porefillingice_thickness_prev_iter ! ngrid x nslope: Thickness of the ice table at the previous iteration [m]
 real, dimension(:),       allocatable :: delta_h2o_icetablesublim           ! ngrid x Total mass of the H2O that has sublimated / condenses from the ice table [kg]
@@ -496 +497 @@
 allocate(watersoil_density_PEM_timeseries(ngrid,nsoilmx_PEM,nslope,timelen))
 allocate(delta_co2_adsorbded(ngrid))
+allocate(co2ice_has_disappeared(ngrid,nslope))
 allocate(porefillingice_thickness_prev_iter(ngrid,nslope))
 allocate(delta_h2o_icetablesublim(ngrid))
@@ -564 +566 @@
 ini_co2ice_sublim = .false.
 ini_h2oice_sublim = .false.
+co2ice_has_disappeared(:,:) = .false.
 do i = 1,ngrid
     Total_surface = Total_surface + cell_area(i)
@@ -789 +792 @@
     do ig = 1,ngrid
         do islope = 1,nslope
-            if (co2_ice_ini(ig,islope) > 0. .and. co2_ice(ig,islope) < 1.e-10) then ! co2 ice disappeared, look for closest point without co2ice
+            if (co2_ice_ini(ig,islope) > 0. .and. co2_ice(ig,islope) < 1.e-10 .and. (.not.(co2ice_has_disappeared(ig,islope)))) then ! co2 ice disappeared, look for closest point without co2ice
+                co2ice_has_disappeared(ig,islope) = .true.
                 if (latitude_deg(ig) > 0) then
                     do ig_loop = ig,ngrid
@@ -813 +817 @@
                     enddo
                 endif
-                co2_ice_ini(ig,islope) = 0.
                 if ((co2_ice(ig,islope) < 1.e-10) .and. (h2o_ice(ig,islope) > frost_albedo_threshold)) then
                     albedo(ig,1,islope) = albedo_h2o_frost
@@ -832 +835 @@
                 enddo
                 tsurf_ave(ig,islope) = ave/timelen
-                ! set the surface albedo to be the ice albedo
-!                if (latitude_deg(ig) > 0) then
-!                   icap = 1
-!                else
-!                   icap = 2
-!                endif
-!                albedo(ig,1,islope) = albedice(icap)
-!                albedo(ig,2,islope) = albedice(icap)
-!                emis(ig,islope) = emisice(icap)
             endif
         enddo