Initialisation on height and velocity

    This study needs to have special initialisation on water height because most of cells should be dry, moreover a an initialisation had been done on velocity in order to accelerate convergence.

Adaptation on condin subroutine that have been done in princi.f :

C  ZONE A MODIFIER

        IF(LNG.EQ.1) WRITE(LU,10)
        IF(LNG.EQ.2) WRITE(LU,11)

10      FORMAT(1X,'CONDIN : CONDIN MODIFIE')
11      FORMAT(1X,'CONDIN : WITH SPECIAL INITIAL CONDITIONS'
     *         ,/,'         YOU HAVE TO MODIFY CONDIN')

      CALL OV( 'X=C     ' , U , Y , Z , 0.D0 , NPOIN )
      CALL OV( 'X=C     ' , V , Y , Z , 0.D0 , NPOIN )

 DO i=1,NPOIN
   IF ((Y(i).ge.993).AND.(Y(i).le.1007)) then
      H(i)=1.
      U(i)=0.87
    ELSE
       H(i)=0.
    ENDIF
ENDDO

C  FIN DE LA ZONE A MODIFIER

init visualisation (middle of minor bed above, and throught the minor bed) :


Initialisation on Strickler

    This adaptation have been made with the subroutine strche in princi.f :

C-----------------------------------------------------------------------
C
C  STRICKLER
C

       CALL OV( 'X=C     ' , CHESTR , Y , Y , 20.d0 , NPOIN )

    DO j=1,NPOIN

           IF ((Y(j).ge.993).AND.(Y(j).le.1007)) then
                 CHESTR(j)=30.D0
           ENDIF

    ENDDO
C
C-----------------------------------------------------------------------
C


Flow definition

    This definition have been made with the subroutine Q in princi.f (not already tested) :

C-----------------------------------------------------------------------
C
c      Q = DEBIT(I)
C
C-----------------------------------------------------------------------
C
      if (i.eq.1) then
           if (T.le.1200.) then
                 Q = DEBIT(I)
           elseif (T.le.8400.) then
                 Q = (T-1200.)*0.346+10
           elseif (T.le.15600.) then
                 Q = -1*(T-8400.)*0.346+2491.2
           else
                 Q = DEBIT(I)
           endif
      elseif (i.eq.2)
         (somme des H des noeuds de sortie)/(nb de noeuds de sortie) pour connaitre la hauteur d'eau en sortie
         tester sur h  si on est dans le lit mineur Q(h)
         sinon lit majeur Q(h)
      endif