Index: LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/aeropt_5wv_rrtm.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/aeropt_5wv_rrtm.F90	(revision 5158)
+++ LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/aeropt_5wv_rrtm.F90	(revision 5159)
@@ -1,3 +1,3 @@
-!
+
 ! $Id: aeropt_5wv_rrtm.F90 3288 2018-03-19 20:58:31Z oboucher $
 !
@@ -16,5 +16,5 @@
   USE YOMCST, ONLY: RD,RG
 
-  !
+
   !    Yves Balkanski le 12 avril 2006
   !    Celine Deandreis
@@ -23,37 +23,37 @@
   !    Olivier Boucher mars 2014 pour adaptation RRTM
   !    
-  !
+
   ! Refractive indices for seasalt come from Shettle and Fenn (1979)
-  !
+
   ! Refractive indices from water come from Hale and Querry (1973)
-  !
+
   ! Refractive indices from Ammonium Sulfate Toon and Pollack (1976)
-  !
+
   ! Refractive indices for Dust, internal mixture of minerals coated with 1.5% hematite 
   ! by Volume (Balkanski et al., 2006)
-  !
+
   ! Refractive indices for POM: Kinne (pers. Communication 
-  !
+
   ! Refractive index for BC from Shettle and Fenn (1979)
-  !
+
   ! Shettle, E. P., & Fenn, R. W. (1979), Models for the aerosols of the lower atmosphere and 
   ! the effects of humidity variations on their optical properties, U.S. Air Force Geophysics 
   ! Laboratory Rept. AFGL-TR-79-0214, Hanscomb Air Force Base, MA.
-  !
+
   ! Hale, G. M. and M. R. Querry, Optical constants of water in the 200-nm to 200-m 
   ! wavelength region, Appl. Opt., 12, 555-563, 1973.
-  !
+
   ! Toon, O. B. and J. B. Pollack, The optical constants of several atmospheric aerosol species:
   ! Ammonium sulfate, aluminum oxide, and sodium chloride, J. Geohys. Res., 81, 5733-5748,
   ! 1976.
-  !
+
   ! Balkanski, Y., M. Schulz, T. Claquin And O. Boucher, Reevaluation of mineral aerosol 
   ! radiative forcings suggests a better agreement with satellite and AERONET data, Atmospheric 
   ! Chemistry and Physics Discussions., 6, pp 8383-8419, 2006.
-  !
+
   IMPLICIT NONE
-  !
+
   ! Input arguments:
-  !
+
   REAL, DIMENSION(klon,klev), INTENT(IN)   :: pdel
   REAL, DIMENSION(klon,klev,naero_tot), INTENT(IN) :: m_allaer
@@ -63,14 +63,14 @@
   REAL, DIMENSION(klon,klev), INTENT(IN)   :: pplay
   REAL, DIMENSION(klon,klev), INTENT(IN)   :: t_seri
-  !
+
   ! Output arguments:
-  !
+
   REAL, DIMENSION(klon), INTENT(OUT)                      :: ai      ! POLDER aerosol index 
   REAL, DIMENSION(klon,nwave,naero_tot), INTENT(OUT)      :: tausum
   REAL, DIMENSION(klon,naero_tot), INTENT(OUT)            :: drytausum
   REAL, DIMENSION(klon,klev,nwave,naero_tot), INTENT(OUT) :: tau
-  !
+
   ! Local
-  !
+
   INTEGER, PARAMETER :: las = nwave_sw
   LOGICAL :: soluble
@@ -110,7 +110,7 @@
   REAL :: abs_aeri_5wv(las,naero_insoluble)           ! Abs. coeff. ** m2/g 
 
-  !
+
   ! BC internal mixture
-  !
+
   INTEGER, PARAMETER ::  nbclassbc = 6  ! Added by Rong Wang/OB for the 5 fractions
                                         ! of BC in the soluble mode:
@@ -123,7 +123,7 @@
   REAL :: abs_MG_5wv(nbre_RH,las,nbclassbc)
 
-  !
+
   ! Proprietes optiques
-  !
+
   REAL :: fact_RH(nbre_RH), BC_massfra
   INTEGER :: n, classbc
@@ -196,5 +196,5 @@
    ! Nitrate insoluble
   0.726, 0.753, 0.780, 0.797, 0.811 /
-!
+
  DATA abs_aers_5wv/ &
    ! absorption BC Accumulation Soluble (AS)
@@ -303,5 +303,5 @@
    4.505,  4.505,  4.505,  4.505,  4.520,  4.444,  4.356,  4.243,  4.089,  3.997,  3.912,  4.179, &
    4.295,  4.295,  4.295,  4.295,  4.307,  4.239,  4.157,  4.045,  3.876,  3.757,  3.602,  3.569  /
-!
+
    DATA abs_MG_5wv/ &
  !--BC content=0.001
@@ -405,5 +405,5 @@
   !      interpolate from Sext to retrieve Sext_at_gridpoint_per_species
   !      compute optical_thickness_at_gridpoint_per_species
-  !
+
   ! Calculations that need to be done since we are not in the subroutines INCA
   !      
Index: LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/aeropt_6bands_rrtm.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/aeropt_6bands_rrtm.F90	(revision 5158)
+++ LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/aeropt_6bands_rrtm.F90	(revision 5159)
@@ -1,5 +1,5 @@
-!
+
 ! $Id: aeropt_6bands_rrtm.F90 4165 2022-05-26 19:56:53Z oboucher $
-!
+
 SUBROUTINE AEROPT_6BANDS_RRTM ( &
      pdel, m_allaer, RHcl, &
@@ -20,9 +20,9 @@
   !    Olivier Boucher f�vrier 2014 pour passage � RRTM
   !    a partir des propri�t�s optiques fournies par Yves Balkanski
-  !
+
   IMPLICIT NONE
   !!
   ! Input arguments:
-  !
+
   REAL, DIMENSION(klon,klev),     INTENT(IN)  :: pdel
   REAL, DIMENSION(klon,klev,naero_tot),   INTENT(IN)  :: m_allaer
@@ -33,15 +33,15 @@
   REAL, DIMENSION(klon,klev),     INTENT(IN)  :: zrho
   LOGICAL,                        INTENT(IN)  :: ok_volcan ! volcanic diags
-  !
+
   ! Output arguments:
   ! 2= total aerosols 
   ! 1= natural aerosols
-  !
+
   REAL, DIMENSION(klon,klev,2,nbands_sw_rrtm), INTENT(OUT) :: tau_allaer ! epaisseur optique aerosol
   REAL, DIMENSION(klon,klev,2,nbands_sw_rrtm), INTENT(OUT) :: piz_allaer ! single scattering albedo aerosol
   REAL, DIMENSION(klon,klev,2,nbands_sw_rrtm), INTENT(OUT) :: cg_allaer  ! asymmetry parameter aerosol
-  !
+
   ! Local
-  !
+
   LOGICAL :: soluble
   INTEGER :: i, k,n, inu, m
@@ -77,7 +77,7 @@
   REAL, DIMENSION(klon,klev,id_ASBCM_phy:id_ASBCM_phy,nbands_sw_rrtm) :: piz_ae_pi
   REAL, DIMENSION(klon,klev,id_ASBCM_phy:id_ASBCM_phy,nbands_sw_rrtm) :: cg_ae_pi
-  !
+
   ! Proprietes optiques
-  !
+
   REAL:: alpha_aers_6bands(nbre_RH,nbands_sw_rrtm,naero_soluble)   !--unit m2/g
   REAL:: alpha_aeri_6bands(nbands_sw_rrtm,naero_insoluble)         !--unit m2/g
@@ -88,5 +88,5 @@
   ! 
   ! BC internal mixture
-  !
+
   INTEGER, PARAMETER ::  nbclassbc = 6  ! Added by Rong Wang/OB for the 5 fractions
                                         ! of BC in the soluble mode:
@@ -99,9 +99,9 @@
   REAL :: cg_MG_6bands(nbre_RH,nbands_sw_rrtm,nbclassbc)
   REAL :: piz_MG_6bands(nbre_RH,nbands_sw_rrtm,nbclassbc)
-  !
+
   INTEGER :: aerindex, classbc, classbc_pi
   REAL :: tmp_var, tmp_var_pi, BC_massfra, BC_massfra_pi
   CHARACTER*20 :: modname
-  !
+
   REAL, PARAMETER :: tau_min = 1.e-7
 
@@ -526,5 +526,5 @@
   ENDIF
 
-  !
+
   ! loop over modes, use of precalculated nmd and corresponding sigma
   !    loop over wavelengths
@@ -796,9 +796,9 @@
 
          IF (.NOT. ok_volcan) THEN
-!
+
 !--this is the default case
 !--in this case, index 1 of tau_allaer contains natural aerosols only
 !--because the objective is to perform the double radiation call with and without anthropogenic aerosols 
-!
+
            tau_allaer(i,k,1,inu)=tau_ae_pi(i,k,id_ASSO4M_phy,inu)+tau_ae_pi(i,k,id_CSSO4M_phy,inu)+ &
                                  tau_ae_pi(i,k,id_ASBCM_phy,inu)+tau_ae_pi(i,k,id_AIBCM_phy,inu)+   &
@@ -834,12 +834,12 @@
                                  (tau_allaer(i,k,1,inu)*piz_allaer(i,k,1,inu))
            cg_allaer(i,k,1,inu)=MIN(MAX(cg_allaer(i,k,1,inu),0.0),1.0)
-!
+
          ELSE
-!
+
 !--this is the case for VOLMIP
 !--in this case index 1 of tau_allaer contains all (natural+anthropogenic) aerosols (same as index 2 above)
 !--but stratospheric aerosols will not be added in rrtm/readaerosolstrato2 as
 !--the objective is to have the double radiation call with and without stratospheric aerosols
-!
+
            tau_allaer(i,k,1,inu)=tau_allaer(i,k,2,inu)
             
@@ -847,5 +847,5 @@
             
            cg_allaer(i,k,1,inu) =cg_allaer(i,k,2,inu)
-!
+
          ENDIF
         ENDDO
Index: LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/aeropt_lw_rrtm.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/aeropt_lw_rrtm.F90	(revision 5158)
+++ LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/aeropt_lw_rrtm.F90	(revision 5159)
@@ -1,7 +1,7 @@
-!
+
 ! aeropt_lw_rrtm.F90 2014-05-13 C. Kleinschmitt
 !                    2016-05-03 O. Boucher 
 !                    2016-12-17 O. Boucher
-!
+
 ! This routine feeds aerosol LW properties to RRTM
 ! we only consider absorption (not scattering)
@@ -20,17 +20,17 @@
   IMPLICIT NONE
 
-  !
+
   ! Input arguments:
-  !
+
   LOGICAL, INTENT(IN)                              :: ok_alw
   INTEGER, INTENT(IN)                              :: flag_aerosol
   REAL, DIMENSION(klon,klev), INTENT(IN)           :: pdel, zrho
   REAL, DIMENSION(klon,klev,naero_tot), INTENT(IN) :: m_allaer, m_allaer_pi
-  !
+
   INTEGER inu, i, k
   REAL :: zdh(klon,klev)
   REAL :: tmp_var, tmp_var_pi
   CHARACTER*20 modname
-  !
+
   !--absorption coefficient for CIDUST
   REAL:: alpha_abs_CIDUST_16bands(nbands_lw_rrtm)   !--unit m2/g 
@@ -38,7 +38,7 @@
   0.001, 0.003, 0.005, 0.006, 0.012, 0.030, 0.148, 0.098, &
   0.017, 0.053, 0.031, 0.008, 0.010, 0.011, 0.013, 0.015  /
-  !
+
   modname='aeropt_lw_rrtm'
-  !
+
   IF (NLW.NE.nbands_lw_rrtm) THEN
     CALL abort_physic(modname,'Erreur NLW doit etre egal a 16 pour cette routine',1)
@@ -46,17 +46,17 @@
   ! 
   IF (ok_alw) THEN                                   !--aerosol LW effects
-   !
+
    IF (flag_aerosol.EQ.5.OR.flag_aerosol.EQ.6.OR.flag_aerosol.EQ.7) THEN  !-Dust
-    !
+
     zdh(:,:)=pdel(:,:)/(RG*zrho(:,:))      ! m
-    !
+
     DO k=1, klev
       DO i=1, klon
-         !
+
          tmp_var   =m_allaer(i,k,id_CIDUSTM_phy)   /1.e6*zdh(i,k)  !--g/m2
          tmp_var_pi=m_allaer_pi(i,k,id_CIDUSTM_phy)/1.e6*zdh(i,k)  !--g/m2
-         !
+
          DO inu=1, NLW
-           !
+
            !--total aerosol
            tau_aero_lw_rrtm(i,k,2,inu) = MAX(1.e-15,tmp_var*alpha_abs_CIDUST_16bands(inu))
@@ -64,16 +64,16 @@
 !           tau_aero_lw_rrtm(:,:,1,inu) = MAX(1.e-15,tmp_var_pi*alpha_abs_CIDUST_16bands(inu))
            tau_aero_lw_rrtm(i,k,1,inu) = 1.e-15  !--test
-           !
+
          ENDDO
       ENDDO
-      !
+
     ENDDO
     ! 
    ENDIF
-   !
+
   ELSE !--no aerosol LW effects
-    !
+
     tau_aero_lw_rrtm = 1.e-15
   ENDIF
-  !
+
 END SUBROUTINE AEROPT_LW_RRTM
Index: LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/dates.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/dates.F90	(revision 5158)
+++ LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/dates.F90	(revision 5159)
@@ -1,7 +1,7 @@
 subroutine dates_demo
 ! --------------------------------------------------------------
-!
+
 ! Conseils a l'utilisateur:
-!
+
 ! 1. VOUS COMPILEZ LES ENTIERS EN 32 BITS:
 ! Utilisez alors les routines
@@ -17,5 +17,5 @@
 ! les parametres subsequents assurant que seuls des entiers
 ! representables en 32 bits y soient utilises.
-!
+
 ! 2. VOUS COMPILEZ LES ENTIERS EN 64 BITS:
 ! Vous pouvez alors utiliser toutes les routines ci-dessus
@@ -30,14 +30,14 @@
 ! - amqhmsree_vers_dj: Conversion date gr�gorienne (en un seul r�el) > date julienne.
 ! - dj_vers_amqhmsree: Conversion date julienne > date gr�gorienne (en un seul r�el).
-!
-! --------------------------------------------------------------
-!
+
+! --------------------------------------------------------------
+
 ! D�finition des dates employ�es ci-dessous:
-!
+
 ! Date julienne DJ:
 !       Elle est compos�e d'un r�el.
 !       R1: Ce r�el cro�t de 1 tous les jours,
 !               et vaut 2451545.0 le 1er janvier 2000 � 12 UTC.
-!
+
 ! Date gr�gorienne "en clair" AMQHMS:
 !       Elle est compos�e de 5 entiers et d'un r�el.
@@ -80,9 +80,9 @@
 character*200 clzue,clze,clech
 character *(*) cdtit
-!
+
 !-------------------------------------------------
 ! Date de validit�.
 !-------------------------------------------------
-!
+
 zs=0.
 zsssss=psssss/3600.
@@ -95,33 +95,33 @@
 call dj_vers_amqhms(zdj,ianv,imov,iquv,ihev,imiv,zsv) ! date gr�gorienne de validit�.
 if(pstati < 3600.) then
-!
+
 !-------------------------------------------------
 ! Ech�ance en minutes.
 !-------------------------------------------------
-!
+
     zech=pstati/60. ; clzue='mn'
 elseif(pstati < 259200.) then
-!
+
 !-------------------------------------------------
 ! Ech�ance en heures.
 !-------------------------------------------------
-!
+
     zech=pstati/3600. ; clzue='h'
 else
-!
+
 !-------------------------------------------------
 ! Ech�ance en jours.
 !-------------------------------------------------
-!
+
     zech=pstati/86400. ; clzue='j'
 endif
-!
+
 ! Affichage de l'echeance avec deux chiffres apres la virgule.
-!
+
 write(clze,fmt='(f9.2)') zech
-!
+
 ! Si l'echeance est voisine d'un entier a mieux que 10**-2 pres,
 ! on l'affiche au format entier.
-!
+
 if(clze(len_trim(clze)-2:len_trim(clze)) == '.00') then
     clze=clze(1:len_trim(clze)-3)
@@ -130,27 +130,27 @@
 ilze=len_trim(clze)
 clech=clze(1:ilze)//clzue
-!
+
 !-------------------------------------------------
 ! Titre 3, de type
 ! BASE 2000.01.15 00:00 +72H VALID 2000.01.18 15:00.
 !-------------------------------------------------
-!
+
 if(imi == 0 .and. imiv == 0) then
-!
+
 !-------------------------------------------------
 ! Les minutes de base et validit� sont nulles.
 ! On ne les affiche pas.
 !-------------------------------------------------
-!
+
     write(cdtit,fmt='(a,i2,a,i2.2,a,i4.4,a,i2.2,3a,i2,a,i2.2,a,i4.4,a,i2.2,a)')&
     &'BASE ',kqu,'.',kmo,'.',kan,' ',ihe,'h UTC + ',clech(1:len_trim(clech))&
     &,', VALID ',iquv,'.',imov,'.',ianv,' ',ihev,'h UTC'
 else
-!
+
 !-------------------------------------------------
 ! Les minutes de base ou validit� sont non nulles.
 ! On les affiche.
 !-------------------------------------------------
-!
+
     write(cdtit,fmt='(a,i2,a,i2.2,a,i4.4,a,i2.2,a,i2.2,3a,i2,a,i2.2,a,i4.4,a,i2.2,a,i2.2,a)')&
     &'BASE ',kqu,'.',kmo,'.',kan,' ',ihe,':',imi,' UTC + ',clech(1:len_trim(clech))&
@@ -207,40 +207,40 @@
 character*3 cljour(0:6)
 data cljour/'Dim','Lun','Mar','Mer','Jeu','Ven','Sam'/
-!
+
 !-------------------------------------------------
 ! Date courante � la f90.
 !-------------------------------------------------
-!
+
 clgol1=' '
 clgol2=' '
 clgol3=' '
 call date_and_time(clgol1,clgol2,clgol3,idatat)
-!
+
 !-------------------------------------------------
 ! clgol1 est du type "AAAAMMQQ".
 !-------------------------------------------------
-!
+
 read(clgol1,fmt='(i4,2i2)') kaaaa,kmm,kqq
-!
+
 !-------------------------------------------------
 ! clgol2 est du type "HHMMSS.SSS".
 !-------------------------------------------------
-!
+
 read(clgol2,fmt='(2i2)') khh,kmi
 read(clgol2(5:),fmt=*) zs
 kss=nint(zs)
 read(clgol1,fmt='(i8)') iaaaammqq
-!
+
 !-------------------------------------------------
 ! Jour de la semaine.
 !-------------------------------------------------
-!
+
 kjs=ijoursem(iaaaammqq)
 cdjs=cljour(kjs)
-!
+
 !-------------------------------------------------
 ! Date totale.
 !-------------------------------------------------
-!
+
 write(cddt,fmt='(i4.4,a,2(i2.2,a),2a,i2.2,a,i2.2,a,i2.2)') &
 &kaaaa,'_',kmm,'_',kqq,'_',cdjs,'_',khh,':',kmi,':',kss
@@ -253,8 +253,8 @@
 ! -------
 ! 1999-08-17, J.M. Piriou.
-!
+
 ! Modifications:
 ! --------------
-!
+
 ! --------------------------------------------------------------------------
 ! En entree:
@@ -284,16 +284,16 @@
 idate1=20000101
 idate2=kaaaa*10000+kmm*100+kqq
-!
+
 !-------------------------------------------------
 ! Nombre de jours �coul�s entre la date
 ! d'entr�e � 0h UTC et le 1er janvier 2000 � 0h UTC.
 !-------------------------------------------------
-!
+
 call ecartdj(idate1,idate2,iecart)
-!
+
 !-------------------------------------------------
 ! Date julienne.
 !-------------------------------------------------
-!
+
 pdj=2451545.0- 0.5 +real(iecart)+real(khh)/24. &
 & +real(kmn)/1440.+ps/86400.
@@ -306,8 +306,8 @@
 ! -------
 ! 94-10-31, J.M. Piriou.
-!
+
 ! Modifications:
 ! --------------
-!
+
 ! --------------------------------------------------------------------------
 ! En entree:
@@ -326,5 +326,5 @@
 ! En sortie:
 ! kdat2 date finale.
-!
+
 ! --------------------------------------------------------------------------
 ! Exemple: call DAPLUS(19940503,1,456,ires) fournira
@@ -373,5 +373,5 @@
 ! Cette routine est utilisable avec des entiers 32 bits ou 64 bits.
 ! +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
-!
+
 ! -------------------------------------------------
 ! Date d'arrivee au jour pres.
@@ -409,5 +409,5 @@
 ! si l'ecart entre les deux dates est inferieur a 2**31 secondes,
 ! soit 68 ans!...
-!
+
 ! Au-dela de cette duree, les entiers doivent etre 64 bits.
 ! +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
@@ -429,31 +429,31 @@
 INTEGER(KIND=4) :: KEC
 character*(*) cd1,cd2
-!
+
 ! -------------------------------------------------
 ! On lit les dates sur des entiers.
 ! -------------------------------------------------
-!
+
 read(cd1,fmt='(i8,3i2)') iamq1,ih1,im1,is1
-!
+
 ! -------------------------------------------------
 ! Calculs d'ecarts et de leur partition
 ! en multiples de 86400 et sous-multiples.
 ! -------------------------------------------------
-!
+
 isec=ih1*3600+im1*60+is1 ! nombre de secondes ecoulees depuis cd10h.
 idelta=kec+isec ! nombre de secondes entre cd10h et cd2.
 ireste=modulo(idelta,86400) ! nombre de secondes entre cd20h et cd2.
 iecjours=(idelta-ireste)/86400 ! nombre de jours entre cd10h et cd20h.
-!
+
 ! -------------------------------------------------
 ! Date d'arrivee au jour pres.
 ! -------------------------------------------------
-!
+
 call daplus(iamq1,1,iecjours,iamq2)
-!
+
 ! -------------------------------------------------
 ! Date finale a la seconde pres.
 ! -------------------------------------------------
-!
+
 ih2=ireste/3600
 ireste=ireste-3600*ih2
@@ -470,8 +470,8 @@
 ! -------
 ! 1999-08-17, J.M. Piriou.
-!
+
 ! Modifications:
 ! --------------
-!
+
 ! --------------------------------------------------------------------------
 ! En entree:
@@ -485,5 +485,5 @@
 ! ps    seconde
 ! --------------------------------------------------------------------------
-!
+
 !-------------------------------------------------
 ! Nombre de jours entre le 1er janvier 2000 � 0 UTC
@@ -508,16 +508,16 @@
 REAL(KIND=8) :: ZFRAC
 zecart=pdj-2451544.5
-!
+
 !-------------------------------------------------
 ! Nombre entier de jours.
 !-------------------------------------------------
-!
+
 zfrac=modulo(zecart, 1._8 )
 iecart=nint(zecart-zfrac)
-!
+
 !-------------------------------------------------
 ! Date gr�gorienne associ�e.
 !-------------------------------------------------
-!
+
 idate1=20000101
 call daplusj(idate1,iecart,idate2)
@@ -526,9 +526,9 @@
 kmm=mod(knouv,100)
 kaaaa=knouv/100
-!
+
 !-------------------------------------------------
 ! Calcul de des heure, minute et seconde.
 !-------------------------------------------------
-!
+
 zfrac=(zecart-real(iecart))*24.
 khh=int(zfrac)
@@ -544,8 +544,8 @@
 ! -------
 ! 2002-11, J.M. Piriou.
-!
+
 ! Modifications:
 ! --------------
-!
+
 ! --------------------------------------------------------------------------
 ! En entree:
@@ -561,15 +561,15 @@
 REAL(KIND=8) :: ZS
 INTEGER(KIND=4) :: iaaaa,imm,iqq,ihh,imn
-!
+
 !-------------------------------------------------
 ! Conversion gr�gorien julien; cible 5 entiers et un r�el.
 !-------------------------------------------------
-!
+
 call dj_vers_amqhms(pdj,iaaaa,imm,iqq,ihh,imn,zs)
-!
+
 !-------------------------------------------------
 ! On passe de ces 5 entiers et un r�el � un seul r�el.
 !-------------------------------------------------
-!
+
 pgrer=real(iaaaa)*10000.+real(imm)*100. &
 & + real(iqq)+real(ihh)/100. &
@@ -583,8 +583,8 @@
 ! -------
 ! 2002-11, J.M. Piriou.
-!
+
 ! Modifications:
 ! --------------
-!
+
 ! --------------------------------------------------------------------------
 ! En entree:
@@ -600,10 +600,10 @@
 REAL(KIND=8) :: ZS,zloc
 INTEGER(KIND=4) :: iaaaa,imm,iqq,ihh,imn,iloc
-!
+
 !-------------------------------------------------
 ! On passe de cette date gr�gorienne donn�e
 ! comme un seul r�el � 5 entiers et un r�el.
 !-------------------------------------------------
-!
+
 iloc=int(pgrer)
 iqq=mod(iloc,100)
@@ -617,9 +617,9 @@
 imn=mod(iloc,100)
 ihh=iloc/100
-!
+
 !-------------------------------------------------
 ! Conversion gr�gorien julien; cible 5 entiers et un r�el.
 !-------------------------------------------------
-!
+
 call amqhms_vers_dj(iaaaa,imm,iqq,ihh,imn,zs,pdj)
 end
@@ -631,8 +631,8 @@
 ! -------
 ! 97-01-09, J.M. Piriou.
-!
+
 ! Modifications:
 ! --------------
-!
+
 ! --------------------------------------------------------------------------
 ! En entree: kopt option de precision sur les dates:
@@ -703,8 +703,8 @@
 ! si l'ecart entre les deux dates est inferieur a 2**31 jours,
 ! soit 5879489 ans!...
-!
+
 ! Au-dela de cette duree, les entiers doivent etre 64 bits.
 ! +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
-!
+
 ! -------------------------------------------------
 ! Ecart entre les deux dates au jour pres.
@@ -742,5 +742,5 @@
 ! si l'ecart entre les deux dates est inferieur a 2**31 secondes,
 ! soit 68 ans!...
-!
+
 ! Au-dela de cette duree, les entiers doivent etre 64 bits.
 ! +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
@@ -759,22 +759,22 @@
 INTEGER(KIND=4) :: KECQ
 character*(*) cd1,cd2
-!
+
 ! -------------------------------------------------
 ! On lit les dates sur des entiers.
 ! -------------------------------------------------
-!
+
 read(cd1,fmt='(i8,3i2)') iamq1,ih1,im1,is1
 read(cd2,fmt='(i8,3i2)') iamq2,ih2,im2,is2
-!
+
 ! -------------------------------------------------
 ! Ecart entre les deux dates au jour pres.
 ! -------------------------------------------------
-!
+
 call ecartd(iamq1,iamq2,1,kecq)
-!
+
 ! -------------------------------------------------
 ! Ecart en secondes.
 ! -------------------------------------------------
-!
+
 kec=kecq*86400+(ih2-ih1)*3600+(im2-im1)*60+is2-is1
 end
@@ -786,8 +786,8 @@
 ! -------
 ! 92-05-27, J.M. Piriou.
-!
+
 ! Modifications:
 ! --------------
-!
+
 ! --------------------------------------------------------------------------
 ! En entree: kopt option de precision sur les dates:
@@ -837,8 +837,8 @@
 INTEGER(KIND=4) :: KOPT
 data idebm/0,31,59,90,120,151,181,212,243,273,304,334/
-!
+
 ! --------------------------------------------------------------------------
 ! **      1. Calcul du nb de jours separant ki2 du 1er janv 1900
-!
+
 ! *       1.1 Extraction des quantieme, mois et annee
 if(kopt == 1) then
@@ -909,5 +909,5 @@
 ! --------------------------------------------------------------------------
 ! **      2. Calcul du nb de jours separant ii1 du 1er janv 1900
-!
+
 ! *       2.1 Extraction des quantieme, mois et annee
 ii1=19000101
@@ -952,8 +952,8 @@
 ! -------
 ! 92-05-27, J.M. Piriou.
-!
+
 ! Modifications:
 ! --------------
-!
+
 ! --------------------------------------------------------------------------
 ! En entree: kopt option de precision sur les dates:
@@ -1004,5 +1004,5 @@
 ! --------------------------------------------------------------------------
 ! **   On determine la date approximative d'arrivee en annees decimales
-!
+
 if(kopt == 1) then
   ! Date de type AAAAMMQQ
@@ -1026,5 +1026,5 @@
 ! --------------------------------------------------------------------------
 ! **   On determine la date en clair ii2p associee a la date decimale
-!
+
 iaaaa=int(zarrdec)
 zarrdec=12.*(zarrdec-real(iaaaa))
@@ -1035,5 +1035,5 @@
 ! --------------------------------------------------------------------------
 ! **   On calcule le nombre de jours separant 19000101 de ii2p
-!
+
 call gregod(ii2p,1,igii2p)
 imod=mod(kgre,iconv)
@@ -1042,5 +1042,5 @@
 ! --------------------------------------------------------------------------
 ! **   On avance de iec jours par rapport a ii2p
-!
+
 ! *       L'annee est-elle bissextile?
 ! Une annee est bissextile ssi elle est
@@ -1071,5 +1071,5 @@
 ! --------------------------------------------------------------------------
 ! **   On met en forme la date finale
-!
+
 idat=iqq+imm*100+iaaaa*10000
 if(kopt == 2) then
@@ -1101,8 +1101,8 @@
 ! -------
 ! 94-10-31, J.M. Piriou.
-!
+
 ! Modifications:
 ! --------------
-!
+
 ! --------------------------------------------------------------------------
 ! En entree:
@@ -1133,8 +1133,8 @@
 ! -------
 ! 92-05-27, J.M. Piriou.
-!
+
 ! Modifications:
 ! --------------
-!
+
 ! --------------------------------------------------------------------------
 
Index: LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/eq_regions_mod.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/eq_regions_mod.F90	(revision 5158)
+++ LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/eq_regions_mod.F90	(revision 5159)
@@ -1,4 +1,4 @@
 module eq_regions_mod
-!
+
 !     Purpose.
 !     --------
@@ -7,5 +7,5 @@
 !           equal area and small diameter.
 !           the type.
-!
+
 !     Background.
 !     -----------
@@ -21,8 +21,8 @@
 !     points in an IFS gaussian grid and provide an optimal (i.e. exact) 
 !     distribution of grid points over regions.
-!
+
 !     The following copyright notice for the eq_regions package is included from 
 !     the original MatLab release.
-!
+
 !     +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
 !     + Release 1.10 2005-06-26                                                 +
@@ -50,13 +50,13 @@
 !     +                                                                         +
 !     +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
-!
+
 !     Author.
 !     -------
 !        George Mozdzynski *ECMWF*
-!
+
 !     Modifications.
 !     --------------
 !        Original : 2006-04-15
-!
+
 !--------------------------------------------------------------------------------
 
@@ -85,8 +85,8 @@
 
 subroutine eq_regions(N)
-!
+
 ! eq_regions uses the zonal equal area sphere partitioning algorithm to partition
 ! the surface of a sphere into N regions of equal area and small diameter.
-!
+
 integer(kind=jpim),intent(in) :: N
 integer(kind=jpim) :: n_collars,j
@@ -100,7 +100,7 @@
 if( N == 1 )then
 
-  !
+
   ! We have only one region, which must be the whole sphere.
-  !
+
   n_regions(1)=1
   n_regions_ns=1
@@ -108,17 +108,17 @@
 else
 
-  !
+
   ! Given N, determine c_polar
   ! the colatitude of the North polar spherical cap.
-  !
+
   c_polar = polar_colat(N)
-  !
+
   ! Given N, determine the ideal angle for spherical collars.
   ! Based on N, this ideal angle, and c_polar,
   ! determine n_collars, the number of collars between the polar caps.
-  !
+
   n_collars = num_collars(N,c_polar,ideal_collar_angle(N))
   n_regions_ns=n_collars+2
-  !
+
   ! Given N, c_polar and n_collars, determine r_regions,
   ! a list of the ideal real number of regions in each collar,
@@ -130,5 +130,5 @@
   allocate(r_regions(n_collars+2))
   call ideal_region_list(N,c_polar,n_collars,r_regions)
-  !
+
   ! Given N and r_regions, determine n_regions, a list of the natural number 
   ! of regions in each collar and the polar caps.
@@ -138,5 +138,5 @@
   ! n_regions[n_collars+2] is 1.
   ! The sum of n_regions is N.
-  !
+
   call round_to_naturals(N,n_collars,r_regions)
   deallocate(r_regions)
@@ -160,10 +160,10 @@
 
 function num_collars(N,c_polar,a_ideal) result(num_c)
-!
+
 !NUM_COLLARS The number of collars between the polar caps
-!
+
 ! Given N, an ideal angle, and c_polar,
 ! determine n_collars, the number of collars between the polar caps.
-!
+
 integer(kind=jpim),intent(in) :: N
 real(kind=jprb),intent(in) :: a_ideal,c_polar
@@ -180,9 +180,9 @@
 
 subroutine ideal_region_list(N,c_polar,n_collars,r_regions)
-!
+
 !IDEAL_REGION_LIST The ideal real number of regions in each zone
-!
+
 ! List the ideal real number of regions in each collar, plus the polar caps.
-!
+
 ! Given N, c_polar and n_collars, determine r_regions, a list of the ideal real 
 ! number of regions in each collar, plus the polar caps.
@@ -191,5 +191,5 @@
 ! r_regions[n_collars+2] is 1.
 ! The sum of r_regions is N.
-!
+
 integer(kind=jpim),intent(in) :: N,n_collars
 real(kind=jprb),intent(in) :: c_polar
@@ -201,8 +201,8 @@
 r_regions(1) = 1.0_jprb
 if( n_collars > 0 )then
-  !
+
   ! Based on n_collars and c_polar, determine a_fitting,
   ! the collar angle such that n_collars collars fit between the polar caps.
-  !
+
   a_fitting = (pi-2.0_jprb*c_polar)/float(n_collars)
   ideal_region_area = area_of_ideal_region(N)
@@ -218,10 +218,10 @@
 
 function ideal_collar_angle(N) result(ideal)
-!
+
 ! IDEAL_COLLAR_ANGLE The ideal angle for spherical collars of an EQ partition
-!
+
 ! IDEAL_COLLAR_ANGLE(N) sets ANGLE to the ideal angle for the
 ! spherical collars of an EQ partition of the unit sphere S^2 into N regions.
-!
+
 integer(kind=jpim),intent(in) :: N
 real(kind=jprb) :: ideal
@@ -231,7 +231,7 @@
 
 subroutine round_to_naturals(N,n_collars,r_regions)
-!
+
 ! ROUND_TO_NATURALS Round off a given list of numbers of regions
-!
+
 ! Given N and r_regions, determine n_regions, a list of the natural number 
 ! of regions in each collar and the polar caps.
@@ -241,5 +241,5 @@
 ! n_regions[n_collars+2] is 1.
 ! The sum of n_regions is N.
-!
+
 integer(kind=jpim),intent(in) :: N,n_collars
 real(kind=jprb),intent(in) :: r_regions(n_collars+2)
@@ -256,7 +256,7 @@
 
 function polar_colat(N) result(polar_c)
-!
+
 ! Given N, determine the colatitude of the North polar spherical cap.
-!
+
 integer(kind=jpim),intent(in) :: N
 real(kind=jprb) :: area
@@ -272,8 +272,8 @@
 
 function area_of_ideal_region(N) result(area)
-!
+
 ! AREA_OF_IDEAL_REGION(N) sets AREA to be the area of one of N equal
 ! area regions on S^2, that is 1/N times AREA_OF_SPHERE.
-!
+
 integer(kind=jpim),intent(in) :: N
 real(kind=jprb) :: area_of_sphere
@@ -285,8 +285,8 @@
 
 function sradius_of_cap(area) result(sradius)
-!
+
 ! SRADIUS_OF_CAP(AREA) returns the spherical radius of
 ! an S^2 spherical cap of area AREA.
-!
+
 real(kind=jprb),intent(in) :: area
 real(kind=jprb) :: sradius
@@ -296,11 +296,11 @@
 
 function area_of_collar(a_top, a_bot) result(area)
-!
+
 ! AREA_OF_COLLAR Area of spherical collar
-!
+
 ! AREA_OF_COLLAR(A_TOP, A_BOT) sets AREA to be the area of an S^2 spherical 
 ! collar specified by A_TOP, A_BOT, where A_TOP is top (smaller) spherical radius,
 ! A_BOT is bottom (larger) spherical radius.
-!
+
 real(kind=jprb),intent(in) :: a_top,a_bot
 real(kind=jprb) area
@@ -310,10 +310,10 @@
 
 function area_of_cap(s_cap) result(area)
-!
+
 ! AREA_OF_CAP Area of spherical cap
-!
+
 ! AREA_OF_CAP(S_CAP) sets AREA to be the area of an S^2 spherical
 ! cap of spherical radius S_CAP.
-!
+
 real(kind=jprb),intent(in) :: s_cap
 real(kind=jprb) area
Index: LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/gfl_subs.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/gfl_subs.F90	(revision 5158)
+++ LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/gfl_subs.F90	(revision 5159)
@@ -582,16 +582,16 @@
 
 !SUBROUTINE DEACT_CLOUD_GFL  ! commente par MPL 10.12.08 (et REACT_CLOUD_GFL)
-!
+
 !**** *DEACT_CLOUD_GFL* Deactivate prognostic cloud variables
-!
+
 !     ------------------------------------------------------------------
-!
+
 !INTEGER(KIND=JPIM) :: JGFL
 !REAL(KIND=JPRB) :: ZHOOK_HANDLE
-!
+
 !#include "suslb.intfb.h"
-!
+
 !IF (LHOOK) CALL DR_HOOK('GFL_SUBS:DEACT_CLOUD_GFL',0,ZHOOK_HANDLE)
-!
+
 !IF (.NOT.L_CLD_DEACT .AND. &
 ! & (YL%LACTIVE .OR. YI%LACTIVE .OR. &
@@ -611,5 +611,5 @@
 !    IF (YA%LGP) YGFL%NUMGPFLDS=YGFL%NUMGPFLDS-1
 !    IF (YCPF%LGP) YGFL%NUMGPFLDS=YGFL%NUMGPFLDS-1
-!
+
 !    IF (YL%LT5) YGFL%NUMFLDS5=YGFL%NUMFLDS5-1
 !    IF (YI%LT5) YGFL%NUMFLDS5=YGFL%NUMFLDS5-1
@@ -618,5 +618,5 @@
 !    IF (YA%LT5) YGFL%NUMFLDS5=YGFL%NUMFLDS5-1
 !    IF (YCPF%LT5) YGFL%NUMFLDS5=YGFL%NUMFLDS5-1
-!
+
 !    CALL FALSIFY_GFLC(YL)
 !    CALL FALSIFY_GFLC(YI)
@@ -643,19 +643,19 @@
 ! ENDDO
 ! CALL SUSLB
-!
+
 ! L_CLD_DEACT=.TRUE.
 ! WRITE(NULOUT,*)' CLOUD FIELDS DE-ACTIVATAD, YGFL%NUMGPFLDS=', &
 ! & YGFL%NUMGPFLDS,' YGFL%NUMFLDS_SL1=', YGFL%NUMFLDS_SL1
 !ENDIF
-!
+
 !IF (LHOOK) CALL DR_HOOK('GFL_SUBS:DEACT_CLOUD_GFL',1,ZHOOK_HANDLE)
-!
+
 !END SUBROUTINE DEACT_CLOUD_GFL
-!
+
 !!=========================================================================
-!
+
 !SUBROUTINE REACT_CLOUD_GFL
 !!**** *REACT_CLOUD_GFL* Reactivate prognostic cloud variables
-!
+
 !INTEGER(KIND=JPIM) :: JGFL
 !REAL(KIND=JPRB) :: ZHOOK_HANDLE
@@ -664,5 +664,5 @@
 !!     ------------------------------------------------------------------
 !IF (LHOOK) CALL DR_HOOK('GFL_SUBS:REACT_CLOUD_GFL',0,ZHOOK_HANDLE)
-!
+
 !IF (L_CLD_DEACT) THEN
 !  LLGPL = YL%LGP
@@ -673,9 +673,9 @@
 !  CALL COPY_GFLC_GFLC(YI,YI_SAVE)
 !  CALL COPY_GFLC_GFLC(YA,YA_SAVE)
-!
+
 !  IF (.NOT. LLGPL .AND. YL%LGP) YGFL%NUMGPFLDS=YGFL%NUMGPFLDS+1
 !  IF (.NOT. LLGPI .AND. YI%LGP) YGFL%NUMGPFLDS=YGFL%NUMGPFLDS+1
 !  IF (.NOT. LLGPA .AND. YA%LGP) YGFL%NUMGPFLDS=YGFL%NUMGPFLDS+1
-!
+
 !  YGFL%NUMFLDS_SL1 = 0
 !  DO JGFL=1,YGFL%NUMFLDS
@@ -688,12 +688,12 @@
 !  ENDDO
 !  CALL SUSLB
-!
+
 !  L_CLD_DEACT=.FALSE.
 !  WRITE(NULOUT,*)' CLOUD FIELDS RE-ACTIVATAD, YGFL%NUMGPFLDS=', &
 !  & YGFL%NUMGPFLDS,' YGFL%NUMFLDS_SL1=', YGFL%NUMFLDS_SL1
 !ENDIF
-!
+
 !IF (LHOOK) CALL DR_HOOK('GFL_SUBS:REACT_CLOUD_GFL',1,ZHOOK_HANDLE)
-!
+
 !!     ------------------------------------------------------------------
 !END SUBROUTINE REACT_CLOUD_GFL
Index: LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/gppref.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/gppref.F90	(revision 5158)
+++ LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/gppref.F90	(revision 5159)
@@ -21,5 +21,5 @@
 !                              PRESH(KPROMA,0:KFLEV) - HALF LEVEL PRESSURE
 !                              PRESF(KPROMA,KFLEV)   - FULL LEVEL PRESSURE
-!
+
 !        Implicit arguments :  NONE.
 !        --------------------
@@ -105,10 +105,10 @@
 ! assumption that the top level input for pressure is 0 hPa.
 ! This restriction is only necessary in the case of use of NDLNPR=1.
-!
+
 ! LVERTFE : .T./.F. Finite element/conventional vertical discretisation.
 ! NDLNPR  : NDLNPR=0: conventional formulation of delta, i.e. ln(P(l)/P(l-1)).
 !           NDLNPR=1: formulation of delta used in non hydrostatic model,
 ! LAPRXPK : way of computing full-levels pressures in primitive equation
-!
+
 LVERTFE=.TRUE.    !!!!! A REVOIR (MPL) comment faut-il vraiment calculer PRESF ?
 
Index: LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/lwu.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/lwu.F90	(revision 5158)
+++ LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/lwu.F90	(revision 5159)
@@ -1,5 +1,5 @@
-!
+
 ! $Id$
-!
+
 SUBROUTINE LWU &
         & (KIDIA, KFDIA, KLON, KLEV, &
Index: LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/namphy0.h
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/namphy0.h	(revision 5158)
+++ LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/namphy0.h	(revision 5159)
@@ -36,5 +36,5 @@
 !     For ACCVIMPGY 
   &,ALFX,TCTC,TVFC,GAMAP1,RKDN,VVN,VVX,FENTRT,HCMIN,FQLIC,FNEBC,FEVAPC &
-!
+
   &,RDPHIC,GWBFAUT,RWBF1,RWBF2,RAUITN,RAUITX,RAUIUSTE &
   &,RSMDNEBX,RSMDTX,NSMTPA,NSMTPB
Index: LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/parkind1.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/parkind1.F90	(revision 5158)
+++ LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/parkind1.F90	(revision 5159)
@@ -1,12 +1,12 @@
 MODULE PARKIND1
-!
+
 !     *** Define usual kinds for strong typing ***
-!
+
 IMPLICIT NONE
 SAVE
-!
+
 !     Integer Kinds
 !     -------------
-!
+
 INTEGER, PARAMETER :: JPIT = SELECTED_INT_KIND(2)
 INTEGER, PARAMETER :: JPIS = SELECTED_INT_KIND(4)
@@ -22,12 +22,12 @@
 #endif
 
-!
+
 !     Real Kinds
 !     ----------
-!
+
 INTEGER, PARAMETER :: JPRT = SELECTED_REAL_KIND(2,1)
 INTEGER, PARAMETER :: JPRS = SELECTED_REAL_KIND(4,2)
 INTEGER, PARAMETER :: JPRM = SELECTED_REAL_KIND(6,37)
 INTEGER, PARAMETER :: JPRB = SELECTED_REAL_KIND(13,300)
-!
+
 END MODULE PARKIND1
Index: LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/parkind2.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/parkind2.F90	(revision 5158)
+++ LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/parkind2.F90	(revision 5159)
@@ -1,17 +1,17 @@
 MODULE PARKIND2
-!
+
 !     *** Define huge kinds for strong typing ***
-!
+
 IMPLICIT NONE
 SAVE
-!
+
 !     Integer Kinds
 !     -------------
-!
+
 INTEGER, PARAMETER :: JPIH = SELECTED_INT_KIND(18)
-!
+
 !     Real Kinds
 !     ----------
-!
+
 #ifdef REALHUGE
 INTEGER, PARAMETER :: JPRH = SELECTED_REAL_KIND(31,291)
@@ -19,4 +19,4 @@
 INTEGER, PARAMETER :: JPRH = SELECTED_REAL_KIND(13,300)
 #endif
-!
+
 END MODULE PARKIND2
Index: LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/radlsw.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/radlsw.F90	(revision 5158)
+++ LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/radlsw.F90	(revision 5159)
@@ -1048,5 +1048,5 @@
 
 !     ------------------------------------------------------------------
-!
+
 !*         2.7    DIFFUSIVITY FACTOR OR SATELLITE VIEWING ANGLE
 !                 ---------------------------------------------
Index: LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/read_rsun_rrtm.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/read_rsun_rrtm.F90	(revision 5158)
+++ LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/read_rsun_rrtm.F90	(revision 5159)
@@ -3,5 +3,5 @@
 !****************************************************************************************
 ! This routine will read the solar constant fraction per band
-!
+
 ! Olivier Boucher with inputs from Marion Marchand 
 !****************************************************************************************
Index: LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/readaerosol_optic_rrtm.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/readaerosol_optic_rrtm.F90	(revision 5158)
+++ LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/readaerosol_optic_rrtm.F90	(revision 5159)
@@ -1,4 +1,4 @@
 ! $Id: readaerosol_optic_rrtm.F90 4124 2022-04-08 14:47:04Z dcugnet $
-!
+
 SUBROUTINE readaerosol_optic_rrtm(debut, aerosol_couple, ok_alw, ok_volcan, &
      flag_aerosol, flag_bc_internal_mixture, itap, rjourvrai, &
@@ -99,11 +99,11 @@
   !    
   !****************************************************************************************
-  !
-  !
+
+
   IF (aerosol_couple) THEN   !--we get aerosols from tr_seri array from INCA
-     !
+
      !--copy fields from INCA tr_seri 
      !--convert to ug m-3 unit for consistency with offline fields
-     !
+
      itr = 0
      DO iq = 1,nqtot
@@ -142,5 +142,5 @@
      nitrcoarse(:,:)   =   tr_seri(:,:,id_CSNO3M)                        *zrho(:,:)*1.e9  ! CSNO3M
      nitrinscoarse(:,:)=   tr_seri(:,:,id_CINO3M)                        *zrho(:,:)*1.e9  ! CINO3M
-     !
+
      bcsol_pi(:,:)        =   0.0 ! ASBCM pre-ind
      pomsol_pi(:,:)       =   0.0 ! ASPOMM pre-ind
@@ -156,7 +156,7 @@
      nitrcoarse_pi(:,:)   =   0.0 ! CSNO3M pre-ind
      nitrinscoarse_pi(:,:)=   0.0 ! CINO3M
-     !
+
   ELSE !--not aerosol_couple
-     !
+
      ! Read and interpolate sulfate
      IF ( flag_aerosol .EQ. 1 .OR. flag_aerosol .EQ. 6 .OR. flag_aerosol .EQ. 7 ) THEN 
@@ -219,5 +219,5 @@
         loaddust=0.
      ENDIF
-     !
+
      ! Read and interpolate asno3m, csno3m, cino3m
      IF (flag_aerosol .EQ. 6 .OR. flag_aerosol .EQ. 7 ) THEN 
@@ -237,5 +237,5 @@
         loadno3(:)=0.0
      ENDIF
-     !
+
      ! CSSO4M is set to 0 as not reliable
      sulfcoarse(:,:)      =   0.0 ! CSSO4M (=SO4) + CSMSAM (=MSA) 
@@ -244,7 +244,7 @@
   ENDIF !--not aerosol_couple
 
-  !
+
   ! Store all aerosols in one variable
-  !
+
   m_allaer(:,:,id_ASBCM_phy)  = bcsol(:,:)        ! ASBCM
   m_allaer(:,:,id_ASPOMM_phy) = pomsol(:,:)       ! ASPOMM
@@ -278,5 +278,5 @@
   m_allaer_pi(:,:,id_STRAT_phy)  = 0.0
 
-  !
+
   ! Calculate the total mass of all soluble aersosols
   ! to be revisited for AR6
@@ -286,5 +286,5 @@
   !****************************************************************************************
   ! 2) Calculate optical properties for the aerosols
-  !
+
   !****************************************************************************************
   DO k = 1, klev
Index: LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/readaerosolstrato1_rrtm.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/readaerosolstrato1_rrtm.F90	(revision 5158)
+++ LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/readaerosolstrato1_rrtm.F90	(revision 5159)
@@ -1,3 +1,3 @@
-!
+
 ! $Id: readaerosolstrato1_rrtm.F90 2526 2016-05-26 22:13:40Z oboucher $
 !
@@ -57,14 +57,14 @@
     DATA cg_sw_strat   /0.6997170, 0.6810035, 0.7403592, 0.7562674, 0.6676504, 0.3478689/
     DATA piz_sw_strat  /0.9999998, 0.9999998, 1.000000000, 0.9999958, 0.9977155, 0.4510679/
-!
+
 !--diagnostics AOD in the SW
 ! alpha_sw_strat_wave is *not* normalised by the 550 nm extinction coefficient
     REAL, DIMENSION(nwave_sw) :: alpha_sw_strat_wave
     DATA alpha_sw_strat_wave/3.708007,4.125824,4.136584,3.887478,3.507738/
-!
+
 !--diagnostics AOD in the LW at 10 um (not normalised by the 550 nm ext coefficient
     REAL :: alpha_lw_strat_wave(nwave_lw)
     DATA alpha_lw_strat_wave/0.2746812/
-!
+
     REAL, DIMENSION(nbands_lw_rrtm) :: alpha_lw_abs_rrtm
     DATA alpha_lw_abs_rrtm/   8.8340312E-02, 6.9856711E-02, 6.2652975E-02, 5.7188231E-02, &
@@ -178,9 +178,9 @@
     
     IF (is_mpi_root.AND.is_omp_root) THEN
-!
+
     DEALLOCATE(tauaerstrat)
     DEALLOCATE(tauaerstrat_mois)
     DEALLOCATE(tauaerstrat_mois_glo)
-!
+
     ENDIF !--is_mpi_root and is_omp_root
 
Index: LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/readaerosolstrato2_rrtm.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/readaerosolstrato2_rrtm.F90	(revision 5158)
+++ LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/readaerosolstrato2_rrtm.F90	(revision 5159)
@@ -1,5 +1,5 @@
-!
+
 ! $Id: readaerosolstrato2_rrtm.F90 2526 2016-05-26 22:13:40Z oboucher $
-!
+
 SUBROUTINE readaerosolstrato2_rrtm(debut, ok_volcan)
 
@@ -341,5 +341,5 @@
 
     IF (.NOT. ok_volcan) THEN
-!
+
 !--this is the default case 
 !--stratospheric aerosols are added to both index 2 and 1 for double radiation calls
@@ -367,7 +367,7 @@
     ENDWHERE
     ENDDO
-!
+
     ELSE
-!
+
 !--this is the VOLMIP case
 !--stratospheric aerosols are only added to index 2 in this case 
@@ -410,7 +410,7 @@
       ENDWHERE
     ENDDO
-!
+
     ELSE
-!
+
 !--this is the VOLMIP case
     DO band=1, NLW
@@ -418,5 +418,5 @@
 !--and we copy index 2 in index 1 because we want the same dust aerosol LW properties as above
       tau_aero_lw_rrtm(:,:,1,band)  = tau_aero_lw_rrtm(:,:,2,band)
-!
+
       WHERE (stratomask.GT.0.999999)
 !--stratospheric aerosols are only added to index 2
Index: LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/recmwf_aero.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/recmwf_aero.F90	(revision 5158)
+++ LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/recmwf_aero.F90	(revision 5159)
@@ -1,5 +1,5 @@
-!
+
 ! $Id: recmwf_aero.F90 4875 2024-03-26 10:29:06Z lguez $
-!
+
 !OPTIONS XOPT(NOEVAL)
 SUBROUTINE RECMWF_AERO (KST, KEND, KPROMA, KTDIA , KLEV,&
@@ -464,5 +464,5 @@
 
   !*         4.1     CALL TO ACTUAL RADIATION SCHEME
-  !
+
   !----now we make multiple calls to the radiation according to which 
   !----aerosol flags are on
Index: LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/rrtm_ecrt_140gp.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/rrtm_ecrt_140gp.F90	(revision 5158)
+++ LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/rrtm_ecrt_140gp.F90	(revision 5159)
@@ -1,5 +1,5 @@
-!
+
 ! $Id$
-!
+
 !****************** SUBROUTINE RRTM_ECRT_140GP **************************
 
Index: LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/sdl_module.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/sdl_module.F90	(revision 5158)
+++ LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/sdl_module.F90	(revision 5159)
@@ -43,5 +43,5 @@
 !CHARACTER(LEN=*), PARAMETER :: CLNECMSG = '*** Calling NEC traceback ***'
 !#endif
-!
+
 !IF (PRESENT(KTID)) THEN
 !  ITID = KTID
@@ -49,5 +49,5 @@
 !  ITID = OML_MY_THREAD()
 !ENDIF
-!
+
 !IF (LHOOK) THEN
 !  IPRINT_OPTION = 2
@@ -55,5 +55,5 @@
 !  CALL C_DRHOOK_PRINT(0, ITID, IPRINT_OPTION, ILEVEL) ! from drhook.c
 !ENDIF
-!
+
 !#ifdef VPP
 !  CALL ERRTRA
@@ -117,9 +117,9 @@
 !MPL 4.12.08
 !#ifdef VPP
-!
+
 !CALL VPP_ABORT()
-!
+
 !#else
-!
+
 !IRETURN_CODE=1
 !CALL MPI_ABORT(KCOMM,IRETURN_CODE,IERROR)
Index: LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/srtm_cmbgb16.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/srtm_cmbgb16.F90	(revision 5158)
+++ LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/srtm_cmbgb16.F90	(revision 5159)
@@ -1,9 +1,9 @@
 SUBROUTINE SRTM_CMBGB16
 
-!
+
 !  Original version:       Michael J. Iacono; July, 1998
 !  Revision for RRTM_SW:   Michael J. Iacono; November, 2002
 !  Revision for RRTMG_SW:  Michael J. Iacono; December, 2003
-!
+
 !  The subroutines CMBGB16->CMBGB29 input the absorption coefficient
 !  data for each band, which are defined for 16 g-points and 14 spectral
@@ -12,7 +12,7 @@
 !  function data in array SFLUXREF are combined without weighting.  All
 !  g-point reduced data are put into new arrays for use in RRTMG_SW.
-!
+
 !  BAND 16:  2600-3250 cm-1 (low key- H2O,CH4; high key - CH4)
-!
+
 !-----------------------------------------------------------------------
 
Index: LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/srtm_srtm_224gp.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/srtm_srtm_224gp.F90	(revision 5158)
+++ LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/srtm_srtm_224gp.F90	(revision 5159)
@@ -1,5 +1,5 @@
-!
+
 ! $Id$
-!
+
 SUBROUTINE SRTM_SRTM_224GP &
         & (KIDIA, KFDIA, KLON, KLEV, KSW, KOVLP, &
Index: LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/su_aerop.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/su_aerop.F90	(revision 5158)
+++ LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/su_aerop.F90	(revision 5159)
@@ -71,5 +71,5 @@
 ! For sea-salt (_SS), 3 bins are considered      (0.03, 0.50, 5.0, 20.)
 ! For desert dust (_DD), 3 bins are considered     (0.03, 0.55, 0.9, 20.)
-!
+
 ! IF BIN LIMITS ARE CHANGED, MAKE SURE THAT THE RELEVANT SEDIMENTATION SPEEDS ARE 
 ! RECOMPUTED ACCORDINGLY
Index: LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/su_aerp.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/su_aerp.F90	(revision 5158)
+++ LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/su_aerp.F90	(revision 5159)
@@ -107,11 +107,11 @@
 !RVDPOCE = (/   0.1_JPRB,   1.2_JPRB,   0.1_JPRB,   0.1_JPRB,   0.1_JPRB,   0.1_JPRB &
 !          &,   0.1_JPRB,   1.2_JPRB,   1.2_JPRB,   1.2_JPRB,   1.5_JPRB,   1.5_JPRB /)
-!
+
 !RVDPSIC = (/   0.1_JPRB,   1.2_JPRB,   0.1_JPRB,   0.1_JPRB,   0.1_JPRB,   0.1_JPRB &
 !          &,   0.1_JPRB,   1.2_JPRB,   1.2_JPRB,   1.2_JPRB,   1.5_JPRB,   1.5_JPRB /)
-!
+
 !RVDPLND = (/   0.1_JPRB,   1.2_JPRB,   0.1_JPRB,   0.1_JPRB,   0.1_JPRB,   0.1_JPRB &
 !          &,   0.1_JPRB,   1.2_JPRB,   1.2_JPRB,   1.2_JPRB,   1.5_JPRB,   1.5_JPRB /)
-!
+
 !RVDPLIC = (/   0.1_JPRB,   1.2_JPRB,   0.1_JPRB,   0.1_JPRB,   0.1_JPRB,   0.1_JPRB &
 !          &,   0.1_JPRB,   1.2_JPRB,   1.2_JPRB,   1.2_JPRB,   1.5_JPRB,   1.5_JPRB /)
Index: LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/suclopn.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/suclopn.F90	(revision 5158)
+++ LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/suclopn.F90	(revision 5159)
@@ -1002,5 +1002,5 @@
 !    RASWCE(JNU)=ZASWCE14(JNU)
 !    RASWCF(JNU)=ZASWCF14(JNU)*1.E-03_JPRB
-!
+
 !    REBCUA(JNU)=ZEBCUA14(JNU)
 !    REBCUB(JNU)=ZEBCUB14(JNU)
@@ -1009,5 +1009,5 @@
 !    REBCUE(JNU)=ZEBCUE14(JNU)
 !    REBCUF(JNU)=ZEBCUF14(JNU)
-!
+
 !    RYFWCA(JNU)=ZYFWCA14(JNU)
 !    RYFWCB(JNU)=ZYFWCB14(JNU)
@@ -1016,5 +1016,5 @@
 !    RYFWCE(JNU)=ZYFWCE14(JNU)
 !    RYFWCF(JNU)=ZYFWCF14(JNU)
-!
+
 !    RSUSHE(JNU)=ZSUSHE14(JNU)*1.E-02_JPRB
 !    RSUSHF(JNU)=ZSUSHF14(JNU)*1.E-02_JPRB
@@ -1053,5 +1053,5 @@
 !!    RFUDD2(JNU)=ZFUDD214(JNU)
 !!    RFUDD3(JNU)=ZFUDD314(JNU)
-!
+
 !    PRINT *,'SUCLOPN: 14-SPECTRAL INTERVALS --> RRTM_SW'     
 !  ENDDO
Index: LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/suecrad.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/suecrad.F90	(revision 5158)
+++ LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/suecrad.F90	(revision 5159)
@@ -1,5 +1,5 @@
-!
+
 ! $Id: suecrad.F90 4251 2022-09-20 00:22:43Z fhourdin $
-!
+
 SUBROUTINE SUECRAD (KULOUT, KLEV, PETAH)
 
Index: LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/sugfl.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/sugfl.F90	(revision 5158)
+++ LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/sugfl.F90	(revision 5159)
@@ -104,5 +104,5 @@
 
 !      1.   CASE LFPART2=F
-!
+
 !      1.1  Initial settings.
 
Index: LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/suphec.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/suphec.F90	(revision 5158)
+++ LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/suphec.F90	(revision 5159)
@@ -129,7 +129,7 @@
 
 IF (LHOOK) CALL DR_HOOK('SUPHEC',0,ZHOOK_HANDLE)
-!
+
   IF (OK_BAD_ECMWF_THERMO) THEN
-!
+
      ! Modify constants defined in suphel.F90 and set RVTMP2 to 0.
      ! CALL GSTATS(1811,0) ! MPL 28.11.08
@@ -169,5 +169,5 @@
      ! Keep constants defined in suphel.F90
      RTICE=RTT-23._JPRB
-!
+
   ENDIF  ! (OK_BAD_ECMWF_THERMO)
 
Index: LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/suphy.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/suphy.F90	(revision 5158)
+++ LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/suphy.F90	(revision 5159)
@@ -34,5 +34,5 @@
 !     SUPHYFL
 !     SUHLPH
-!
+
 !     Reference.
 !     ----------
@@ -99,5 +99,5 @@
 print *,'---- SUPHY: avant SUPHMF'
 CALL SUPHMF(KULOUT)
-!
+
 print *,'---- SUPHY: avant SUGFL'
 !SUGFL: Set up unified_treatment grid-point fields
@@ -117,5 +117,5 @@
 ! Commente par MPL 20.11.08
 !CALL SUHLPH(KULOUT)    
-!
+
 !     ------------------------------------------------------------------
 
Index: LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/suphy0.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/suphy0.F90	(revision 5158)
+++ LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/suphy0.F90	(revision 5159)
@@ -341,5 +341,5 @@
 RWBF1=300._JPRB
 RWBF2=4._JPRB
-!
+
 RAUITN=233.15_JPRB
 RAUITX=263.15_JPRB
Index: LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/surdi.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/surdi.F90	(revision 5158)
+++ LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/surdi.F90	(revision 5159)
@@ -1,5 +1,5 @@
-!
+
 ! $Id$
-!
+
 SUBROUTINE SURDI
 
Index: LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/surhcri.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/surhcri.F90	(revision 5158)
+++ LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/surhcri.F90	(revision 5159)
@@ -2,40 +2,40 @@
       SUBROUTINE SURHCRI(KULOUT)
 !-----------------------------------------------------------------------
-!
+
 !**** *SURHCRI * - COMPUTATION OF THE CRTITICAL RELATIVE HUMIDITY
 !                  PROFILE FOR SMITH'S CONDENSATION SCHEME.
-!
+
 !**   Interface.
 !     ----------
 !        *CALL* *SURHCRI*
-!
+
 !-----------------------------------------------------------------------
-!
+
 ! -   ARGUMENTS D'ENTREE./INPUT ARGUMENTS.
 !     ------------------------------------
-!
+
 !-----------------------------------------------------------------------
-!
+
 ! -   ARGUMENTS IMPLICITES.
 !     ---------------------
-!
+
 ! COMMON /YOMPHY/
 ! COMMON /YOMPHY0/
-!
+
 !*
 !-----------------------------------------------------------------------
-!
+
 !     Auteur.
 !     -------
 !         05-03, Luc Gerard (from Ph. Lopez acrhcri)
-!
+
 !     Modifications.
 !     --------------
 !         06-10, nettoyage - R. Brozkova
 !-----------------------------------------------------------------------
-!
+
 USE PARKIND1  ,ONLY : JPIM     ,JPRB
 USE YOMHOOK   ,ONLY : LHOOK,   DR_HOOK
-!
+
 USE YOMPHY0   , ONLY : RHCRIT1, RHCRIT2, RETAMIN, GRHCMOD, RHCRI ,NRHCRI 
 USE YOMDIM    , ONLY : NFLEVG
Index: LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/susrtm.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/susrtm.F90	(revision 5158)
+++ LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/susrtm.F90	(revision 5159)
@@ -40,5 +40,5 @@
 ! minimize the effect on the resulting fluxes and cooling rates, and
 ! caution should be used if the mapping is modified.
-!
+
 ! JPGPT   The total number of new g-points (NGPT)
 ! NGC     The number of new g-points in each band
Index: LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/yoeaeratm.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/yoeaeratm.F90	(revision 5158)
+++ LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/yoeaeratm.F90	(revision 5159)
@@ -27,7 +27,7 @@
 ! RMFMIN     : minimum mass flux for convective aerosol transport
 ! RMASSE     : Molar mass: N.B.: either g/mol or Avogadro number
-!
+
 ! REPSCAER   : security on aerosol concentration: always >= 1.E-15
-!
+
 ! LAERCLIMG  : .T. to start prognostic aerosols with geographical monthly 
 !                  mean climatology
Index: LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/yoewcou.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/yoewcou.F90	(revision 5158)
+++ LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/yoewcou.F90	(revision 5159)
@@ -39,8 +39,8 @@
 !     *RNORTW*   REAL       NORTH BOUNDARY OF THE WAVE MODEL.
 !     *RDEGREW*  REAL       RESOLUTION OF THE WAVE MODEL (DEGREES).
-!
+
 !     *MASK_WAVE_IN*  INTEGER  COMMS MASK FOR INPUT TO WAVE MODEL
 !     *MASK_WAVE_OUT* INTEGER  COMMS MASK FOR OUTPUT FROM WAVE MODEL
-!
+
 !     *LWVIN_MASK_NOT_SET* LOGICAL indicates whether mask_wave_in
 !                           has been updated on the first call to the 
@@ -64,5 +64,5 @@
 !     *MWVIN_SENDIND* INTEGER global indexes of data on remote tasks that
 !                           the local task needs
-!
+
 !     *MWVIN_RECVOFF* INTEGER nproc sized array containing offsets into 
 !                           the MWVIN_RECVBUF array
Index: LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/yomct0.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/yomct0.F90	(revision 5158)
+++ LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/yomct0.F90	(revision 5159)
@@ -463,5 +463,5 @@
 LOGICAL ::  LPC_OLD
 LOGICAL ::  LPC_NESC
-!
+
 ! * FORCING
 LOGICAL ::  LSFORC
Index: LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/yomgrb.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/yomgrb.F90	(revision 5158)
+++ LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/yomgrb.F90	(revision 5159)
@@ -280,5 +280,5 @@
 ! NGRBCO2B         - 210068 CO2 - biosphere flux
 ! NGRBCO2A         - 210069 CO2 - anthropogenic emissions
-!
+
 !---------------------------------------------------
 ! NGRBGRG(JPGRG)   - 210121 GRG1: Nitrogen dioxide
Index: LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/yomphy0.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/yomphy0.F90	(revision 5158)
+++ LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/yomphy0.F90	(revision 5159)
@@ -352,5 +352,5 @@
 !    TFVS     : VITESSE DE CHUTE DES PRECIPITATIONS SOLIDES.
 !    GRHCMOD  : MODULATION IN CRITICAL RELATIVE HUMIDITY COMPUTATION.
-!
+
 ! Pseudo prognostic TKE scheme
 !    NUPTKE   : TUNABLE VALUE FOR PSEUDO TKE SCHEME FOLLOWING
@@ -374,14 +374,14 @@
 
 !   GCVALMX   : MAXIMUM ACCEPTABLE VALUE FOR TOTAL MESH FRACTION
-!
+
 !              Ascent properties
 !      ECMNPI : ECMNP for ice, the original ECMNP being kept for liquid.
 !      GFRIC  : INVERSE OF CHARACTERISTIC TIME for ICE CONDENSATION in ud.
-!
+
 !              Squeezing:
 !   GCVSQDN   : threshold value of sigma_d*q_cd*dp
 !   GCVSQDR   : fraction of the max to consider for squeezing
 !   GCVSQDCX  : maximum acceptable compression (<1)
-!
+
 !              Downdraught:
 !   GDDWPF    : Influence of rain fall velocity on downdraught
@@ -389,5 +389,5 @@
 !   GDDBETA   : DOWNDRAUGHT EXPLICIT DETRAINMENT COEFFICIENT
 !   TENTRD    : DOWNDRAUGHT ENTRAINMENT RATE (S^2/M^2)
-!
+
 !              Intensive Precipitation:
 !   GRRMINA   : MINIMUM REALISTIC PRECIPITATING MESH FRACTION
@@ -671,5 +671,5 @@
 !     For ALARO-0 :
 !     ------------------------------------------------------------------
-!
+
 ! RDPHIC   : REFERENCE GEOPOTENTIAL FOR CLOUDINESS ADJUSTMENT.
 ! GWBFAUT  : GAIN FOR THE WEGENER BERGERON FINDEISEN PROCESS IN ACAUTO
@@ -677,5 +677,5 @@
 ! RWBF2    : SECOND TUNING CONSTANT FOR BERGERON FINDEISEN PROCESS
 !               IN ACPLUIE_PROG
-!
+
 ! RAUITN   : TEMPERATURE LEVELS FOR RQICRMIN
 ! RAUITX   : TEMPERATURE LEVELS FOR RQICRMAX
@@ -687,7 +687,7 @@
 ! NSMTPA   : NUMBER OF LEVELS BELOW Tt LEVEL FOR LSMTPS SMOOTHING
 ! NSMTPB   : NUMBER OF LEVELS ABOVE Tt LEVEL FOR LSMTPS SMOOTHING
-!
+
 ! ------------------------------------------------------------------
-!
+
 REAL(KIND=JPRB) :: RDPHIC
 REAL(KIND=JPRB) :: GWBFAUT
Index: LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/yomphy3.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/yomphy3.F90	(revision 5158)
+++ LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/yomphy3.F90	(revision 5159)
@@ -88,5 +88,5 @@
 !               : AT THE UPPER CASE OF PADE FUNCTIONS FOR GASES.
 ! Parameters for cloud model:
-!
+
 !   Notations:
 !     g      - asymmetry factor            (unscaled)
@@ -98,10 +98,10 @@
 !     iwc    - ice water content
 !     lwc    - liquid water content
-!
+
 !   First index of FCM arrays (FCM = Fitting parameters for Cloud Model)
 !   denotes spectral band:
 !     1      - solar
 !     2      - thermal
-!
+
 !   FCM_DEL_A(2)    : Critical value of delta0 for computation of c_abs.
 !   FCM_DEL_D(2)    : Critical value of delta0 for computation of c_scat.
