Changeset 2159 for LMDZ5/trunk/libf/phylmd/phys_state_var_mod.F90

Timestamp:

Nov 27, 2014, 4:48:31 PM (10 years ago)

Author:

jyg

Message:

1/ Splitting of the boundary layer : the climbing down and up of Pbl_surface is
split between the off-wake and wake regions ; the thermal scheme is applied
only to the off-wake region.
2/ Elimination of wake_scal and calwake_scal.

File:

• LMDZ5/trunk/libf/phylmd/phys_state_var_mod.F90 (modified) (6 diffs)

LMDZ5/trunk/libf/phylmd/phys_state_var_mod.F90

@@ -66 +66 @@
       REAL, ALLOCATABLE, SAVE :: coefm(:,:,:) ! Kz momentum
 !$OMP THREADPRIVATE(pbl_tke, coefh,coefm)
+!nrlmd<
+      REAL, ALLOCATABLE, SAVE :: delta_tsurf(:,:) ! Surface temperature difference inside-outside cold pool
+!$OMP THREADPRIVATE(delta_tsurf)
+!>nrlmd
       REAL, ALLOCATABLE, SAVE :: zmax0(:), f0(:) ! 
 !$OMP THREADPRIVATE(zmax0,f0)
@@ -230 +234 @@
 !$OMP THREADPRIVATE(dq_wake)
 !
+!jyg<
+! variables related to the spitting of the PBL between wake and 
+! off-wake regions.
+! wake_delta_pbl_TKE : difference TKE_w - TKE_x
+      REAL,ALLOCATABLE,SAVE :: wake_delta_pbl_TKE(:,:,:)
+!$OMP THREADPRIVATE(wake_delta_pbl_TKE)
+!>jyg
+!
 ! pfrac_impa : Produits des coefs lessivage impaction
 ! pfrac_nucl : Produits des coefs lessivage nucleation
@@ -406 +418 @@
       ALLOCATE(ratqs(klon,klev))
       ALLOCATE(pbl_tke(klon,klev+1,nbsrf+1))
+!nrlmd<
+      ALLOCATE(delta_tsurf(klon,nbsrf))
+!>nrlmd
       ALLOCATE(coefh(klon,klev+1,nbsrf+1))
       ALLOCATE(coefm(klon,klev+1,nbsrf+1))
@@ -475 +490 @@
       ALLOCATE(wake_pe(klon), wake_fip(klon))
       ALLOCATE(dt_wake(klon,klev), dq_wake(klon,klev))
+!jyg<
+      ALLOCATE(wake_delta_pbl_TKE(klon,klev+1,nbsrf))
+!>jyg
       ALLOCATE(pfrac_impa(klon,klev), pfrac_nucl(klon,klev))
       ALLOCATE(pfrac_1nucl(klon,klev))
@@ -551 +569 @@
       deallocate(        tr_ancien)                           !RomP
       deallocate(ratqs, pbl_tke,coefh,coefm)
+!nrlmd<
+      deallocate(delta_tsurf)
+!>nrlmd
       deallocate(zmax0, f0)
       deallocate(sig1, w01)
@@ -601 +622 @@
       deallocate(wake_Cstar, wake_s, wake_pe, wake_fip)
       deallocate(dt_wake, dq_wake)
+!jyg<
+      deallocate(wake_delta_pbl_TKE)
+!>jyg
       deallocate(pfrac_impa, pfrac_nucl)
       deallocate(pfrac_1nucl)