Changeset 704 for LMDZ4/branches/V3_test/libf/phylmd/calltherm.F

Timestamp:

Aug 17, 2006, 5:41:51 PM (18 years ago)

Author:

Laurent Fairhead

Message:

Inclusion des modifs de Y. Meurdesoif pour la version V3
LF

File:

• LMDZ4/branches/V3_test/libf/phylmd/calltherm.F (modified) (3 diffs)

LMDZ4/branches/V3_test/libf/phylmd/calltherm.F

@@ -8 +8 @@
      s      ,fm_therm,entr_therm)
 
+      USE dimphy
       implicit none
-#include "dimensions.h"
-#include "dimphy.h"
+cym#include "dimensions.h"
+cym#include "dimphy.h"
 #include "thermcell.h"
 
@@ -35 +36 @@
       REAL d_u_the(klon,klev),d_v_the(klon,klev)
 c
-      real zfm_therm(klon,klev+1),zentr_therm(klon,klev),zdt
-      save zentr_therm,zfm_therm
-
+cym      real zfm_therm(klon,klev+1),zentr_therm(klon,klev),zdt
+cym      save zentr_therm,zfm_therm
+      real,allocatable,dimension(:,:),save :: zfm_therm,zentr_therm
+c$OMP THREADPRIVATE(zfm_therm,zentr_therm)
+      real :: zdt
+      logical,save :: first=.true.
+c$OMP THREADPRIVATE(first)
       integer i,k
 
@@ -47 +52 @@
        print*,'avant isplit ',nsplit_thermals
 
-
+cym initialisation dynamique
+         if (first) then
+           allocate(zfm_therm(klon,klev+1),zentr_therm(klon,klev))
+           zfm_therm(:,:)=0.
+           zentr_therm(:,:)=0.
+           first=.false.
+         endif
+         
          fm_therm(:,:)=0.
          entr_therm(:,:)=0.