Changeset 3512 for trunk/LMDZ.COMMON/libf/evolution/NS_fast_subs_mars.F90

Timestamp:

Nov 12, 2024, 6:35:33 PM (11 days ago)

Author:

jbclement

Message:

PEM:
Few corrections related to r3498 (time step from integer to real) and r3493 (Norbert Schorghofer's subroutines for dynamic ice table) in order to make the code work properly.
JBC

File:

• trunk/LMDZ.COMMON/libf/evolution/NS_fast_subs_mars.F90 (modified) (1 diff)

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

@@ -116 +116 @@
      call icechanges(nz,z(:),typeF,avdrho(k),avdrhoP(k),ypp(:), &
           & Diff,porosity,icefrac,bigstep,zdepthT(k),porefill(:,k),typeG)
-
+     if (typeP>2) then
      if (mode2 .and. porefill(typeP,k)>=1. .and. porefill(typeP-1,k)==0.) then  ! nothing changed
         porefill(typeP-1,k)=1.  ! paste a layer
    !     write(34,*) '# mode 2 growth occurred',typeP,typeF,typeT
+     endif
      endif