Commit 8d70e43d authored by Ignacio Pisso's avatar Ignacio Pisso
Browse files

Update SPECIES directory

parent f4a43397
&SPECIES_PARAMS
PSPECIES="BC", ! Tracer name
PDECAY=-999.99, ! Species half life
PWETA_GAS=-0.9E-9, ! Below-cloud scavenging (gases) - A (weta_gas)
PWETB_GAS=-9.9, ! Below-cloud scavenging (gases) - B (wetb_gas)
PCRAIN_AERO=1.0, ! Below-cloud scavenging (particles) - Crain (crain_aero)
PCSNOW_AERO=1.0, ! Below-cloud scavenging (particles) - Csnow (csnow_aero)
PCCN_AERO=0.9, ! In-cloud scavenging (particles) - CCNeff (ccn_aero)
PIN_AERO=0.1, ! In-cloud scavenging (particles) - INeff (in_aero)
PDENSITY=2000.0, ! Dry deposition (particles) - rho
PDQUER=4.0E-07, ! Dry deposition (particles) - dquer
PDSIGMA=0.3, ! Dry deposition (particles) - dsig
PDRYVEL=-9.9, ! Alternative: dry deposition velocity
PRELDIFF=-9.9, ! Dry deposition (gases) - D
PHENRY=-0.9E-9, ! Dry deposition (gases) - Henrys const.
PF0=-9, ! Dry deposition (gases) - f0 (reactivity)
PWEIGHTMOLAR=-9.9, ! molweight
POHCCONST=-0.9E-9, ! OH Reaction rate - C [cm^3/molecule/sec]
POHDCONST=-9.9, ! OH Reaction rate - D [K]
POHNCONST=2.0, ! OH Reaction rate - C [cm^3/molecule/sec]
/
!**********************************************************************
! Copyright 1998,1999,2000,2001,2002,2005,2007,2008,2009,2010 *
! Andreas Stohl, Petra Seibert, A. Frank, Gerhard Wotawa, *
! Caroline Forster, Sabine Eckhardt, John Burkhart, Harald Sodemann *
! *
! This file is part of FLEXPART. *
! *
! FLEXPART is free software: you can redistribute it and/or modify *
! it under the terms of the GNU General Public License as published by*
! the Free Software Foundation, either version 3 of the License, or *
! (at your option) any later version. *
! *
! FLEXPART is distributed in the hope that it will be useful, *
! but WITHOUT ANY WARRANTY; without even the implied warranty of *
! MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the *
! GNU General Public License for more details. *
! *
! You should have received a copy of the GNU General Public License *
! along with FLEXPART. If not, see <http://www.gnu.org/licenses/>. *
!**********************************************************************
implicit none
character(len=11) :: speciesfn
character(len=3) :: aspec
character(len=16) :: pspecies
real :: pdecay, pweta_gas, pwetb_gas, preldiff, phenry, pf0, pdensity, pdquer
real :: pdsigma, pdryvel, pweightmolar, pohcconst, pohdconst, pohnconst
real :: pcrain_aero, pcsnow_aero, pccn_aero, pin_aero
integer :: readerror, unitspecies, specnumber
! declare namelist
namelist /species_params/ &
pspecies, pdecay, pweta_gas, pwetb_gas, &
pcrain_aero, pcsnow_aero, pccn_aero, pin_aero, &
preldiff, phenry, pf0, pdensity, pdquer, &
pdsigma, pdryvel, pweightmolar, pohcconst, pohdconst, pohnconst
unitspecies=4
write(*,*) ' Species | | WetDep(gas) | DryDep(gas) |WetDep(below-C)| WetDep(in-C)|'//&
' DryDepo(particles) Altern| Radioact. | OH Reaction |'
write(*,*) ' Name |molwght| A B | D H f0 | Crain Csnow | ccn in |' //&
' rho dquer dsig vd | Halflife[s]| C** D[T] N*** |'
write(*,*) '--------------|-------|-----------------|------------------|---------------|-------------|'//&
'-------------------------------|------------|----------------------|'
! write(*,*) ' Specie | Radioact. | WetDep(gas) |WetDep(below-C)| WetDep(in-C)| DryDepo(gas) |'//&
! ' DryDepo(particles) Altern| | OH Reaction |'
! write(*,*) ' Name | Halflife[s]| A B | Crain Csnow | ccn in | D H f0 |' //&
! ' rho dquer dsig vd |molwght| C** D[T] N*** |'
! write(*,*) '--------------|------------|-----------------|---------------|-------------|-------------------|'//&
! '-------------------------------|-------|----------------------|'
do specnumber=1,100
write (aspec,'(i0.3)') specnumber
speciesfn='SPECIES_'//aspec
! write(*,*) 'Processing: ',speciesfn
pspecies="" ! read failure indicator value
pdecay=-9.9
pweta_gas=-0.9E-09
pwetb_gas=0.0
pcrain_aero=-9.9
pcsnow_aero=-9.9
pccn_aero=-9.9
pin_aero=-9.9
preldiff=-9.9
phenry=0.0
pf0=0.0
pdensity=-0.9E09
pdquer=0.0
pdsigma=0.0
pdryvel=-9.99
pohcconst=-9.9
pohdconst=-9.9
pohnconst=2.0
pweightmolar=-9.9
! Open the SPECIES file and read species names and properties
!************************************************************
open(unitspecies,file=speciesfn,status='old',form='formatted',err=998)
read(unitspecies,species_params,err=998)
close(unitspecies)
write(*,45) specnumber,' ',pspecies,'|',pweightmolar,'|',pweta_gas,' ',pwetb_gas,'|', &
preldiff,' ',phenry,' ',pf0,'|', &
pcrain_aero,' ',pcsnow_aero,'|',pccn_aero,' ',pin_aero,'|', &
pdensity,pdquer,pdsigma,pdryvel,'|',pdecay,'|',pohcconst,pohdconst,pohnconst,'|'
45 format(i3,a1,a11,a1,f7.1,a1,e8.1,a1,f8.2,a1, &
f4.1,a1,e8.1,a1,f4.1,a1, &
f7.1,a1,f7.1,a1,f6.1,a1,f6.1,a1, &
e8.1,e9.1,f7.1,f7.2,a1,f12.1,a1,e8.1,f7.1,f7.1,a1)
! write(*,45) specnumber,' ',pspecies,'|',pdecay,'|',pweta_gas,' ',pwetb_gas,'|',pcrain_aero,' ', &
! pcsnow_aero,'|',pccn_aero,' ',pin_aero,'|',preldiff,' ',phenry,' ',pf0,'|', &
! pdensity,pdquer,pdsigma,pdryvel,'|',pweightmolar,'|',pohcconst,pohdconst,pohnconst,'|'
!5 format(i3,a1,a11,a1,f12.1,a1,e8.1,a1,f8.2,a1,f7.1,a1,f7.1,a1,f6.1,a1,f6.1,a1,f5.1,a1,e8.1,a1,f4.1,a1, &
! e8.1,e9.1,f7.1,f7.2,a1,f7.1,a1,e8.1,f7.1,f7.1,a1)
998 continue
enddo
write(*,*) '** unit [cm^3/molec/s] (in FLEXPART version 9.2 and below this had unit [cm3/s], note the unit is now changed!)'
write(*,*) '*** no unit'
end
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment