Esope source files

   1 : C ELINIT    SOURCE    CHAT      05/01/12    23:33:16     5004           
   2 : C
   3 : C
   4 :       SUBROUTINE ELINIT(ICC)
   5 : C     **********************************************************************
   6 : C     OBJET ELINIT : INITIALISE LES CONSTANTES DES COMMON 
   7 : C
   8 : C     EN ENTREE :
   9 : C      ICC : LE CODE DE CALCUL
  10 : C     REMARQUE  :
  11 : C       L'INITIALISATION DES CONSTANTES EST REALISE PAR PROCEDURE
  12 : C       PLUTOT QUE PAR UN BLOCK DATA POUR DES RAISONS DE PORTABILITE.
  13 : C       ELINI DOIT ETRE APPELEE AU DEBUT DE CHAQUE PROGRAMME
  14 : C     **********************************************************************
  15 : C      NBCODE   : NOMBRE D'ELEMENTS 
  16 : C      ISCODE(I): CODE STANDARDISE POUR LES ELEMENTS DE CODE I
  17 : C      ISNOM(I) : NOM POUR LES ELEMENTS DE CODE I
  18 : C      ISNBNO(I): NOMBRE DE NOEUDS POUR LES ELEMENTS DE CODE I
  19 : C
  20 : C
  21 : C     **********************************************************************
  22 :       IMPLICIT INTEGER(I-N)
  23 :       INTEGER ICC
  24 :       INTEGER ICODE,I,J
  25 : C
  26 :       INTEGER   MAXEL,MAXTYP,MAXORD
  27 :       PARAMETER   (MAXEL = 21, MAXTYP = 8, MAXORD = 2)
  28 :       INTEGER NBCODE,ISCODE,ISNBNO,ISTYO,ISTY,ISTO,ISIDE
  29 :       CHARACTER*5 ISNOM  
  30 :       COMMON /ELCODE/ ISCODE(MAXEL)
  31 :       COMMON /ELNBNO/ ISNBNO(MAXEL)
  32 :       COMMON /ELTYPO/ ISTYO(MAXTYP,MAXORD)
  33 :       COMMON /ELTYPC/ ISTY(MAXEL)
  34 :       COMMON /ELORDR/ ISTO(MAXEL)
  35 :       COMMON /ELIDE/  ISIDE(MAXEL)
  36 :       COMMON /ELNOM/  ISNOM(MAXEL)
  37 : C
  38 :       INTEGER     NBCODL,NBCARL
  39 :       PARAMETER   (NBCODL = 21, NBCARL = 5)
  40 :       INTEGER     ISTYPO(MAXTYP,MAXORD),ISTYPE(NBCODL)
  41 :       INTEGER     ISDIDE(NBCODL),ISORDR(NBCODL),ISNBNL(NBCODL)
  42 : C
  43 : C     ---- NOMBRE DE NOEUDS POUR CHAQUE ELEMENT -----------------------
  44 :       DATA  ISNBNL /2,3,3,6,4,6,8,9,4,10,6,9,12,15,18,8,12,16,18,20,27 /
  45 : C
  46 : C     ---- NUMERO DE L'ELEMENT DE TYPE DONNE (TRIANGLE, TETRA...)
  47 : C          POUR UN ORDRE DONNE (LINEAIRE ,QUADRATYQUE)
  48 :       DATA ISTYPO / 0, 1, 3, 5, 9,11,-1,16,
  49 :      >             -1, 2, 4, 7,10,14,-1,20 /
  50 : C
  51 : C     ---- TYPE POUR CHAQUE ELEMENT -----------------------
  52 :       DATA  ISTYPE /2,2,3,3,4,4,4,4,5,5,6,6,6,6,6,8,8,8,8,8,8 /
  53 : C
  54 : C     ---- ORDRE POUR CHAQUE ELEMENT -----------------------
  55 :       DATA  ISORDR /1,2,1,2,1,-1,2,3,1,2,1,-1,-1,2,-1,1,-1,-1,-1,2,4 /
  56 : C
  57 : C     ---- DIMENSION POUR CHAQUE ELEMENT -----------------------
  58 :       DATA  ISDIDE /1,1,2,2,2,2,2,2,3,3,3,3,3,3,3,3,3,3,3,3,3 /
  59 : C
  60 :       NBCODE = NBCODL
  61 :       DO 10 ICODE=1,NBCODE
  62 :         ISNBNO(ICODE) = ISNBNL(ICODE)
  63 :         ISTY(ICODE)   = ISTYPE(ICODE)
  64 :         ISTO(ICODE)   = ISORDR(ICODE)
  65 :         ISIDE(ICODE)  = ISDIDE(ICODE)
  66 :    10 CONTINUE
  67 :       DO 30 I=1,MAXORD
  68 :       DO 20 J=1,MAXTYP
  69 :        ISTYO(J,I) = ISTYPO(J,I)
  70 :    20 CONTINUE
  71 :    30 CONTINUE
  72 : C   
  73 : C              ===============
  74 : C     --------- POUR VIPLEF3D ------------
  75 : C              ===============
  75 :       CALL EL2000(ISCODE,ISNOM,NBCODE)
  76 : C
  77 :  9999 END
  78 :  
  79 :