Changes in src/readreleases.f90 [8ee24a5:c8fc724] in flexpart.git

File:

• src/readreleases.f90 (modified) (8 diffs)

src/readreleases.f90

@@ -218 +218 @@
   rewind(unitreleases)
 
+  if (nspec.gt.maxspec) goto 994
+
   ! allocate arrays of matching size for number of species (namelist output)
   deallocate(mass)
@@ -274 +276 @@
   do i=1,maxspec
     DRYDEPSPEC(i)=.false.
+    WETDEPSPEC(i)=.false.
   end do
 
@@ -367 +370 @@
          &(dquer(i).gt.0. .and. (crain_aero(i) .gt. 0. .or. csnow_aero(i).gt.0.)))  then
       WETDEP=.true.
+      WETDEPSPEC(i)=.true.
       if (lroot) then
         write (*,*) '  Below-cloud scavenging: ON'
@@ -378 +382 @@
     if (dquer(i).gt.0..and.(ccn_aero(i).gt.0. .or. in_aero(i).gt.0.))  then
       WETDEP=.true.
+      WETDEPSPEC(i)=.true.
       if (lroot) then
         write (*,*) '  In-cloud scavenging: ON'
@@ -397 +402 @@
     endif
 
-  end do
+  end do ! end loop over species
 
   if (WETDEP.or.DRYDEP) DEP=.true.
@@ -515 +520 @@
     stop
   endif
-
-  ! If FLEXPART is run for backward deposition force zpoint
-  !*********************************************************************
-  if (WETBKDEP) then
-      zpoint1(numpoint)=0.
-      zpoint2(numpoint)=20000.
-      kindz(numpoint)=1
-  endif
-  if (DRYBKDEP) then
-      zpoint1(numpoint)=0.
-      zpoint2(numpoint)=2.*href
-      kindz(numpoint)=1
-  endif
-
 
   ! Check whether x coordinates of release point are within model domain
@@ -588 +579 @@
   endif
 
+
   ! Determine the release rate (particles per second) and total number
   ! of particles released during the simulation
@@ -647 +639 @@
   endif
 
+  if (lroot) then
+    write(*,FMT='(A,ES14.7)') ' Total mass released:', sum(xmass(1:numpoint,1:nspec))
+  end if
+
   return