Merge branch 'dev' into release-10

options/SPECIES/SPECIES_015

 &SPECIES_PARAMS
- PSPECIES=""I-131",           ! Tracer name
+ PSPECIES="I-131",            ! Tracer name
  PDECAY=691200.0,             ! 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)

options/SPECIES/SPECIES_018

 &SPECIES_PARAMS
- PSPECIES=""Ru-106",          ! Tracer name
+ PSPECIES="Ru-106",          ! Tracer name
  PDECAY=31536000.0,           ! 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)

src/get_wetscav.f90

@@ -296,7 +296,8 @@ subroutine get_wetscav(itime,ltsample,loutnext,jpart,ks,grfraction,inc_count,blc
           else if (act_temp .ge. 273.) then
             liq_frac=1
           else
-            liq_frac =((act_temp-273.)/(273.-253.))**2.
+! sec, bugfix after FLEXPART paper review, liq_frac was 1-liq_frac
+            liq_frac =((act_temp-253.)/(273.-253.))**2.
           end if
 ! ZHG: Calculate the aerosol partition based on cloud phase and Ai and Bi
           frac_act = liq_frac*ccn_aero(ks) +(1-liq_frac)*in_aero(ks)