condim_susp.f

Go to the documentation of this file.

!                   **********************
                    SUBROUTINE condim_susp
!                   **********************
!
     &(cs,cs0,nsicla)
!
!***********************************************************************
! SISYPHE   V6P1                                   21/07/2011
!***********************************************************************
!
!brief    INITIALISES THE SUSPENDED SEDIMENT CONCENTRATION
!+               (CONDIM_SISYPHE.F IS READ EVEN IF CHARR=NO).
!
!history  M. GONZALES DE LINARES
!+        2004
!+        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
!
!history  C.VILLARET (EDF-LNHE), P.TASSI (EDF-LNHE)
!+        19/07/2011
!+        V6P1
!+  Name of variables
!+
!
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!| CS             |<->| SUSPENDED SEDIMENT CONCENTRATION
!| CS0            |-->| INITIAL CONCENTRATIONS (CONSTANT VALUES)/CLASS
!| NSICLA         |-->| NUMBER OF SIZE CLASSES FOR BED MATERIALS
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!
      USE bief
!
      USE declarations_special
      IMPLICIT NONE
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      INTEGER, INTENT(IN)           :: NSICLA
      DOUBLE PRECISION,INTENT(IN)   :: CS0(nsicla)
!
      TYPE(bief_obj), INTENT(INOUT) :: CS
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      INTEGER I
!
!-----------------------------------------------------------------------
!
!  --------------------------------------------------------------
!  INITIALISES THE ARRAYS THAT HAVE NOT BEEN READ IN THE RESULTS FILE:
!  --------------------------------------------------------------
!
      IF(nsicla.GT.0) THEN
        DO i=1,nsicla
          CALL os('X=C     ',x=cs%ADR(i)%P,c=cs0(i))
        ENDDO
      ENDIF
!
!-----------------------------------------------------------------------
!
      RETURN
      END