Changeset 2524 for LMDZ5

Timestamp:

May 26, 2016, 7:39:44 AM (9 years ago)

Author:

oboucher

Message:

Enabling the possibility to prescribe daily values for solaire and RSUN
Only works for RRTM new radiation scheme
This is activated by ok_suntime_rrtm
This does not have any impact on REPROBUS

Location:

LMDZ5/trunk/libf/phylmd

Files:

• clesphys.h (modified) (3 diffs)
• conf_phys_m.F90 (modified) (4 diffs)
• physiq_mod.F90 (modified) (5 diffs)
• rrtm/read_rsun_rrtm.F90 (added)

LMDZ5/trunk/libf/phylmd/clesphys.h

@@ -15 +15 @@
        INTEGER nbapp_rad, iflag_con,iflag_ener_conserv
        REAL co2_ppm, co2_ppm0, solaire
+       LOGICAL ok_suntime_rrtm
        REAL(kind=8) RCO2, RCH4, RN2O, RCFC11, RCFC12  
        REAL(kind=8) RCO2_act, RCH4_act, RN2O_act, RCFC11_act, RCFC12_act  
@@ -110 +111 @@
      &     , top_height                                                 &
      &     , cycle_diurne, soil_model, new_oliq                         &
-     &     ,ok_orodr, ok_orolf, ok_limitvrai, nbapp_rad, iflag_con      &
+     &     , ok_orodr, ok_orolf, ok_limitvrai, nbapp_rad, iflag_con     &
      &     , iflag_ener_conserv                                         &
      &     , ok_4xCO2atm                                                & 
+     &     , ok_suntime_rrtm                                            & 
      &     , overlap                                                    &
      &     , ok_kzmin                                                   &
@@ -129 +131 @@
        save /clesphys/
 !$OMP THREADPRIVATE(/clesphys/)
- 

LMDZ5/trunk/libf/phylmd/conf_phys_m.F90

@@ -148 +148 @@
 
     REAL,SAVE :: R_ecc_omp,R_peri_omp,R_incl_omp,solaire_omp
+    LOGICAL,SAVE :: ok_suntime_rrtm_omp
     REAL,SAVE :: co2_ppm_omp, RCO2_omp, co2_ppm_per_omp, RCO2_per_omp
     REAL,SAVE :: CH4_ppb_omp, RCH4_omp, CH4_ppb_per_omp, RCH4_per_omp
@@ -515 +516 @@
     call getin('solaire', solaire_omp)
     !
+    !Config Key  = ok_sun_time
+    !Config Desc = oui ou non variabilite solaire
+    !Config Def  = .false.
+    !Config Help =
+    !
+    !
+    !valeur AMIP II
+    ok_suntime_rrtm_omp = .false.
+    call getin('ok_suntime_rrtm',ok_suntime_rrtm_omp)
+    !
     !Config Key  = co2_ppm
     !Config Desc = concentration du gaz carbonique en ppmv
@@ -1967 +1978 @@
     R_incl = R_incl_omp
     solaire = solaire_omp
+    ok_suntime_rrtm = ok_suntime_rrtm_omp
     co2_ppm = co2_ppm_omp
     RCO2 = RCO2_omp
@@ -2281 +2293 @@
     write(lunout,*)' Inclinaison =',R_incl
     write(lunout,*)' Constante solaire =',solaire
+    write(lunout,*)' ok_suntime_rrtm =',ok_suntime_rrtm
     write(lunout,*)' co2_ppm =',co2_ppm
     write(lunout,*)' RCO2_act = ',RCO2_act

LMDZ5/trunk/libf/phylmd/physiq_mod.F90

@@ -69 +69 @@
 #ifdef CPP_RRTM
     USE YOERAD, ONLY : NRADLP
+    USE YOESW, ONLY : RSUN
 #endif
     USE ioipsl_getin_p_mod, ONLY : getin_p
@@ -432 +433 @@
     REAL d_t_adjwk(klon,klev)                !jyg
     REAL d_q_adjwk(klon,klev)                !jyg
-
-    REAL,SAVE :: oliqmax=999.    ! Seuil en kg/kg pour le contenu max
-    !$OMP THREADPRIVATE(oliqmax) ! en eau liquide dans l'atmosphère.
-
     LOGICAL,SAVE :: ok_adjwk=.FALSE.
     !$OMP THREADPRIVATE(ok_adjwk) 
@@ -1054 +1051 @@
        CALL getin_p('random_notrig_max',random_notrig_max)
        CALL getin_p('ok_adjwk',ok_adjwk)
-       CALL getin_p('oliqmax',oliqmax)
     ENDIF
 
@@ -2884 +2880 @@
     CALL add_phys_tend(du0,dv0,d_t_lsc,d_q_lsc,d_ql_lsc,d_qi_lsc,paprs, &
          'lsc',abortphy)
-
-   !---------------------------------------------------------------------------
-   ! Elimination d'eau liquide > oliqmax
-   !---------------------------------------------------------------------------
-    rain_num(:)=0.
-    DO k = 1, klev  
-       DO i = 1, klon
-          IF (ql_seri(i,k)>oliqmax) THEN
-             rain_num(i)=rain_num(i)+(ql_seri(i,k)-oliqmax)*zmasse(i,k)/pdtphys
-             ql_seri(i,k)=oliqmax
-          ENDIF
-       ENDDO
-    ENDDO
-    rain_lsc(:)=rain_lsc(:)+rain_num(:)
-
+    !---------------------------------------------------------------------------
     DO k = 1, klev
        DO i = 1, klon
@@ -3429 +3411 @@
        !enddo
        !albedo SB <<<
+
+!--if ok_suntime_rrtm we use ancillay data for RSUN 
+!--previous values are therefore overwritten
+!--this is needed for CMIP6 runs
+!--and only possible for new radiation scheme
+       IF (iflag_rrtm.EQ.1.AND.ok_suntime_rrtm) THEN 
+         CALL read_rsun_rrtm(debut)
+       ENDIF
+
+       if (mydebug) then
+          call writefield_phy('u_seri',u_seri,nbp_lev)
+          call writefield_phy('v_seri',v_seri,nbp_lev)
+          call writefield_phy('t_seri',t_seri,nbp_lev)
+          call writefield_phy('q_seri',q_seri,nbp_lev)
+       endif
 
        !