--- trunk/Sources/dyn3d/inithist.f	2017/01/12 12:31:31	212
+++ trunk/Sources/dyn3d/inithist.f	2018/01/23 15:49:10	253
@@ -8,85 +8,59 @@
 
   subroutine inithist(tstep, nq, t_ops, t_wrt)
 
-    ! From inithist.F, version 1.1.1.1 2004/05/19 12:53:05
+    ! From inithist.F, version 1.1.1.1, 2004/05/19 12:53:05
     ! L. Fairhead, LMD, 03/99
 
-    ! Routine d'initialisation des écritures des fichiers histoires
-    ! LMDZ au format IOIPSL.
+    ! Routine d'initialisation des écritures des fichiers histoires au
+    ! format IOIPSL.
 
     USE dimens_m, ONLY: jjm, llm
     USE disvert_m, ONLY: presnivs
     use dynetat0_m, only: day_ref, annee_ref, rlatu, rlatv, rlonu, rlonv
-    USE histbeg_totreg_m, ONLY : histbeg_totreg
-    USE histdef_m, ONLY : histdef
-    USE histend_m, ONLY : histend
-    USE histvert_m, ONLY : histvert
+    use histbeg_totreg_m, only: histbeg_totreg
+    USE histdef_m, ONLY: histdef
+    USE histend_m, ONLY: histend
+    USE histvert_m, ONLY: histvert
     USE iniadvtrac_m, ONLY: ttext
     USE nr_util, ONLY: pi
     USE paramet_m, ONLY: iip1, jjp1
     USE temps, ONLY: itau_dyn
-    USE ymds2ju_m, ONLY: ymds2ju
+    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
 
-    ! Variables locales:
-    real zjulian
+    ! Local:
+    real julian
     integer iq
-    real rlong(iip1, jjp1), rlat(iip1, jjp1)
     integer uhoriid, vhoriid, thoriid, zvertiid
-    integer ii, jj
 
     !-----------------------------------------------------------------------
 
-    CALL ymds2ju(annee_ref, 1, day_ref, 0.0, zjulian)
+    print *, "Call sequence information: inithist"
+    CALL ymds2ju(annee_ref, 1, day_ref, 0., julian)
 
-    do jj = 1, jjp1
-       do ii = 1, iip1
-          rlong(ii, jj) = rlonu(ii) * 180. / pi
-          rlat(ii, jj) = rlatu(jj) * 180. / pi
-       enddo
-    enddo
-
-    call histbeg_totreg("dyn_histu.nc", rlong(:,1), rlat(1,:), 1, iip1, 1, &
-         jjp1, itau_dyn, zjulian, tstep, uhoriid, histuid)
-
-    do jj = 1, jjm
-       do ii = 1, iip1
-          rlong(ii, jj) = rlonv(ii) * 180. / pi
-          rlat(ii, jj) = rlatv(jj) * 180. / pi
-       enddo
-    enddo
+    call histbeg_totreg("dyn_histu.nc", rlonu * 180. / pi, rlatu * 180. / pi, &
+         1, iip1, 1, jjp1, itau_dyn, julian, tstep, 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', rlong(:, 1), rlat(1, :jjm), 1, iip1, &
-         1, jjm, itau_dyn, zjulian, tstep, vhoriid, histvid)
-
-    do jj = 1, jjp1
-       do ii = 1, iip1
-          rlong(ii, jj) = rlonv(ii) * 180. / pi
-          rlat(ii, jj) = rlatu(jj) * 180. / pi
-       enddo
-    enddo
-
-    call histbeg_totreg("dyn_hist.nc", rlong(:, 1), rlat(1, :), 1, iip1, 1, &
-         jjp1, itau_dyn, zjulian, tstep, thoriid, histid)
-
-    ! Appel a histvert pour la grille verticale
+    call histbeg_totreg('dyn_histv.nc', rlonv * 180. / pi, rlatv * 180. / pi, &
+         1, iip1, 1, jjm, itau_dyn, julian, tstep, vhoriid, histvid)
 
-    call histvert(histid, 'presnivs', 'Niveaux pression','mb', presnivs/100., &
-         zvertiid,'down')
-    call histvert(histvid, 'presnivs', 'Niveaux pression','mb', &
-         presnivs/100., zvertiid,'down')
-    call histvert(histuid, 'presnivs', 'Niveaux pression','mb', &
-         presnivs/100., zvertiid,'down')
+    call histbeg_totreg("dyn_hist.nc", rlonv * 180. / pi, rlatu * 180. / pi, &
+         1, iip1, 1, jjp1, itau_dyn, julian, tstep, thoriid, histid)
 
-    ! Appels a histdef pour la definition des variables a sauvegarder
+    call histvert(histid, 'presnivs', 'Niveaux pression', 'mb', presnivs/100., &
+         zvertiid, 'down')
+    call histvert(histvid, 'presnivs', 'Niveaux pression', 'mb', &
+         presnivs/100., zvertiid, 'down')
+    call histvert(histuid, 'presnivs', 'Niveaux pression', 'mb', &
+         presnivs/100., zvertiid, 'down')
 
     call histdef(histuid, 'u', 'vent u', 'm/s', iip1, jjp1, uhoriid, llm, 1, &
          llm, zvertiid, 'inst(X)', t_ops, t_wrt)
@@ -98,7 +72,7 @@
          llm, 1, llm, zvertiid, 'inst(X)', t_ops, t_wrt)
 
     ! Traceurs
-    DO iq=1, nq
+    DO iq = 1, nq
        call histdef(histid, ttext(iq), ttext(iq), '-', iip1, jjp1, thoriid, &
             llm, 1, llm, zvertiid, 'inst(X)', t_ops, t_wrt)
     enddo