Changeset 2516

Timestamp:

May 19, 2016, 8:58:23 PM (9 years ago)

Author:

fhourdin

Message:

Introduction d'un seuil oliqmax sur l'eau condensée liquide maximum
conservée dans l'atmosphère.
Par défaut, oliqmax=999 et la correction est inactive.
Pour des valeurs plus faibles, l'eau en sus est convertie en
précipitation ajoutée à rain_lsc (plul dans les sorties).

Cette pluie additionelle est également sauvegardée dans rain_num

et sortie dans plun.

Cette modification est une tentative pour résoudre certains plantages

en haute résolution.

Location:

LMDZ5/trunk/libf/phylmd

Files:

• phys_local_var_mod.F90 (modified) (3 diffs)
• phys_output_ctrlout_mod.F90 (modified) (1 diff)
• phys_output_write_mod.F90 (modified) (3 diffs)
• physiq_mod.F90 (modified) (3 diffs)

LMDZ5/trunk/libf/phylmd/phys_local_var_mod.F90

@@ -281 +281 @@
       REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: zxfluxlat, zxtsol, snow_lsc, zxfqfonte
 !$OMP THREADPRIVATE(zxfluxlat, zxtsol, snow_lsc, zxfqfonte)
-      REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: zxqsurf, rain_lsc
-!$OMP THREADPRIVATE(zxqsurf, rain_lsc)
+      REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: zxqsurf, rain_lsc, rain_num
+!$OMP THREADPRIVATE(zxqsurf, rain_lsc, rain_num)
 !
 !jyg+nrlmd<
@@ -541 +541 @@
       ALLOCATE(zxtsol(klon), snow_lsc(klon), zxfqfonte(klon), zxqsurf(klon))
       ALLOCATE(rain_lsc(klon))
+      ALLOCATE(rain_num(klon))
 !
       ALLOCATE(sens_x(klon), sens_w(klon))
@@ -746 +747 @@
       DEALLOCATE(zxtsol, snow_lsc, zxfqfonte, zxqsurf)
       DEALLOCATE(rain_lsc)
+      DEALLOCATE(rain_num)
 !
       DEALLOCATE(sens_x, sens_w)

LMDZ5/trunk/libf/phylmd/phys_output_ctrlout_mod.F90

@@ -154 +154 @@
   TYPE(ctrl_out), SAVE :: o_plul = ctrl_out((/ 1, 1, 1, 10, 10, 10, 11, 11, 11 /), &
     'plul', 'Large-scale Precip.', 'kg/(s*m2)', (/ ('', i=1, 9) /))
+  TYPE(ctrl_out), SAVE :: o_plun = ctrl_out((/ 1, 1, 1, 10, 10, 10, 11, 11, 11 /), &
+    'plun', 'Numerical Precip.', 'kg/(s*m2)', (/ ('', i=1, 9) /))
   TYPE(ctrl_out), SAVE :: o_pluc = ctrl_out((/ 1, 1, 1, 10, 5, 10, 11, 11, 11 /), &
     'pluc', 'Convective Precip.', 'kg/(s*m2)', (/ ('', i=1, 9) /))

LMDZ5/trunk/libf/phylmd/phys_output_write_mod.F90

@@ -36 +36 @@
          o_wind10m, o_wind10max, o_gusts, o_sicf, &
          o_psol, o_mass, o_qsurf, o_qsol, &
-         o_precip, o_ndayrain, o_plul, o_pluc, &
+         o_precip, o_ndayrain, o_plul, o_pluc, o_plun, &
          o_snow, o_msnow, o_fsnow, o_evap, o_ep,o_epmax_diag, & ! epmax_cape
          o_tops, o_tops0, o_topl, o_topl0, &
@@ -195 +195 @@
          t2m_min_mon, t2m_max_mon, evap, &
          zu10m, zv10m, zq2m, zustar, zxqsurf, &
-         rain_lsc, snow_lsc, bils, sens, fder, &
+         rain_lsc, rain_num, snow_lsc, bils, sens, fder, &
          zxffonte, zxfqcalving, zxfqfonte, fluxu, &
          fluxv, zxsnow, qsnow, snowhgt, to_ice, &
@@ -441 +441 @@
        ENDIF
        CALL histwrite_phy(o_plul, zx_tmp_fi2d)
+       CALL histwrite_phy(o_plun, rain_num)
 
        IF (vars_defined) THEN

LMDZ5/trunk/libf/phylmd/physiq_mod.F90

@@ -432 +432 @@
     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) 
@@ -1050 +1054 @@
        CALL getin_p('random_notrig_max',random_notrig_max)
        CALL getin_p('ok_adjwk',ok_adjwk)
+       CALL getin_p('oliqmax',oliqmax)
     ENDIF
 
@@ -2879 +2884 @@
     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