Esope source files

   1 : C INTCN0    SOURCE    PV        09/03/12    21:25:18     6325           
   2 : 
   3 : C=======================================================================
   4 : C=                            I N T C N 0                              =
   5 : *                             -----------                              =
   6 : *
   7 : * FONCTION:
   8 : * ---------
   9 : *     INTEGRATION NUMERIQUE DANS UN DOMAINE UNIDIMENSIONNEL DU PRODUIT:
  10 : *                     COEF. * TRANSPOSEE( N )
  11 : *     COEF. : GRANDEUR PHYSIQUE REPRESENTEE PAR UN CHAMELEM
  12 : *     N     : FONCTIONS DE FORME DE L'ELEMENT MAILLANT LE DOMAINE
  13 : *             CONSIDERE
  14 : *
  15 : * PARAMETRES:   (E)=ENTREE   (S)=SORTIE   (+ = CONTENU DANS UN COMMUN)
  16 : * -----------
  17 : *     IPVANO  (E)  POINTEUR SUR UN SEGMENT MELVAL CONTENANT LES
  18 : *                  VALEURS NODALES DE NOTRE COEFFICIENT
  19 : *     IPGEOM  (E)  POINTEUR SUR UN OBJET MAILLAGE ELEMENTAIRE
  20 : *                  DU DOMAINE D'INTEGRATION
  21 : *     IPINTE  (E)  POINTEUR SUR UN SEGMENT MINTE CONTENANT LES
  22 : *                  CARACTERISTIQUES D'INTEGRATION DES FACES
  23 : *     IPVAEQ  (S)  POINTEUR SUR UN SEGMENT MELVAL CONTENANT LES VALEURS
  24 : *                  NODALES EQUIVALENTES (APRES INTEGRATION NUMERIQUE)
  25 : *
  26 : * VARIABLES:
  27 : * ----------
  28 : *     XE(3,NBPTEL)  =  CORDONNEES DES ELEMENTS DANS LE REPERE GLOBAL
  29 : *     FACE          =  LONGUEUR DE DE L'ELEMENT CONSIDERE
  30 : *     A             =  TABLEAU DE TRAVAIL
  31 : *
  32 : * REMARQUES:
  33 : * ----------
  34 : *     L'UTILISATION DE CE S-P PRESUPPOSE UN PRE ET POST-TRAITEMENT
  35 : *     DES SEGMENTS MELVAL PASSES EN TANT QUE PARAMETRES
  36 : C=======================================================================
  37 : 
  38 :       SUBROUTINE INTCN0 (IPVANO,IPGEOM,IPINTE,IPVAEQ)
  39 : 
  40 :       IMPLICIT INTEGER(I-N)
  41 :       IMPLICIT REAL*8 (A-H,O-Z)
  42 : 
  43 : -INC CCOPTIO
  44 : -INC CCREEL
  45 : -INC SMCHAML
  46 : -INC SMELEME
  47 : -INC SMINTE
  48 : -INC SMCOORD
  49 : 
  50 :       PARAMETER (X2Pi=6.283185307179586476925286766559D0)
  51 :       PARAMETER (X4Pi=12.566370614359172953850573533118D0)
  52 : 
  53 :       DIMENSION XE(3,1)
  54 : 
  55 : *  OPTION INDISPONIBLE
  56 :       IF (IDIM.NE.1) THEN
  56 :         CALL ERREUR(19)
  57 :         RETURN
  58 :       ENDIF
  59 : 
  60 : *  ON RECUPERE LES CARACTERISTIQUES D'INTEGRATION DES FACES
  61 :       MINTE=IPINTE
  62 :       SEGACT,MINTE
  63 :       NBPGAU=POIGAU(/1)
  64 : 
  65 : * ON RECUPERE UN DES MAILLAGES ELEMENTAIRES DE L'ENVELOPPE
  66 :       MELEME=IPGEOM
  67 :       SEGACT,MELEME
  68 :       NBPTEL=NUM(/1)
  69 :       NEL=NUM(/2)
  70 : 
  71 : C*OF  IF ((NBPGAU.NE.1).OR.(NBPTEL.NE.1)) THEN
  72 : C*OF    WRITE(6,*) 'ERREUR FATALE : INTCN0'
  73 : C*OF    GOTO 10
  74 : C*OF  ENDIF
  75 : 
  76 : *  ON RECUPERE LES VALEURS NODALES DU COEFFICIENT
  77 :       MELVA1=IPVANO
  78 :       SEGACT,MELVA1
  79 :       NBPTE1=MELVA1.VELCHE(/1)
  80 :       NEL1=MELVA1.VELCHE(/2)
  81 : C*OF  INMIN=MIN(NBPTE1,1)
  82 : 
  83 : *  INITIALISATION DU MELVAL QUI CONTIENDRA LES VALEURS EQUIVALENTES
  84 :       N1PTEL=NBPTEL
  85 :       N1EL=NEL
  86 :       N2PTEL=0
  87 :       N2EL=0
  88 :       SEGINI,MELVAL
  89 : 
  90 : *     BOUCLE SUR LES ELEMENTS
  91 :       DO iElt=1,NEL
  92 : 
  93 : *  COORDONNEES DES NOEUDS DE L'ELEMENT DANS LE REPERE GLOBAL
  93 :         CALL DOXE(XCOOR,IDIM,NBPTEL,NUM,iElt,XE)
  94 : 
  95 : *  CAS DES ELEMENTS AXISYMETRIQUES et spheriques
  96 :         IF (IFOUR.GE.12.AND.IFOUR.LE.14) THEN
  97 :           A11=X2Pi*XE(1,1)
  98 :         ELSE IF (IFOUR.EQ.15) THEN
  99 :           RR=XE(1,1)
 100 :           A11=X4Pi*RR*RR
 101 :         ELSE
 102 :           A11=1.
 103 :         ENDIF
 104 : C*OF    IEMIN=MIN(NEL1,iElt)
 105 : C*OF    VELCHE(1,iElt)=MELVA1.VELCHE(INMIN,IEMIN)*A11
 106 :         VELCHE(1,iElt)=MELVA1.VELCHE(1,iElt)*A11
 107 :       ENDDO
 108 : 
 109 :       SEGDES,MELVAL,MELVA1
 110 :  10   SEGDES,MELEME,MINTE
 111 :       IPVAEQ=MELVAL
 112 : 
 113 :       RETURN
 114 :       END
 115 : 
 116 :  
 117 :  
 118 :