--- trunk/libf/bibio/initdynav.f90	2011/01/25 15:11:05	39
+++ trunk/libf/bibio/initdynav.f90	2011/02/22 13:49:36	40
@@ -9,7 +9,9 @@
   subroutine initdynav(day0, anne0, tstep, nq, fileid, t_ops, t_wrt)
 
     ! From initdynav.F, version 1.1.1.1, 2004/05/19 12:53:05
-    !   L. Fairhead, LMD
+    ! L. Fairhead, LMD
+    ! Routine d'initialisation des écritures des fichiers histoires LMDZ
+    ! au format IOIPSL. Initialisation du fichier histoire moyenne.
 
     use calendar, ONLY: ymds2ju
     USE comvert, ONLY : nivsigs
@@ -21,43 +23,24 @@
     USE paramet_m, ONLY : iip1, jjp1
     USE temps, ONLY : itau_dyn
 
-    !   Routine d'initialisation des ecritures des fichiers histoires LMDZ
-    !   au format IOIPSL. Initialisation du fichier histoire moyenne.
+    integer, intent(in):: day0, anne0 ! date de reference
+    real, intent(in):: tstep ! frequence d'ecriture
+    real, intent(in):: t_ops ! frequence de l'operation pour IOIPSL
+    real, intent(in):: t_wrt ! frequence d'ecriture sur le fichier
+    integer, intent(out):: fileid ! ID du fichier netcdf cree
+    integer, intent(in):: nq ! nombre de traceurs
 
-    !   Appels succesifs des routines: histbeg
-    !                                  histhori
-    !                                  histver
-    !                                  histdef
-    !                                  histend
-
-    !   Entree:
-    !      day0,anne0: date de reference
-    !      tstep : frequence d'ecriture
-    !      t_ops: frequence de l'operation pour IOIPSL
-    !      t_wrt: frequence d'ecriture sur le fichier
-    !      nq: nombre de traceurs
-
-    !   Sortie:
-    !      fileid: ID du fichier netcdf cree
-
-    !   Arguments
-    integer day0, anne0
-    real, intent(in):: tstep, t_ops, t_wrt
-    integer fileid
-    integer nq
+    ! Variables locales
     integer thoriid, zvertiid
-
-    !   Variables locales
-
     real zjulian
     integer iq
-    real rlong(iip1,jjp1), rlat(iip1,jjp1)
-    integer ii,jj
+    real rlong(iip1, jjp1), rlat(iip1, jjp1)
+    integer ii, jj
     integer zan, dayref
 
     !----------------------------------------------------
 
-    !  Appel a histbeg: creation du fichier netcdf et initialisations diverses
+    print *, "Call sequence information: initdynav"
 
     zan = anne0
     dayref = day0
@@ -65,80 +48,45 @@
 
     do jj = 1, jjp1
        do ii = 1, iip1
-          rlong(ii,jj) = rlonv(ii) * 180. / pi
-          rlat(ii,jj)  = rlatu(jj) * 180. / pi
+          rlong(ii, jj) = rlonv(ii) * 180. / pi
+          rlat(ii, jj) = rlatu(jj) * 180. / pi
        enddo
     enddo
 
-    call histbeg_totreg('dyn_hist_ave.nc', rlong(:,1), rlat(1,:), &
-         1, iip1, 1, jjp1, &
-         itau_dyn, zjulian, tstep, thoriid, fileid)
-
-
-    !  Appel a histvert pour la grille verticale
-
-    call histvert(fileid, 'sigss', 'Niveaux sigma','Pa', &
-         llm, nivsigs, zvertiid)
-
-    !  Appels a histdef pour la definition des variables a sauvegarder
-
-    !  Vents U
+    call histbeg_totreg('dyn_hist_ave.nc', rlong(:, 1), rlat(1, :), 1, iip1, &
+         1, jjp1, itau_dyn, zjulian, tstep, thoriid, fileid)
+    call histvert(fileid, 'sigss', 'Niveaux sigma', 'Pa', llm, nivsigs, &
+         zvertiid)
 
-    write(6,*)'inithistave',tstep
     call histdef(fileid, 'u', 'vents u scalaires moyennes', &
          'm/s', iip1, jjp1, thoriid, llm, 1, llm, zvertiid, &
          'ave(X)', t_ops, t_wrt)
-
-
-    !  Vents V
-
     call histdef(fileid, 'v', 'vents v scalaires moyennes', &
          'm/s', iip1, jjp1, thoriid, llm, 1, llm, zvertiid, &
          'ave(X)', t_ops, t_wrt)
-
-
-    !  Temperature
-
     call histdef(fileid, 'temp', 'temperature moyennee', 'K', &
          iip1, jjp1, thoriid, llm, 1, llm, zvertiid, &
          'ave(X)', t_ops, t_wrt)
-
-    !  Temperature potentielle
-
     call histdef(fileid, 'theta', 'temperature potentielle', 'K', &
          iip1, jjp1, thoriid, llm, 1, llm, zvertiid, &
          'ave(X)', t_ops, t_wrt)
-
-
-
-    !  Geopotentiel
-
     call histdef(fileid, 'phi', 'geopotentiel moyenne', '-', &
          iip1, jjp1, thoriid, llm, 1, llm, zvertiid, &
          'ave(X)', t_ops, t_wrt)
 
-    !  Traceurs
-
-    DO iq=1,nq
+    ! Traceurs
+    DO iq = 1, nq
        call histdef(fileid, ttext(iq), ttext(iq), '-', &
             iip1, jjp1, thoriid, llm, 1, llm, zvertiid, &
             'ave(X)', t_ops, t_wrt)
     enddo
 
-    !  Masse
-
     call histdef(fileid, 'masse', 'masse', 'kg', &
          iip1, jjp1, thoriid, 1, 1, 1, -99, &
          'ave(X)', t_ops, t_wrt)
-
-    !  Pression au sol
-
     call histdef(fileid, 'ps', 'pression naturelle au sol', 'Pa', &
          iip1, jjp1, thoriid, 1, 1, 1, -99, &
          'ave(X)', t_ops, t_wrt)
-
-    !  Pression au sol
-
     call histdef(fileid, 'phis', 'geopotentiel au sol', '-', &
          iip1, jjp1, thoriid, 1, 1, 1, -99, &
          'ave(X)', t_ops, t_wrt)