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