Changeset 1114 for LMDZ4/branches/LMDZ4-dev/libf/dyn3dpar/calfis_p.F

Timestamp:

Mar 3, 2009, 5:40:26 PM (16 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/dyn3dpar/calfis_p.F (modified) (17 diffs)

LMDZ4/branches/LMDZ4-dev/libf/dyn3dpar/calfis_p.F

@@ -4 +4 @@
 C
 C
-      SUBROUTINE calfis_p(nq,
-     $                  lafin,
+      SUBROUTINE calfis_p(lafin,
      $                  rdayvrai,
      $                  heure,
@@ -40 +39 @@
       USE Times
       USE IOPHY
+      USE infotrac
+
       IMPLICIT NONE
 c=======================================================================
@@ -98 +99 @@
 #include "paramet.h"
 #include "temps.h"
-#include "advtrac.h"
-
-      INTEGER ngridmx,nq
+
+      INTEGER ngridmx
       PARAMETER( ngridmx = 2+(jjm-1)*iim - 1/jjm   )
 
@@ -119 +119 @@
       REAL pteta(iip1,jjp1,llm)
       REAL pmasse(iip1,jjp1,llm)
-      REAL pq(iip1,jjp1,llm,nqmx)
+      REAL pq(iip1,jjp1,llm,nqtot)
       REAL pphis(iip1,jjp1)
       REAL pphi(iip1,jjp1,llm)
@@ -126 +126 @@
       REAL pducov(iip1,jjp1,llm)
       REAL pdteta(iip1,jjp1,llm)
-      REAL pdq(iip1,jjp1,llm,nqmx)
+      REAL pdq(iip1,jjp1,llm,nqtot)
 c
       REAL pps(iip1,jjp1)
@@ -135 +135 @@
       REAL pdufi(iip1,jjp1,llm)
       REAL pdhfi(iip1,jjp1,llm)
-      REAL pdqfi(iip1,jjp1,llm,nqmx)
+      REAL pdqfi(iip1,jjp1,llm,nqtot)
       REAL pdpsfi(iip1,jjp1)
 
@@ -253 +253 @@
       ALLOCATE(zphi(klon,llm),zphis(klon))
       ALLOCATE(zufi(klon,llm), zvfi(klon,llm))
-      ALLOCATE(ztfi(klon,llm),zqfi(klon,llm,nqmx))
+      ALLOCATE(ztfi(klon,llm),zqfi(klon,llm,nqtot))
       ALLOCATE(pcvgu(klon,llm), pcvgv(klon,llm))
       ALLOCATE(pcvgt(klon,llm), pcvgq(klon,llm,2))
 c      ALLOCATE(pvervel(klon,llm))
       ALLOCATE(zdufi(klon,llm),zdvfi(klon,llm))
-      ALLOCATE(zdtfi(klon,llm),zdqfi(klon,llm,nqmx))
+      ALLOCATE(zdtfi(klon,llm),zdqfi(klon,llm,nqtot))
       ALLOCATE(zdpsrf(klon))
       ALLOCATE(zdufi2(klon+iim,llm),zdvfi2(klon+iim,llm))
@@ -335 +335 @@
 c
 
-      DO iq=1,nq
+      DO iq=1,nqtot
          iiq=niadv(iq)
 c$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
@@ -527 +527 @@
 cc$OMP  PARALLEL DEFAULT(NONE) 
 cc$OMP+ PRIVATE(i,l,offset,iq)
-cc$OMP+ SHARED(klon_omp_nb,nq,klon_omp_begin,
+cc$OMP+ SHARED(klon_omp_nb,nqtot,klon_omp_begin,
 cc$OMP+        debut,lafin,rdayvrai,heure,dtphys,zplev,zplay,
 cc$OMP+        zphi,zphis,presnivs,clesphy0,zufi,zvfi,ztfi,
@@ -549 +549 @@
         allocate(zvfi_omp(klon,llm))
         allocate(ztfi_omp(klon,llm))
-        allocate(zqfi_omp(klon,llm,nq))
+        allocate(zqfi_omp(klon,llm,nqtot))
 c        allocate(pvervel_omp(klon,llm))
         allocate(zdufi_omp(klon,llm))
         allocate(zdvfi_omp(klon,llm))
         allocate(zdtfi_omp(klon,llm))
-        allocate(zdqfi_omp(klon,llm,nq))
+        allocate(zdqfi_omp(klon,llm,nqtot))
         allocate(zdpsrf_omp(klon))
         allocate(flxwfi_omp(klon,llm))
@@ -609 +609 @@
       enddo
         
-      do iq=1,nq
+      do iq=1,nqtot
         do l=1,llm
           do i=1,klon
@@ -641 +641 @@
       enddo
         
-      do iq=1,nq
+      do iq=1,nqtot
         do l=1,llm
           do i=1,klon
@@ -664 +664 @@
       CALL physiq (klon,
      .             llm,
-     .             nq,
      .             debut,
      .             lafin,
@@ -743 +742 @@
       enddo
         
-      do iq=1,nq
+      do iq=1,nqtot
         do l=1,llm
           do i=1,klon
@@ -775 +774 @@
       enddo
         
-      do iq=1,nq
+      do iq=1,nqtot
         do l=1,llm
           do i=1,klon
@@ -937 +936 @@
 c   ---------------------
 
-      DO iq=1,nqmx
+      DO iq=1,nqtot
 c$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
          DO l=1,llm
@@ -976 +975 @@
 C
 
-      DO iq=1,nq
+      DO iq=1,nqtot
          iiq=niadv(iq)
 c$OMP DO SCHEDULE(STATIC,OMP_CHUNK)