Index: trunk/LMDZ.MARS/libf/phymars/physiq_mod.F
===================================================================
--- trunk/LMDZ.MARS/libf/phymars/physiq_mod.F	(revision 3146)
+++ trunk/LMDZ.MARS/libf/phymars/physiq_mod.F	(revision 3147)
@@ -2665,15 +2665,20 @@
 c        -------------------------------------------------------------------
 
-       if (water) then
          do ig=1,ngrid
 	   if(mu0(ig).le.0.01) then
 	    fluxsurf_dir_dn_sw(ig) = 0.
 	   else
-            fluxsurf_dir_dn_sw(ig) = flux_1AU/dist_sol/dist_sol*mu0(ig)*
+            if (water) then
+             ! both water and dust contribute
+             fluxsurf_dir_dn_sw(ig) = flux_1AU/dist_sol/dist_sol*mu0(ig)*
      &                    exp(-(tau(ig,iaer_dust_doubleq)+
      &                          tau(ig,iaer_h2o_ice))/mu0(ig))
-           endif
+            else
+             ! only dust contributes
+             fluxsurf_dir_dn_sw(ig) = flux_1AU/dist_sol/dist_sol*mu0(ig)*
+     &                    exp(-(tau(ig,iaer_dust_doubleq))/mu0(ig))
+            endif ! of if (water)
+           endif ! of if(mu0(ig).le.0.01)
          enddo
-       endif
 
            ! Density-scaled opacities
@@ -3286,4 +3291,5 @@
         call write_output("local_time","Local time",
      &                   'sol',local_time(:))
+        if (water) then
             if (.not.activice) then
                CALL write_output('tauTESap',
@@ -3295,4 +3301,5 @@
      &                         '',taucloudtes(:))
              endif
+        endif ! of if (water)
 #else
       !!! this is to ensure correct initialisation of mesoscale model