Index: trunk/LMDZ.MARS/libf/phymars/callradite_mod.F
===================================================================
--- trunk/LMDZ.MARS/libf/phymars/callradite_mod.F	(revision 2678)
+++ trunk/LMDZ.MARS/libf/phymars/callradite_mod.F	(revision 2685)
@@ -7,6 +7,7 @@
       SUBROUTINE callradite(icount,ngrid,nlayer,nq,zday,ls,pq,albedo,
      $     emis,mu0,pplev,pplay,pt,tsurf,fract,dist_sol,igout,
-     $     dtlw,dtsw,fluxsurf_lw,fluxsurf_sw,fluxtop_lw,
-     $     fluxtop_sw,tau_pref_scenario,tau_pref_gcm,
+     $     dtlw,dtsw,fluxsurf_lw,fluxsurf_dn_sw,fluxsurf_up_sw,
+     $     fluxtop_lw,fluxtop_dn_sw,fluxtop_up_sw,
+     &     tau_pref_scenario,tau_pref_gcm,
      &     tau,aerosol,dsodust,tauscaling,dust_rad_adjust,IRtoVIScoef,
      $     taucloudtes,rdust,rice,nuice,riceco2,nuiceco2,co2ice,
@@ -141,23 +142,4 @@
 c   nueffrad(ngrid,nlayer,naerkind) Aerosol effective variance
 
-c
-c  output:
-c  -------
-c dtlw (ngrid,nlayer)       longwave (IR) heating rate (K/s)
-c dtsw(ngrid,nlayer)        shortwave (Solar) heating rate (K/s)
-c fluxsurf_lw(ngrid)        surface downward flux tota LW (thermal IR) (W.m-2)
-c fluxsurf_sw(ngrid,1)      surface downward flux SW for solar band#1 (W.m-2)
-c fluxsurf_sw(ngrid,2)      surface downward flux SW for solar band#2 (W.m-2)
-c
-c fluxtop_lw(ngrid)         outgoing upward flux tota LW (thermal IR) (W.m-2)
-c fluxtop_sw(ngrid,1)       outgoing upward flux SW for solar band#1 (W.m-2)
-c fluxtop_sw(ngrid,2)       outgoing upward flux SW for solar band#2 (W.m-2)
-
-c   tau          Column total visible dust optical depth at each point
-c   aerosol(ngrid,nlayer,naerkind)    aerosol extinction optical depth
-c                         at reference wavelength "longrefvis" set
-c                         in dimradmars_h , in each layer, for one of
-c                         the "naerkind" kind of aerosol optical
-c                         properties.
 c=======================================================================
 c
@@ -189,17 +171,31 @@
       REAL,INTENT(IN) :: tsurf(ngrid)
       REAL,INTENT(IN) :: dist_sol,mu0(ngrid),fract(ngrid)
-      REAL,INTENT(OUT) :: dtlw(ngrid,nlayer),dtsw(ngrid,nlayer)
-      REAL,INTENT(OUT) :: fluxsurf_lw(ngrid), fluxtop_lw(ngrid)
-      REAL,INTENT(OUT) :: fluxsurf_sw(ngrid,2), fluxtop_sw(ngrid,2)
+      REAL,INTENT(OUT) :: dtlw(ngrid,nlayer) ! longwave (IR) heating rate (K/s)
+      REAL,INTENT(OUT) :: dtsw(ngrid,nlayer) ! shortwave (Solar) heating rate (K/s)
+      REAL,INTENT(OUT) :: fluxsurf_lw(ngrid) ! total LW (thermal IR) downward flux 
+                                             ! (W.m-2) at the surface
+      REAL,INTENT(OUT) :: fluxtop_lw(ngrid) ! outgoing total LW (thermal IR)
+                                     ! upward flux (W.m-2) at the top of the atm.
+      REAL,INTENT(OUT) :: fluxsurf_dn_sw(ngrid,2) ! surface downward SW flux for 
+                                      ! solar bands #1 and #2 (W.m-2)
+      REAL,INTENT(OUT) :: fluxsurf_up_sw(ngrid,2) ! surface upward SW flux for 
+                                      ! solar bands #1 and #2 (W.m-2)
+      REAL,INTENT(OUT) :: fluxtop_dn_sw(ngrid,2) ! incoming downward SW flux for 
+                                 ! solar bands #1 and #2 (W.m-2) at top of atm.
+      REAL,INTENT(OUT) :: fluxtop_up_sw(ngrid,2) ! outgoing upward SW flux for 
+                                 ! solar bands #1 and #2 (W.m-2) at top of atm.
       REAL,INTENT(OUT) :: tau_pref_scenario(ngrid) ! prescribed dust column
                           ! visible opacity at odpref from scenario
       REAL,INTENT(OUT) :: tau_pref_gcm(ngrid) ! computed dust column
                           ! visible opacity at odpref in the GCM
-      REAL,INTENT(OUT) :: tau(ngrid,naerkind)
+      REAL,INTENT(OUT) :: tau(ngrid,naerkind) ! Column visible optical depth
+                          ! for each aerosol
       REAL,INTENT(OUT) :: taucloudtes(ngrid)! Cloud opacity at infrared
                                !   reference wavelength using
                                !   Qabs instead of Qext
                                !   (direct comparison with TES)
-      REAL,INTENT(OUT) :: aerosol(ngrid,nlayer,naerkind)
+      REAL,INTENT(OUT) :: aerosol(ngrid,nlayer,naerkind) ! aerosol extinction
+                          ! optical depth at reference wavelength "longrefvis",
+                          ! set in dimradmars_h, for each kind of aerosol
       REAL,INTENT(INOUT) :: dsodust(ngrid,nlayer)
       REAL,INTENT(OUT) :: rdust(ngrid,nlayer)  ! Dust geometric mean radius (m)
@@ -585,9 +581,18 @@
         enddo
 
+        ! copy SW fluxes at surface and TOA
         do ig = 1,nd
-          fluxsurf_sw(ig0+ig,1) = zfluxd_sw(ig,1,1)
-          fluxsurf_sw(ig0+ig,2) = zfluxd_sw(ig,1,2)
-          fluxtop_sw(ig0+ig,1) = zfluxu_sw(ig,nlaylte+1,1)
-          fluxtop_sw(ig0+ig,2) = zfluxu_sw(ig,nlaylte+1,2)
+          ! surface downward SW flux
+          fluxsurf_dn_sw(ig0+ig,1) = zfluxd_sw(ig,1,1)
+          fluxsurf_dn_sw(ig0+ig,2) = zfluxd_sw(ig,1,2)
+          ! surface upward SW flux
+          fluxsurf_up_sw(ig0+ig,1) = zfluxu_sw(ig,1,1)
+          fluxsurf_up_sw(ig0+ig,2) = zfluxu_sw(ig,1,2)
+          ! downward SW flux at top of atmosphere
+          fluxtop_dn_sw(ig0+ig,1) = zfluxd_sw(ig,nlaylte+1,1)
+          fluxtop_dn_sw(ig0+ig,2) = zfluxd_sw(ig,nlaylte+1,2)          
+          ! upward SW flux at top of atmosphere
+          fluxtop_up_sw(ig0+ig,1) = zfluxu_sw(ig,nlaylte+1,1)
+          fluxtop_up_sw(ig0+ig,2) = zfluxu_sw(ig,nlaylte+1,2)
         enddo
 
@@ -613,5 +618,5 @@
             PRINT*,albedo(igout,1),emis(igout),mu0(igout),
      s           fract(igout), fluxsurf_lw(igout),
-     $     fluxsurf_sw(igout,1)+fluxsurf_sw(igout,2)
+     $     fluxsurf_dn_sw(igout,1)+fluxsurf_dn_sw(igout,2)
             PRINT*,'Tlay Tlev Play Plev dT/dt SW dT/dt LW (K/s)'
             PRINT*,'daysec',daysec
