Changeset 3563 for trunk/LMDZ.COMMON/libf/evolution/NS_dyn_ss_ice_m.F90

Timestamp:

Dec 17, 2024, 5:52:24 PM (5 days ago)

Author:

jbclement

Message:

PEM:
Correction of Norbert Schorghofer's code due to missing initialization and bad shape array as subroutine argument + some cleanings.
JBC

File:

• trunk/LMDZ.COMMON/libf/evolution/NS_dyn_ss_ice_m.F90 (modified) (2 diffs)

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

@@ -37 +37 @@
   real(8)  ssi_depth_in, ssi_depth, T1
   real(8), dimension(NP) :: zdepthF, zdepthE, zdepthT, zdepthG
-  real(8), dimension(NMAX,NP) :: porefill, porefill_in
+  real(8), dimension(nz,NP) :: porefill, porefill_in
   real(8), dimension(nz) :: Tb
   real(8), dimension(NP) :: Tmean1, Tmean3, avrho1
@@ -152 +152 @@
   !print *,'Zt0=  ',ZdepthT
      call icelayer_mars(timestep,nz,NP,thIn,rhoc,z,porosity,pfrost,Tb,zdepthF, &
-          & zdepthE,porefill(1:nz,:),Tmean1,Tmean3,zdepthG, &
+          & zdepthE,porefill,Tmean1,Tmean3,zdepthG, &
           & albedo,p0,icefrac,zdepthT,avrho1, &
           & avrho1prescribed)