Index: trunk/LMDZ.MARS/libf/phymars/initracer.F90
===================================================================
--- trunk/LMDZ.MARS/libf/phymars/initracer.F90	(revision 4044)
+++ trunk/LMDZ.MARS/libf/phymars/initracer.F90	(revision 4052)
@@ -14,4 +14,5 @@
                           callnlte, nltemodel, dustinjection,    &
                           co2clouds, co2useh2o, meteo_flux
+use ioipsl_getin_p_mod, only : getin_p
 
 implicit none
@@ -618,4 +619,6 @@
     ! '''''''''''
     reff_lift = 3.0e-6 !3.e-6 !Effective radius of lifted dust (m)
+    call getin_p("reff_lift_dust",reff_lift)
+    
     alpha_devil(igcm_dust_mass)=9.e-9   !  dust devil lift mass coeff
     ! alpha_lift(igcm_dust_mass)=3.0e-15  !  Lifted mass coeff
@@ -629,5 +632,5 @@
     alpha_lift(igcm_dust_mass)=1.e-6 !1.e-6 !Lifted mass coeff
     if (dustinjection >= 1) then
-        reff_lift = 3.0e-6 ! Effective radius of lifted dust (m)
+        !reff_lift = 3.0e-6 ! Effective radius of lifted dust (m)! already set
         alpha_lift(igcm_dust_mass)=(4/3.)*reff_lift*rho_dust/2.4
     endif
@@ -649,5 +652,5 @@
     ! normal dust
     if (rdstorm) then
-        reff_storm=3.e-6 ! reff_lift !3.e-6
+        reff_storm=reff_lift !3.e-6
         r0_storm=reff_storm/ref_r0
         rho_q(igcm_stormdust_mass)=rho_dust
@@ -672,6 +675,6 @@
         rho_q(igcm_topdust_mass)=rho_dust
         rho_q(igcm_topdust_number)=rho_dust
-        radius(igcm_topdust_mass) = 3.e-6
-        radius(igcm_topdust_number) = 3.e-6
+        radius(igcm_topdust_mass) = reff_lift !3.e-6
+        radius(igcm_topdust_number) = reff_lift !3.e-6
     end if !(topflows)
     ! ----------------------------------------------------------------------