init_zero.f

Go to the documentation of this file.

!                   ********************
                    SUBROUTINE init_zero
!                   ********************
!
!
!***********************************************************************
! SISYPHE   V6P1                                   21/07/2011
!***********************************************************************
!
!brief    INITIALISES VARIABLES.
!
!history  JMH
!+        04/06/2008
!+
!+   INITIALISATION OF MASDEP
!
!history  C. VILLARET
!+
!+        V5P9
!+
!
!history  N.DURAND (HRW), S.E.BOURBAN (HRW)
!+        13/07/2010
!+        V6P0
!+   Translation of French comments within the FORTRAN sources into
!+   English comments
!
!history  N.DURAND (HRW), S.E.BOURBAN (HRW)
!+        21/08/2010
!+        V6P0
!+   Creation of DOXYGEN tags for automated documentation and
!+   cross-referencing of the FORTRAN sources
!
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!
      USE declarations_sisyphe
      USE bief
!
      USE declarations_special
      IMPLICIT NONE
!
!-----------------------------------------------------------------------
!
      INTEGER I
!
!-----------------------------------------------------------------------
!
!========================================================================
!                         INITIALISES
! =======================================================================
!
!---- THE WORKING ARRAYS
!
      IF(npriv > 0) CALL os ('X=0     ', x=prive)
      CALL os('X=0     ', x=t12   )
      CALL os('X=0     ', x=coefpn)
!
!---- THE SEDIMENT VARIABLES :
!
      CALL os('X=0     ', x=qs)
      CALL os('X=0     ', x=qsx)
      CALL os('X=0     ', x=qsy)
      CALL os('X=0     ', x=qsclxc )
      CALL os('X=0     ', x=qsclyc )
      CALL os('X=0     ', x=qsclxs )
      CALL os('X=0     ', x=qsclys )
!
! PROVISIONAL INITIALISATION FOR FIRST OUTPUT IN RESULTS FILE
!
      CALL os('X=0     ', x=qscl )  ! BLOCK OF SIZE NSICLA
      CALL os('X=0     ', x=qs_s )
      CALL os('X=0     ', x=qsxs )
      CALL os('X=0     ', x=qsys )
      CALL os('X=0     ', x=qs_c )
      CALL os('X=0     ', x=qsxc )
      CALL os('X=0     ', x=qsyc )
!
! PROBABLY USEFUL ONLY IF(CHARR) AND WITH FINITE ELEMENTS
!
      CALL os('X=0     ', x=zfcl_c )
!
!---- THE DEPOSITION MASSES FOR EVERY CLASS IN SUSPENSION
!
      DO i=1,nsicla
        masdep(i)=0.d0
      ENDDO
!
      CALL os('X=0     ', x=e    )
      CALL os('X=0     ', x=esomt)
      CALL os('X=0     ', x=cs   )
!
!---- THE HYDRODYNAMIC VARIABLES :
!
      CALL os('X=0     ', x=qu )
      CALL os('X=0     ', x=qv )
      CALL os('X=0     ', x=u2d )
      CALL os('X=0     ', x=v2d )
      CALL os('X=0     ', x=hn )
      CALL os('X=0     ', x=q  )
      CALL os('X=0     ', x=tob)
!
!---- THE WAVE PARAMETERS IF NEED BE
!
!     ALL INITIALISATIONS OF THE WAVES ARE TO BE REMOVED
!     WHEN ALL CHECKS WILL BE DONE
!     SEE BEDLOAD_BAILARD, DIBWAT, BIJKER AND SOULSBY
!
!
!     FW=0.3 CORRESPONDS TO NO WAVES, 0 WOULD DO A LOG(0)
      CALL os('X=C     ', x=fw ,c=0.3d0   )   !
      CALL os('X=0     ', x=hw    )   !
      CALL os('X=0     ', x=tw    )   !
      CALL os('X=C     ', x=thetaw, c=90.d0)  !
      CALL os('X=0     ', x=uw    )   !
      CALL os('X=0     ', x=tobw)     !
!
!-----------------------------------------------------------------------
!
      RETURN
      END