Merge remote-tracking branch 'sabine/dev' into dev

options/SPECIES/SPECIES.README

@@ -52,6 +52,7 @@ WET DEPOSITION
 Below-cloud scavenging
 - Gases:
 	A and B are the parameters for below-cloud scavenging defined by Eq. scavenging coefficient = A x precipitation-rate^B
+        If you migrate from FLEXPART version 9.2 and below this is equal to A and B as it was.
 - Particles: 
 	A "rain collection efficiency"  
 	B "snow collection efficiency" are both defined by Equations 46-47.
@@ -61,6 +62,7 @@ Below-cloud scavenging
 In-cloud scavenging
 - Gases:
     Is switched on if Ai or Bi are given greater than 0
+    If you migrate from FLEXPART version 9.2 and below this is equal to Ai and Bi should be set to 1
 - Particles:
 	Ai ("Cloud Condensation Nuclei (CCN) efficiency") 
 	Bi ("Ice Nuclei (IN) efficiency")                  are both from Eq. 51. 	

options/SPECIES/SPECIES_031

+****************************************************************************
+*                                                                          *
+*     Input file for the Lagrangian particle dispersion model FLEXPART     *
+*               Definition file of chemical species/radionuclides          *
+*                                                                          *
+****************************************************************************
+PCB28              Tracer name
+-999.9             Species half life
+3.5E-05            Below-cloud scavenging (gases and particles) - A
+0.62               Below-cloud scavenging (gases and particles)- B
+1.0                In-cloud scavenging (gases and particles) - Ai
+1.0                In-cloud scavenging (gases and particles) - Bi
+ 4.3               Dry deposition (gases) - D
+3.1e-00            Dry deposition (gases) - Henrys const.
+0.1                Dry deposition (gases) - f0 (reactivity)
+-9.9E09            Dry deposition (particles) - rho
+                   Dry deposition (particles) - dquer
+                   Dry deposition (particles) - dsig
+-9.99              Alternative: dry deposition velocity
+257.54             molweight
+1.07E-11           OH Reaction rate - C [cm^3/molecule/sec]
+1203.0             OH reaction rate - D [K]
+2.0                OH Reaction rate - N (no unit)
+-9                 number of associated specias (neg. none)
+-99.99             KOA - organic matter air partitioning

options/SPECIES/SPECIES_034

+****************************************************************************
+*                                                                          *
+*     Input file for the Lagrangian particle dispersion model FLEXPART     *
+*               Definition file of chemical species/radionuclides          *
+*                                                                          *
+****************************************************************************
+G-HCH              Tracer name
+-999.9             Species half life
+3.5E-05            Below-cloud scavenging (gases and particles) - A
+0.62               Below-cloud scavenging (gases and particles)- B
+1.0                In-cloud scavenging (gases and particles) - Ai
+1.0                In-cloud scavenging (gases and particles) - Bi
+4.3                Dry deposition (gases) - D
+4.22E+02           Dry deposition (gases) - Henrys const.
+0.1                Dry deposition (gases) - f0 (reactivity)
+-9.9E09            Dry deposition (particles) - rho
+                   Dry deposition (particles) - dquer
+                   Dry deposition (particles) - dsig
+-9.99              Alternative: dry deposition velocity
+290.83             molweight
+6.21E-11           OH Reaction rate - C [cm^3/molecule/sec]
+1203.0             OH reaction rate - D [K]
+2.0                OH Reaction rate - N (no unit)
+-9                 number of associated specias (neg. none)
+-99.99             KOA - organic matter air partitioning

src/readOHfield.f90

@@ -71,6 +71,9 @@ subroutine readOHfield
     thefile=trim(ohfields_path)//'OH_FIELDS/'//'geos-chem.OH.2005'//mm//'01.nc'
     ierr=nf_open(trim(thefile),NF_NOWRITE,nid)
     if(ierr.ne.0) then
+      write (*,*) 'The OH field at: '//thefile// ' could not be opened'
+      write (*,*) 'please copy the OH fields there, or change the path in the'
+      write (*,*) 'COMMAND namelist!'
       write(*,*) nf_strerror(ierr)
       stop
     endif