Ignore:
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:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ6/trunk/libf/phylmd/radlwsw_m.F90

    r3914 r3918  
    13211321!       longitude(rad), sin(latitude), PMASQ_ ???
    13221322      & ZGELAM, ZGEMU, zmasq, &
     1323!       pression et temp aux milieux
    13231324      & pplay_i, t_i, &
    1324 !       PTEMPERATURE_H ?, &
     1325!       PTEMPERATURE_H ?,
    13251326      & paprs_i, ZTH_i, q_i, qsat_i, &
     1327!       Gas
    13261328      & ZCO2, ZCH4, ZN2O, ZNO2, ZCFC11, ZCFC12, ZHCFC22, ZCCL4, ZO3_DP_i, &
     1329!       
    13271330      & cldfra_i, flwc_i, fiwc_i, ZQ_RAIN, ZQ_SNOW, & 
    13281331      & ref_liq_i, ref_ice_i, &
     
    13951398      DO k=0,klev
    13961399         DO i=1,klon
    1397          ZEMTD(i,k+1)  = ZEMTD_i(i,k+1)
    1398          ZEMTU(i,k+1)  = ZEMTU_i(i,k+1)
    1399          ZTRSO(i,k+1)  = ZTRSO_i(i,k+1)
    1400          ZTH(i,k+1)    = ZTH_i(i,k+1)
     1400         ZEMTD(i,k+1)  = ZEMTD_i(i,klev+1-k)
     1401         ZEMTU(i,k+1)  = ZEMTU_i(i,klev+1-k)
     1402         ZTRSO(i,k+1)  = ZTRSO_i(i,klev+1-k)
     1403!         ZTH(i,k+1)    = ZTH_i(i,klev+1-k)
    14011404! AI ATTENTION
    14021405          ZLWFT(i,k+1)  = ZLWFT_i(i,klev+1-k)
     
    14041407          ZSWFT0_i(i,k+1) = ZSWFT0_i(i,klev+1-k)*fract(i)
    14051408!
    1406          ZFLUP(i,k+1)  = ZFLUX_i(i,1,k+1)
    1407          ZFLDN(i,k+1)  = ZFLUX_i(i,2,k+1)
    1408          ZFLUP0(i,k+1) = ZFLUC_i(i,1,k+1)
    1409          ZFLDN0(i,k+1) = ZFLUC_i(i,2,k+1)
    1410          ZFSDN(i,k+1)  = ZFSDWN_i(i,k+1)*fract(i)
    1411          ZFSDN0(i,k+1) = ZFCDWN_i(i,k+1)*fract(i)
    1412          ZFSDNC0(i,k+1)= ZFCCDWN_i(i,k+1)*fract(i)
    1413          ZFSUP (i,k+1) = ZFSUP_i(i,k+1)*fract(i)
    1414          ZFSUP0(i,k+1) = ZFCUP_i(i,k+1)*fract(i)
    1415          ZFSUPC0(i,k+1)= ZFCCUP_i(i,k+1)*fract(i)
    1416          ZFLDNC0(i,k+1)= ZFLCCDWN_i(i,k+1)
    1417          ZFLUPC0(i,k+1)= ZFLCCUP_i(i,k+1)
     1409         ZFLUP(i,k+1)  = ZFLUX_i(i,1,klev+1-k)
     1410         ZFLDN(i,k+1)  = ZFLUX_i(i,2,klev+1-k)
     1411         ZFLUP0(i,k+1) = ZFLUC_i(i,1,klev+1-k)
     1412         ZFLDN0(i,k+1) = ZFLUC_i(i,2,klev+1-k)
     1413         ZFSDN(i,k+1)  = ZFSDWN_i(i,klev+1-k)*fract(i)
     1414         ZFSDN0(i,k+1) = ZFCDWN_i(i,klev+1-k)*fract(i)
     1415         ZFSDNC0(i,k+1)= ZFCCDWN_i(i,klev+1-k)*fract(i)
     1416         ZFSUP (i,k+1) = ZFSUP_i(i,klev+1-k)*fract(i)
     1417         ZFSUP0(i,k+1) = ZFCUP_i(i,klev+1-k)*fract(i)
     1418         ZFSUPC0(i,k+1)= ZFCCUP_i(i,klev+1-k)*fract(i)
     1419         ZFLDNC0(i,k+1)= ZFLCCDWN_i(i,klev+1-k)
     1420         ZFLUPC0(i,k+1)= ZFLCCUP_i(i,klev+1-k)
    14181421         IF (ok_volcan) THEN
    1419             ZSWADAERO(i,k+1)=ZSWADAERO(i,k+1)*fract(i) !--NL
     1422            ZSWADAERO(i,k+1)=ZSWADAERO(i,klev+1-k)*fract(i) !--NL
    14201423         ENDIF
    14211424         
     
    15931596 ENDDO ! j = 1, nb_gr
    15941597
     1598!IF (lldebug) THEN
     1599! print*,'================== Sortie de radlw ================='
     1600! print*,'cool = ', cool
     1601! print*,'heat = ', heat
     1602! print*,'topsw = ', topsw
     1603! print*,'toplw = ', toplw
     1604! print*,'sollw = ', sollw
     1605! print*,'solsw = ', solsw
     1606! print*,'lwdn = ', lwdn
     1607! print*,'lwup = ', lwup
     1608! print*,'swdn = ', swdn
     1609! print*,'swup =', swup
     1610!ENDIF
     1611
    15951612END SUBROUTINE radlwsw
    15961613
Note: See TracChangeset for help on using the changeset viewer.