inistb.f

Go to the documentation of this file.

!                       *****************
                        SUBROUTINE inistb
!                       *****************
!
     &(npoin1,typele,maille,precis,ngeo,nsec2,nsec11,nsec12)
!
!***********************************************************************
! PROGICIEL : STBTEL V5.2            09/08/89    J.C. GALLAND
!***********************************************************************
!
!   FONCTION  : RECHERCHE LES NOMBRES TOTAUX DE NOEUDS ET D'ELEMENTS DU
!               MAILLAGE DANS LE FICHIER UNIVERSEL DE SUPERTAB
!
!-----------------------------------------------------------------------
!                             ARGUMENTS
! .________________.____.______________________________________________
! |      NOM       |MODE|                   ROLE
! |________________|____|______________________________________________
! | NPOIN1         |<-- | NOMBRE REEL DE POINTS DU MAILLAGE
! |                |    | (NPOIN REPRESENTE L'INDICE MAX DES NOEUDS CAR
! |                |    | SUPERTAB LAISSE DES TROUS DANS LA NUMEROTATION
! | TYPELE         |<-- | TYPE D'ELEMENTS
! | MAILLE         |<-- | NOM DU MAILLEUR
! | PRECIS         |<-- | FORMAT DE LECTURE DES COORDONNEES DES NOEUDS
! |________________|____|______________________________________________
! | COMMON:        |    |
! |  GEO:          |    |
! |    MESH        |<-- | TYPE DES ELEMENTS DU MAILLAGE
! |    NDP         | -->| NOMBRE DE NOEUDS PAR ELEMENTS
! |    NPOIN       |<-- | NOMBRE TOTAL DE NOEUDS DU MAILLAGE
! |    NELEM       |<-- | NOMBRE TOTAL D'ELEMENTS DU MAILLAGE
! |    NPMAX       | -->| DIMENSION EFFECTIVE DES TABLEAUX X ET Y
! |                |    | (NPMAX = NPOIN + 0.1*NELEM)
! |    NELMAX      | -->| DIMENSION EFFECTIVE DES TABLEAUX CONCERNANT
! |                |    | LES ELEMENTS (NELMAX = NELEM + 0.2*NELEM)
! |  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 TRIGRID
! |  SECT:         |    |
! |    NSEC11      |--> | INDICATEUR DU SECTEUR CONTENANT LES NOEUDS
! |                |    | (LECTURE EN SIMPLE PRECISION)
! |    NSEC12      |--> | INDICATEUR DU SECTEUR CONTENANT LES NOEUDS
! |                |    | (LECTURE EN DOUBLE PRECISION)
! |    NSEC2       |--> | INDICATEUR DU SECTEUR CONTENANT LES ELEMENTS
! |    NSEC3       |--> | INDICATEUR DU SECTEUR CONTENANT LE TITRE
! |________________|____|______________________________________________
! MODE : -->(DONNEE NON MODIFIEE), <--(RESULTAT), <-->(DONNEE MODIFIEE)
!-----------------------------------------------------------------------
! APPELE PAR : HOMERE
! APPEL DE : -
!***********************************************************************
!
      USE declarations_special
      USE declarations_stbtel, ONLY: nelem,mesh,ndp,npoin
      IMPLICIT NONE
!
      INTEGER, INTENT(INOUT) :: NPOIN1
      CHARACTER(LEN=11), INTENT(INOUT) :: TYPELE
      CHARACTER(LEN=9), INTENT(IN) ::  MAILLE
      CHARACTER(LEN=6), INTENT(INOUT) ::  PRECIS
      INTEGER, INTENT(IN) :: NSEC11 , NSEC12 , NGEO, NSEC2
!
      DOUBLE PRECISION X1
      INTEGER NPOIN2
      INTEGER N1, NSEC
      INTEGER N2
      INTEGER INDI11 , INDI12 , INDIC2
!
      CHARACTER(LEN=2)  MOINS1
      CHARACTER(LEN=4)  BLANC*4
!
! COMMON
!
!
!=======================================================================
! INITIALISATION
!=======================================================================
!
      rewind ngeo
      npoin1  = 0
      npoin2  = 0
      nelem   = 0
      indi11  = 0
      indi12  = 0
      indic2  = 0
!
!=======================================================================
! LECTURE SEQUENTIELLE DU FICHIER ET RECHERCHE DES INDICATEURS
! NSEC1 ET NSEC2
!=======================================================================
!
 10   READ(ngeo,1000,err=110,end=120) blanc,moins1
      IF (moins1.NE.'-1'.OR.blanc.NE.'    ') GOTO 10
 1000 FORMAT(a4,a2)
!
 20   READ(ngeo,2000,err=110,end=120) nsec
      IF (nsec.EQ.-1) THEN
        GOTO 20
!
!=======================================================================
! NOMBRE DE POINTS
!=======================================================================
!
! LECTURE EN SIMPLE PRECISION
!
      ELSE IF (nsec.EQ.nsec11) THEN
        indi11 = 1
!
 30     CONTINUE
        READ(ngeo,2000,err=110,end=120) nsec
!
        IF (nsec.NE.-1) THEN
          npoin1 = npoin1+1
          npoin2 = max0(nsec,npoin1)
          GOTO 30
        ELSE
          GOTO 50
        ENDIF
!
! LECTURE EN DOUBLE PRECISION
!
      ELSE IF (nsec.EQ.nsec12) THEN
        indi12 = 1
!
 31     CONTINUE
        READ(ngeo,2000,err=110,end=120) nsec
!
        IF (nsec.NE.-1) THEN
          npoin1 = npoin1+1
          npoin2 = max0(nsec,npoin1)
        ELSE
          GOTO 50
        ENDIF
!
        READ(ngeo,4000,err=110,end=120) x1
!
        GOTO 31
!
!=======================================================================
! NOMBRE ET TYPE D'ELEMENTS
!=======================================================================
!
      ELSE IF (nsec.EQ.nsec2) THEN
        indic2 = 1
        IF (maille.EQ.'SUPERTAB4') THEN
!  LECTURE AU FORMAT SUPERTAB VERSION 4
          READ(ngeo,3000,err=110,end=120) n1,n2,mesh
        ELSE
!  LECTURE AU FORMAT SUPERTAB VERSION 6
          READ(ngeo,3000,err=110,end=120) n1,mesh,n2
        ENDIF
        nelem = 1
 40     READ(ngeo,2000,err=110,end=120) nsec
        IF (nsec.NE.-1) THEN
          nelem = nelem+1
          GOTO 40
        ELSE
          GOTO 50
        ENDIF
      ENDIF
!
 50   IF ((indi11.EQ.1.OR.indi12.EQ.1).AND.indic2.EQ.1) THEN
        GOTO 60
      ELSE
        GOTO 10
      ENDIF
!
 110  CONTINUE
      WRITE(lu,4100)
      CALL plante(1)
      stop
!
 120  CONTINUE
      IF (indi11.NE.1.AND.indi12.NE.1) WRITE(lu,4200)
      IF (indi12.NE.1) WRITE(lu,4300)
      stop
!
 60   CONTINUE
!
!=======================================================================
! AFFECTATION DES VALEURS LUES AUX VARIABLES CONCERNEES
!=======================================================================
!
      IF (indi11.EQ.1) precis='SIMPLE'
      IF (indi12.EQ.1) precis='DOUBLE'
!
      nelem =nelem / 2
!
      npoin = npoin2
!
!=======================================================================
! MISE DES VALEURS DE MESH AU STANDARD TELEMAC
!=======================================================================
!
      IF (mesh.EQ.94) THEN
        mesh = 2
        ndp  = 4
        typele = 'QUADRANGLES'
      ELSEIF (mesh.EQ.91) THEN
        mesh = 3
        ndp  = 3
        typele = 'TRIANGLES  '
      ELSE
        WRITE(lu,3140) mesh
 3140   FORMAT(' INISTB : TYPE OF MESH NOT AVAILABLE IN TELEMAC,
     &           MESH = ',i4)
      ENDIF
!
!=======================================================================
! IMPRESSION DES RESULTATS
!=======================================================================
!
 4100 FORMAT(//,'****************************************',/,
     &          'ERROR IN READING UNIVERSAL FILE (INISTB)',/,
     &          '****************************************')
 4200 FORMAT(//,'********************************************',/,
     &          'END OF THE UNIVERSAL FILE : NO NODE (INISTB)',/,
     &          '********************************************')
 4300 FORMAT(//,'***********************************************',/,
     &          'END OF THE UNIVERSAL FILE : NO ELEMENT (INISTB)',/,
     &          '***********************************************')
 2000 FORMAT(i10)
 3000 FORMAT(3i10)
 4000 FORMAT(d25.16)
!
      RETURN
      END