Esope source files

   1 : C MURLOC    SOURCE    CHAT      05/01/13    01:55:49     5004           
   2 : C MURLOC    SOURCE    PP        27/04/04   21:25:46     4501
   3 :       SUBROUTINE MURLOC(XE,NBNO,LHOOK,LRE,BPSS,XH,BGENE)
   4 : C=======================================================================
   5 : C
   6 : C    -TEST DE PLANEITE DES FACES DE L'ELEMENT JOI2
   7 : C    -CALCUL DE LA MATRICE DE PASSAGE BPSS
   8 : C    -CALCUL DES COORDONNEES LOCALES XEL
   9 : C          ROUTINE FORTRAN PUR
  10 : C          CODE S. FELIX SEPTEMBRE 92
  11 : C=======================================================================
  12 : C  INPUT
  13 : C     XE     = COORDONNEES  DE L ELEMENT
  14 : C     NBNO   = NOMBRE DE NOEUDS DE L'ELEMENT
  15 : C  OUTPUT
  16 : C     BPSS   = MATRICE DE PASAGE REPERE GLOBAL/REPERE LOCAL
  17 : C     XH     = HAUTEUR DE L'ELEMENT
  18 : C     BGENE  = MATRICE B
  19 : C
  20 : C=======================================================================
  21 :       IMPLICIT INTEGER(I-N)
  22 :       IMPLICIT REAL*8(A-H,O-Z)
  23 :       DIMENSION XE(3,*),BPSS(3,3),BGENE(LHOOK,LRE)
  24 :       DIMENSION S1(2),SN(2)
  25 :       PARAMETER(ZERO=0.D0, UN=1.D0)
  26 : C
  27 : C----------     CALCUL DE LA MATRICE DE  PASSAGE
  28 : C
  29 :        S1(1) = XE(1,3) - XE(1,1)
  30 :        S1(2) = XE(2,3) - XE(2,1)
  31 : C
  32 :        XH = SQRT((S1(1)**2) + (S1(2)**2))
  33 :        S1(1) = S1(1) / XH
  34 :        S1(2) = S1(2) / XH
  35 : C
  36 :        SN(1) =-S1(2)
  37 :        SN(2) = S1(1)
  38 : C
  39 :        DO I=1,2
  40 :           BPSS(1,I) = S1(I)
  41 :           BPSS(2,I) = SN(I)
  42 :        ENDDO
  43 : C
  44 : C----------     CALCUL DE BGENE
  45 : C
  46 :        DO I=1,LHOOK
  47 :          DO J=1,LRE
  48 :            BGENE(I,J)=ZERO
  49 :          ENDDO
  50 :        ENDDO
  51 : C
  52 :        BGENE(1,1)=-BPSS(1,1)
  53 :        BGENE(1,2)=-BPSS(1,2)
  54 :        BGENE(1,7)= UN
  55 : C
  56 :        BGENE(2,3)=-UN
  57 :        BGENE(2,8)= UN
  58 : C
  59 :        BGENE(3,1)=-BPSS(2,1) 
  60 :        BGENE(3,2)=-BPSS(2,2) 
  61 :        BGENE(3,4)= BPSS(2,1) 
  62 :        BGENE(3,5)= BPSS(2,2) 
  63 :        BGENE(3,8)= XH
  64 : C
  65 :        BGENE(4,4)= BPSS(1,1) 
  66 :        BGENE(4,5)= BPSS(1,2) 
  67 :        BGENE(4,7)=-UN        
  68 : C
  69 :        BGENE(5,6)= UN
  70 :        BGENE(5,8)=-UN
  71 : C
  72 :       RETURN
  73 :       END
  74 : 
  75 :  
  76 :