Changeset 462f74b in flexpart.git for src/readcommand.f90

Timestamp:

Jul 12, 2016, 11:02:42 AM (8 years ago)

Author:

Sabine Eckhardt <sabine@…>

Branches:

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

Children:

842074e

Parents:

f28aa0a

Message:

first version of the backward scavenging

File:

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

src/readcommand.f90

@@ -339 +339 @@
   endif
 
+ if ((ldirect.eq.-1).and.(iout.eq.6)) then
+      if ((ind_receptor .eq. 1) .and.  (ind_source .eq. 1)) then
+         write(*,*) ' #### FLEXPART SCAVENGING DEPOSIT BACKWARD MODE    #### '
+         SCAVDEP=.true.
+         allocate(xscav_frac1(maxpart,maxspec))
+      else
+        write(*,*) '#### FLEXPART MODEL ERROR! FILE COMMAND:     ####'
+        write(*,*) '#### FOR SCAVDEP MODE ind_source and         ####'
+        write(*,*) '#### ind_receptor have to be 1 !             ####'
+        stop
+      endif
+    else
+      SCAVDEP=.false.
+    endif
+
   ! Check input dates
   !******************
@@ -384 +399 @@
   !**********************************************************************
 
-  if ((iout.lt.1).or.(iout.gt.5)) then
+  if ((iout.lt.1).or.(iout.gt.6)) then
     write(*,*) ' #### FLEXPART MODEL ERROR! FILE COMMAND:     #### '
-    write(*,*) ' #### IOUT MUST BE 1, 2, 3, 4, OR 5 FOR       #### '
+    write(*,*) ' #### IOUT MUST BE 1, 2, 3, 4, 5 OR 6 FOR     #### '
     write(*,*) ' #### STANDARD FLEXPART OUTPUT OR  9 - 13     #### '
     write(*,*) ' #### FOR NETCDF OUTPUT                       #### '