Changeset 5589 for LMDZ6/branches/contrails/libf/phylmd/radlwsw_m.F90

Timestamp:

Mar 26, 2025, 6:05:40 PM (2 months ago)

Author:

aborella

Message:

Multiple changes:

• added new radiative diagnostics for contrails
• added ok_rad_contrail option to allow for a double call of RRTM (w/ and w/o contrails)
• transformed resuspension of snow into ice sedimentation (poprecip)
• some modifications in poprecip in line with the ones from EV
• cleaned sublimation of ice clouds in lmdz_lscp_condensation, option ok_ice_supersat
• aviation emissions can now be read with IOIPSL (in lon/lat mode)

File:

• LMDZ6/branches/contrails/libf/phylmd/radlwsw_m.F90 (modified) (7 diffs)

LMDZ6/branches/contrails/libf/phylmd/radlwsw_m.F90

@@ -45 +45 @@
        ZLWFT0_i, ZFLDN0, ZFLUP0, &
        ZSWFT0_i, ZFSDN0, ZFSUP0, &
-       cloud_cover_sw)
+       cloud_cover_sw, &
+       !--AB contrails radiative effects
+       cldfra_nocont, fiwc_nocont, ref_ice_nocont, &
+       topsw_nocont, solsw_nocont, toplw_nocont, sollw_nocont)
 
     ! Modules necessaires
@@ -279 +282 @@
     REAL, DIMENSION(kdlon,kflev+1), INTENT(out) :: ZSWFT0_i
     REAL, DIMENSION(kdlon,kflev+1), INTENT(out) :: ZLWFT0_i
+    !--AB contrails radiative effects
+    REAL, DIMENSION(klon,klev), INTENT(in)  :: cldfra_nocont
+    REAL, DIMENSION(klon,klev), INTENT(in)  :: fiwc_nocont
+    REAL, DIMENSION(klon,klev), INTENT(in)  :: ref_ice_nocont
+    REAL, DIMENSION(klon),      INTENT(out) :: topsw_nocont
+    REAL, DIMENSION(klon),      INTENT(out) :: solsw_nocont
+    REAL, DIMENSION(klon),      INTENT(out) :: toplw_nocont
+    REAL, DIMENSION(klon),      INTENT(out) :: sollw_nocont
 
     ! Local variables
@@ -467 +478 @@
     REAL(KIND=8) ZFLCCDWN_i (klon,klev+1)
     REAL(KIND=8) ZFLCCUP_i (klon,klev+1)
+    !--AB contrails radiative effects
+    REAL(KIND=8) cldfra_nocont_i(klon,klev)
+    REAL(KIND=8) fiwc_nocont_i(klon,klev)
+    REAL(KIND=8) ref_ice_nocont_i(klon,klev)
+    REAL(KIND=8) ZTOPSWNOCONT(klon)
+    REAL(KIND=8) ZSOLSWNOCONT(klon)
+    REAL(KIND=8) ZTOPLWNOCONT(klon)
+    REAL(KIND=8) ZSOLLWNOCONT(klon)
     ! 3 lignes suivantes a activer pour CCMVAL (MPL 20100412)
     !      REAL(KIND=8) RSUN(3,2)
@@ -900 +919 @@
              ref_liq_pi_i(1:klon,k) =ref_liq_pi(1:klon,klev+1-k)
              ref_ice_pi_i(1:klon,k) =ref_ice_pi(1:klon,klev+1-k)
+             IF (ok_rad_contrail) THEN
+               !--AB contrails radiative effects
+               cldfra_nocont_i(1:klon,k)  = cldfra_nocont(1:klon,klev+1-k)
+               fiwc_nocont_i(1:klon,k)    = fiwc_nocont(1:klon,klev+1-k)
+               ref_ice_nocont_i(1:klon,k) = ref_ice_nocont(1:klon,klev+1-k)
+             ENDIF
           ENDDO
           DO k=1,kflev
@@ -980 +1005 @@
                ZLWADAERO, & !--NL
                volmip_solsw, flag_volc_surfstrat, & !--VOLMIP
-               ok_ade, ok_aie, ok_volcan, flag_aerosol,flag_aerosol_strat, flag_aer_feedback) ! flags aerosols
+               ok_ade, ok_aie, ok_volcan, flag_aerosol,flag_aerosol_strat, flag_aer_feedback, & ! flags aerosols
+               !--AB contrails radiative effect
+               ok_rad_contrail, cldfra_nocont_i, fiwc_nocont_i, ref_ice_nocont_i, &
+               ZTOPSWNOCONT, ZSOLSWNOCONT, ZTOPLWNOCONT, ZSOLLWNOCONT)
 
           !--OB diagnostics
@@ -1080 +1108 @@
           ZSOLSWCF_AERO(:,2)=ZSOLSWCF_AERO(:,2)*fract(:)
           ZSOLSWCF_AERO(:,3)=ZSOLSWCF_AERO(:,3)*fract(:)
+
+          IF (ok_rad_contrail) THEN
+            !--AB contrails radiative effect
+            ZTOPSWNOCONT(:) = ZTOPSWNOCONT(:) * fract(:)
+            ZSOLSWNOCONT(:) = ZSOLSWNOCONT(:) * fract(:)
+          ENDIF
 
           ! ---------
@@ -1655 +1689 @@
           ENDDO
        ENDIF
+       !--AB radiative effect from contrails
+       IF (ok_rad_contrail) THEN
+         DO i = 1, kdlon
+           topsw_nocont(iof+i) = ZTOPSWNOCONT(i)
+           solsw_nocont(iof+i) = ZSOLSWNOCONT(i)
+           toplw_nocont(iof+i) = ZTOPLWNOCONT(i)
+           sollw_nocont(iof+i) = ZSOLLWNOCONT(i)
+         ENDDO
+       ENDIF
        DO k = 1, kflev
           DO i = 1, kdlon