Changeset 940 for LMDZ4/trunk/libf/phylmd/concvl.F

Timestamp:

Apr 7, 2008, 4:33:30 PM (16 years ago)

Author:

Laurent Fairhead

Message:

On remplace le fichier include dimphy.h par le module dimphy.F90i pour etre
coherent avec le partout
LF

File:

• LMDZ4/trunk/libf/phylmd/concvl.F (modified) (9 diffs)

LMDZ4/trunk/libf/phylmd/concvl.F

@@ +1 @@
+
 !
 ! $Header$
@@ -23 +24 @@
 *
 c
-c      USE dimphy
+      USE dimphy
       IMPLICIT none
 c======================================================================
@@ -66 +67 @@
 c
 #include "dimensions.h"
-#include "dimphy.h"
+cccccc#include "dimphy.h"
 c
       integer NTRAC
@@ -122 +123 @@
        INTEGER i,k,itra
        REAL qs(klon,klev),qs_wake(klon,klev)
-       REAL cbmf(klon)
-       SAVE cbmf
-!       REAL cbmflast(klon)
+cLF       REAL cbmf(klon)
+cLF       SAVE cbmf
+       REAL, SAVE, ALLOCATABLE :: cbmf(:)
+c$OMP THREADPRIVATE(cbmf)!       
+       REAL cbmflast(klon)
        INTEGER ifrst
        SAVE ifrst
@@ -133 +136 @@
 C     Variables supplementaires liees au bilan d'energie
 c      Real paire(klon)
-      Real ql(klon,klev)
+cLF      Real ql(klon,klev)
 c      Save paire
-      Save ql
-      Real t1(klon,klev),q1(klon,klev)
-      Save t1,q1
+cLF      Save ql
+cLF      Real t1(klon,klev),q1(klon,klev)
+cLF      Save t1,q1
 c      Data paire /1./
+       REAL, SAVE, ALLOCATABLE :: ql(:,:), q1(:,:), t1(:,:)
+c$OMP THREADPRIVATE(ql, q1, t1)
 c
 C     Variables liees au bilan d'energie et d'enthalpi
@@ -162 +167 @@
       REAL ZRCPD
 c-jld ec_conser
+cLF
+      INTEGER nloc 
+      logical, save :: first=.true.
 c
 #include "YOMCST.h"
@@ -168 +176 @@
 #include "FCTTRE.h"
 c
+      if (first) then
+c Allocate some variables LF 04/2008
+c
+        allocate(cbmf(klon))
+        allocate(ql(klon,klev))
+        allocate(t1(klon,klev))
+        allocate(q1(klon,klev))
+      endif
 
 c    Copy T into Tconv
@@ -189 +205 @@
       snow(:)=0
       
-      IF (ifrst .EQ. 0) THEN
-         ifrst = 1
+c      IF (ifrst .EQ. 0) THEN
+c         ifrst = 1
+       if (first) then
+         first=.false.
 c
 C===========================================================================
@@ -313 +331 @@
       else
 
-      CALL cva_driver(klon,klev,klev+1,ntra,
+cLF   necessary for gathered fields
+      nloc=klon
+      CALL cva_driver(klon,klev,klev+1,ntra,nloc,
      $              iflag_con,iflag_mix,iflag_clos,dtime,
      :              t,q,qs,t_wake,q_wake,qs_wake,u,v,tra,