Changeset 3048 for LMDZ6/trunk/libf/phylmd/physiq_mod.F90

Timestamp:

Oct 27, 2017, 6:14:12 PM (7 years ago)

Author:

musat

Message:

Set up lwoff configuration : add ok_lwoff=y to activate it.
Add solbnd CMIP6 output.
Enable 3D outputs of SW/LW total sky and clear sky.
IM

File:

• LMDZ6/trunk/libf/phylmd/physiq_mod.F90 (modified) (8 diffs)

LMDZ6/trunk/libf/phylmd/physiq_mod.F90

@@ -1105 +1105 @@
     REAL, dimension(klon, klev) :: cldfrarad   ! fraction nuageuse
 
+    !lwoff=y : offset LW CRE for radiation code and other schemes
+    REAL, SAVE :: betalwoff 
+    !OMP THREADPRIVATE(betalwoff)
+!
     INTEGER :: nbtr_tmp ! Number of tracer inside concvl
     REAL, dimension(klon,klev) :: sh_in ! Specific humidity entering in phytrac
@@ -1788 +1792 @@
        WRITE(*,*)'lat2_beta=',lat2_beta
        WRITE(*,*)'mskocean_beta=',mskocean_beta
+
+      !lwoff=y : offset LW CRE for radiation code and other schemes
+      !lwoff=y : betalwoff=1.
+      betalwoff=0.
+      IF (ok_lwoff) THEN
+         betalwoff=1.
+      ENDIF
+      WRITE(*,*)'ok_lwoff=',ok_lwoff
+      ! 
+      !lwoff=y to begin only sollw and sollwdown are set up to CS values
+      sollw = sollw + betalwoff * (sollw0 - sollw)
+      sollwdown(:)= sollwdown(:) + betalwoff *(-1.*ZFLDN0(:,1) - &
+                    sollwdown(:))
     ENDIF
     !
@@ -3603 +3620 @@
           ! global
           !
+!IM 251017 begin
+                print*,'physiq betaCRF global zdtime=',zdtime
+!IM 251017 end
           DO k=1, klev
              DO i=1, klon
@@ -3717 +3737 @@
              print *,' ->radlwsw, number 1 '
           ENDIF
+
           !
           CALL radlwsw &
@@ -3733 +3754 @@
                tau_aero_lw_rrtm, & 
                cldtaupirad,new_aod, &
+!              zqsat, flwcrad, fiwcrad, &
                zqsat, flwc, fiwc, &
                ref_liq, ref_ice, ref_liq_pi, ref_ice_pi, &
@@ -3755 +3777 @@
                ZSWFT0_i, ZFSDN0, ZFSUP0)
 
+          !lwoff=y, betalwoff=1. : offset LW CRE for radiation code and other
+          !schemes
+          toplw = toplw + betalwoff * (toplw0 - toplw)
+          sollw = sollw + betalwoff * (sollw0 - sollw)
+          lwdn = lwdn + betalwoff * (lwdn0 - lwdn)
+          lwup = lwup + betalwoff * (lwup0 - lwup)
+          sollwdown(:)= sollwdown(:) + betalwoff *(-1.*ZFLDN0(:,1) - &
+                        sollwdown(:))
+          cool = cool + betalwoff * (cool0 - cool)
+ 
 #ifndef CPP_XIOS
           !--OB 30/05/2016 modified 21/10/2016
@@ -3805 +3837 @@
                      tau_aero_lw_rrtm, &
                      cldtaupi,new_aod, &
+!                    zqsat, flwcrad, fiwcrad, &
                      zqsat, flwc, fiwc, &
                      ref_liq, ref_ice, ref_liq_pi, ref_ice_pi, &
@@ -3826 +3859 @@
                      ZLWFT0_i, ZFLDN0, ZFLUP0, &
                      ZSWFT0_i, ZFSDN0, ZFSUP0)
-          endif
-          !
+          endif !ok_4xCO2atm
        ENDIF ! aerosol_couple
        itaprad = 0