Changeset 4296

Timestamp:

Oct 11, 2022, 10:18:49 AM (2 years ago)

Author:

Ehouarn Millour

Message:

Add "clean" initializations of coefh and coefm which are defined on klev+1 vertical levels but mostly used/computed on klev throughout the code; unitialized topmost values cause problems in outputs when in debug mode.
EM

File:

• LMDZ6/trunk/libf/phylmd/phys_state_var_mod.F90 (modified) (1 diff)

LMDZ6/trunk/libf/phylmd/phys_state_var_mod.F90

@@ -560 +560 @@
       ALLOCATE(coefh(klon,klev+1,nbsrf+1))
       ALLOCATE(coefm(klon,klev+1,nbsrf+1))
+      ! initialize cleanly coefh,coefm
+      ! (most of the time in the code these are assumed to be on klev levels)
+      coefh(:,:,:)=0
+      coefm(:,:,:)=0
       ALLOCATE(zmax0(klon), f0(klon))
       ALLOCATE(sig1(klon,klev), w01(klon,klev))