Changeset 3918 for LMDZ6/trunk/libf/phylmd/radlwsw_m.F90

Timestamp:

May 28, 2021, 3:39:30 PM (3 years ago)

Author:

idelkadi

Message:

Correction in the case of use of the ECRAD code: Reversal of flows on the vertical

File:

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

LMDZ6/trunk/libf/phylmd/radlwsw_m.F90

@@ -1321 +1321 @@
 !       longitude(rad), sin(latitude), PMASQ_ ???
       & ZGELAM, ZGEMU, zmasq, &
+!       pression et temp aux milieux
       & pplay_i, t_i, &
-!       PTEMPERATURE_H ?, &
+!       PTEMPERATURE_H ?, 
       & paprs_i, ZTH_i, q_i, qsat_i, & 
+!       Gas
       & ZCO2, ZCH4, ZN2O, ZNO2, ZCFC11, ZCFC12, ZHCFC22, ZCCL4, ZO3_DP_i, &
+!       
       & cldfra_i, flwc_i, fiwc_i, ZQ_RAIN, ZQ_SNOW, &  
       & ref_liq_i, ref_ice_i, &
@@ -1395 +1398 @@
       DO k=0,klev
          DO i=1,klon
-         ZEMTD(i,k+1)  = ZEMTD_i(i,k+1)
-         ZEMTU(i,k+1)  = ZEMTU_i(i,k+1)
-         ZTRSO(i,k+1)  = ZTRSO_i(i,k+1)
-         ZTH(i,k+1)    = ZTH_i(i,k+1)
+         ZEMTD(i,k+1)  = ZEMTD_i(i,klev+1-k)
+         ZEMTU(i,k+1)  = ZEMTU_i(i,klev+1-k)
+         ZTRSO(i,k+1)  = ZTRSO_i(i,klev+1-k)
+!         ZTH(i,k+1)    = ZTH_i(i,klev+1-k)
 ! AI ATTENTION
           ZLWFT(i,k+1)  = ZLWFT_i(i,klev+1-k)
@@ -1404 +1407 @@
           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)
-         ZFLUP0(i,k+1) = ZFLUC_i(i,1,k+1)
-         ZFLDN0(i,k+1) = ZFLUC_i(i,2,k+1)
-         ZFSDN(i,k+1)  = ZFSDWN_i(i,k+1)*fract(i)
-         ZFSDN0(i,k+1) = ZFCDWN_i(i,k+1)*fract(i)
-         ZFSDNC0(i,k+1)= ZFCCDWN_i(i,k+1)*fract(i)
-         ZFSUP (i,k+1) = ZFSUP_i(i,k+1)*fract(i)
-         ZFSUP0(i,k+1) = ZFCUP_i(i,k+1)*fract(i)
-         ZFSUPC0(i,k+1)= ZFCCUP_i(i,k+1)*fract(i)
-         ZFLDNC0(i,k+1)= ZFLCCDWN_i(i,k+1)
-         ZFLUPC0(i,k+1)= ZFLCCUP_i(i,k+1)
+         ZFLUP(i,k+1)  = ZFLUX_i(i,1,klev+1-k)
+         ZFLDN(i,k+1)  = ZFLUX_i(i,2,klev+1-k)
+         ZFLUP0(i,k+1) = ZFLUC_i(i,1,klev+1-k)
+         ZFLDN0(i,k+1) = ZFLUC_i(i,2,klev+1-k)
+         ZFSDN(i,k+1)  = ZFSDWN_i(i,klev+1-k)*fract(i)
+         ZFSDN0(i,k+1) = ZFCDWN_i(i,klev+1-k)*fract(i)
+         ZFSDNC0(i,k+1)= ZFCCDWN_i(i,klev+1-k)*fract(i)
+         ZFSUP (i,k+1) = ZFSUP_i(i,klev+1-k)*fract(i)
+         ZFSUP0(i,k+1) = ZFCUP_i(i,klev+1-k)*fract(i)
+         ZFSUPC0(i,k+1)= ZFCCUP_i(i,klev+1-k)*fract(i)
+         ZFLDNC0(i,k+1)= ZFLCCDWN_i(i,klev+1-k)
+         ZFLUPC0(i,k+1)= ZFLCCUP_i(i,klev+1-k)
          IF (ok_volcan) THEN
-            ZSWADAERO(i,k+1)=ZSWADAERO(i,k+1)*fract(i) !--NL
+            ZSWADAERO(i,k+1)=ZSWADAERO(i,klev+1-k)*fract(i) !--NL
          ENDIF
          
@@ -1593 +1596 @@
  ENDDO ! j = 1, nb_gr
 
+!IF (lldebug) THEN
+! print*,'================== Sortie de radlw ================='
+! print*,'cool = ', cool
+! print*,'heat = ', heat
+! print*,'topsw = ', topsw
+! print*,'toplw = ', toplw
+! print*,'sollw = ', sollw
+! print*,'solsw = ', solsw
+! print*,'lwdn = ', lwdn
+! print*,'lwup = ', lwup
+! print*,'swdn = ', swdn
+! print*,'swup =', swup
+!ENDIF
+
 END SUBROUTINE radlwsw