Sources Esope | Cast3M

q4crsp
C Q4CRSP    SOURCE    CHAT      06/03/29    21:30:52     5360
C
      SUBROUTINE Q4CRSP(NB1,NB2,ITRNOE,NBNMAX,NBE,NBN,iarr)
C     **********************************************************************
C     OBJET Q4CRSP : CREATION D'UN MAILLAGE GRILLE (CAS SIMPLE)
C     OBJET          MEME NOMBRE DE NOEUDS SUR LES COTES OPPOSES.
C
C     EN ENTREE   :
C       NB1,NB2 : NOMBRE DE NOEUDS SUR UNE LIGNE, SUR UNE COLONNE
C       ITRNOE : LE TABLEAU DES MAILLES (A REMPLIR)
C                 SA TAILLE DOIT ETRE EGALE A (NB1 * NB2 * NBNMAX)
C       NBNMAX  : NOMBRE MAXIMUM DE NOEUDS PAR ELEMENTS (4 OU +)
C
C     EN SORTIE   :
C       ITRNOE : TABLEAU DE DEFINITION DES ELEMENTS (NUMERO DES NOEUDS)
C       NBE     : NOMBRE D'ELEMENTS
C       NBN     : NOMBRE DE NOEUDS
C       iarr    : CODE D'ERREUR
C                 -1 SI NB1, NB2 OU NBNMAX EST INCORRECT
C     **********************************************************************
      IMPLICIT INTEGER(I-N)
      INTEGER NB1,NB2
      INTEGER ITRNOE(*),NBNMAX,NBE,NBN,iarr
C
      INTEGER I,J
C
      iarr = -1
      NBE = 0
      NBN = 0
      IF( NBNMAX.LT.4 )GOTO 999
      IF((NB1.LE.0).OR.(NB2.LE.0))GOTO 999
      iarr = 0
      NBE = NB1*NB2
      NBN = (NB1+1)*(NB2+1)
      DO 20 I=1,NB1
        DO 10 J=1,NB2
          ITRNOE(((I-1)*NB2+(J-1))*NBNMAX+1 ) =
     >      J+((NB2+1)*(I-1))
          ITRNOE(((I-1)*NB2+(J-1))*NBNMAX+2 ) =
     >      J+((NB2+1)*(I-1)) + 1
          ITRNOE(((I-1)*NB2+(J-1))*NBNMAX+3 ) =
     >      J+((NB2+1)* I   ) + 1
          ITRNOE(((I-1)*NB2+(J-1))*NBNMAX+4 ) =
     >      J+((NB2+1)* I   )
   10   CONTINUE
   20  CONTINUE
  999 END