Ignore:
Timestamp:
Jul 28, 2025, 7:23:15 PM (7 days ago)
Author:
aborella
Message:

Merge with trunk r5789

Location:
LMDZ6/branches/contrails
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • LMDZ6/branches/contrails

  • LMDZ6/branches/contrails/libf/phylmd/rrtm/rrtm_rrtm_140gp.F90

    r2146 r5791  
    4141 & PCLDF , PTAUCLD,&
    4242 & PTAU_LW,&
    43  & PEMIT , PFLUX , PFLUC, PTCLEAR &
    44  &
     43 & PEMIT , PFLUX , PFLUC, PTCLEAR, &
     44 & PTOAG, PTOACG) !FC
    4545
    4646! *** This program is the driver for RRTM, the AER rapid model. 
     
    8888REAL(KIND=JPRB)   ,INTENT(OUT)   :: PFLUC(KLON,2,KLEV+1) ! LW clear sky flux (1=up, 2=down)
    8989REAL(KIND=JPRB)   ,INTENT(OUT)   :: PTCLEAR(KLON) ! clear-sky fraction of column
     90
     91REAL(KIND=JPRB)   ,INTENT(OUT)   :: PTOAG(KLON,JPGPT) ! full-sky TOA G !FC
     92REAL(KIND=JPRB)   ,INTENT(OUT)   :: PTOACG(KLON,JPGPT) ! clear-sky TOA G !FC
     93
    9094INTEGER(KIND=JPIM) :: ICLDLYR(JPLAY)        ! Cloud indicator
    9195REAL(KIND=JPRB) :: Z_CLDFRAC(JPLAY)           ! Cloud fraction
     
    116120REAL(KIND=JPRB) :: Z_TOTUFLUX(0:JPLAY)
    117121
    118 INTEGER(KIND=JPIM) :: i, icld, iplon, I_K
     122REAL(KIND=JPRB) :: Z_TOAG(JPGPT) !FC
     123REAL(KIND=JPRB) :: Z_TOACG(JPGPT) !FC
     124
     125
     126INTEGER(KIND=JPIM) :: i, icld, iplon, I_K ,JI !FC
    119127INTEGER(KIND=JPIM) :: ISTART
    120128INTEGER(KIND=JPIM) :: IEND
     
    253261!  Clear radiative transfer is done for clear layers and cloudy radiative
    254262!  transfer is done for cloudy layers as identified by icldlyr.
     263!FC
    255264
    256265  CALL RRTM_RTRN1A_140GP (KLEV,ISTART,IEND,ICLDLYR,Z_CLDFRAC,Z_TAUCLD,Z_ABSS1,&
    257266   & Z_OD,Z_TAUSF1,Z_CLFNET,Z_CLHTR,Z_FNET,Z_HTR,Z_TOTDFLUC,Z_TOTDFLUX,Z_TOTUFLUC,Z_TOTUFLUX,&
     267   & Z_TOAG, Z_TOACG,&   
    258268   & Z_TAVEL,Z_PZ,Z_TZ,Z_TBOUND,Z_PFRAC,Z_SEMISS,Z_SEMISLW,IREFLECT) 
    259269
     
    270280    PFLUX(iplon,2,i+1) = -Z_TOTDFLUX(i)*Z_FLUXFAC
    271281  ENDDO
     282!FC
     283  DO JI = 1, JPGPT
     284   PTOAG(iplon,JI) = Z_TOAG(JI)*Z_FLUXFAC !FC faire attention aux niveaux (i+1 et i)
     285   PTOACG(iplon,JI) = Z_TOACG(JI)*Z_FLUXFAC
     286  ENDDO
     287!FC
    272288ENDDO
    273289
Note: See TracChangeset for help on using the changeset viewer.