--- trunk/phylmd/Radlwsw/sw.f90	2019/11/26 20:47:46	345
+++ trunk/phylmd/Radlwsw/sw.f90	2019/12/09 20:15:29	346
@@ -5,9 +5,9 @@
 contains
 
   SUBROUTINE SW(PSCT, PRMU0, PFRAC, PPMB, PDP, PPSOL, PALBD, PALBP, PTAVE, &
-       PWV, PQS, POZON, PCLDSW, PTAU, POMEGA, PCG, PHEAT, PHEAT0, PALBPLA, &
-       PTOPSW, PSOLSW, PTOPSW0, PSOLSW0, ZFSUP, ZFSDN, ZFSUP0, ZFSDN0, &
-       PTOPSWAD, PSOLSWAD, ok_ade)
+       PWV, PQS, POZON, PCLDSW, PTAU, POMEGA, PCG, PHEAT, PHEAT0, PTOPSW, &
+       PSOLSW, PTOPSW0, PSOLSW0, ZFSUP, ZFSDN, ZFSUP0, ZFSDN0, PTOPSWAD, &
+       PSOLSWAD, ok_ade)
 
     ! Purpose.
     ! This routine computes the shortwave radiation fluxes in two
@@ -38,7 +38,7 @@
 
     ! ARGUMENTS:
 
-    DOUBLE PRECISION PSCT ! constante solaire (valeur conseillee: 1370)
+    DOUBLE PRECISION PSCT ! constante solaire
     DOUBLE PRECISION PRMU0(KDLON) ! COSINE OF ZENITHAL ANGLE
     DOUBLE PRECISION PFRAC(KDLON) ! fraction de la journee
     DOUBLE PRECISION PPMB(KDLON, KFLEV+1) ! HALF-LEVEL PRESSURE (MB)
@@ -56,7 +56,6 @@
     DOUBLE PRECISION PCG(KDLON, 2, KFLEV) ! ASYMETRY FACTOR
     DOUBLE PRECISION PHEAT(KDLON, KFLEV) ! SHORTWAVE HEATING (K/DAY)
     DOUBLE PRECISION PHEAT0(KDLON, KFLEV)! SHORTWAVE HEATING (K/DAY) clear-sky
-    DOUBLE PRECISION PALBPLA(KDLON) ! PLANETARY ALBEDO
     DOUBLE PRECISION PTOPSW(KDLON) ! SHORTWAVE FLUX AT T.O.A.
     DOUBLE PRECISION PSOLSW(KDLON) ! SHORTWAVE FLUX AT SURFACE
     DOUBLE PRECISION PTOPSW0(KDLON) ! SHORTWAVE FLUX AT T.O.A. (CLEAR-SKY)
@@ -201,8 +200,6 @@
        ENDDO
     ENDDO
     DO i = 1, KDLON
-       PALBPLA(i) = ZFSUP(i, KFLEV+1)/(ZFSDN(i, KFLEV+1)+1.0e-20)
-
        PSOLSW(i) = ZFSDN(i, 1) - ZFSUP(i, 1)
        PTOPSW(i) = ZFSDN(i, KFLEV+1) - ZFSUP(i, KFLEV+1)