--- mo_EVA_orig.f90	2017-11-02 14:17:29.000000000 +0100
+++ mo_EVA.f90	2019-10-06 21:43:06.000000000 +0200
@@ -29,8 +29,9 @@
    !----------------------------------------------------------------------------------------
    ! Input files
 
-   CHARACTER(len=50)    :: eruption_list_filename = 'Eruption_list_GMD_1960_2015.nc'
-   CHARACTER(len=50)    :: parameter_set_filename = 'EVAv1_parameter_set_v1.1.nc'
+!   CHARACTER(len=50)    :: eruption_list_filename = 'Eruption_list_GMD_1960_2015.nc'
+   CHARACTER(len=50)    :: eruption_list_filename = 'eruptions.nc'
+   CHARACTER(len=50)    :: parameter_set_filename = 'EVAv1_parameter_set_v1.0.nc'
    CHARACTER(len=50)    :: Lookuptable_filename   = 'eva_Mie_lookuptables.nc'
 
    ! --------------------------------------------------------------------------------------
@@ -173,7 +174,9 @@
      real :: aod, SO4_lat
      real :: v_shape(nz), v_shape_norm(nz)
 
-     nonlin_threshold = (so4_to_aod_nonlin/so4_to_aod_linear)**(3.)
+!     nonlin_threshold = (so4_to_aod_nonlin/so4_to_aod_linear)**(3.)
+!     jan:
+     nonlin_threshold = (so4_to_aod_nonlin/so4_to_aod_linear)**(3./(1.0+3.0*0.1))
 
       ! find index of lat0 in lat
       do i=1,nlat
@@ -189,7 +192,9 @@
       if ( SO4_lat < nonlin_threshold ) then
          aod=so4_to_aod_linear*SO4_lat
       else
-         aod=so4_to_aod_nonlin*(SO4_lat**(2.0/3.0))
+!         aod=so4_to_aod_nonlin*(SO4_lat**(2.0/3.0))
+!         jan
+         aod=so4_to_aod_nonlin*(SO4_lat**(2.0/3.0-0.1))
       end if
 
       v_shape=shape_EQ(lat0_ind,:)*SO4(2) &