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 :