Esope source files

   1 : C GENJLF    SOURCE    CHAT      05/01/13    00:16:57     5004           
   2 :       SUBROUTINE GENJLF(P,SZ,NELM,FLA, IFA,NPTO, NUM,SZONE,NUMEL)
   3 : C--------------------------------------------------------------------
   4 : C
   5 : C
   6 : C     PP 9/97 
   7 : C     Pierre Pegon/JRC Ispra
   8 : C--------------------------------------------------------------------
   9 :       IMPLICIT INTEGER(I-N)
  10 :       IMPLICIT REAL*8 (A-H,O-Z)
  11 : C
  12 :       INTEGER P,SZ,FLA ,SZONE
  13 :       DIMENSION P(NPTO,1),SZ(1),NELM(1),FLA(1)
  14 :       DIMENSION NUM(1)
  15 :       DIMENSION N14TO4(4,6),N16TO4(4,3),N25TO4(4)
  16 :       DIMENSION N16TO3(3,2),N23TO3(3,4),N25TO3(3,4)
  17 :       DATA N14TO4/1,2,6,5, 2,3,7,6, 3,4,8,7, 4,1,5,8, 1,2,3,4, 5,6,7,8/
  18 :       DATA N16TO4/1,2,5,4, 2,3,6,5, 3,1,4,6/
  19 :       DATA N25TO4/1,2,3,4/
  20 :       DATA N16TO3/1,2,3, 4,5,6/
  21 :       DATA N23TO3/1,2,4, 2,3,4, 3,1,4, 1,2,3/
  22 :       DATA N25TO3/1,2,5, 2,3,5, 3,4,5, 4,1,5/
  23 : C
  24 :       RETURN
  25 : C
  26 :       ENTRY GF1424(P,SZ,NELM,FLA, IFA,NPTO, NUM,SZONE,NUMEL)
  27 :       DO IE1=1,6
  28 :         IFA=IFA+1
  29 :         SZ  (IFA)=SZONE
  30 :         NELM(IFA)=NUMEL
  31 :         FLA (IFA)=0
  32 :         DO IE2=1,4
  33 :           P(IE2,IFA)=NUM(N14TO4(IE2,IE1))
  34 :           FLA(IFA)=FLA(IFA)+P(IE2,IFA)
  35 :         ENDDO
  36 :       ENDDO
  37 :       RETURN
  38 : C
  39 :       ENTRY GF1624(P,SZ,NELM,FLA, IFA,NPTO, NUM,SZONE,NUMEL)
  40 :       DO IE1=1,3
  41 :         IFA=IFA+1
  42 :         SZ  (IFA)=SZONE
  43 :         NELM(IFA)=NUMEL
  44 :         FLA (IFA)=0
  45 :         DO IE2=1,4
  46 :           P(IE2,IFA)=NUM(N16TO4(IE2,IE1))
  47 :           FLA(IFA)=FLA(IFA)+P(IE2,IFA)
  48 :         ENDDO
  49 :       ENDDO
  50 :       RETURN
  51 : C
  52 :       ENTRY GF2524(P,SZ,NELM,FLA, IFA,NPTO, NUM,SZONE,NUMEL)
  53 :       IFA=IFA+1
  54 :       SZ  (IFA)=SZONE
  55 :       NELM(IFA)=NUMEL
  56 :       FLA (IFA)=0
  57 :       DO IE2=1,4
  58 :         P(IE2,IFA)=NUM(N25TO4(IE2))
  59 :         FLA(IFA)=FLA(IFA)+P(IE2,IFA)
  60 :       ENDDO
  61 :       RETURN
  62 : C
  63 :       ENTRY GF1623(P,SZ,NELM,FLA, IFA,NPTO, NUM,SZONE,NUMEL)
  64 :       DO IE1=1,2
  65 :         IFA=IFA+1
  66 :         SZ  (IFA)=SZONE
  67 :         NELM(IFA)=NUMEL
  68 :         FLA (IFA)=0
  69 :         DO IE2=1,3
  70 :           P(IE2,IFA)=NUM(N16TO3(IE2,IE1))
  71 :           FLA(IFA)=FLA(IFA)+P(IE2,IFA)
  72 :         ENDDO
  73 :       ENDDO
  74 :       RETURN
  75 : C
  76 :       ENTRY GF2323(P,SZ,NELM,FLA, IFA,NPTO, NUM,SZONE,NUMEL)
  77 :       DO IE1=1,4
  78 :         IFA=IFA+1
  79 :         SZ  (IFA)=SZONE
  80 :         NELM(IFA)=NUMEL
  81 :         FLA (IFA)=0
  82 :         DO IE2=1,3
  83 :           P(IE2,IFA)=NUM(N23TO3(IE2,IE1))
  84 :           FLA(IFA)=FLA(IFA)+P(IE2,IFA)
  85 :         ENDDO
  86 :       ENDDO
  87 :       RETURN
  88 : C
  89 :       ENTRY GF2523(P,SZ,NELM,FLA, IFA,NPTO, NUM,SZONE,NUMEL)
  90 :       DO IE1=1,4
  91 :         IFA=IFA+1
  92 :         SZ  (IFA)=SZONE
  93 :         NELM(IFA)=NUMEL
  94 :         FLA (IFA)=0
  95 :         DO IE2=1,3
  96 :           P(IE2,IFA)=NUM(N25TO3(IE2,IE1))
  97 :           FLA(IFA)=FLA(IFA)+P(IE2,IFA)
  98 :         ENDDO
  99 :       ENDDO
 100 :       RETURN
 101 : C
 102 :       END
 103 :  
 104 :