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