Changes in src/readreleases.f90 [c8fc724:8ee24a5] 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)
@@ -276 +274 @@
   do i=1,maxspec
     DRYDEPSPEC(i)=.false.
-    WETDEPSPEC(i)=.false.
   end do
 
@@ -370 +367 @@
          &(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'
@@ -382 +378 @@
     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'
@@ -402 +397 @@
     endif
 
-  end do ! end loop over species
+  end do
 
   if (WETDEP.or.DRYDEP) DEP=.true.
@@ -520 +515 @@
     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
@@ -579 +588 @@
   endif
 
-
   ! Determine the release rate (particles per second) and total number
   ! of particles released during the simulation
@@ -639 +647 @@
   endif
 
-  if (lroot) then
-    write(*,FMT='(A,ES14.7)') ' Total mass released:', sum(xmass(1:numpoint,1:nspec))
-  end if
-
   return