Changeset 3208 for LMDZ6/trunk/libf/phylmd/phys_state_var_mod.F90

Timestamp:

Feb 16, 2018, 12:42:18 PM (6 years ago)

Author:

jyg

Message:

Implementation of a first crude model of the
dynamic of wake population.

File:

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

LMDZ6/trunk/libf/phylmd/phys_state_var_mod.F90

@@ -247 +247 @@
 ! wake_deltaq : ecart d'humidite avec la zone non perturbee
 ! wake_s      : fraction surfacique occupee par la poche froide
+! awake_dens  : number of active wakes per unit area
 ! wake_dens   : number of wakes per unit area
 ! wake_occ    : occurence of wakes (= 1 if wakes occur, =0 otherwise)
@@ -258 +259 @@
       REAL,ALLOCATABLE,SAVE :: wake_s(:)
 !$OMP THREADPRIVATE(wake_s)
-      REAL,ALLOCATABLE,SAVE :: wake_dens(:)
-!$OMP THREADPRIVATE(wake_dens)
+      REAL,ALLOCATABLE,SAVE :: awake_dens(:), wake_dens(:)
+!$OMP THREADPRIVATE(awake_dens, wake_dens)
       REAL,ALLOCATABLE,SAVE :: wake_Cstar(:)
 !$OMP THREADPRIVATE(wake_Cstar)
@@ -539 +540 @@
       ALLOCATE(wght_th(klon,klev))
       ALLOCATE(wake_deltat(klon,klev), wake_deltaq(klon,klev))
-      ALLOCATE(wake_s(klon), wake_dens(klon))
+      ALLOCATE(wake_s(klon), awake_dens(klon), wake_dens(klon))
       ALLOCATE(wake_Cstar(klon))
       ALLOCATE(wake_pe(klon), wake_fip(klon))
@@ -682 +683 @@
       deallocate(lalim_conv, wght_th)
       deallocate(wake_deltat, wake_deltaq)
-      deallocate(wake_s, wake_dens)
+      deallocate(wake_s, awake_dens, wake_dens)
       deallocate(wake_Cstar, wake_pe, wake_fip)
 !jyg<