Index: trunk/MESOSCALE/LMD_MM_MARS/SRC/LESnophys_/modif/module_initialize_les.F
===================================================================
--- trunk/MESOSCALE/LMD_MM_MARS/SRC/LESnophys_/modif/module_initialize_les.F	(revision 158)
+++ trunk/MESOSCALE/LMD_MM_MARS/SRC/LESnophys_/modif/module_initialize_les.F	(revision 159)
@@ -164,5 +164,9 @@
 !  stretch_grid = .true.
 !  FOR LES, set stretch to false
+#ifdef NOPHYS
+   stretch_grid = .true.
+#else
    stretch_grid = .false.
+#endif
    delt = 3.
 !   z_scale = .50
@@ -233,4 +237,5 @@
 
    IF (stretch_grid) THEN ! exponential stretch for eta (nearly constant dz)
+     print *, '*** STRETCH GRID ***' 
      DO k=1, kde
       grid%znw(k) = (exp(-(k-1)/float(kde-1)/z_scale) - exp(-1./z_scale))/ &
Index: trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/mars_lmd_new_storm/libf/phymars/meso_physiq.F
===================================================================
--- trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/mars_lmd_new_storm/libf/phymars/meso_physiq.F	(revision 158)
+++ trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/mars_lmd_new_storm/libf/phymars/meso_physiq.F	(revision 159)
@@ -211,4 +211,5 @@
       logical tracerdyn
       CHARACTER (len=20) :: wtnom(nqmx) ! tracer name
+      REAL toto
 
 c   outputs:
@@ -725,4 +726,26 @@
      $     zdtlw,zdtsw,fluxsurf_lw,fluxsurf_sw,fluxtop_lw,fluxtop_sw,
      &     tauref,tau,aerosol,ccn,rdust,rice,nuice,zdqnorm,dsodust)
+
+!!!!!!!!!!!!!!!!!!!!!! TEST TEST
+!         DO l=2,nlayer
+!            DO ig=1,ngrid
+!               toto = abs(rdust(ig,l) - rdust(ig,l-1))
+!               if ( toto .gt. 0.5*rdust(ig,l) ) then
+!                    write (*,*) "LIMITING RADIUS", ig, l
+!                    write (*,*) rdust(ig,l-1), rdust(ig,l)
+!                   rdust(ig,l) = rdust(ig,l-1)
+!               endif         
+!            ENDDO
+!         ENDDO
+!
+!         DO l=2,nlayer
+!            DO ig=2,ngrid
+!               toto = abs(rdust(ig-1,l) - rdust(ig,l))/rdust(ig,l)
+!               if ( toto .gt. 0.45 ) then
+!                write (*,*) "RDUST ", ig, l, rdust(ig-1,l), rdust(ig,l)
+!                rdust(ig,l) = rdust(ig,l-1)
+!               endif         
+!            ENDDO
+!         ENDDO