Changeset 5903 for LMDZ6/trunk/libf/phylmd/rrtm

Timestamp:

Dec 3, 2025, 6:01:05 PM (3 months ago)

Author:

fcheruy

Message:

correction les flux spectraux

File:

• LMDZ6/trunk/libf/phylmd/rrtm/recmwf_aero.F90 (modified) (18 diffs)

LMDZ6/trunk/libf/phylmd/rrtm/recmwf_aero.F90

@@ -349 +349 @@
   REAL(KIND=JPRB) ::  LWUP0_AERO(KPROMA,KLEV+1,5)
   REAL(KIND=JPRB) ::  LWDN0_AERO(KPROMA,KLEV+1,5)
+  !--LW SPEC FC) :
+  REAL(KIND=JPRB) ::  PTOAB_AERO(KPROMA,NLW,5)
+  REAL(KIND=JPRB) ::  PTOACB_AERO(KPROMA,NLW,5)
+
 
 #include "radlsw.intfb.h"
@@ -374 +378 @@
   LWUP0_AERO(:,:,:)=0.
   LWDN0_AERO(:,:,:)=0.
+
+  !FC
+  PTOAB_AERO (:,:,:) =0.
+  PTOACB_AERO (:,:,:)=0.
 
   PTAU_ZERO(:,:,:) =1.e-15
@@ -473 +481 @@
   !----aerosol flags are on
 
+!  write(*,*) ' FC flag_aerosol flag_aerosol_strat '
+!  write(*,*) flag_aerosol, flag_aerosol_strat 
+!  write(*,*) ' FC ok_ade ok_aie swaero_diag '
+!  write(*,*) ok_ade, ok_aie ,swaero_diag 
+
   IF (flag_aerosol .GT. 0 .OR. flag_aerosol_strat) THEN
 
@@ -501 +514 @@
              & PFSDN , PFSUP  , PFSCDN  , PFSCUP ,&
              & PTOAB, PTOACB  )   !FC
+!        write(*,*) 'FC cas 1'
 
         !* SAVE VARIABLES IN INTERIM VARIABLES A LA SW_AEROAR4
@@ -514 +528 @@
         LWUP_AERO(:,:,1) = PFLUX(:,1,:)
         LWDN_AERO(:,:,1) = PFLUX(:,2,:)
+
+        PTOAB_AERO(:,:,1)  = PTOAB(:,:) 
+        PTOACB_AERO(:,:,1)  = PTOACB(:,:) 
 
      ENDIF
@@ -543 +560 @@
              & PTOAB,  PTOACB )  !FC
 
+!        write(*,*) 'FC cas 2 ok_ade ok_aie swaero_diag', ok_ade , ok_aie , swaero_diag
         !* SAVE VARIABLES IN INTERIM VARIABLES A LA SW_AEROAR4
         ZFSUP0_AERO(:,:,2) = PFSCUP(:,:)
@@ -555 +573 @@
         LWUP_AERO(:,:,2) = PFLUX(:,1,:)
         LWDN_AERO(:,:,2) = PFLUX(:,2,:)
+
+        PTOAB_AERO(:,:,2)  = PTOAB(:,:) 
+        PTOACB_AERO(:,:,2)  = PTOACB(:,:) 
 
      ENDIF ! ok_aie      
@@ -584 +605 @@
              & PTOAB,  PTOACB )   !FC
 
+!        write(*,*) 'FC cas 3 ok_ade ok_aie swaero_diag' ,ok_ade , ok_aie , swaero_diag
         !* SAVE VARIABLES IN INTERIM VARIABLES A LA SW_AEROAR4
         ZFSUP0_AERO(:,:,3) = PFSCUP(:,:)
@@ -596 +618 @@
         LWUP_AERO(:,:,3) = PFLUX(:,1,:)
         LWDN_AERO(:,:,3) = PFLUX(:,2,:)
+
+        PTOAB_AERO(:,:,3)  = PTOAB(:,:) 
+        PTOACB_AERO(:,:,3)  = PTOACB(:,:) 
 
      ENDIF !-end ok_ade
@@ -624 +649 @@
              & PTOAB, PTOACB )   !FC
 
+!        write(*,*) 'FC cas 4 ok_ade ok_aie ', ok_ade , ok_aie 
         !* SAVE VARIABLES IN INTERIM VARIABLES A LA SW_AEROAR4
         ZFSUP0_AERO(:,:,4) = PFSCUP(:,:)
@@ -636 +662 @@
         LWUP_AERO(:,:,4) = PFLUX(:,1,:)
         LWDN_AERO(:,:,4) = PFLUX(:,2,:)
+
+        PTOAB_AERO(:,:,4)  = PTOAB(:,:) 
+        PTOACB_AERO(:,:,4)  = PTOACB(:,:) 
 
      ENDIF ! ok_ade .and. ok_aie
@@ -665 +694 @@
           & PFSDN , PFSUP  , PFSCDN  , PFSCUP , &
           & PTOAB, PTOACB )   !FC
+!        write(*,*) 'FC cas 5'
+!        write(*,*) 'FC cas 5 AEROSOLFEEDBACK_ACTIVE  flag_aerosol  swaerofree_diag', AEROSOLFEEDBACK_ACTIVE , flag_aerosol , swaerofree_diag
 
      !* SAVE VARIABLES IN INTERIM VARIABLES A LA SW_AEROAR4
@@ -678 +709 @@
      LWUP_AERO(:,:,5) = PFLUX(:,1,:)
      LWDN_AERO(:,:,5) = PFLUX(:,2,:)
+
+     PTOAB_AERO(:,:,5)  = PTOAB(:,:) 
+     PTOACB_AERO(:,:,5)  = PTOACB(:,:) 
 
   ENDIF ! .not. AEROSOLFEEDBACK_ACTIVE
@@ -754 +788 @@
         PFLUC(:,1,:) =  LWUP0_AERO(:,:,4)
         PFLUC(:,2,:) =  LWDN0_AERO(:,:,4)    
+
+        PTOAB(:,:)  = PTOAB_AERO(:,:,4)
+        PTOACB(:,:)  = PTOACB_AERO(:,:,4)
+!        write(*,*) 'je remplis 4'
      ENDIF
 
@@ -766 +804 @@
         PFLUC(:,1,:) =  LWUP0_AERO(:,:,3)
         PFLUC(:,2,:) =  LWDN0_AERO(:,:,3) 
+
+        PTOAB(:,:)  = PTOAB_AERO(:,:,3)
+        PTOACB(:,:)  = PTOACB_AERO(:,:,3)
+!        write(*,*) 'je remplis 3'
      ENDIF
 
@@ -778 +820 @@
         PFLUC(:,1,:) =  LWUP0_AERO(:,:,2)
         PFLUC(:,2,:) =  LWDN0_AERO(:,:,2) 
+
+        PTOAB(:,:)  = PTOAB_AERO(:,:,2)
+        PTOACB(:,:)  = PTOACB_AERO(:,:,2)
+!        write(*,*) 'je remplis 2'
      ENDiF
 
@@ -790 +836 @@
         PFLUC(:,1,:) =  LWUP0_AERO(:,:,1)
         PFLUC(:,2,:) =  LWDN0_AERO(:,:,1) 
+
+        PTOAB(:,:)  = PTOAB_AERO(:,:,1)
+        PTOACB(:,:)  = PTOACB_AERO(:,:,1)
+!        write(*,*) 'je remplis 1'
      ENDIF
 
      ! The following allows to compute the forcing diagostics without
      ! letting the aerosol forcing act on the meteorology
-     ! SEE logic above
 
   ELSE  !--not AEROSOLFEEDBACK_ACTIVE
@@ -808 +857 @@
      PFLUC(:,2,:) =  LWDN0_AERO(:,:,5) 
 
+     PTOAB(:,:)  = PTOAB_AERO(:,:,5)
+     PTOACB(:,:)  = PTOACB_AERO(:,:,5)
+!        write(*,*) 'je remplis 5'
   ENDIF