Changeset 4881 for LMDZ6/trunk/libf/phylmd/phys_local_var_mod.F90

Timestamp:

Mar 28, 2024, 6:06:40 PM (7 weeks ago)

Author:

evignon

Message:

extraction plus propre de la dissipation de TKE

File:

• LMDZ6/trunk/libf/phylmd/phys_local_var_mod.F90 (modified) (3 diffs)

LMDZ6/trunk/libf/phylmd/phys_local_var_mod.F90

@@ -22 +22 @@
       REAL, SAVE, ALLOCATABLE :: d_rneb_dyn(:,:)
       !$OMP THREADPRIVATE(d_rneb_dyn)
-      REAL, SAVE, ALLOCATABLE :: l_mixmin(:,:,:),l_mix(:,:,:),tke_dissip(:,:,:),wprime(:,:,:)
-      !$OMP THREADPRIVATE(l_mixmin, l_mix, tke_dissip,wprime)
+      REAL, SAVE, ALLOCATABLE :: l_mixmin(:,:,:),l_mix(:,:,:),wprime(:,:,:)
+      !$OMP THREADPRIVATE(l_mixmin, l_mix, wprime)
+      REAL, SAVE, ALLOCATABLE :: pbl_eps(:,:,:)
+      !$OMP THREADPRIVATE(pbl_eps)
       REAL, SAVE, ALLOCATABLE :: tr_seri(:,:,:)
       !$OMP THREADPRIVATE(tr_seri)
@@ -666 +668 @@
       ALLOCATE(t_seri(klon,klev),q_seri(klon,klev),ql_seri(klon,klev),qs_seri(klon,klev), qbs_seri(klon,klev))
       ALLOCATE(u_seri(klon,klev),v_seri(klon,klev))
-      ALLOCATE(l_mixmin(klon,klev+1,nbsrf),l_mix(klon,klev+1,nbsrf),tke_dissip(klon,klev+1,nbsrf),wprime(klon,klev+1,nbsrf))
-      l_mix(:,:,:)=0.;l_mixmin(:,:,:)=0.;tke_dissip(:,:,:)=0.;wprime(:,:,:)=0. ! doit etre initialse car pas toujours remplis
-
+      ALLOCATE(l_mixmin(klon,klev+1,nbsrf),l_mix(klon,klev+1,nbsrf),wprime(klon,klev+1,nbsrf))
+      ALLOCATE(pbl_eps(klon,klev+1,nbsrf+1))
+      l_mix(:,:,:)=0.;l_mixmin(:,:,:)=0.;wprime(:,:,:)=0. ! doit etre initialse car pas toujours remplis
       ALLOCATE(rhcl(klon,klev))
       ALLOCATE(tr_seri(klon,klev,nbtr))
@@ -1019 +1021 @@
       DEALLOCATE(t_seri,q_seri,ql_seri,qs_seri, qbs_seri)
       DEALLOCATE(u_seri,v_seri)
-      DEALLOCATE(l_mixmin,l_mix, tke_dissip,wprime)
-
+      DEALLOCATE(l_mixmin,l_mix,wprime)
+      DEALLOCATE(pbl_eps)
       DEALLOCATE(rhcl)
       DEALLOCATE(tr_seri)