Index: trunk/LMDZ.MARS/libf/phymars/physiq_mod.F
===================================================================
--- trunk/LMDZ.MARS/libf/phymars/physiq_mod.F	(revision 2079)
+++ trunk/LMDZ.MARS/libf/phymars/physiq_mod.F	(revision 2080)
@@ -3068,8 +3068,8 @@
            if (rdstorm) then
              call writediagfi(1,'aerosol_dust','opacity of env. dust',''
-     &                        ,1,aerosol(:,:,igcm_dust_mass))
+     &                        ,1,aerosol(:,:,1))
              call writediagfi(1,'aerosol_stormdust',
      &                         'opacity of storm dust',''
-     &                        ,1,aerosol(:,:,igcm_stormdust_mass))
+     &                        ,1,aerosol(:,:,2))
              call WRITEDIAGFI(ngrid,'dqsdifdustq','diffusion',
      &                       'kg.m-2.s-1',0,zdqsdif(1,igcm_dust_mass))
Index: trunk/LMDZ.MARS/libf/phymars/vdifc_mod.F
===================================================================
--- trunk/LMDZ.MARS/libf/phymars/vdifc_mod.F	(revision 2079)
+++ trunk/LMDZ.MARS/libf/phymars/vdifc_mod.F	(revision 2080)
@@ -548,5 +548,5 @@
 
 ! calcul de zc et zd pour la couche top en prenant en compte le terme
-! de variation de masse (on fait une boucle pour que ça converge)
+! de variation de masse (on fait une boucle pour que \E7a converge)
 
 ! Identification des points de grilles qui ont besoin de la correction
@@ -740,8 +740,10 @@
                   endif ! if (submicron)
      		ELSE ! outside dust injection time frame
-                 	pdqsdif(ig,igcm_dust_mass)= 0.
-                        pdqsdif(ig,igcm_dust_number)= 0.
+                  pdqsdif(ig,igcm_dust_mass)= 0.
+                  pdqsdif(ig,igcm_dust_number)= 0.
+                  if (rdstorm) then      
                  	pdqsdif(ig,igcm_stormdust_mass)= 0.
                         pdqsdif(ig,igcm_stormdust_number)= 0.
+                  end if
      		ENDIF