Changeset 1114 for LMDZ4/branches/LMDZ4-dev/libf/phylmd/physiq.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/physiq.F (modified) (16 diffs)

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

@@ -2 +2 @@
 c#define IO_DEBUG
 
-      SUBROUTINE physiq (nlon,nlev,nqmax,
+      SUBROUTINE physiq (nlon,nlev,
      .            debut,lafin,rjourvrai,gmtime,pdtphys,
      .            paprs,pplay,pphi,pphis,presnivs,clesphy0,
@@ -15 +15 @@
       USE write_field_phy
       USE dimphy
+      USE infotrac
       USE mod_grid_phy_lmdz
       USE mod_phys_lmdz_para
@@ -50 +51 @@
 c nlon----input-I-nombre de points horizontaux
 c nlev----input-I-nombre de couches verticales
-c nqmax---input-I-nombre de traceurs (y compris vapeur d'eau) = 1
 c debut---input-L-variable logique indiquant le premier passage
 c lafin---input-L-variable logique indiquant le dernier passage
@@ -87 +87 @@
 #include "clesphys.h"
 #include "control.h"
-!#include "logic.h"
 #include "temps.h"
-cym#include "comgeomphy.h"
-#include "advtrac.h"
 #include "iniprint.h"
 #include "thermcell.h"
@@ -184 +181 @@
       INTEGER nlon
       INTEGER nlev
-      INTEGER nqmax
       REAL rjourvrai
       REAL gmtime
@@ -200 +196 @@
       REAL v(klon,klev)
       REAL t(klon,klev),theta(klon,klev)
-      REAL qx(klon,klev,nqmax)
+      REAL qx(klon,klev,nqtot)
       REAL flxmass_w(klon,klev)
       REAL omega(klon,klev) ! vitesse verticale en Pa/s
@@ -206 +202 @@
       REAL d_v(klon,klev)
       REAL d_t(klon,klev)
-      REAL d_qx(klon,klev,nqmax)
+      REAL d_qx(klon,klev,nqtot)
       REAL d_ps(klon)
       real da(klon,klev),phi(klon,klev,klev),mp(klon,klev)
@@ -1110 +1106 @@
          write(lunout,*) 'DEBUT DE PHYSIQ !!!!!!!!!!!!!!!!!!!!'
          write(lunout,*)
-     s 'nlon,nlev,nqmax,debut,lafin,rjourvrai,gmtime,pdtphys'
+     s 'nlon,nlev,nqtot,debut,lafin,rjourvrai,gmtime,pdtphys'
          write(lunout,*)
-     s  nlon,nlev,nqmax,debut,lafin,rjourvrai,gmtime,pdtphys 
+     s  nlon,nlev,nqtot,debut,lafin,rjourvrai,gmtime,pdtphys 
 
          write(lunout,*) 'papers, play, phi, u, v, t, omega'
@@ -1170 +1166 @@
       END IF 
       ok_sync=.TRUE.
-      IF (nqmax .LT. 2) THEN
-         abort_message = 'eaux vapeur et liquide sont indispensables'
-         CALL abort_gcm (modname,abort_message,1)
-      ENDIF
+
       IF (debut) THEN
          CALL suphel ! initialiser constantes et parametres phys.
@@ -1500 +1493 @@
 
 c$OMP MASTER
-       call phys_output_open(jjmp1,nqmax,nlevSTD,clevSTD,nbteta, 
+       call phys_output_open(jjmp1,nlevSTD,clevSTD,nbteta, 
      &                        ctetaSTD,dtime,presnivs,ok_veget,
      &                        type_ocean,iflag_pbl,ok_mensuel,ok_journe,
@@ -1554 +1547 @@
      $                   calday,
      $                   klon,
-     $                   nqmax,
+     $                   nqtot,
      $                   pdtphys,
      $                   annee_ref,
@@ -1610 +1603 @@
       ENDDO
       ENDDO
-      DO iq = 1, nqmax
+      DO iq = 1, nqtot
       DO k = 1, klev
       DO i = 1, klon
@@ -1633 +1626 @@
       ENDDO
       ENDDO
-      IF (nqmax.GE.3) THEN
-      DO iq = 3, nqmax
+      IF (nqtot.GE.3) THEN
+      DO iq = 3, nqtot
       DO  k = 1, klev
       DO  i = 1, klon
@@ -3163 +3156 @@
      I                   debut,
      I                   lafin,
-     I                   nqmax-2,
      I                   nlon,
      I                   nlev,
@@ -3386 +3378 @@
       ENDDO
 c
-      IF (nqmax.GE.3) THEN
-      DO iq = 3, nqmax
+      IF (nqtot.GE.3) THEN
+      DO iq = 3, nqtot
       DO  k = 1, klev
       DO  i = 1, klon
@@ -3421 +3413 @@
       write(lunout,*) 'FIN DE PHYSIQ !!!!!!!!!!!!!!!!!!!!'
       write(lunout,*)
-     s 'nlon,nlev,nqmax,debut,lafin,rjourvrai,gmtime,pdtphys pct tlos'
+     s 'nlon,nlev,nqtot,debut,lafin,rjourvrai,gmtime,pdtphys pct tlos'
       write(lunout,*)
-     s  nlon,nlev,nqmax,debut,lafin,rjourvrai,gmtime,pdtphys,
+     s  nlon,nlev,nqtot,debut,lafin,rjourvrai,gmtime,pdtphys,
      s  pctsrf(igout,is_ter), pctsrf(igout,is_lic),pctsrf(igout,is_oce),
      s  pctsrf(igout,is_sic)