Changeset 1146 for LMDZ4/trunk/libf/dyn3dpar/advtrac_p.F

Timestamp:

Apr 9, 2009, 12:11:35 PM (15 years ago)

Author:

Laurent Fairhead

Message:

Réintegration dans le tronc des modifications issues de la branche LMDZ-dev
comprises entre la révision 1074 et 1145
Validation: une simulation de 1 jour en séquentiel sur PC donne les mêmes
résultats entre la trunk et la dev
LF

Location:

LMDZ4/trunk

Files:

• . (modified) (1 prop)
• libf/dyn3dpar/advtrac_p.F (modified) (9 diffs)

LMDZ4/trunk/libf/dyn3dpar/advtrac_p.F

@@ -22 +22 @@
       USE Vampir
       USE times
+      USE infotrac
       IMPLICIT NONE
 c
@@ -35 +36 @@
 #include "ener.h"
 #include "description.h"
-#include "advtrac.h"
 
 c-------------------------------------------------------------------
@@ -46 +46 @@
       INTEGER iapptrac
       REAL pbaru(ip1jmp1,llm),pbarv(ip1jm,llm)
-      REAL q(ip1jmp1,llm,nqmx),masse(ip1jmp1,llm)
+      REAL q(ip1jmp1,llm,nqtot),masse(ip1jmp1,llm)
       REAL p( ip1jmp1,llmp1 ),teta(ip1jmp1,llm)
       REAL pk(ip1jmp1,llm)
@@ -59 +59 @@
       REAL,SAVE::pbarug(ip1jmp1,llm),pbarvg(ip1jm,llm),wg(ip1jmp1,llm) 
       REAL (kind=kind(1.d0)) :: t_initial, t_final, tps_cpu
-      real cpuadv(nqmx)
-      common/cpuadv/cpuadv
-
       INTEGER iadvtr
       INTEGER ij,l,iq,iiq
@@ -76 +73 @@
       REAL psppm(iim,jjp1) ! pression  au sol
       REAL unatppm(iim,jjp1,llm),vnatppm(iim,jjp1,llm)
-      REAL qppm(iim*jjp1,llm,nqmx)
+      REAL qppm(iim*jjp1,llm,nqtot)
       REAL fluxwppm(iim,jjp1,llm)
       REAL apppm(llmp1), bpppm(llmp1)
@@ -88 +85 @@
       REAL,SAVE :: teta_tmp(ip1jmp1,llm)
       REAL,SAVE :: pk_tmp(ip1jmp1,llm)
-      
+
       ijb_u=ij_begin
       ije_u=ij_end
@@ -196 +193 @@
       call Register_SwapFieldHallo(pk_tmp,pk_tmp,ip1jmp1,llm,
      *                             jj_Nb_vanleer,1,1,Request_vanleer)
-      do j=1,nqmx
+      do j=1,nqtot
         call Register_SwapFieldHallo(q(1,1,j),q(1,1,j),ip1jmp1,llm,
      *                             jj_nb_vanleer,0,0,Request_vanleer)
@@ -279 +276 @@
 c     Appel des sous programmes d'advection
 c-----------------------------------------------------------
-      do iq=1,nqmx
+      do iq=1,nqtot
 c        call clock(t_initial)
         if(iadv(iq) == 0) cycle 
@@ -479 +476 @@
 c$OMP END MASTER
 
-        do j=1,nqmx
+        do j=1,nqtot
           call Register_SwapFieldHallo(q(1,1,j),q(1,1,j),ip1jmp1,llm,
      *                             jj_nb_caldyn,0,0,Request_vanleer)