Changeset 3179 for LMDZ6/trunk/libf/phylmd/phys_local_var_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/phys_local_var_mod.F90 (modified) (5 diffs)

LMDZ6/trunk/libf/phylmd/phys_local_var_mod.F90

@@ -347 +347 @@
 !!!$OMP THREADPRIVATE(q_x, q_w)
 !>jyg
-!!! Sorties ferret
-      REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: dtvdf_x, dtvdf_w
-!$OMP THREADPRIVATE(dtvdf_x, dtvdf_w)
-      REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: dqvdf_x, dqvdf_w
-!$OMP THREADPRIVATE(dqvdf_x, dqvdf_w)
 ! Variables suppl\E9mentaires dans physiq.F relative au splitting de la surface
       REAL,ALLOCATABLE,SAVE,DIMENSION(:,:,:) :: pbl_tke_input
@@ -579 +574 @@
       ALLOCATE(plul_st(klon),plul_th(klon))
       ALLOCATE(d_t_vdf(klon,klev),d_q_vdf(klon,klev),d_t_diss(klon,klev))
-!nrlmd+jyg<
-      ALLOCATE(d_t_vdf_w(klon,klev),d_q_vdf_w(klon,klev))
-      ALLOCATE(d_t_vdf_x(klon,klev),d_q_vdf_x(klon,klev))
-!>nrlmd+jyg
       ALLOCATE(d_u_vdf(klon,klev),d_v_vdf(klon,klev))
       ALLOCATE(d_t_oli(klon,klev),d_t_oro(klon,klev))
@@ -724 +715 @@
 !!      ALLOCATE(q_x(klon,klev), q_w(klon,klev))
 !>jyg
-      ALLOCATE(dtvdf_x(klon,klev), dtvdf_w(klon,klev))
-      ALLOCATE(dqvdf_x(klon,klev), dqvdf_w(klon,klev))
+      ALLOCATE(d_t_vdf_x(klon,klev), d_t_vdf_w(klon,klev))
+      ALLOCATE(d_q_vdf_x(klon,klev), d_q_vdf_w(klon,klev))
       ALLOCATE(pbl_tke_input(klon,klev+1,nbsrf))
       ALLOCATE(t_therm(klon,klev), q_therm(klon,klev),u_therm(klon,klev), v_therm(klon,klev))
@@ -881 +872 @@
       DEALLOCATE(plul_st,plul_th)
       DEALLOCATE(d_t_vdf,d_q_vdf,d_t_diss)
-!nrlmd+jyg<
-      DEALLOCATE(d_t_vdf_w,d_q_vdf_w)
-      DEALLOCATE(d_t_vdf_x,d_q_vdf_x)
-!>nrlmd+jyg
       DEALLOCATE(d_u_vdf,d_v_vdf)
       DEALLOCATE(d_t_oli,d_t_oro)
@@ -1016 +1003 @@
 !!      DEALLOCATE(q_x, q_w)
 !>jyg
-      DEALLOCATE(dtvdf_x, dtvdf_w)
-      DEALLOCATE(dqvdf_x, dqvdf_w)
+      DEALLOCATE(d_t_vdf_x, d_t_vdf_w)
+      DEALLOCATE(d_q_vdf_x, d_q_vdf_w)
       DEALLOCATE(pbl_tke_input)
       DEALLOCATE(t_therm, q_therm, u_therm, v_therm)