Changeset 2255

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

Location:

trunk/LMDZ.MARS

Files:

• README (modified) (1 diff)
• libf/phymars/topmons_mod.F90 (modified) (2 diffs)

trunk/LMDZ.MARS/README

@@ -2875 +2875 @@
 == 06/03/2020 == AB
 Bug fix following r2248 in aeropacity_mod and topmons_mod : since dsodust, dsords and dsotop are diagnostic physiq_mod variables, we don't want them to be reinitialized at each call of aeropacity_mod and topmons_mod, but we initialize them once and for all at the beginning of physiq_mod instead.
+
+== 10/03/2020 == AB
+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)

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