Changeset 2255 for trunk/LMDZ.MARS/libf/phymars

Timestamp:

Mar 10, 2020, 5:01:24 PM (5 years ago)

Author:

abierjon

Message:

Mars GCM:
Bug fix in topmons_mod : the array rhobarz(ig,l) (density at interlayer levels) had only nlayer allocated instead of nlayer+1 (caused issues in calculating the topdust detrainment at the top of the model)
AB

File:

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

trunk/LMDZ.MARS/libf/phymars/topmons_mod.F90

@@ -156 +156 @@
 !     Entrainment
       REAL entr(ngrid)                                             ! entrainment flux
-      REAL rho(ngrid,nlayer),rhobarz(ngrid,nlayer)                 ! density, density at the levels
+      REAL rho(ngrid,nlayer),rhobarz(ngrid,nlayer+1)                 ! density, density at the levels
       REAL masse(ngrid,nlayer)                                     ! air mass
       REAL masse_pbl(ngrid)                                        ! total air mass within the PBL
@@ -458 +458 @@
           rhobarz(:,l)=pplev(:,l)/(r*0.5*(pt(:,l)+pt(:,l-1)))
         ENDDO
+        rhobarz(:,nlayer+1)=0. !top of the model
         !! Mass computation
         DO l=1,nlayer