Index: /LMDZ4/branches/LMDZ4-dev/libf/phylmd/physiq.F
===================================================================
--- /LMDZ4/branches/LMDZ4-dev/libf/phylmd/physiq.F	(revision 1158)
+++ /LMDZ4/branches/LMDZ4-dev/libf/phylmd/physiq.F	(revision 1159)
@@ -2941,5 +2941,5 @@
 
          CALL radlwsw_aero
-     e        (dist, rmu0, fract, solaire,
+     e        (dist, rmu0, fract, 
      e        paprs, pplay,zxtsol,albsol1, albsol2, 
      e        t_seri,q_seri,wo,
@@ -2948,4 +2948,5 @@
      e        tau_aero, piz_aero, cg_aero,
      e        cldtaupi,new_aod,
+     e        zqsat, flwc, fiwc,
      s        heat,heat0,cool,cool0,radsol,albpla,
      s        topsw,toplw,solsw,sollw,
Index: /LMDZ4/branches/LMDZ4-dev/libf/phylmd/radlwsw_aero.F90
===================================================================
--- /LMDZ4/branches/LMDZ4-dev/libf/phylmd/radlwsw_aero.F90	(revision 1158)
+++ /LMDZ4/branches/LMDZ4-dev/libf/phylmd/radlwsw_aero.F90	(revision 1159)
@@ -1,5 +1,5 @@
 SUBROUTINE radlwsw_aero( &
-   dist, rmu0, fract, solaire, &
-   paprs, pplay,tsol,albedo, alblw, &
+   dist, rmu0, fract, &
+   paprs, pplay,tsol,alb1, alb2, &
    t,q,wo,&
    cldfra, cldemi, cldtaupd,&
@@ -7,4 +7,5 @@
    tau_aero, piz_aero, cg_aero,&
    cldtaupi, new_aod, &
+   qsat, flwc, fiwc, &
    heat,heat0,cool,cool0,radsol,albpla,&
    topsw,toplw,solsw,sollw,&
@@ -32,9 +33,9 @@
   ! fract----input-R- duree d'ensoleillement normalisee
   ! co2_ppm--input-R- concentration du gaz carbonique (en ppm)
-  ! solaire--input-R- constante solaire (W/m**2)
   ! paprs----input-R- pression a inter-couche (Pa)
   ! pplay----input-R- pression au milieu de couche (Pa)
   ! tsol-----input-R- temperature du sol (en K)
-  ! albedo---input-R- albedo du sol (entre 0 et 1)
+  ! alb1-----input-R- albedo du sol(entre 0 et 1) dans l'interval visible 
+  ! alb2-----input-R- albedo du sol(entre 0 et 1) dans l'interval proche infra-rouge   
   ! t--------input-R- temperature (K)
   ! q--------input-R- vapeur d'eau (en kg/kg)
@@ -92,11 +93,12 @@
   include "YOETHF.h"
   include "YOMCST.h"
+  include "clesphys.h"
+  include "iniprint.h"
 
 ! Input arguments
-  REAL,    INTENT(in)  :: solaire
   REAL,    INTENT(in)  :: dist
   REAL,    INTENT(in)  :: rmu0(KLON), fract(KLON)
   REAL,    INTENT(in)  :: paprs(KLON,KLEV+1), pplay(KLON,KLEV)
-  REAL,    INTENT(in)  :: albedo(KLON), alblw(KLON), tsol(KLON)
+  REAL,    INTENT(in)  :: alb1(KLON), alb2(KLON), tsol(KLON)
   REAL,    INTENT(in)  :: t(KLON,KLEV), q(KLON,KLEV), wo(KLON,KLEV)
   LOGICAL, INTENT(in)  :: ok_ade, ok_aie                                 ! switches whether to use aerosol direct (indirect) effects or not
@@ -107,4 +109,7 @@
   REAL,    INTENT(in)  :: cldtaupi(KLON,KLEV)                            ! cloud optical thickness for pre-industrial aerosol concentrations
   LOGICAL, INTENT(in)  :: new_aod                                        ! flag pour retrouver les resultats exacts de l'AR4 dans le cas ou l'on ne travaille qu'avec les sulfates
+  REAL,    INTENT(in)  :: qsat(klon,klev) ! Variable pour iflag_rrtm=1
+  REAL,    INTENT(in)  :: flwc(klon,klev) ! Variable pour iflag_rrtm=1
+  REAL,    INTENT(in)  :: fiwc(klon,klev) ! Variable pour iflag_rrtm=1
 
 ! Output arguments
@@ -207,8 +212,8 @@
       zfract(i) = fract(iof+i)
       zrmu0(i) = rmu0(iof+i)
-      PALBD(i,1) = albedo(iof+i)
-      PALBD(i,2) = alblw(iof+i)
-      PALBP(i,1) = albedo(iof+i)
-      PALBP(i,2) = alblw(iof+i)
+      PALBD(i,1) = alb1(iof+i)
+      PALBD(i,2) = alb2(iof+i)
+      PALBP(i,1) = alb1(iof+i)
+      PALBP(i,2) = alb2(iof+i)
       PEMIS(i) = 1.0 
       PVIEW(i) = 1.66
@@ -281,54 +286,66 @@
       ENDDO
     ENDDO
-    !
-    !======================================================================
-    CALL LW_LMDAR4(&
-       PPMB, PDP,&
-       PPSOL,PDT0,PEMIS,&
-       PTL, PTAVE, PWV, POZON, PAER,&
-       PCLDLD,PCLDLU,&
-       PVIEW,&
-       zcool, zcool0,&
-       ztoplw,zsollw,ztoplw0,zsollw0,&
-       zsollwdown,&
-       ZFLUP, ZFLDN, ZFLUP0,ZFLDN0)
-
-
-    IF (.NOT. new_aod) THEN 
-       ! use old version
-        CALL SW_LMDAR4(PSCT, zrmu0, zfract,&
-           PPMB, PDP, &
-           PPSOL, PALBD, PALBP,&
-           PTAVE, PWV, PQS, POZON, PAER,&
-           PCLDSW, PTAU, POMEGA, PCG,&
-           zheat, zheat0,&
-           zalbpla,ztopsw,zsolsw,ztopsw0,zsolsw0,&
-           ZFSUP,ZFSDN,ZFSUP0,ZFSDN0,&
-           tau_aero(:,:,5,:), piz_aero(:,:,5,:), cg_aero(:,:,5,:),& 
-           PTAUA, POMEGAA,&
-           ztopswadaero,zsolswadaero,&
-           ztopswaiaero,zsolswaiaero,& 
-           ok_ade, ok_aie) 
-    ELSE
-
-        CALL SW_AEROAR4(PSCT, zrmu0, zfract,&
-           PPMB, PDP,&
-           PPSOL, PALBD, PALBP,&
-           PTAVE, PWV, PQS, POZON, PAER,&
-           PCLDSW, PTAU, POMEGA, PCG,&
-           zheat, zheat0,&
-           zalbpla,ztopsw,zsolsw,ztopsw0,zsolsw0,&
-           ZFSUP,ZFSDN,ZFSUP0,ZFSDN0,&
-           tauaero, pizaero, cgaero, &
-           PTAUA, POMEGAA,&
-           ztopswadaero,zsolswadaero,&
-           ztopswad0aero,zsolswad0aero,&
-           ztopswaiaero,zsolswaiaero, & 
-           ztopsw_aero,ztopsw0_aero,&
-           zsolsw_aero,zsolsw0_aero,&
-           ok_ade, ok_aie) 
-    
-    ENDIF
-    !======================================================================
+
+!
+!===== iflag_rrtm ================================================
+!      
+    IF (iflag_rrtm == 0) THEN
+       ! Old radiation scheme, used for AR4 runs
+       CALL LW_LMDAR4(&
+            PPMB, PDP,&
+            PPSOL,PDT0,PEMIS,&
+            PTL, PTAVE, PWV, POZON, PAER,&
+            PCLDLD,PCLDLU,&
+            PVIEW,&
+            zcool, zcool0,&
+            ztoplw,zsollw,ztoplw0,zsollw0,&
+            zsollwdown,&
+            ZFLUP, ZFLDN, ZFLUP0,ZFLDN0)
+
+
+       IF (.NOT. new_aod) THEN 
+          ! use old version
+          CALL SW_LMDAR4(PSCT, zrmu0, zfract,&
+               PPMB, PDP, &
+               PPSOL, PALBD, PALBP,&
+               PTAVE, PWV, PQS, POZON, PAER,&
+               PCLDSW, PTAU, POMEGA, PCG,&
+               zheat, zheat0,&
+               zalbpla,ztopsw,zsolsw,ztopsw0,zsolsw0,&
+               ZFSUP,ZFSDN,ZFSUP0,ZFSDN0,&
+               tau_aero(:,:,5,:), piz_aero(:,:,5,:), cg_aero(:,:,5,:),& 
+               PTAUA, POMEGAA,&
+               ztopswadaero,zsolswadaero,&
+               ztopswaiaero,zsolswaiaero,& 
+               ok_ade, ok_aie) 
+          
+       ELSE ! new_aod=T         
+          CALL SW_AEROAR4(PSCT, zrmu0, zfract,&
+               PPMB, PDP,&
+               PPSOL, PALBD, PALBP,&
+               PTAVE, PWV, PQS, POZON, PAER,&
+               PCLDSW, PTAU, POMEGA, PCG,&
+               zheat, zheat0,&
+               zalbpla,ztopsw,zsolsw,ztopsw0,zsolsw0,&
+               ZFSUP,ZFSDN,ZFSUP0,ZFSDN0,&
+               tauaero, pizaero, cgaero, &
+               PTAUA, POMEGAA,&
+               ztopswadaero,zsolswadaero,&
+               ztopswad0aero,zsolswad0aero,&
+               ztopswaiaero,zsolswaiaero, & 
+               ztopsw_aero,ztopsw0_aero,&
+               zsolsw_aero,zsolsw0_aero,&
+               ok_ade, ok_aie) 
+          
+       ENDIF
+
+    ELSE  
+!===== iflag_rrtm=1, on passe dans SW via RECMWFL ===============
+       WRITE(lunout,*) "Option iflag_rrtm=T ne fonctionne pas encore !!!"
+       CALL abort_gcm('radlwsw_aero','iflag_rrtm=T not valid',1) 
+
+    ENDIF ! iflag_rrtm
+!======================================================================
+
     DO i = 1, kdlon
       radsol(iof+i) = zsolsw(i) + zsollw(i)
@@ -405,9 +422,9 @@
       ENDDO
     ENDDO
-    !
-  ENDDO
-
-
-ENDSUBROUTINE radlwsw_aero
-
-
+
+ ENDDO ! j = 1, nb_gr
+
+
+END SUBROUTINE radlwsw_aero
+
+
Index: /LMDZ4/branches/LMDZ4-dev/libf/phylmd/sw_aeroAR4.F90
===================================================================
--- /LMDZ4/branches/LMDZ4-dev/libf/phylmd/sw_aeroAR4.F90	(revision 1158)
+++ /LMDZ4/branches/LMDZ4-dev/libf/phylmd/sw_aeroAR4.F90	(revision 1159)
@@ -114,13 +114,15 @@
   PARAMETER(swpas=1)
 
-  INTEGER itapsw
-  LOGICAL appel1er
-  DATA itapsw /0/
-  DATA appel1er /.TRUE./
-  SAVE itapsw,appel1er
+  INTEGER, SAVE :: itapsw = 0
+  !$OMP THREADPRIVATE(itapsw)
+  LOGICAL, SAVE :: appel1er = .TRUE.
   !$OMP THREADPRIVATE(appel1er)
-  !$OMP THREADPRIVATE(itapsw)
+  LOGICAL, SAVE :: initialized = .FALSE.
+  !$OMP THREADPRIVATE(initialized)
+
   !jq-Introduced for aerosol forcings
-  REAL*8 flag_aer
+  REAL*8, SAVE :: flag_aer
+  !$OMP THREADPRIVATE(flag_aer)
+
   LOGICAL ok_ade, ok_aie    ! use aerosol forcings or not?
   REAL*8 tauaero(kdlon,kflev,9,2)  ! aerosol optical properties
@@ -162,12 +164,4 @@
   REAL*8,ALLOCATABLE,SAVE ::  ZFSDN0_AERO(:,:,:)
   !$OMP THREADPRIVATE(ZFSDN0_AERO)
-
-  LOGICAL initialized
-  !rv
-  SAVE flag_aer
-  !$OMP THREADPRIVATE(flag_aer)
-  DATA initialized/.FALSE./
-  SAVE initialized
-  !$OMP THREADPRIVATE(initialized)
 
 
@@ -196,5 +190,4 @@
      ZFSDN0_AERO(:,:,:)=0.
   ENDIF
-  !rv
 
 
