Index: options/SPECIES/specoverview.f90
===================================================================
--- options/SPECIES/specoverview.f90	(revision 95a45d34ce18935ea257ed298f98e77418e8f2e1)
+++ options/SPECIES/specoverview.f90	(revision dba422187e10650e160348f0a0271adfac6b8b88)
@@ -19,4 +19,5 @@
 ! along with FLEXPART.  If not, see <http://www.gnu.org/licenses/>.   *
 !**********************************************************************
+program specoverview
   implicit none
 
@@ -27,89 +28,84 @@
   real :: pdsigma, pdryvel, pweightmolar, pohcconst, pohdconst, pohnconst
   real :: pcrain_aero, pcsnow_aero, pccn_aero, pin_aero
-  integer :: readerror, unitspecies, specnumber
+  real :: parea_dow(7), parea_hour(24), ppoint_dow(7), ppoint_hour(24)
+  integer :: ierr, unitspecies, specnumber
 
-! declare namelist
+  ! 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 
+       pdsigma, pdryvel, pweightmolar, pohcconst, pohdconst, pohnconst, &
+       parea_dow, parea_hour, ppoint_dow, ppoint_hour
 
   unitspecies=4
 
-  write(*,*) '    Species   |       |   WetDep(gas)   |    DryDep(gas)   |WetDep(below-C)| WetDep(in-C)|'//& 
-             '     DryDepo(particles)  Altern| Radioact.  |     OH Reaction      |'
+  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(*,*) '    Name      |molwght| A          B    | D    H        f0 | Crain   Csnow |  ccn    in  |'// &
+       '   rho    dquer    dsig    vd  | Halflife[s]|   C**     D[T]  N*** |'
 
-  write(*,*) '--------------|-------|-----------------|------------------|---------------|-------------|'//&
-             '-------------------------------|------------|----------------------|'
+  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
 
-  do specnumber=1,100
-  
-  write (aspec,'(i0.3)') specnumber
-  speciesfn='SPECIES_'//aspec
- 
-! write(*,*) 'Processing: ',speciesfn
+     write (aspec,'(i0.3)') specnumber
+     speciesfn='SPECIES_'//aspec
 
-  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
+     ! write(*,*) 'Processing: ',speciesfn
 
-! 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,'|'
+     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
 
-45 format(i3,a1,a11,a1,f7.1,a1,e8.1,a1,f8.2,a1, &
+     ! Open the SPECIES file and read species names and properties
+     !************************************************************
+     open(unitspecies,file=speciesfn,status='old',form='formatted',iostat=ierr)
+     if (ierr /= 0) then
+        ! Report missing files below 40
+        if (specnumber <= 40) then
+           write(*,'(i3,2(1x,a))') specnumber, 'no such file:', trim(speciesfn)
+        end if
+        cycle
+     end if
+     read(unitspecies,species_params)  ! to skip errors add: ,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,'|'
+998  continue
+  enddo
 
-!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'
 
-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'
+  print*,'rho: density'
 
- print*,'rho: density'
-
-end
+end program specoverview
