Changeset 18adf60 in flexpart.git for src/timemanager.f90

Timestamp:

May 9, 2016, 9:35:01 AM (8 years ago)

Author:

Espen Sollum ATMOS <eso@…>

Branches:

master, 10.4.1_pesei, GFS_025, bugfixes+enhancements, dev, release-10, release-10.4.1, scaling-bug, univie

Children:

3b80e98, 93786a1

Parents:

62e65c7

Message:

Bug found in the MQUASILAG=1 option, xmass array accessed with out-of-bounds index. Quick fix implemented (see changes in timemanager and advance).

File:

• src/timemanager.f90 (modified) (2 diffs)

src/timemanager.f90

@@ -589 +589 @@
             endif
 
-
-            if (mdomainfill.eq.0) then
+! Skip check on mass fraction when npoint represents particle number
+            if (mdomainfill.eq.0.and.mquasilag.eq.0) then
               if (xmass(npoint(j),ks).gt.0.) &
                    xmassfract=max(xmassfract,real(npart(npoint(j)))* &
@@ -605 +605 @@
 !ZHG 2015
             else
-              xmassfract=1.
-            endif
+              xmassfract=1.0
+            end if
           end do