--- trunk/Sources/dyn3d/inithist.f	2018/01/23 15:49:10	253
+++ trunk/dyn3d/inithist.f	2018/09/06 13:19:51	302
@@ -6,7 +6,7 @@
 
 contains
 
-  subroutine inithist(tstep, nq, t_ops, t_wrt)
+  subroutine inithist(t_ops, t_wrt)
 
     ! From inithist.F, version 1.1.1.1, 2004/05/19 12:53:05
     ! L. Fairhead, LMD, 03/99
@@ -14,7 +14,8 @@
     ! Routine d'initialisation des écritures des fichiers histoires au
     ! format IOIPSL.
 
-    USE dimens_m, ONLY: jjm, llm
+    use comconst, only: dtvr
+    USE dimensions, ONLY: jjm, llm, nqmx
     USE disvert_m, ONLY: presnivs
     use dynetat0_m, only: day_ref, annee_ref, rlatu, rlatv, rlonu, rlonv
     use histbeg_totreg_m, only: histbeg_totreg
@@ -27,8 +28,6 @@
     USE temps, ONLY: itau_dyn
     use ymds2ju_m, ONLY: ymds2ju
 
-    real, intent(in):: tstep ! durée du pas de temps en secondes
-    integer, intent(in):: nq ! nombre de traceurs
     real, intent(in):: t_ops ! fréquence de l'opération pour IOIPSL
     real, intent(in):: t_wrt ! fréquence d'écriture sur le fichier
 
@@ -43,17 +42,17 @@
     CALL ymds2ju(annee_ref, 1, day_ref, 0., julian)
 
     call histbeg_totreg("dyn_histu.nc", rlonu * 180. / pi, rlatu * 180. / pi, &
-         1, iip1, 1, jjp1, itau_dyn, julian, tstep, uhoriid, histuid)
+         1, iip1, 1, jjp1, itau_dyn, julian, dtvr, uhoriid, histuid)
 
     ! Creation du fichier histoire pour la grille en V (oblige pour l'instant,
     ! IOIPSL ne permet pas de grilles avec des nombres de point differents dans
     ! un meme fichier)
 
     call histbeg_totreg('dyn_histv.nc', rlonv * 180. / pi, rlatv * 180. / pi, &
-         1, iip1, 1, jjm, itau_dyn, julian, tstep, vhoriid, histvid)
+         1, iip1, 1, jjm, itau_dyn, julian, dtvr, vhoriid, histvid)
 
     call histbeg_totreg("dyn_hist.nc", rlonv * 180. / pi, rlatu * 180. / pi, &
-         1, iip1, 1, jjp1, itau_dyn, julian, tstep, thoriid, histid)
+         1, iip1, 1, jjp1, itau_dyn, julian, dtvr, thoriid, histid)
 
     call histvert(histid, 'presnivs', 'Niveaux pression', 'mb', presnivs/100., &
          zvertiid, 'down')
@@ -66,13 +65,15 @@
          llm, zvertiid, 'inst(X)', t_ops, t_wrt)
     call histdef(histvid, 'v', 'vent v', 'm/s', iip1, jjm, vhoriid, llm, 1, &
          llm, zvertiid, 'inst(X)', t_ops, t_wrt)
-    call histdef(histid, 'teta', 'temperature potentielle', '-', iip1, jjp1, &
+    call histdef(histid, 'temp', 'temperature', 'K', iip1, jjp1, &
          thoriid, llm, 1, llm, zvertiid, 'inst(X)', t_ops, t_wrt)
-    call histdef(histid, 'phi', 'geopotentiel', '-', iip1, jjp1, thoriid, &
+    call histdef(histid, 'theta', 'temperature potentielle', 'K', iip1, jjp1, &
+         thoriid, llm, 1, llm, zvertiid, 'inst(X)', t_ops, t_wrt)
+    call histdef(histid, 'phi', 'geopotential', 'm2 s-2', iip1, jjp1, thoriid, &
          llm, 1, llm, zvertiid, 'inst(X)', t_ops, t_wrt)
 
     ! Traceurs
-    DO iq = 1, nq
+    DO iq = 1, nqmx
        call histdef(histid, ttext(iq), ttext(iq), '-', iip1, jjp1, thoriid, &
             llm, 1, llm, zvertiid, 'inst(X)', t_ops, t_wrt)
     enddo
@@ -81,8 +82,6 @@
          llm, zvertiid, 'inst(X)', t_ops, t_wrt)
     call histdef(histid, 'ps', 'pression naturelle au sol', 'Pa', iip1, jjp1, &
          thoriid, 1, 1, 1, -99, 'inst(X)', t_ops, t_wrt)
-    call histdef(histid, 'phis', 'geopotentiel au sol', '-', iip1, jjp1, &
-         thoriid, 1, 1, 1, -99, 'inst(X)', t_ops, t_wrt)
 
     call histend(histid)
     call histend(histuid)