Changeset 2524 for LMDZ5/trunk/libf/phylmd/physiq_mod.F90

Timestamp:

May 26, 2016, 7:39:44 AM (8 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

File:

• LMDZ5/trunk/libf/phylmd/physiq_mod.F90 (modified) (5 diffs)

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
 
        !