Changeset 54cbd6c in flexpart.git for src/readcommand.f90


Ignore:
Timestamp:
Sep 30, 2016, 11:01:54 AM (7 years ago)
Author:
Sabine <sabine.eckhardt@…>
Branches:
master, 10.4.1_pesei, GFS_025, bugfixes+enhancements, dev, release-10, release-10.4.1, scaling-bug, univie
Children:
9669e1e
Parents:
dced13c
Message:

all f90 files for dry/wet backward mode

File:
1 edited

Legend:

Unmodified
Added
Removed
  • src/readcommand.f90

    r462f74b r54cbd6c  
    296296  !Af          1 = mass units
    297297  !Af          2 = mass mixing ratio units
     298  !Af          3 = wet deposition in outputfield
     299  !Af          4 = dry deposition in outputfield
    298300
    299301  if ( ldirect .eq. 1 ) then  ! FWD-Run
     
    318320     endif
    319321  !Af set release-switch
    320      if (ind_receptor .eq. 1) then !mass
     322     WETBKDEP=.false.
     323     DRYBKDEP=.false.
     324     select case (ind_receptor)
     325     case (1)  !  1 .. concentration at receptor
    321326        ind_rel = 1
    322      else ! mass mix
     327     case (2)  !  2 .. mixing ratio at receptor
    323328        ind_rel = 0
    324      endif
     329     case (3)  ! 3 .. wet deposition in outputfield
     330        ind_rel = 3
     331         write(*,*) ' #### FLEXPART WET DEPOSITION BACKWARD MODE    #### '
     332         WETBKDEP=.true.
     333         allocate(xscav_frac1(maxpart,maxspec))
     334     case (4)  ! 4 .. dry deposition in outputfield
     335         ind_rel = 4
     336         write(*,*) ' #### FLEXPART DRY DEPOSITION BACKWARD MODE    #### '
     337         DRYBKDEP=.true.
     338         allocate(xscav_frac1(maxpart,maxspec))
     339     end select
    325340  endif
    326341
     
    338353    stop
    339354  endif
    340 
    341  if ((ldirect.eq.-1).and.(iout.eq.6)) then
    342       if ((ind_receptor .eq. 1) .and.  (ind_source .eq. 1)) then
    343          write(*,*) ' #### FLEXPART SCAVENGING DEPOSIT BACKWARD MODE    #### '
    344          SCAVDEP=.true.
    345          allocate(xscav_frac1(maxpart,maxspec))
    346       else
    347         write(*,*) '#### FLEXPART MODEL ERROR! FILE COMMAND:     ####'
    348         write(*,*) '#### FOR SCAVDEP MODE ind_source and         ####'
    349         write(*,*) '#### ind_receptor have to be 1 !             ####'
    350         stop
    351       endif
    352     else
    353       SCAVDEP=.false.
    354     endif
    355355
    356356  ! Check input dates
     
    401401  if ((iout.lt.1).or.(iout.gt.6)) then
    402402    write(*,*) ' #### FLEXPART MODEL ERROR! FILE COMMAND:     #### '
    403     write(*,*) ' #### IOUT MUST BE 1, 2, 3, 4, 5 OR 6 FOR     #### '
     403    write(*,*) ' #### IOUT MUST BE 1, 2, 3, 4 OR 5 FOR        #### '
    404404    write(*,*) ' #### STANDARD FLEXPART OUTPUT OR  9 - 13     #### '
    405405    write(*,*) ' #### FOR NETCDF OUTPUT                       #### '
Note: See TracChangeset for help on using the changeset viewer.
hosted by ZAMG