--- trunk/libf/bibio/initfluxsto.f	2008/03/31 13:58:05	9
+++ trunk/libf/bibio/initfluxsto.f	2010/04/01 09:07:28	30
@@ -2,11 +2,9 @@
 ! $Header: /home/cvsroot/LMDZ4/libf/bibio/initfluxsto.F,v 1.1.1.1 2004/05/19 12:53:05 lmdzadmin Exp $
 !
       subroutine initfluxsto
-     .  (infile,tstep,t_ops,t_wrt,nq,
+     .  (tstep,t_ops,t_wrt,nq,
      .                    fileid,filevid,filedid)
 
-       USE IOIPSL
-
 C
 C   Routine d'initialisation des ecritures des fichiers histoires LMDZ
 C   au format IOIPSL
@@ -19,7 +17,6 @@
 C
 C   Entree:
 C
-C      infile: nom du fichier histoire a creer
 C      day0,anne0: date de reference
 C      tstep: duree du pas de temps en seconde
 C      t_ops: frequence de l'operation pour IOIPSL
@@ -35,6 +32,8 @@
 C =====================================================================
 C
 C   Declarations
+       USE calendar
+       use histcom
        use dimens_m
       use paramet_m
       use comconst
@@ -48,21 +47,19 @@
 
 C   Arguments
 C
-      character*(*) infile
       integer*4 itau
-      real tstep, t_ops, t_wrt
+      real, intent(in):: tstep
+      real t_ops, t_wrt
       integer fileid, filevid,filedid
       integer nq,ndex(1)
       real nivd(1)
 
 C   Variables locales
-C
-      integer tau0
       real zjulian
       character*3 str
       character*10 ctrac
       integer iq
-      real rlong(iip1,jjp1), rlat(iip1,jjp1),rl(1,1)
+      real rlong(iip1,jjp1), rlat(iip1,jjp1)
       integer uhoriid, vhoriid, thoriid, zvertiid,dhoriid,dvertiid
       integer ii,jj
       integer zan, idayref
@@ -81,7 +78,6 @@
       zan = annee_ref
       idayref = day_ref
       CALL ymds2ju(zan, 1, idayref, 0.0, zjulian)
-      tau0 = itau_dyn
 	
 	do jj = 1, jjp1
         do ii = 1, iip1
@@ -90,9 +86,9 @@
         enddo
       enddo
  
-      call histbeg_totreg(infile, iip1, rlong(:,1), jjp1, rlat(1,:),
+      call histbeg_totreg('fluxstoke', rlong(:,1), rlat(1,:),
      .             1, iip1, 1, jjp1,
-     .             tau0, zjulian, tstep, uhoriid, fileid)
+     .             itau_dyn, zjulian, tstep, uhoriid, fileid)
 C
 C  Creation du fichier histoire pour la grille en V (oblige pour l'instant,
 C  IOIPSL ne permet pas de grilles avec des nombres de point differents dans 
@@ -106,14 +102,13 @@
         enddo
       enddo
 
-      call histbeg_totreg('fluxstokev.nc', iip1, rlong(:,1), jjm,
-     .             rlat(1,:),1, iip1, 1, jjm,
-     .             tau0, zjulian, tstep, vhoriid, filevid)
+      call histbeg_totreg('fluxstokev.nc', rlong(:,1),
+     .             rlat(1,:jjm),1, iip1, 1, jjm,
+     .             itau_dyn, zjulian, tstep, vhoriid, filevid)
 	
-	rl(1,1) = 1.	
-      call histbeg_regular('defstoke.nc', 1, rl, 1, rl,
+      call histbeg_totreg('defstoke.nc', (/1./), (/1./),
      .             1, 1, 1, 1,
-     .             tau0, zjulian, tstep, dhoriid, filedid)
+     .             itau_dyn, zjulian, tstep, dhoriid, filedid)
 
 C
 C  Appel a histhori pour rajouter les autres grilles horizontales
@@ -125,7 +120,7 @@
         enddo
       enddo
 
-      call histhori(fileid, iip1, rlong, jjp1, rlat, 'scalar',
+      call histhori_regular(fileid, iip1, rlong, jjp1, rlat, 'scalar',
      .              'Grille points scalaires', thoriid)
 	
 C
@@ -148,23 +143,23 @@
 C  Appels a histdef pour la definition des variables a sauvegarder
 	
 	CALL histdef(fileid, "phis", "Surface geop. height", "-",
-     .                iip1,jjp1,thoriid, 1,1,1, -99, 32,
+     .                iip1,jjp1,thoriid, 1,1,1, -99,
      .                "once", t_ops, t_wrt)
 
          CALL histdef(fileid, "aire", "Grid area", "-",
-     .                iip1,jjp1,thoriid, 1,1,1, -99, 32,
+     .                iip1,jjp1,thoriid, 1,1,1, -99,
      .                "once", t_ops, t_wrt)
 	
 	CALL histdef(filedid, "dtvr", "tps dyn", "s",
-     .                1,1,dhoriid, 1,1,1, -99, 32,
+     .                1,1,dhoriid, 1,1,1, -99,
      .                "once", t_ops, t_wrt)
         
          CALL histdef(filedid, "istdyn", "tps stock", "s",
-     .                1,1,dhoriid, 1,1,1, -99, 32,
+     .                1,1,dhoriid, 1,1,1, -99,
      .                "once", t_ops, t_wrt)
          
          CALL histdef(filedid, "istphy", "tps stock phy", "s",
-     .                1,1,dhoriid, 1,1,1, -99, 32,
+     .                1,1,dhoriid, 1,1,1, -99,
      .                "once", t_ops, t_wrt)
 
 
@@ -173,33 +168,33 @@
 C
       call histdef(fileid, 'masse', 'Masse', 'kg',
      .             iip1, jjp1, thoriid, llm, 1, llm, zvertiid,
-     .             32, 'inst(X)', t_ops, t_wrt)
+     .             'inst(X)', t_ops, t_wrt)
 C
 C  Pbaru 
 C
       call histdef(fileid, 'pbaru', 'flx de masse zonal', 'kg m/s',
      .             iip1, jjp1, uhoriid, llm, 1, llm, zvertiid,
-     .             32, 'inst(X)', t_ops, t_wrt)
+     .             'inst(X)', t_ops, t_wrt)
 
 C
 C  Pbarv 
 C
       call histdef(filevid, 'pbarv', 'flx de masse mer', 'kg m/s',
      .             iip1, jjm, vhoriid, llm, 1, llm, zvertiid,
-     .             32, 'inst(X)', t_ops, t_wrt)
+     .             'inst(X)', t_ops, t_wrt)
 C
 C  w 
 C
       call histdef(fileid, 'w', 'flx de masse vert', 'kg m/s',
      .             iip1, jjp1, thoriid, llm, 1, llm, zvertiid,
-     .             32, 'inst(X)', t_ops, t_wrt)
+     .             'inst(X)', t_ops, t_wrt)
 
 C
 C  Temperature potentielle
 C
       call histdef(fileid, 'teta', 'temperature potentielle', '-',
      .             iip1, jjp1, thoriid, llm, 1, llm, zvertiid,
-     .             32, 'inst(X)', t_ops, t_wrt)
+     .             'inst(X)', t_ops, t_wrt)
 C
 
 C
@@ -207,7 +202,7 @@
 C
       call histdef(fileid, 'phi', 'geopotentiel instantane', '-',
      .             iip1, jjp1, thoriid, llm, 1, llm, zvertiid,
-     .             32, 'inst(X)', t_ops, t_wrt)
+     .             'inst(X)', t_ops, t_wrt)
 C
 C  Fin
 C