--- trunk/libf/dyn3d/bilan_dyn.f	2008/02/27 13:16:39	3
+++ trunk/libf/dyn3d/bilan_dyn.f	2008/10/15 16:19:57	20
@@ -178,19 +178,17 @@
 c   ncum est la frequence de stokage en pas de temps
         ncum=dt_cum/dt_app
         if (abs(ncum*dt_app-dt_cum).gt.1.e-5*dt_app) then
-           WRITE(lunout,*)
+           print *,
      .            'Pb : le pas de cumule doit etre multiple du pas'
-           WRITE(lunout,*)'dt_app=',dt_app
-           WRITE(lunout,*)'dt_cum=',dt_cum
+           print *,'dt_app=',dt_app
+           print *,'dt_cum=',dt_cum
            stop
         endif
 
         if (i_sortie.eq.1) then
          file='dynzon'
-         call inigrads(ifile
-     s  ,(/0./),180./pi,0.,0.,rlatv,-90.,90.,180./pi
-     s  ,presnivs,1.
-     s  ,dt_cum,file,'dyn_zon ')
+         call inigrads(ifile ,(/0./),180./pi,0.,0.,rlatv,-90.,90.,
+     $        180./pi ,presnivs,1. ,dt_cum,file,'dyn_zon ')
         endif
 
         nom(itemp)='T'
@@ -223,7 +221,7 @@
       rlong=0.
       rlatg=rlatv*180./pi
        
-      call histbeg_totreg(infile, 1, rlong(:1), jjm, rlatg,
+      call histbeg_totreg(infile, rlong(:1), rlatg,
      .             1, 1, 1, jjm,
      .             tau0, zjulian, dt_cum, thoriid, fileid)
 
@@ -253,18 +251,18 @@
       do iQ=1,nQ
          do itr=1,ntr
       IF (prt_level > 5)
-     . WRITE(lunout,*)'var ',itr,iQ
+     . print *,'var ',itr,iQ
      .      ,znom(itr,iQ),znoml(itr,iQ),zunites(itr,iQ)
             call histdef(fileid,znom(itr,iQ),znoml(itr,iQ),
      .        zunites(itr,iQ),1,jjm,thoriid,llm,1,llm,zvertiid,
-     .        32,'ave(X)',dt_cum,dt_cum)
+     .        'ave(X)',dt_cum,dt_cum)
          enddo
 c   Declarations pour les fonctions de courant
 c      print*,'2HISTDEF'
           call histdef(fileid,'psi'//nom(iQ)
      .      ,'stream fn. '//znoml(itot,iQ),
      .      zunites(itot,iQ),1,jjm,thoriid,llm,1,llm,zvertiid,
-     .      32,'ave(X)',dt_cum,dt_cum)
+     .      'ave(X)',dt_cum,dt_cum)
       enddo
 
 
@@ -272,15 +270,15 @@
 c      print*,'3HISTDEF'
       call histdef(fileid, 'masse', 'masse',
      .             'kg', 1, jjm, thoriid, llm, 1, llm, zvertiid,
-     .             32, 'ave(X)', dt_cum, dt_cum)
+     .             'ave(X)', dt_cum, dt_cum)
       call histdef(fileid, 'v', 'v',
      .             'm/s', 1, jjm, thoriid, llm, 1, llm, zvertiid,
-     .             32, 'ave(X)', dt_cum, dt_cum)
+     .             'ave(X)', dt_cum, dt_cum)
 c   Declarations pour les fonctions de courant
 c      print*,'4HISTDEF'
           call histdef(fileid,'psi','stream fn. MMC ','mega t/s',
      .      1,jjm,thoriid,llm,1,llm,zvertiid,
-     .      32,'ave(X)',dt_cum,dt_cum)
+     .      'ave(X)',dt_cum,dt_cum)
 
 
 c   Declaration des champs 1D de transport en latitude
@@ -289,7 +287,7 @@
          do itr=2,ntr
             call histdef(fileid,'a'//znom(itr,iQ),znoml(itr,iQ),
      .        zunites(itr,iQ),1,jjm,thoriid,1,1,1,-99,
-     .        32,'ave(X)',dt_cum,dt_cum)
+     .        'ave(X)',dt_cum,dt_cum)
          enddo
       enddo
 
@@ -340,7 +338,7 @@
       endif
 
       IF (prt_level > 5)
-     . WRITE(lunout,*)'dans bilan_dyn ',icum,'->',icum+1
+     . print *,'dans bilan_dyn ',icum,'->',icum+1
       icum=icum+1
 
 c   accumulation des flux de masse horizontaux
@@ -406,7 +404,7 @@
          enddo
       enddo
       IF (prt_level > 5)
-     . WRITE(lunout,*)'Apres les calculs fait a chaque pas'
+     . print *,'Apres les calculs fait a chaque pas'
 c=====================================================================
 c   PAS DE TEMPS D'ECRITURE
 c=====================================================================
@@ -414,7 +412,7 @@
 c=====================================================================
 
       IF (prt_level > 5)
-     . WRITE(lunout,*)'Pas d ecriture'
+     . print *,'Pas d ecriture'
 
 c   Normalisation
       do iQ=1,nQ
@@ -534,19 +532,15 @@
       if (i_sortie.eq.1) then
       do iQ=1,nQ
          do itr=1,ntr
-            call histwrite(fileid,znom(itr,iQ),itau,zvQ(:,:,itr,iQ)
-     s      ,jjm*llm,ndex3d)
+            call histwrite(fileid,znom(itr,iQ),itau,zvQ(:,:,itr,iQ))
          enddo
-         call histwrite(fileid,'psi'//nom(iQ),itau,psiQ(:,1:llm,iQ)
-     s      ,jjm*llm,ndex3d)
+         call histwrite(fileid,'psi'//nom(iQ),itau,psiQ(:,1:llm,iQ))
       enddo
 
-      call histwrite(fileid,'masse',itau,zmasse
-     s   ,jjm*llm,ndex3d)
-      call histwrite(fileid,'v',itau,zv
-     s   ,jjm*llm,ndex3d)
+      call histwrite(fileid,'masse',itau,zmasse)
+      call histwrite(fileid,'v',itau,zv)
       psi=psi*1.e-9
-      call histwrite(fileid,'psi',itau,psi(:,1:llm),jjm*llm,ndex3d)
+      call histwrite(fileid,'psi',itau,psi(:,1:llm))
 
       endif
 
@@ -566,8 +560,7 @@
                zavQ(:,itr,iQ)=zavQ(:,itr,iQ)+zvQ(:,l,itr,iQ)*zmasse(:,l)
             enddo
             zavQ(:,itr,iQ)=zavQ(:,itr,iQ)/zamasse(:)
-            call histwrite(fileid,'a'//znom(itr,iQ),itau,zavQ(:,itr,iQ)
-     s      ,jjm*llm,ndex3d)
+            call histwrite(fileid,'a'//znom(itr,iQ),itau,zavQ(:,itr,iQ))
          enddo
       enddo