Changeset 634 for LMDZ4/branches/LMDZ4_par_0/libf/phylmd/calltherm.F

Timestamp:

May 4, 2005, 5:11:29 PM (19 years ago)

Author:

Laurent Fairhead

Message:

Modifications faites à la physique pour la rendre parallele YM
Une branche de travail LMDZ4_par_0 a été créée provisoirement afin de tester
les modifs pleinement avant leurs inclusions dans le tronc principal
LF

File:

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

LMDZ4/branches/LMDZ4_par_0/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
+      real :: zdt
+      logical,save :: first=.true.
       integer i,k
 
@@ -47 +50 @@
        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.