INITIAL CONDITIONS

We have used special initial conditions to simulate the avalanche triggering-off. In fact, we used previous work on Rupture de Barrage made in MCIP by former students of the MFN option because the wet snow avalanche can be compared with this sort of flow.

We had to enter the initial conditions in the condin.f file : the "water" height is forced to 5 meters in the first half of the avalanche channel and equals to zero everywhere else.  The modified part of the file is written with bold font in the following listing.

 *****************                       SUBROUTINE CONDIN  C                          ***************** C      *(U,V,H,T,AK,EP,VISC,AUX3,X,Y,ZF,NPOIN,TRAC,TRAC0,ITURB,PROPNU,      * TEMPS,AKEP,CDTINI,COTINI,HAUTIN,NBOR,NPTFR,PRIVE) C C*********************************************************************** C  TELEMAC 2D VERSION 4.0    17/08/94    J-M HERVOUET (LNH) 30 87 80 18 C C*********************************************************************** C C     FONCTION  : INITIALISATION DES TABLEAUX DES GRANDEURS PHYSIQUES  C C   INITIALISATION DE H , LA HAUTEUR D'EAU  C       IF(INCLUS(CDTINI,'COTE NULLE').OR.      *   INCLUS(CDTINI,'ZERO ELEVATION') ) THEN         COTE = 0.D0         CALL OV( 'X=C     ' , H , Y  , Z , COTE , NPOIN )         CALL OV( 'X=X-Y   ' , H , ZF , Z , C    , NPOIN )       ELSEIF(INCLUS(CDTINI,'COTE CONSTANTE').OR.      *       INCLUS(CDTINI,'CONSTANT ELEVATION') ) THEN         COTE = COTINI         CALL OV( 'X=C     ' , H , Y  , Z , COTE , NPOIN )         CALL OV( 'X=X-Y   ' , H , ZF , Z , C    , NPOIN )       ELSEIF(INCLUS(CDTINI,'HAUTEUR NULLE').OR.      *       INCLUS(CDTINI,'ZERO DEPTH') ) THEN         CALL OV( 'X=C     ' , H , Y  , Z , 0.D0 , NPOIN )       ELSEIF(INCLUS(CDTINI,'HAUTEUR CONSTANTE').OR.      *       INCLUS(CDTINI,'CONSTANT DEPTH') ) THEN         CALL OV( 'X=C     ' , H , Y  , Z , HAUTIN , NPOIN )       ELSEIF(INCLUS(CDTINI,'PARTICULIERES').OR.      *       INCLUS(CDTINI,'SPECIAL')        ) THEN         DO I=1,NPOIN          IF (Y(I).LT.375.) THEN             H(I)=0.0          ELSE             H(I)=5.          ENDIF       ENDDO           ELSE         IF(LNG.EQ.1) WRITE(LU,20) CDTINI         IF(LNG.EQ.2) WRITE(LU,21) CDTINI 20      FORMAT(1X,'CONDIN : CONDITION INITIALE INCONNUE :',/,A72) 21      FORMAT(1X,'CONDIN : UNKNOWN INITIAL CONDITION :',/,A72)         CALL PLANTE(0)         STOP         ENDIF Visualisation of the initial condition on water height in the domain :