Ignore:
Timestamp:
Apr 4, 2014, 2:51:02 PM (11 years ago)
Author:
Laurent Fairhead
Message:

Nouvelle version qui inclut les effets des aérosols et propose les mêmes diagnostics des effets
directs et indirects que l'ancienne version du rayonnement.
OB


New RRTM version that includes the effects of aerosols and outputs the same direct and indirect effects
diagnostics as the old version
OB

Location:
LMDZ5/trunk/libf/phylmd/rrtm
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • LMDZ5/trunk/libf/phylmd/rrtm/rrtm_taumol1.F90

    r1990 r2003  
    189189
    190190IF (LHOOK) CALL DR_HOOK('RRTM_TAUMOL1',0,ZHOOK_HANDLE)
     191!--ajout OB
     192IF (K_LAYTROP.GT.100) THEN
     193PRINT *,'ATTENTION KLAY_TROP > 100 PROBLEME ARRAY DANS RRTM ON ARRETE'
     194STOP
     195!--fin ajout OB
     196ENDIF
    191197DO I_LAY = 1, K_LAYTROP
    192198  IND0(I_LAY) = ((K_JP(I_LAY)-1)*5+(K_JT(I_LAY)-1))*NSPA(1) + 1
  • LMDZ5/trunk/libf/phylmd/rrtm/swclr.F90

    r1990 r2003  
    55 & PRK   , PRMU0 , PTAUAZ, PTRA1 , PTRA2 , PTRCLR, &
    66!++MODIFCODE
    7   & LRDUST,PPIZA_DST, PCGA_DST, PTAUREL_DST )
     7  & LDDUST,PPIZA_DST, PCGA_DST, PTAU_DST )
    88!--MODIFCODE
    99
     
    5858!        Y.Seity 05-10-10 : add add 3 optional arg. for dust SW properties
    5959!        Y.Seity 06-09-09 : add modset from O.Thouron (MesoNH) under NOVLP tests
     60!        O.Boucher fev.2014: modification sur les aerosols pour utiliser les variables DST
    6061!     ------------------------------------------------------------------
    6162
     
    7172
    7273IMPLICIT NONE
    73 
    74 include "clesphys.h"
     74INCLUDE "clesphys.h"
    7575
    7676INTEGER(KIND=JPIM),INTENT(IN)    :: KLON
     
    8686REAL(KIND=JPRB)   ,INTENT(IN)    :: PSEC(KLON)
    8787!++MODIFCODE
    88 LOGICAL           ,INTENT(IN)    :: LRDUST                   ! flag for DUST
     88LOGICAL           ,INTENT(IN)    :: LDDUST                   ! flag for DUST
    8989REAL(KIND=JPRB)   ,INTENT(IN)    :: PPIZA_DST(KLON,KLEV)
    9090REAL(KIND=JPRB)   ,INTENT(IN)    :: PCGA_DST(KLON,KLEV)
    91 REAL(KIND=JPRB)   ,INTENT(IN)    :: PTAUREL_DST(KLON,KLEV)
     91REAL(KIND=JPRB)   ,INTENT(IN)    :: PTAU_DST(KLON,KLEV)
    9292!--MODIFCODE
    9393REAL(KIND=JPRB)   ,INTENT(OUT)   :: PCGAZ(KLON,KLEV)
     
    157157
    158158!++MODIFCODE 
     159!--OB on fait passer les aerosols LMDZ dans la variable DST
    159160  IF(NOVLP < 5)THEN !ECMWF VERSION
    160   DO JAE=1,6
     161!  DO JAE=1,6
    161162      DO JL = KIDIA,KFDIA
    162         PTAUAZ(JL,JK)=PTAUAZ(JL,JK)+PAER(JL,JAE,IKL)*RTAUA(KNU,JAE)
    163         PPIZAZ(JL,JK)=PPIZAZ(JL,JK)+PAER(JL,JAE,IKL)&
    164          & * RTAUA(KNU,JAE)*RPIZA(KNU,JAE) 
    165         PCGAZ(JL,JK) =  PCGAZ(JL,JK) +PAER(JL,JAE,IKL)&
    166          & * RTAUA(KNU,JAE)*RPIZA(KNU,JAE)*RCGA(KNU,JAE) 
     163!        PTAUAZ(JL,JK)=PTAUAZ(JL,JK)+PAER(JL,JAE,IKL)*RTAUA(KNU,JAE)
     164        PTAUAZ(JL,JK)=PTAU_DST(JL,IKL)
     165!        PPIZAZ(JL,JK)=PPIZAZ(JL,JK)+PAER(JL,JAE,IKL)&
     166!         & * RTAUA(KNU,JAE)*RPIZA(KNU,JAE) 
     167        PPIZAZ(JL,JK)=PTAU_DST(JL,IKL)*PPIZA_DST(JL,IKL)
     168!        PCGAZ(JL,JK) =  PCGAZ(JL,JK) +PAER(JL,JAE,IKL)&
     169!         & * RTAUA(KNU,JAE)*RPIZA(KNU,JAE)*RCGA(KNU,JAE) 
     170        PCGAZ(JL,JK)=PTAU_DST(JL,IKL)*PPIZA_DST(JL,IKL)*PCGA_DST(JL,IKL)
    167171      ENDDO
    168     ENDDO
     172!    ENDDO
    169173  ELSE ! MESONH VERSION
    170      DO JAE=1,6
     174!--OB on utilise directement les aerosols LMDZ
     175!     DO JAE=1,6
    171176        DO JL = KIDIA,KFDIA
    172177           !Special optical properties for dust
    173            IF (LRDUST.AND.(JAE==3)) THEN
     178!           IF (LDDUST.AND.(JAE==3)) THEN
    174179           !Ponderation of aerosol optical properties:first step
    175180           !ti
    176             PTAUAZ(JL,JK)=PTAUAZ(JL,JK) + PAER(JL,JAE,IKL) * PTAUREL_DST(JL,IKL)
     181!            PTAUAZ(JL,JK)=PTAUAZ(JL,JK) + PAER(JL,JAE,IKL) * PTAUREL_DST(JL,IKL)
     182            PTAUAZ(JL,JK)= PTAU_DST(JL,IKL)
    177183           !wi*ti
    178              PPIZAZ(JL,JK)=PPIZAZ(JL,JK) + PAER(JL,JAE,IKL)  &
    179                    & *PTAUREL_DST(JL,IKL)*PPIZA_DST(JL,IKL)
     184!             PPIZAZ(JL,JK)=PPIZAZ(JL,JK) + PAER(JL,JAE,IKL)  &
     185!                   & *PTAUREL_DST(JL,IKL)*PPIZA_DST(JL,IKL)
     186             PPIZAZ(JL,JK)=PTAU_DST(JL,IKL)*PPIZA_DST(JL,IKL)
    180187           !wi*ti*gi
    181              PCGAZ(JL,JK) = PCGAZ(JL,JK) + PAER(JL,JAE,IKL) &
    182                 &  *PTAUREL_DST(JL,IKL)*PPIZA_DST(JL,IKL)*PCGA_DST(JL,IKL)
     188!             PCGAZ(JL,JK) = PCGAZ(JL,JK) + PAER(JL,JAE,IKL) &
     189!                &  *PTAUREL_DST(JL,IKL)*PPIZA_DST(JL,IKL)*PCGA_DST(JL,IKL)
     190             PCGAZ(JL,JK) = PTAU_DST(JL,IKL)*PPIZA_DST(JL,IKL)*PCGA_DST(JL,IKL)
    183191           !wi*ti*(gi**2)
    184              ZFACOA_NEW(JL,JK)= ZFACOA_NEW(JL,JK)+PAER(JL, JAE, IKL)&
    185                 & *PTAUREL_DST(JL,IKL) *PPIZA_DST(JL,IKL)*PCGA_DST(JL,IKL)*&
     192!             ZFACOA_NEW(JL,JK)= ZFACOA_NEW(JL,JK)+PAER(JL, JAE, IKL)&
     193!                & *PTAUREL_DST(JL,IKL) *PPIZA_DST(JL,IKL)*PCGA_DST(JL,IKL)*&
     194!                & PCGA_DST(JL,IKL)
     195             ZFACOA_NEW(JL,JK)= ZFACOA_NEW(JL,JK)+&
     196                & PTAU_DST(JL,IKL) *PPIZA_DST(JL,IKL)*PCGA_DST(JL,IKL)*&
    186197                & PCGA_DST(JL,IKL)
    187            ELSE
     198!           ELSE
    188199           !Ponderation of aerosol optical properties:first step
    189200           !ti
    190              PTAUAZ(JL,JK)=PTAUAZ(JL,JK)+PAER(JL, JAE, IKL)*RTAUA(KNU,JAE)
     201!             PTAUAZ(JL,JK)=PTAUAZ(JL,JK)+PAER(JL, JAE, IKL)*RTAUA(KNU,JAE)
    191202           !wi*ti
    192              PPIZAZ(JL,JK)=PPIZAZ(JL,JK)+PAER(JL, JAE, IKL)&
    193                 &* RTAUA(KNU,JAE)*RPIZA(KNU,JAE)
     203!             PPIZAZ(JL,JK)=PPIZAZ(JL,JK)+PAER(JL, JAE, IKL)&
     204!                &* RTAUA(KNU,JAE)*RPIZA(KNU,JAE)
    194205           !wi*ti*gi
    195              PCGAZ(JL,JK) =  PCGAZ(JL,JK) +PAER(JL, JAE, IKL)&
    196                 &* RTAUA(KNU,JAE)*RPIZA(KNU,JAE)*RCGA(KNU,JAE)
     206!             PCGAZ(JL,JK) =  PCGAZ(JL,JK) +PAER(JL, JAE, IKL)&
     207!                &* RTAUA(KNU,JAE)*RPIZA(KNU,JAE)*RCGA(KNU,JAE)
    197208           !wi*ti*(gi**2)
    198              ZFACOA_NEW(JL,JK)= ZFACOA_NEW(JL,JK)+PAER(JL, JAE, IKL)&
    199                 &* RTAUA(KNU,JAE)*RPIZA(KNU,JAE)*RCGA(KNU,JAE)*RCGA(KNU,JAE)
    200            ENDIF
     209!             ZFACOA_NEW(JL,JK)= ZFACOA_NEW(JL,JK)+PAER(JL, JAE, IKL)&
     210!                &* RTAUA(KNU,JAE)*RPIZA(KNU,JAE)*RCGA(KNU,JAE)*RCGA(KNU,JAE)
     211!           ENDIF
    201212        ENDDO
    202      ENDDO
     213!     ENDDO
    203214  ENDIF
    204215!--MODIFCODE 
Note: See TracChangeset for help on using the changeset viewer.