--- trunk/libf/dyn3d/init_dynzon.f90	2012/01/30 12:54:02	57
+++ trunk/dyn3d/init_dynzon.f	2014/03/05 14:57:53	82
@@ -15,16 +15,17 @@
 
     ! From LMDZ4/libf/dyn3d/bilan_dyn.F, version 1.5 2005/03/16 10:12:17
 
-    USE conf_gcm_m, ONLY : day_step, iperiod, periodav
-    USE histcom, ONLY : histbeg_totreg, histdef, histend, histvert
-    USE calendar, ONLY : ymds2ju
-    USE dimens_m, ONLY : jjm, llm
-    USE comvert, ONLY : presnivs
-    USE comgeom, ONLY : rlatv
-    USE temps, ONLY : annee_ref, day_ref, itau_dyn
-    USE nr_util, ONLY : pi
-
-    ! Arguments:
+    USE conf_gcm_m, ONLY: day_step, iperiod, periodav
+    USE histbeg_totreg_m, ONLY: histbeg_totreg
+    USE histdef_m, ONLY: histdef
+    USE histend_m, ONLY: histend
+    USE histvert_m, ONLY: histvert
+    USE calendar, ONLY: ymds2ju
+    USE dimens_m, ONLY: jjm, llm
+    USE disvert_m, ONLY: presnivs
+    USE comgeom, ONLY: rlatv
+    USE temps, ONLY: annee_ref, day_ref, itau_dyn
+    USE nr_util, ONLY: pi
 
     real, intent(in):: dt_app
 
@@ -34,63 +35,57 @@
     character(len=5), parameter:: unites(nQ) = (/'K    ', 'm2/s2', 'm2/s2', &
          'ang  ', 'm/s  ', 'kg/kg', 'un   '/)
 
-    ! champs de tansport en moyenne zonale
+    ! Champs de tansport en moyenne zonale
     integer itr
-
-    character(len=26) znoml(ntr, nQ)
+    character(len=26) noml(ntr, nQ)
     character(len=12) zunites(ntr, nQ)
-
     character(len=3), parameter:: ctrs(ntr) = (/'   ', 'TOT', 'MMC', 'TRS', &
          'STN'/)
-
     integer iQ
 
     ! Initialisation du fichier contenant les moyennes zonales.
 
-    integer thoriid, zvertiid
-
-    real zjulian
-    integer zan, dayref
-
+    integer horiid, vertiid
+    real julian
+    integer an, dayref
     real rlong(jjm), rlatg(jjm)
 
     !-----------------------------------------------------------------
 
     print *, "Call sequence information: init_dynzon"
 
-    ! initialisation des fichiers
+    ! Initialisation des fichiers
     ! ncum est la frequence de stokage en pas de temps
     ncum = day_step / iperiod * periodav
     dt_cum = ncum * dt_app
 
     ! Initialisation du fichier contenant les moyennes zonales
 
-    zan = annee_ref
+    an = annee_ref
     dayref = day_ref
-    CALL ymds2ju(zan, 1, dayref, 0.0, zjulian)
+    CALL ymds2ju(an, 1, dayref, 0.0, julian)
 
     rlong = 0.
-    rlatg = rlatv*180./pi
+    rlatg = rlatv * 180. / pi
 
     call histbeg_totreg('dynzon', rlong(:1), rlatg, 1, 1, 1, jjm, itau_dyn, &
-         zjulian, dt_cum, thoriid, fileid)
+         julian, dt_cum, horiid, fileid)
 
     ! Appel à histvert pour la grille verticale
 
-    call histvert(fileid, 'presnivs', 'Niveaux sigma', 'mb', llm, presnivs, &
-         zvertiid)
+    call histvert(fileid, 'presnivs', 'Niveaux sigma', 'mb', presnivs, vertiid)
 
     ! Appels à histdef pour la définition des variables à sauvegarder
     do iQ = 1, nQ
        do itr = 1, ntr
           if (itr == 1) then
              znom(itr, iQ) = nom(iQ)
-             znoml(itr, iQ) = nom(iQ)
+             noml(itr, iQ) = nom(iQ)
              zunites(itr, iQ) = unites(iQ)
           else
-             znom(itr, iQ) = ctrs(itr)//'v'//nom(iQ)
-             znoml(itr, iQ) = 'transport : v * '//nom(iQ)//' '//ctrs(itr)
-             zunites(itr, iQ) = 'm/s * '//unites(iQ)
+             znom(itr, iQ) = ctrs(itr) // 'v' // nom(iQ)
+             noml(itr, iQ) = 'transport: v * ' // nom(iQ) // ' ' // ctrs(itr)
+             zunites(itr, iQ) = 'm/s * ' // unites(iQ)
           endif
        enddo
     enddo
@@ -98,34 +93,31 @@
     ! Déclarations des champs avec dimension verticale
     do iQ = 1, nQ
        do itr = 1, ntr
-          call histdef(fileid, znom(itr, iQ), znoml(itr, iQ), &
-               zunites(itr, iQ), 1, jjm, thoriid, llm, 1, llm, zvertiid, &
+          call histdef(fileid, znom(itr, iQ), noml(itr, iQ), &
+               zunites(itr, iQ), 1, jjm, horiid, llm, 1, llm, vertiid, &
                'ave(X)', dt_cum, dt_cum)
        enddo
        ! Déclarations pour les fonctions de courant
-       call histdef(fileid, 'psi'//nom(iQ), 'stream fn. '//znoml(2, iQ), &
-            zunites(2, iQ), 1, jjm, thoriid, llm, 1, llm, zvertiid, &
+       call histdef(fileid, 'psi' // nom(iQ), 'stream fn. ' // noml(2, iQ), &
+            zunites(2, iQ), 1, jjm, horiid, llm, 1, llm, vertiid, &
             'ave(X)', dt_cum, dt_cum)
     enddo
 
     ! Déclarations pour les champs de transport d'air
-    call histdef(fileid, 'masse', 'masse', &
-         'kg', 1, jjm, thoriid, llm, 1, llm, zvertiid, &
-         'ave(X)', dt_cum, dt_cum)
-    call histdef(fileid, 'v', 'v', &
-         'm/s', 1, jjm, thoriid, llm, 1, llm, zvertiid, &
-         'ave(X)', dt_cum, dt_cum)
+    call histdef(fileid, 'masse', 'masse', 'kg', 1, jjm, horiid, llm, 1, &
+         llm, vertiid, 'ave(X)', dt_cum, dt_cum)
+    call histdef(fileid, 'v', 'v', 'm/s', 1, jjm, horiid, llm, 1, llm, &
+         vertiid, 'ave(X)', dt_cum, dt_cum)
     ! Déclarations pour les fonctions de courant
-    call histdef(fileid, 'psi', 'stream fn. MMC ', 'mega t/s', &
-         1, jjm, thoriid, llm, 1, llm, zvertiid, &
-         'ave(X)', dt_cum, dt_cum)
+    call histdef(fileid, 'psi', 'stream fn. MMC ', 'mega t/s', 1, jjm, &
+         horiid, llm, 1, llm, vertiid, 'ave(X)', dt_cum, dt_cum)
 
     ! Déclaration des champs 1D de transport en latitude
     do iQ = 1, nQ
        do itr = 2, ntr
-          call histdef(fileid, 'a'//znom(itr, iQ), znoml(itr, iQ), &
-               zunites(itr, iQ), 1, jjm, thoriid, 1, 1, 1, -99, &
-               'ave(X)', dt_cum, dt_cum)
+          call histdef(fileid, 'a' // znom(itr, iQ), noml(itr, iQ), &
+               zunites(itr, iQ), 1, jjm, horiid, 1, 1, 1, -99, 'ave(X)', &
+               dt_cum, dt_cum)
        enddo
     enddo