Changeset 1795 for LMDZ5/branches/testing/libf/phylmd/phys_state_var_mod.F90

Timestamp:

Jul 18, 2013, 10:20:28 AM (11 years ago)

Author:

Ehouarn Millour

Message:

Version testing basee sur la r1794

---

Testing release based on r1794

Location:

LMDZ5/branches/testing

Files:

• . (modified) (1 prop)
• libf/phylmd/phys_state_var_mod.F90 (modified) (6 diffs)

LMDZ5/branches/testing/libf/phylmd/phys_state_var_mod.F90

@@ -63 +63 @@
 !$OMP THREADPRIVATE(ratqs)
       REAL, ALLOCATABLE, SAVE :: pbl_tke(:,:,:) ! turb kinetic energy
-!$OMP THREADPRIVATE(pbl_tke)
+      REAL, ALLOCATABLE, SAVE :: coefh(:,:,:) ! Kz enthalpie
+      REAL, ALLOCATABLE, SAVE :: coefm(:,:,:) ! Kz momentum
+!$OMP THREADPRIVATE(pbl_tke, coefh,coefm)
       REAL, ALLOCATABLE, SAVE :: zmax0(:), f0(:) ! 
 !$OMP THREADPRIVATE(zmax0,f0)
@@ -359 +361 @@
 !======================================================================
 SUBROUTINE phys_state_var_init(read_climoz)
-use dimphy
+USE dimphy
 USE control_mod
-use aero_mod
-use infotrac, ONLY : nbtr
+USE aero_mod
+USE infotrac, ONLY : nbtr
+USE indice_sol_mod
 IMPLICIT NONE
 
@@ -373 +376 @@
 ! climatology and the daylight climatology
 
-#include "indicesol.h"
       ALLOCATE(rlat(klon), rlon(klon))
       ALLOCATE(pctsrf(klon,nbsrf))
@@ -394 +396 @@
       ALLOCATE(clwcon(klon,klev),rnebcon(klon,klev))
       ALLOCATE(ratqs(klon,klev))
-      ALLOCATE(pbl_tke(klon,klev+1,nbsrf))
+      ALLOCATE(pbl_tke(klon,klev+1,nbsrf+1))
+      ALLOCATE(coefh(klon,klev+1,nbsrf+1))
+      ALLOCATE(coefm(klon,klev+1,nbsrf+1))
       ALLOCATE(zmax0(klon), f0(klon))
       ALLOCATE(ema_work1(klon,klev), ema_work2(klon,klev))
@@ -518 +522 @@
 !======================================================================
 SUBROUTINE phys_state_var_end
-use dimphy
-use control_mod
+USE dimphy
+USE control_mod
+USE indice_sol_mod
 IMPLICIT NONE
-#include "indicesol.h"
 
       deallocate(rlat, rlon, pctsrf, ftsol, falb1, falb2)
@@ -530 +534 @@
       deallocate(        u_ancien, v_ancien                 )
       deallocate(        tr_ancien)                           !RomP
-      deallocate(ratqs, pbl_tke)
+      deallocate(ratqs, pbl_tke,coefh,coefm)
       deallocate(zmax0, f0)
       deallocate(ema_work1, ema_work2)