Changeset 3048 for LMDZ6/trunk/libf

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

Location:

LMDZ6/trunk/libf/phylmd

Files:

• clesphys.h (modified) (2 diffs)
• conf_phys_m.F90 (modified) (5 diffs)
• phys_output_ctrlout_mod.F90 (modified) (1 diff)
• phys_output_mod.F90 (modified) (2 diffs)
• phys_output_write_mod.F90 (modified) (3 diffs)
• physiq_mod.F90 (modified) (8 diffs)

LMDZ6/trunk/libf/phylmd/clesphys.h

@@ -11 +11 @@
        LOGICAL ok_limitvrai
        LOGICAL ok_all_xml
+       LOGICAL ok_lwoff
        INTEGER nbapp_rad, iflag_con, nbapp_cv, nbapp_wk, iflag_ener_conserv
        REAL co2_ppm, co2_ppm0, solaire
@@ -137 +138 @@
      &     , iflag_ice_thermo, ok_gwd_rando, NSW, iflag_albedo          &
      &     , ok_chlorophyll,ok_conserv_q, adjust_tropopause             &
-     &     , ok_daily_climoz, ok_all_xml                                &
+     &     , ok_daily_climoz, ok_all_xml, ok_lwoff                      &
      &     , iflag_phytrac
      

LMDZ6/trunk/libf/phylmd/conf_phys_m.F90

@@ -91 +91 @@
     LOGICAL, SAVE       :: ok_newmicro_omp
     LOGICAL, SAVE       :: ok_all_xml_omp
+    LOGICAL, SAVE       :: ok_lwoff_omp
     LOGICAL, SAVE       :: ok_journe_omp, ok_mensuel_omp, ok_instan_omp, ok_hf_omp        
     LOGICAL, SAVE       :: ok_LES_omp   
@@ -294 +295 @@
     ok_all_xml_omp = .FALSE.
     CALL getin('ok_all_xml', ok_all_xml_omp)
+
+    !
+    !Config Key  = ok_lwoff
+    !Config Desc = inhiber l effet radiatif LW des nuages
+    !Config Def  = .FALSE.
+    !Config Help = 
+    !
+    ok_lwoff_omp = .FALSE.
+    CALL getin('ok_lwoff', ok_lwoff_omp)
     !
 
@@ -2245 +2255 @@
 
     ok_all_xml = ok_all_xml_omp
+    ok_lwoff = ok_lwoff_omp
     ok_newmicro = ok_newmicro_omp
     ok_journe = ok_journe_omp
@@ -2326 +2337 @@
     cvl_corr = cvl_corr_omp
     ok_lic_melt = ok_lic_melt_omp
-    ok_lic_cond = ok_lic_cond_omp
     f_cdrag_ter=f_cdrag_ter_omp
     f_cdrag_oce=f_cdrag_oce_omp
@@ -2504 +2514 @@
     write(lunout,*)' cvl_corr=', cvl_corr
     write(lunout,*)'ok_lic_melt=', ok_lic_melt
-    write(lunout,*)'ok_lic_cond=', ok_lic_cond
     write(lunout,*)'cycle_diurne=',cycle_diurne
     write(lunout,*)'soil_model=',soil_model

LMDZ6/trunk/libf/phylmd/phys_output_ctrlout_mod.F90

@@ -1417 +1417 @@
   TYPE(ctrl_out), SAVE :: o_lcc3dstra = ctrl_out((/ 2,  6, 10, 10, 10, 10, 11, 11, 11, 11/), &
     'lcc3dstra', 'Stratiform cloud liquid fraction', '1', (/ ('', i=1, 10) /))
+  TYPE(ctrl_out), SAVE :: o_solbnd = ctrl_out((/ 2,  6, 10, 10, 10, 10, 11, 11, 11, 11/), &
+    'solbnd', 'Top-of-Atmosphere Solar Insolation for each band', 'W m-2', (/ ('', i=1, 10) /))
   TYPE(ctrl_out), SAVE :: o_stratomask = ctrl_out((/ 2,  6, 10, 10, 10, 10, 11, 11, 11, 11/), &
     'stratomask', 'Stratospheric fraction', '1', (/ ('', i=1, 10) /))

LMDZ6/trunk/libf/phylmd/phys_output_mod.F90

@@ -137 +137 @@
     REAL :: ini_heure
 #endif
+    INTEGER                         :: ISW
+    REAL, DIMENSION(NSW)            :: wl1_sun, wl2_sun !wavelength bounds (in um) for SW
+    REAL, DIMENSION(NSW)            :: wn1_sun, wn2_sun !wavenumber bounds (in m-1) for SW 
+    REAL, DIMENSION(NSW)            :: spectband  !mean wavenumb. of each sp.band
+    REAL, DIMENSION(NSW,2)          :: spbnds_sun !bounds of spectband
 
     WRITE(lunout,*) 'Debut phys_output_mod.F90'
@@ -354 +359 @@
      CALL wxios_add_vaxis("Alt", &
             levmax(iff) - levmin(iff) + 1, pseudoalt)
+
+    IF (NSW.EQ.6) THEN
+!
+!wl1_sun: minimum bound of wavelength (in um)
+!
+      wl1_sun(1)=0.180
+      wl1_sun(2)=0.250
+      wl1_sun(3)=0.440
+      wl1_sun(4)=0.690
+      wl1_sun(5)=1.190
+      wl1_sun(6)=2.380
+!
+!wl2_sun: maximum bound of wavelength (in um)
+!
+      wl2_sun(1)=0.250
+      wl2_sun(2)=0.440
+      wl2_sun(3)=0.690
+      wl2_sun(4)=1.190
+      wl2_sun(5)=2.380
+      wl2_sun(6)=4.000
+!
+    ELSE IF(NSW.EQ.2) THEN
+!
+!wl1_sun: minimum bound of wavelength (in um)
+!
+      wl1_sun(1)=0.250
+      wl1_sun(2)=0.690
+!
+!wl2_sun: maximum bound of wavelength (in um)
+!
+      wl2_sun(1)=0.690
+      wl2_sun(2)=4.000
+    ENDIF
+
+    DO ISW=1, NSW
+     wn1_sun(ISW)=1.e+6/wl1_sun(ISW) 
+     wn2_sun(ISW)=1.e+6/wl2_sun(ISW) 
+     spbnds_sun(ISW,1)=wn2_sun(ISW)
+     spbnds_sun(ISW,2)=wn1_sun(ISW)
+     spectband(ISW)=(wn1_sun(ISW)+wn2_sun(ISW))/2
+    ENDDO
+!
+!!! ajout axe vertical spectband : solar band number
+    CALL wxios_add_vaxis("spectband", NSW, spectband, positif='down')
    ELSE
     ! NMC files

LMDZ6/trunk/libf/phylmd/phys_output_write_mod.F90

@@ -119 +119 @@
          o_swtoacf_zero, o_swsrfcf_zero, &
          o_topswai, o_solswai, o_toplwai, o_sollwai, o_scdnc, &
-         o_cldncl, o_reffclws, o_reffclwc, o_stratomask,&
+         o_cldncl, o_reffclws, o_reffclwc, o_solbnd, o_stratomask,&
          o_cldnvi, o_lcc, o_lcc3d, o_lcc3dcon, &
          o_lcc3dstra, o_reffclwtop, o_ec550aer, &
@@ -389 +389 @@
     REAL, DIMENSION (klon,klev) :: zx_tmp_fi3d, zpt_conv
     REAL, DIMENSION (klon,klev+1) :: zx_tmp_fi3d1
+    REAL, DIMENSION (klon,NSW) :: zx_tmp_fi3dsp
     CHARACTER (LEN=4)              :: bb2
     INTEGER, DIMENSION(nbp_lon*nbp_lat)  :: ndex2d
@@ -1383 +1384 @@
        CALL histwrite_phy(o_zfull,zphi/RG)
 
+!solbnd begin
+#ifdef CPP_RRTM
+      IF (iflag_rrtm.EQ.1) THEN
+       IF (vars_defined) THEN
+        DO ISW=1, NSW
+          zx_tmp_fi3dsp(:,ISW) = swdn(:,klevp1)*swradcorr(:)*RSUN(ISW)
+        ENDDO
+        CALL histwrite_phy(o_solbnd, zx_tmp_fi3dsp)
+       ENDIF
+      ENDIF
+#endif
+
+!solbnd end
+
        IF (flag_aerosol_strat.EQ.2) THEN
          CALL histwrite_phy(o_stratomask, stratomask)

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