Changeset 1114 for LMDZ4/branches/LMDZ4-dev/libf/phylmd/conema3.F

Timestamp:

Mar 3, 2009, 5:40:26 PM (15 years ago)

Author:

jghattas

Message:

Creation du module infotrac:

• contient les variables de advtrac.h
• contient la subroutine iniadvtrac renommer en infotrac_init
• le nombre des traceurs est lu dans tracer.def en dynamique (ou par default ou recu par INCA)
• ce module est utilise dans la dynamique et la physique
• contient aussi la variable nbtr qui avant etait stockee dans dimphy

Le fichier advtrac.h n'existe plus.
La compilation ne prend plus en compte le nombre de traceur.

/JG

File:

• LMDZ4/branches/LMDZ4-dev/libf/phylmd/conema3.F (modified) (5 diffs)

LMDZ4/branches/LMDZ4-dev/libf/phylmd/conema3.F

@@ -10 +10 @@
 
       USE dimphy
+      USE infotrac, ONLY : nbtr
       IMPLICIT none
 c======================================================================
@@ -55 +56 @@
 c
 #include "dimensions.h"
-cym#include "dimphy.h"
 #include "conema3.h"
       INTEGER i, l,m,itra
-      INTEGER ntra,ntrac !number of tracers; if no tracer transport
+      INTEGER ntra       ! if no tracer transport
                          ! is needed, set ntra = 1 (or 0)
-      PARAMETER (ntrac=nqmx-2)
       REAL dtime
 c
@@ -97 +96 @@
       REAL,ALLOCATABLE,SAVE :: em_qs(:) 
 c$OMP THREADPRIVATE(em_qs)  
-cym      REAL em_u(klev), em_v(klev), em_tra(klev,ntrac)
+cym      REAL em_u(klev), em_v(klev), em_tra(klev,nbtr)
       REAL,ALLOCATABLE,SAVE :: em_u(:),em_v(:),em_tra(:,:)
 c$OMP THREADPRIVATE(em_u,em_v,em_tra)      
@@ -111 +110 @@
       REAL,ALLOCATABLE,SAVE :: em_d_t(:),em_d_q(:)
 c$OMP THREADPRIVATE(em_d_t,em_d_q)
-cym      REAL em_d_u(klev), em_d_v(klev), em_d_tra(klev,ntrac)
+cym      REAL em_d_u(klev), em_d_v(klev), em_d_tra(klev,nbtr)
       REAL,ALLOCATABLE,SAVE ::em_d_u(:),em_d_v(:),em_d_tra(:,:)
 c$OMP THREADPRIVATE(em_d_u,em_d_v,em_d_tra)      
@@ -188 +187 @@
         allocate(em_q(klev))
         allocate(em_qs(klev))
-        allocate(em_u(klev), em_v(klev), em_tra(klev,ntrac))
+        allocate(em_u(klev), em_v(klev), em_tra(klev,nbtr))
         allocate(em_ph(klev+1), em_p(klev))
         allocate(em_work1(klev), em_work2(klev))
         allocate(em_d_t(klev), em_d_q(klev))
-        allocate(em_d_u(klev), em_d_v(klev), em_d_tra(klev,ntrac))
+        allocate(em_d_u(klev), em_d_v(klev), em_d_tra(klev,nbtr))
         allocate(em_upwd(klev), em_dnwd(klev), em_dnwdbis(klev))
         allocate(emmip(klev))