Index: trunk/LMDZ.MARS/libf/phymars/dustlift.F
===================================================================
--- trunk/LMDZ.MARS/libf/phymars/dustlift.F	(revision 308)
+++ trunk/LMDZ.MARS/libf/phymars/dustlift.F	(revision 310)
@@ -50,14 +50,22 @@
 !!!! AS: you have to compile with -DMESOSCALE to do so
       REAL alpha
+      REAL r0_lift
       INTEGER ierr
+      REAL ulim
         OPEN(99,file='stress.def',status='old',form='formatted'
      .   ,iostat=ierr)
         !!! no file => default values
         IF(ierr.EQ.0) THEN
-          READ(99,*) stress_seuil
+          READ(99,*) ulim !ulim = sqrt(stress_seuil/rho) avec rho = 0.02. 
+                          !prendre ulim = 1.061 m/s pour avoir stress_seuil = 0.0225 
           READ(99,*) alpha
+          stress_seuil = 0.02 * ulim * ulim
           write(*,*) 'USER-DEFINED threshold: ', stress_seuil, alpha
           CLOSE(99)
-          alpha_lift(1:nq) = alpha
+          alpha_lift(igcm_dust_mass) = alpha
+          r0_lift = radius(igcm_dust_mass) / ref_r0
+          alpha_lift(igcm_dust_number)=r3n_q*
+     &                        alpha_lift(igcm_dust_mass)/r0_lift**3
+          write(*,*) 'set dust number: ', alpha_lift(igcm_dust_number)
         ENDIF
 #endif
Index: trunk/LMDZ.MARS/libf/phymars/vdifc.F
===================================================================
--- trunk/LMDZ.MARS/libf/phymars/vdifc.F	(revision 308)
+++ trunk/LMDZ.MARS/libf/phymars/vdifc.F	(revision 310)
@@ -494,4 +494,5 @@
 c       Dust lifting:
         if (lifting) then
+#ifndef MESOSCALE
            if (doubleq.AND.submicron) then
              do ig=1,ngrid
@@ -522,4 +523,8 @@
      &                    pdqsdif)
            endif !doubleq.AND.submicron
+#else
+            call dustlift(ngrid,nlay,nq,rho,zcdh_true,zcdh,co2ice,
+     &                    pdqsdif)
+#endif
         else
            pdqsdif(1:ngrid,1:nq) = 0.