Changeset 5791 for LMDZ6/branches/contrails/libf/phylmd/rrtm/radlsw.F90

Timestamp:

Jul 28, 2025, 7:23:15 PM (6 days ago)

Author:

aborella

Message:

Merge with trunk r5789

Location:

LMDZ6/branches/contrails

Files:

• . (modified) (1 prop)
• libf/phylmd/rrtm/radlsw.F90 (modified) (9 diffs)

LMDZ6/branches/contrails/libf/phylmd/rrtm/radlsw.F90

@@ -13 +13 @@
  & LRDUST,PPIZA_DST,PCGA_DST,PTAUREL_DST,&
  & PTAU_LW,&
- & PFLUX,PFLUC,PFSDN ,PFSUP , PFSCDN , PFSCUP)
+ & PFLUX,PFLUC,PFSDN ,PFSUP , PFSCDN , PFSCUP,&
+ & PTOAB, PTOACB )   !FC
 
 use write_field_phy
@@ -78 +79 @@
 ! PFLUX  : (KLON,2,KLEV+1)   ; LW total sky flux (1=up, 2=down)
 ! PFLUC  : (KLON,2,KLEV+1)   ; LW clear sky flux (1=up, 2=down)
+!FC PTOAB : (klon,NLW)           ; LW total sky TOA flux bande
+!FC PTOACB : (klon,NLW)          ; LW clear sky TOA flux bande
 ! PFSDN(KLON,KLEV+1)         ; SW total sky flux down
 ! PFSUP(KLON,KLEV+1)         ; SW total sky flux up
@@ -154 +157 @@
 USE YOMLUN_IFSAUX , ONLY : NULOUT
 USE YOMCT3        , ONLY : NSTEP
+!FC
+USE PARRRTM, ONLY : JPGPT, GP_PER_LWBAND, JPBAND
+
 
 IMPLICIT NONE
@@ -218 +224 @@
 REAL(KIND=JPRB)   ,INTENT(OUT)   :: PFLUX(KLON,2,KLEV+1) ! LW total sky flux (1=up, 2=down)
 REAL(KIND=JPRB)   ,INTENT(OUT)   :: PFLUC(KLON,2,KLEV+1) ! LW clear sky flux (1=up, 2=down)
+REAL(KIND=JPRB)   ,INTENT(OUT)   :: PTOAB(klon,NLW)  ! FC LW full sectral bande TOA
+REAL(KIND=JPRB)   ,INTENT(OUT)   :: PTOACB(klon,NLW) ! FC LW clear spectral TOA
 REAL(KIND=JPRB)   ,INTENT(OUT)   :: PFSDN(KLON,KLEV+1)   ! SW total sky flux down
 REAL(KIND=JPRB)   ,INTENT(OUT)   :: PFSUP(KLON,KLEV+1)   ! SW total sky flux up
@@ -270 +278 @@
  & , ZRAINT(KLON)       , ZRES(KLON)&
  & , ZTICE(KLON)        , ZEMIT(KLON),  ZBICFU(KLON)&
- & , ZKICFU(KLON)
+ & , ZKICFU(KLON)&
+ & , PTOAG(KLON,JPGPT),PTOACG(KLON,JPGPT) !FC
+
+
 REAL(KIND=JPRB) :: ZSUDU(KLON)   , ZPARF(KLON)       , ZUVDF(KLON), ZPARCF(KLON)
 INTEGER(KIND=JPIM) :: IKL, JK, JKL, JKLP1, JKP1, JL, JNU, JRTM, JSW, INDLAY
+
+INTEGER(KIND=JPIM) :: JLW , JI , J1, J2, J  !FC
+
 
 REAL(KIND=JPRB) :: ZASYMX, ZDIFFD, ZGI, ZGL, ZGR, ZIWGKG, ZLWGKG,&
@@ -318 +332 @@
   PFLUC(JL,1,KLEV+1) = 0.0_JPRB
   PFLUC(JL,2,KLEV+1) = 0.0_JPRB
+!FC
+  DO JI = 1, JPGPT
+  PTOAG(JL,JI) = 0.0_JPRB
+  PTOACG(JL,JI) = 0.0_JPRB
+  ENDDO
+  DO JI = 1, NLW
+  PTOAB(JL,JI) = 0.0_JPRB
+  PTOACB(JL,JI) = 0.0_JPRB
+  ENDDO
+!FC
+
   ZFSDNN(JL) = 0.0_JPRB
   ZFSDNV(JL) = 0.0_JPRB
@@ -1069 +1094 @@
 IF (.NOT.LPHYLIN) THEN
   IF ( .NOT. LRRTM) THEN
+           !FC ON EST DANS LE CAS OU NOT.RRTM EST TRUE DONC ON NE PASSE PAS LA
+
 
     CALL LW &
@@ -1114 +1141 @@
      & ZCLDSW  , ZTAUCLD,&
      & PTAU_LW,&
-     & ZEMIT , PFLUX , PFLUC , ZTCLEAR )
+     & ZEMIT , PFLUX , PFLUC , ZTCLEAR,&
+     & PTOAG, PTOACG )  !FC 
 !   print *,'RADLSW: apres CALL RRTM_RRTM_140GP'
 
@@ -1218 +1246 @@
 ENDDO
 
+!FC
+  DO JL = KIDIA,KFDIA
+    J1 = 1
+    J2 = 0
+    DO JI = 1, NLW  !NLW=16
+      J2 = J2 + GP_PER_LWBAND(JI)
+      DO J = J1, J2
+        PTOAB(JL,JI) = PTOAB(JL,JI) + PTOAG(JL,J)
+        PTOACB(JL,JI) = PTOACB(JL,JI) + PTOACG(JL,J)
+      ENDDO
+      J1 = J1 + GP_PER_LWBAND(JI)
+    ENDDO
+  ENDDO
+!FC
+
 DO JL = KIDIA,KFDIA
   PFRSOD(JL)=ZFSDWN(JL,1)