Changeset 3250 for trunk/LMDZ.MARS/libf/phymars/dyn1d

Timestamp:

Mar 4, 2024, 11:03:35 AM (16 months ago)

Author:

llange

Message:

Mars PCM

• Fixing bug in initestphys and waterice_tifeedback: when changing the soil

thermal properties to account for the high thermal inertia of subsurface ice,
the wrong layer was considered as 'mixed' (it was the layer after the correct
one which was modified).

LL

File:

• trunk/LMDZ.MARS/libf/phymars/dyn1d/init_testphys1d_mod.F90 (modified) (2 diffs)

trunk/LMDZ.MARS/libf/phymars/dyn1d/init_testphys1d_mod.F90

@@ -678 +678 @@
             do isoil = 1,nsoil
                 if ((ice_depth >= layer(isoil)) .and. (ice_depth < layer(isoil + 1))) then
-                    iref = isoil + 1
+                    iref = isoil 
                     exit
                 endif
@@ -685 +685 @@
             inertiedat(1,:iref - 1) = inertiedat(1,1)
             ! We compute the transition in layer(iref)
-            inertiedat(1,iref) = sqrt((layer(iref) - layer(iref - 1))/(((ice_depth - layer(iref - 1))/inertiedat(1,1)**2) + ((layer(iref) - ice_depth)/inertieice**2)))
+            inertiedat(1,iref) = sqrt((layer(iref+1) - layer(iref))/(((ice_depth - layer(iref))/inertiedat(1,1)**2) + ((layer(iref+1) - ice_depth)/inertieice**2)))
             ! Finally, we compute the underlying ice:
             inertiedat(1,iref + 1:) = inertieice