Changeset 3914

Timestamp:

May 26, 2021, 10:27:28 AM (4 years ago)

Author:

idelkadi

Message:

Correction: calculation of net radiative flux in the case where ECRAD is activated

File:

• LMDZ6/trunk/libf/phylmd/radlwsw_m.F90 (modified) (2 diffs)

LMDZ6/trunk/libf/phylmd/radlwsw_m.F90

@@ -1399 +1399 @@
          ZTRSO(i,k+1)  = ZTRSO_i(i,k+1)
          ZTH(i,k+1)    = ZTH_i(i,k+1)
+! AI ATTENTION
+          ZLWFT(i,k+1)  = ZLWFT_i(i,klev+1-k)
+          ZSWFT(i,k+1)  = ZSWFT_i(i,klev+1-k)*fract(i)
+          ZSWFT0_i(i,k+1) = ZSWFT0_i(i,klev+1-k)*fract(i)
+!
          ZFLUP(i,k+1)  = ZFLUX_i(i,1,k+1)
          ZFLDN(i,k+1)  = ZFLUX_i(i,2,k+1)
@@ -1417 +1422 @@
 !   Nouveau calcul car visiblement ZSWFT et ZSWFC sont nuls dans RRTM cy32
 !   en sortie de radlsw.F90 - MPL 7.01.09
-         ZSWFT(i,k+1)  = (ZFSDWN_i(i,k+1)-ZFSUP_i(i,k+1))*fract(i)
-         ZSWFT0_i(i,k+1) = (ZFCDWN_i(i,k+1)-ZFCUP_i(i,k+1))*fract(i)
-         ZLWFT(i,k+1) =-ZFLUX_i(i,2,k+1)-ZFLUX_i(i,1,k+1)
-         ZLWFT0_i(i,k+1)=-ZFLUC_i(i,2,k+1)-ZFLUC_i(i,1,k+1)
+! AI ATTENTION
+!         ZSWFT(i,k+1)  = (ZFSDWN_i(i,k+1)-ZFSUP_i(i,k+1))*fract(i)
+!         ZSWFT0_i(i,k+1) = (ZFCDWN_i(i,k+1)-ZFCUP_i(i,k+1))*fract(i)
+!         ZLWFT(i,k+1) =-ZFLUX_i(i,2,k+1)-ZFLUX_i(i,1,k+1)
+!         ZLWFT0_i(i,k+1)=-ZFLUC_i(i,2,k+1)-ZFLUC_i(i,1,k+1)
          ENDDO
       ENDDO