--- trunk/libf/bibio/initfluxsto.f90	2011/01/25 15:11:05	39
+++ trunk/Sources/bibio/initfluxsto.f	2015/04/29 15:47:56	134
@@ -23,18 +23,23 @@
 
   !   L. Fairhead, 03/99                                             
 
-  USE calendar
-  USE histcom
-  USE dimens_m
-  USE paramet_m
   USE comconst
-  USE comvert
-  USE logic
   USE comgeom
-  USE serre
-  USE temps, ONLY : annee_ref, day_ref, itau_dyn
-  USE ener
+  use conf_gcm_m
+  USE dimens_m
+  USE disvert_m
+  use dynetat0_m, only: day_ref, annee_ref
+  USE histbeg_totreg_m, ONLY : histbeg_totreg
+  USE histdef_m, ONLY : histdef
+  USE histend_m, ONLY : histend
+  use histhori_regular_m, only: histhori_regular
+  use histsync_m, only: histsync
+  USE histvert_m, ONLY : histvert
   USE nr_util, ONLY : pi
+  USE paramet_m
+  USE serre
+  USE temps, ONLY : itau_dyn
+  use ymds2ju_m, only: ymds2ju
 
   IMPLICIT NONE
 
@@ -48,13 +53,12 @@
 
   !   Variables locales                                                   
   REAL zjulian
-  CHARACTER*3 str
-  CHARACTER*10 ctrac
+  CHARACTER(len=3) str
+  CHARACTER(len=10) ctrac
   INTEGER iq
   REAL rlong(iip1, jjp1), rlat(iip1, jjp1)
   INTEGER uhoriid, vhoriid, thoriid, zvertiid, dhoriid, dvertiid
-  INTEGER ii, jj
-  INTEGER zan, idayref
+  INTEGER ii, jj, l
   LOGICAL ok_sync
 
   !---------------------------------------------------------
@@ -66,9 +70,7 @@
 
   !  Appel a histbeg: creation du fichier netcdf et initialisations diverses
 
-  zan = annee_ref
-  idayref = day_ref
-  CALL ymds2ju(zan, 1, idayref, 0.0, zjulian)
+  CALL ymds2ju(annee_ref, 1, day_ref, 0.0, zjulian)
 
   DO jj = 1, jjp1
      DO ii = 1, iip1
@@ -111,14 +113,14 @@
 
   !  Appel a histvert pour la grille verticale                            
 
-  CALL histvert(fileid, 'sig_s', 'Niveaux sigma', 'sigma_level', llm, &
-       nivsigs, zvertiid)
+  CALL histvert(fileid, 'sig_s', 'Niveaux sigma', 'sigma_level', &
+       (/(real(l), l = 1, llm)/), zvertiid)
   ! Pour le fichier V                                                     
-  CALL histvert(filevid, 'sig_s', 'Niveaux sigma', 'sigma_level', llm, &
-       nivsigs, zvertiid)
+  CALL histvert(filevid, 'sig_s', 'Niveaux sigma', 'sigma_level', &
+       (/(real(l), l = 1, llm)/), zvertiid)
   ! pour le fichier def                                                   
   nivd(1) = 1
-  CALL histvert(filedid, 'sig_s', 'Niveaux sigma', 'sigma_level', 1, nivd, &
+  CALL histvert(filedid, 'sig_s', 'Niveaux sigma', 'sigma_level', nivd, &
        dvertiid)
 
   !  Appels a histdef pour la definition des variables a sauvegarder