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

Timestamp:

Oct 2, 2017, 10:16:06 AM (7 years ago)

Author:

jyg

Message:

(i) Implementation of the conditionning of the Alp
provided by the wakes: when iflag_alp_wk_cond > 0,
Alp_wk is divided by the probability that there is
a gust front within the grid cell.

New subroutine alpale_wk.F90 called by physiq.

(ii) Some changes concerning the initialization of
the variable wake_dens.

File:

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

LMDZ5/trunk/libf/phylmd/physiq_mod.F90

@@ -566 +566 @@
                                                      ! updated within calwake
     !$OMP THREADPRIVATE(iflag_wake_tend)
+    INTEGER,  SAVE               :: iflag_alp_wk_cond=0 ! wake: if =0, then Alp_wk is the average lifting
+                                                        ! power provided by the wakes; else, Alp_wk is the
+                                                        ! lifting power conditionned on the presence of a 
+                                                        ! gust-front in the grid cell.
+    !$OMP THREADPRIVATE(iflag_alp_wk_cond)
     REAL t_w(klon,klev),q_w(klon,klev) ! temperature and moisture profiles in the wake region
     REAL t_x(klon,klev),q_x(klon,klev) ! temperature and moisture profiles in the off-wake region
@@ -591 +596 @@
     !>jyg
 
+    REAL wake_fip_0(klon)           ! Average Front Incoming Power (unconditionned)
     REAL wake_gfl(klon)             ! Gust Front Length
 !!!    REAL wake_dens(klon)         ! moved to phys_state_var_mod
@@ -1226 +1232 @@
     IF (debut) THEN
        CALL suphel ! initialiser constantes et parametres phys.
+       CALL getin_p('iflag_alp_wk_cond', iflag_alp_wk_cond) 
        CALL getin_p('random_notrig_max',random_notrig_max)
        CALL getin_p('ok_adjwk',ok_adjwk) 
@@ -2805 +2812 @@
                wake_deltat, wake_deltaq, wake_s, wake_dens,  &
                wake_dth, wake_h,  &
-               wake_pe, wake_fip, wake_gfl,  &
+!!               wake_pe, wake_fip, wake_gfl,  &
+               wake_pe, wake_fip_0, wake_gfl,  &   !! jyg
                d_t_wake, d_q_wake,  &
                wake_k, t_x, q_x,  &
@@ -2835 +2843 @@
           call prt_enerbil('wake',itap)
        ENDIF   ! (iflag_wake_tend .GT. 0.)
+
+       IF (iflag_alp_wk_cond .GT. 0.) THEN
+
+         CALL alpale_wk(dtime, cell_area, wake_k, wake_s, wake_dens, wake_fip_0, &
+                        wake_fip)
+       ELSE
+         wake_fip(:) = wake_fip_0(:)
+       ENDIF   ! (iflag_alp_wk_cond .GT. 0.)
 
     ENDIF  ! (iflag_wake>=1)