Changeset 1114 for LMDZ4/branches/LMDZ4-dev/libf/dyn3dpar/etat0_netcdf.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/dyn3dpar/etat0_netcdf.F (modified) (5 diffs)

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

@@ -13 +13 @@
       USE phys_state_var_mod
       USE filtreg_mod
+      USE infotrac
       !
       IMPLICIT NONE
@@ -39 +40 @@
       REAL :: vvent(iip1, jjm, llm)
       REAL :: t3d(iip1, jjp1, llm), tpot(iip1, jjp1, llm)
-      REAL :: q3d(iip1, jjp1, llm,nqmx), qsat(iip1, jjp1, llm)
+      REAL :: q3d(iip1, jjp1, llm,nqtot), qsat(iip1, jjp1, llm)
       REAL :: tsol(klon), qsol(klon), sn(klon)
       REAL :: tsolsrf(klon,nbsrf), qsolsrf(klon,nbsrf),snsrf(klon,nbsrf) 
@@ -625 +626 @@
 
 C init pour traceurs
-      call iniadvtrac(nq)
+      call infotrac_init
 C Ecriture
       CALL inidissip( lstardis, nitergdiv, nitergrot, niterh   ,
@@ -649 +650 @@
      *                phi,w, pbaru,pbarv,time+iday-dayref   )
        print*,'sortie caldyn0'     
-      CALL dynredem0("start.nc",dayref,phis,nqmx)
+      CALL dynredem0("start.nc",dayref,phis)
       print*,'sortie dynredem0'
-      CALL dynredem1("start.nc",0.0,vvent,uvent,tpot,q3d,nqmx,masse ,
+      CALL dynredem1("start.nc",0.0,vvent,uvent,tpot,q3d,masse ,
      .                            psol)
       print*,'sortie dynredem1' 
@@ -743 +744 @@
       visu_file='Etat0_visu.nc'
       CALL initdynav(visu_file,dayref,anneeref,time_step,
-     .              t_ops, t_wrt, nqmx, visuid)
-      CALL writedynav(visuid, nqmx, itau,vvent ,
+     .              t_ops, t_wrt, visuid)
+      CALL writedynav(visuid, itau,vvent ,
      .                uvent,tpot,pk,phi,q3d,masse,psol,phis)
       else