create textfile with species overfiew

options/SPECIES/specoverview.f90

+!**********************************************************************
+! 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