initri.f

Go to the documentation of this file.

!                       *****************
                        SUBROUTINE initri
!                       *****************
!
     &( npoin1,typele,ngeo,nfo1)
!
!***********************************************************************
! PROGICIEL : STBTEL V5.2      07.04.1993  P. LANG     (LHF)
!***********************************************************************
!
!     FONCTION  : INITIALISATION DES INFORMATIONS DANS LE CAS DE TRIGRID
!
!-----------------------------------------------------------------------
!                             ARGUMENTS
! .________________.____.______________________________________________
! !      NOM       |MODE!                   ROLE
! !________________|____!______________________________________________
! ! NPOIN1         !<-- ! NOMBRE DE POINTS DU MAILLAGE
! ! TYPELE         !<-- ! TYPE D'ELEMENT DU MAILLAGE (ICI TRIAGLES)
! !________________!____!______________________________________________
! ! COMMON:        !    !
! !  GEO:          !    !
! !    MESH        !<-- ! TYPE DES ELEMENTS DU MAILLAGE
! !    NDP         !<-- ! NOMBRE DE NOEUDS PAR ELEMENTS ( ICI FORCEMENT
! !    NPOIN       !<-- ! NOMBRE TOTAL DE NOEUDS DU MAILLAGE
! !    NELEM       !<-- ! NOMBRE TOTAL D'ELEMENTS DU MAILLAGE
! !  FICH:         !    !
! !    NRES        ! -->! NUMERO DU CANAL DU FICHIER DE SERAFIN
! !    NGEO       ! -->! NUMERO DU CANAL DU FICHIER MAILLEUR
! !    NLIM      ! -->! NUMERO DU CANAL DU FICHIER DYNAM DE TELEMAC
! !    NFO1      ! -->! NUMERO DU CANAL DU FICHIER TRIANGLE DE TRIGRID
! !________________!____!______________________________________________
! MODE : -->(DONNEE NON MODIFIEE), <--(RESULTAT), <-->(DONNEE MODIFIEE)
!-----------------------------------------------------------------------
! APPELE PAR : HOMERE
! APPEL DE :
!***********************************************************************
!
      USE declarations_stbtel, ONLY: nelem,mesh,ndp,npoin
      IMPLICIT NONE
!
      INTEGER, INTENT(IN) :: NGEO, NFO1
      INTEGER, INTENT(INOUT) :: NPOIN1
      CHARACTER(LEN=*), INTENT(INOUT) :: TYPELE
!
      CHARACTER(LEN=1)   ZDUMMY
!
! COMMON
!
!
      rewind(ngeo)
      rewind(nfo1)
      READ (ngeo,*) npoin1
      npoin = npoin1
      nelem = 0
 1    CONTINUE
        READ (nfo1, '(A1)', end=9000) zdummy
        nelem = nelem + 1
      GOTO 1
!
 9000 CONTINUE
      typele = 'TRIANGLES  '
      ndp = 3
      mesh = 3
!
      RETURN
      END