--- trunk/libf/phylmd/phytrac.f90	2008/10/15 16:19:57	20
+++ trunk/libf/phylmd/phytrac.f90	2012/04/20 14:58:43	61
@@ -1,7 +1,5 @@
 module phytrac_m
 
-  ! This module is clean: no C preprocessor directive, no include line.
-
   IMPLICIT none
 
   private
@@ -10,25 +8,24 @@
 contains
 
   SUBROUTINE phytrac(rnpb, itap, lmt_pas, julien, gmtime, firstcal, lafin, &
-       nq_phys, pdtphys, u, v, t_seri, paprs, pplay, pmfu, pmfd, pen_u, &
-       pde_u, pen_d, pde_d, coefh, fm_therm, entr_therm, yu1, yv1, ftsol, &
-       pctsrf, frac_impa, frac_nucl, pphis, pphi, albsol, rh, cldfra, rneb, &
-       diafra, cldliq, itop_con, ibas_con, pmflxr, pmflxs, prfl, psfl, da, &
-       phi, mp, upwd, dnwd, tr_seri, zmasse)
+       nq_phys, pdtphys, u, t_seri, paprs, pplay, pmfu, pmfd, pen_u, pde_u, &
+       pen_d, pde_d, coefh, fm_therm, entr_therm, yu1, yv1, ftsol, pctsrf, &
+       frac_impa, frac_nucl, pphis, albsol, rh, cldfra, rneb, diafra, cldliq, &
+       pmflxr, pmflxs, prfl, psfl, da, phi, mp, upwd, dnwd, tr_seri, zmasse)
 
-    ! From phylmd/phytrac.F, version 1.15 2006/02/21 08:08:30
+    ! From phylmd/phytrac.F, version 1.15 2006/02/21 08:08:30 (SVN revision 679)
 
     ! Authors: Frédéric Hourdin, Abderrahmane Idelkadi, Marie-Alice
     ! Foujols, Olivia
     ! Objet : moniteur général des tendances des traceurs
 
-    ! Remarques :
-    ! 1/ L'appel de "phytrac" se fait avec "nq-2" donc nous avons bien 
+    ! L'appel de "phytrac" se fait avec "nqmx-2" donc nous avons bien
     ! les vrais traceurs (en nombre "nbtr", sans la vapeur d'eau ni l'eau
     ! liquide) dans "phytrac".
-    ! 2/ Le choix du radon et du plomb se fait juste avec un "data" 
-    ! (peu propre).
-    ! Pourrait-on avoir une variable qui indiquerait le type de traceur ?
+
+    ! Modifications pour les traceurs :
+    ! - uniformisation des parametrisations ds phytrac
+    ! - stockage des moyennes des champs necessaires en mode traceur off-line 
 
     use dimens_m, only: llm
     use indicesol, only: nbsrf
@@ -36,23 +33,17 @@
     use clesphys, only: ecrit_tra
     use clesphys2, only: iflag_con
     use abort_gcm_m, only: abort_gcm
-    use YOMCST, only: rg
+    use SUPHEC_M, only: rg
     use ctherm, only: iflag_thermals
     use regr_pr_comb_coefoz_m, only: regr_pr_comb_coefoz
     use phyetat0_m, only: rlat
     use o3_chem_m, only: o3_chem
-    use ini_hist, only: ini_histrac
+    use ini_histrac_m, only: ini_histrac
     use radiornpb_m, only: radiornpb
     use minmaxqfi_m, only: minmaxqfi
-    use numer_rec, only: assert
+    use nr_util, only: assert
     use press_coefoz_m, only: press_coefoz
 
-    ! Arguments:
-
-    !   EN ENTREE:
-
-    !   divers:
-
     logical, intent(in):: rnpb
 
     integer, intent(in):: nq_phys
@@ -61,8 +52,6 @@
     integer, intent(in):: itap  ! number of calls to "physiq"
     integer, intent(in):: lmt_pas ! number of time steps of "physics" per day
     integer, intent(in):: julien !jour julien, 1 <= julien <= 360
-    integer itop_con(klon)
-    integer ibas_con(klon)
     real, intent(in):: gmtime ! heure de la journée en fraction de jour
     real, intent(in):: pdtphys  ! pas d'integration pour la physique (s)
     real, intent(in):: t_seri(klon, llm) ! temperature, in K
@@ -70,8 +59,7 @@
     real, intent(inout):: tr_seri(:, :, :) ! (klon, llm, nbtr)
     ! (mass fractions of tracers, excluding water, at mid-layers)
 
-    real u(klon, llm)
-    real v(klon, llm)
+    real, intent(in):: u(klon, llm)
     real rh(klon, llm)     ! humidite relative
     real cldliq(klon, llm) ! eau liquide nuageuse
     real cldfra(klon, llm) ! fraction nuageuse (tous les nuages)
@@ -88,8 +76,7 @@
     real, intent(in):: pplay(klon, llm)
     ! (pression pour le mileu de chaque couche, en Pa)
 
-    real pphi(klon, llm) ! geopotentiel
-    real pphis(klon)
+    real, intent(in):: pphis(klon)
     logical, intent(in):: firstcal ! first call to "calfis"
     logical, intent(in):: lafin ! fin de la physique
 
@@ -468,12 +455,14 @@
       ! From phylmd/write_histrac.h, version 1.9 2006/02/21 08:08:30
 
       use dimens_m, only: iim, jjm, llm
-      use ioipsl, only: histwrite, histsync
+      use histsync_m, only: histsync
+      use histwrite_m, only: histwrite
       use temps, only: itau_phy
       use iniadvtrac_m, only: tnom
       use comgeomphy, only: airephy
       use dimphy, only: klon
-      use grid_change, only: gr_phy_write_2d, gr_phy_write_3d
+      use grid_change, only: gr_phy_write_2d
+      use gr_phy_write_3d_m, only: gr_phy_write_3d
 
       logical, intent(in):: lessivage
 
@@ -486,7 +475,6 @@
       ! Variables local to the procedure:
       integer it
       integer itau_w   ! pas de temps ecriture
-      REAL zx_tmp_2d(iim, jjm+1), zx_tmp_3d(iim, jjm+1, llm)
       logical, parameter:: ok_sync = .true.
 
       !-----------------------------------------------------
@@ -513,7 +501,7 @@
       ENDDO
 
       CALL histwrite(nid_tra, "pplay", itau_w, gr_phy_write_3d(pplay))
-      CALL histwrite(nid_tra, "t", itau_w, gr_phy_write_3d(t_seri))
+      CALL histwrite(nid_tra, "T", itau_w, gr_phy_write_3d(t_seri))
 
       if (ok_sync) then
          call histsync(nid_tra)