Esope source files

   1 : C REFREE    SOURCE    CHAT      05/01/13    02:48:04     5004           
   2 :       SUBROUTINE REFREE(NDIM,ITY1,NOEL1,IZSH,XREF,XREEL)
   3 : ***********************************************************************
   4 : *** SP 'REFREE' : permet le passage entre coord reference d'un pt de
   5 : *** l'element de reference et coord reelles de ce pt dans l'element reel
   6 : ***
   7 : *** APPELES 1 = 'SHAPE'
   8 : *** APPELES 2 = aucun
   9 : *** 
  10 : *** E = 'NDIM' dimension de l'espace
  11 : ***     'ITY1' entier caracterisant le type de l'element
  12 : ***     'NOEL1' nombre de noeuds de l'element considere
  13 : ***     'IZSH' segment contenant fonctions de forme et de base ainsi
  14 : ***            que coord reelles des noeuds de l'element considere
  15 : ***     'XREF' coordonnees reference du pt considere
  16 : ***
  17 : *** S = 'XREEL' coordonnees reelles du pt considere
  18 : ***********************************************************************
  19 : 
  20 :       IMPLICIT INTEGER(I-N)
  21 :       IMPLICIT REAL*8 (A-H,O-Z)
  22 :       SEGMENT IZSH
  23 :           REAL*8 SHP(6,MNO9),SHY(12,MNO9),XYZL(3,MNO9)
  24 :       ENDSEGMENT
  25 :       DIMENSION XREF(3),XREEL(3)
  26 : ***   calcul des fonctions de forme au pt de calcul de l'element reference
  26 :       CALL SHAPE(XREF(1),XREF(2),XREF(3),ITY1,SHP,IRET)
  27 : ***   passage en coord reelles via les fonctions de forme et valeurs aux noeuds
  28 :       DO 10 I=1,NDIM
  29 :          XREEL(I)=0.D0
  30 :          DO 20 K=1,NOEL1
  31 :             XREEL(I)=XREEL(I)+SHP(1,K)*XYZL(I,K)
  32 :  20      CONTINUE
  33 :  10   CONTINUE
  34 :       RETURN
  35 :       END
  36 :  
  37 :