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

Timestamp:

Apr 3, 2026, 4:34:51 PM (8 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/clim_state_init.F90 (modified) (7 diffs)

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

@@ -205 +205 @@
 use surf_ice,  only: set_is_h2o_perice_PCM, set_co2_perice_PCM
 use soil,      only: set_TI_PCM, set_inertiedat_PCM
+use ice_table, only: set_coef_ssdif_PCM
 use display,   only: print_msg, LVL_NFO
 
@@ -276 +277 @@
 call set_tsoil_PCM(tmp3d)
 
+call get_var_nc('coef_ssdif',tmp2d)
+call set_coef_ssdif_PCM(tmp2d)
+
 deallocate(tmp2d); allocate(tmp2d(ngrid,nsoil_PCM))
 call get_var_nc('inertiedat',tmp2d)
@@ -282 +286 @@
 call get_var_nc('inertiesoil',tmp3d)
 call set_TI_PCM(tmp3d)
-
-! To do?
-!   h2oice_depth
-!   d_coef
-!   lag_co2_ice
 
 ! Close/Deallocate
@@ -398 +397 @@
                     !!! transition
                     delta = depth_breccia
-                    TI(i,index_breccia + 1,islope) = sqrt((layer(index_breccia + 1) - layer(index_breccia))/                       &
-                                                               (((delta - layer(index_breccia))/(TI(i,index_breccia,islope)**2)) + &
-                                                               ((layer(index_breccia + 1) - delta)/(TI_breccia**2))))
+                    TI(i,index_breccia + 1,islope) = sqrt((layer(index_breccia + 1) - layer(index_breccia))/                  &
+                                                          (((delta - layer(index_breccia))/(TI(i,index_breccia,islope)**2)) + &
+                                                          ((layer(index_breccia + 1) - delta)/(TI_breccia**2))))
                     do k = index_breccia + 2,index_bedrock
                         TI(i,k,islope) = TI_breccia
@@ -411 +410 @@
                 !! transition
                 delta = depth_bedrock
-                TI(i,index_bedrock + 1,islope) = sqrt((layer(index_bedrock + 1) - layer(index_bedrock))/                   &
-                                                       (((delta - layer(index_bedrock))/(TI(i,index_bedrock,islope)**2)) + &
-                                                       ((layer(index_bedrock + 1) - delta)/(TI_breccia**2))))
+                TI(i,index_bedrock + 1,islope) = sqrt((layer(index_bedrock + 1) - layer(index_bedrock))/                  &
+                                                      (((delta - layer(index_bedrock))/(TI(i,index_bedrock,islope)**2)) + &
+                                                      ((layer(index_bedrock + 1) - delta)/(TI_breccia**2))))
                 do k = index_bedrock + 2,nsoil
                     TI(i,k,islope) = TI_bedrock
@@ -427 +426 @@
         do i = 1,ngrid
             if (inertiedat(i,index_breccia) < TI_breccia) then
-                inertiedat(i,index_breccia + 1) = sqrt((layer(index_breccia + 1) - layer(index_breccia))/                    &
-                                                        (((delta - layer(index_breccia))/(inertiedat(i,index_breccia)**2)) + &
-                                                        ((layer(index_breccia + 1) - delta)/(TI_breccia**2))))
+                inertiedat(i,index_breccia + 1) = sqrt((layer(index_breccia + 1) - layer(index_breccia))/                   &
+                                                       (((delta - layer(index_breccia))/(inertiedat(i,index_breccia)**2)) + &
+                                                       ((layer(index_breccia + 1) - delta)/(TI_breccia**2))))
                 do k = index_breccia + 2,index_bedrock
                     inertiedat(i,k) = TI_breccia
@@ -442 +441 @@
         delta = depth_bedrock
         do i = 1,ngrid
-            inertiedat(i,index_bedrock + 1) = sqrt((layer(index_bedrock + 1) - layer(index_bedrock))/                    &
-                                                    (((delta - layer(index_bedrock))/(inertiedat(i,index_bedrock)**2)) + &
-                                                    ((layer(index_bedrock + 1) - delta)/(TI_bedrock**2))))
+            inertiedat(i,index_bedrock + 1) = sqrt((layer(index_bedrock + 1) - layer(index_bedrock))/                   &
+                                                   (((delta - layer(index_bedrock))/(inertiedat(i,index_bedrock)**2)) + &
+                                                   ((layer(index_bedrock + 1) - delta)/(TI_bedrock**2))))
         end do