inisim.f

Go to the documentation of this file.

!                       *****************
                        SUBROUTINE inisim
!                       *****************
!
     &(npoin1,typele,inop5,ngeo)
!
!***********************************************************************
! PROGICIEL : STBTEL  V5.2            25/02/92    J.C. GALLAND
!***********************************************************************
!
!   FONCTION  : RECHERCHE LES NOMBRES TOTAUX DE NOEUDS ET D'ELEMENTS DU
!               MAILLAGE, AINSI QUE LA LONGUEUR DU TABLEAU NOP5 DANS LE
!               FICHIER GEOMETRIE DE SIMAIL
!
!-----------------------------------------------------------------------
!                             ARGUMENTS
! .________________.____.______________________________________________
! |      NOM       |MODE|                   ROLE
! |________________|____|______________________________________________
! | NPOIN1         |<-- | NOMBRE REEL DE POINTS DU MAILLAGE
! | TYPELE         |<-- | TYPE D'ELEMENTS
! | IA             |--> | TABLEAU DE TRAVAIL POUR LA LECTURE DE LA SD
! | INOP5          |<-- | DIMENSION DU TABLEAU NOP5 (CONTENANT LES IKLE)
! |                |    | DU FICHIER GEOMETRIE DE SIMAIL
! |________________|____|______________________________________________
! | 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
! |  FICH:         |    |
! |    NRES        |--> | NUMERO DU CANAL DU FICHIER DE SERAFIN
! |    NGEO        |--> | NUMERO DU CANAL DU FICHIER MAILLEUR
! |    NDYNAM      |--> | NUMERO DU CANAL DU FICHIER DYNAM DE TELEMAC
! |    NFO1        |--> | NUMERO DU CANAL DU FICHIER TRIANGLE TRIGRID
! |________________|____|______________________________________________
! 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,INOP5
      INTEGER, INTENT(IN) :: NGEO
      CHARACTER(LEN=11), INTENT(INOUT) :: TYPELE
!
      INTEGER IA(32),LONG,I,J,NTASD,MESHT,MESHQ
!
!=======================================================================
! INITIALISATION
!=======================================================================
!
      rewind ngeo
      npoin  = 0
      npoin1 = 0
      nelem  = 0
!
!=======================================================================
! LECTURE PARTIELLE DU 1ER ENREGISTREMENT DE LA SD
!=======================================================================
!
      READ(ngeo,err=110,end=120) long,(ia(i),i=1,min(long,32))
!
!=======================================================================
! LECTURE PARTIELLE DU TABLEAU NOP0
! RECHERCHE DU NOMBRE DE TABLEAUX ASSOCIES, NTASD
!=======================================================================
!
      READ(ngeo,err=110,end=120) long,(ia(i),i=1,32)
      ntasd = ia(32)
!
!=======================================================================
! LECTURE DU TABLEAU NOP1 ET DES TABLEAUX ASSOCIES
!=======================================================================
!
      IF (ntasd.GT.0) THEN
        READ(ngeo,err=110,end=120) long,(ia(i),i=1,min(long,32))
        DO i=1,ntasd
          READ(ngeo,err=110,end=120) long,(ia(j),j=1,min(long,32))
        ENDDO
      ENDIF
!
!=======================================================================
! LECTURE DU TABLEAU NOP2
! LECTURE DU NOMBRE DE POINTS, DU NOMBRE D'ELEMENTS, DU TYPE D'ELEMENT
! ET DE LA LONGUEUR DU TABLEAU NOP5 (TABLEAU DES IKLE)
! AFFECTATION DES VALEURS LUES AUX VARIABLES CONCERNEES
!=======================================================================
!
      READ(ngeo,err=110,end=120) long,(ia(i),i=1,min(long,32))
      npoin1 = ia(15)
      nelem  = ia(5)
      mesht  = ia(8)
      meshq  = ia(9)
      inop5  = ia(26)
!
      npoin = npoin1
!
!=======================================================================
! MISE DES VALEURS DE MESH AU STANDARD TELEMAC
!=======================================================================
!
      IF (meshq.NE.0) THEN
        mesh = 2
        ndp  = 4
        typele = 'QUADRANGLES'
      ELSEIF (mesht.NE.0) THEN
        mesh = 3
        ndp  = 3
        typele = 'TRIANGLES  '
      ELSE
        WRITE(lu,3100)
        CALL plante(1)
        stop
      ENDIF
!
      GOTO 20
!
 110  CONTINUE
      WRITE(lu,4100)
 120  CONTINUE
      WRITE(lu,4200)
!
20    CONTINUE
!
!=======================================================================
! IMPRESSION DES RESULTATS
!=======================================================================
!
 3100 FORMAT(/,'*************************************************'
     &      ,/,'INISIM : TELEMAC DOESN''T WORK WITH MESHES MIXING '
     &      ,/,'         TRIANGLES AND QUADRILATERALS',
     &      /,'**************************************************')
 4100 FORMAT(//,'*************************************',/,
     &          'INISIM : ERROR IN READING FILE SIMAIL',/,
     &          '*************************************',//)
 4200 FORMAT(//,'*************************************************',/,
     &          'INISIM : ATTEMPT TO READ AFTER END OF FILE SIMAIL',/,
     &          '*************************************************',//)
!
      RETURN
      END