Changeset 3179 for LMDZ6/trunk/libf/phylmd/physiq_mod.F90

Timestamp:

Feb 3, 2018, 9:52:04 AM (6 years ago)

Author:

jyg

Message:

New structure for the representation of vdf
splitting in pbl_surface. Computations are
gathered in two subroutines included in the module
wx_pbl_mod: wx_pbl0_fuse, called before the
sub-surfaces, determines the single column
equivalent to the (w) and (x) columns.
wx_pbl0_split, called after the subsurfaces,
determines the distinct (w) and (x) surface
fluxes.

This is a first version with no surface
temperature difference between (w) and (x) (hence
the index 0).

File:

• LMDZ6/trunk/libf/phylmd/physiq_mod.F90 (modified) (4 diffs)

LMDZ6/trunk/libf/phylmd/physiq_mod.F90

@@ -63 +63 @@
        !
        d_t_vdf,d_q_vdf,d_u_vdf,d_v_vdf,d_t_diss, &
-       d_t_vdf_w,d_q_vdf_w, &
-       d_t_vdf_x,d_q_vdf_x, &
        d_ts, &
        !
@@ -135 +133 @@
        zxfluxlat_x, zxfluxlat_w, &
        !
-       dtvdf_x, dtvdf_w, &
-       dqvdf_x, dqvdf_w, &
+       d_t_vdf_x, d_t_vdf_w, &
+       d_q_vdf_x, d_q_vdf_w, &
        pbl_tke_input, &
        t_therm, q_therm, u_therm, v_therm, &
@@ -2271 +2269 @@
 !!jyg       IF (prt_level .ge. 2 .and. mod(iflag_pbl_split,2) .eq. 1) THEN
        IF (prt_level .ge. 2 .and. mod(iflag_pbl_split,10) .ge. 1) THEN
-          print *,'debut du splitting de la PBL'
+          print *,'debut du splitting de la PBL, wake_s = ', wake_s(:)
+          print *,'debut du splitting de la PBL, wake_deltat = ', wake_deltat(:,1)
+          print *,'debut du splitting de la PBL, wake_deltaq = ', wake_deltaq(:,1)
        ENDIF
        ! !!
@@ -2984 +2984 @@
           call prt_enerbil('wake',itap)
        ENDIF   ! (iflag_wake_tend .GT. 0.)
+       !
+       IF (prt_level .GE. 10) THEN
+         print *,' physiq, after calwake, wake_s: ',wake_s(:)
+         print *,' physiq, after calwake, wake_deltat: ',wake_deltat(:,1)
+         print *,' physiq, after calwake, wake_deltaq: ',wake_deltaq(:,1)
+       ENDIF
 
        IF (iflag_alp_wk_cond .GT. 0.) THEN