Esope source files

   1 : C PRFACF    SOURCE    CHAT      05/01/13    02:26:16     5004           
   2 :         SUBROUTINE PRFACF(MATR)
   3 : C
   4 : C       IMPRIME LES FACTEURS DE FORME     (cf. operateur FFOR)
   5 : C           ou  UNE MATRICE DE RAYONNEMENT(cf. operateur RAYE)
   6 : C
   7 : C
   8 : C     FACTEURS DE FORME
   9 : C           NNBEL1 = NOMBRE DE LIGNES + 1
  10 : C           NBEL2  = NOMBRE DE COLONNES
  11 : C     LFACT(NNBEL1) POINTE SUR LE TABLEAU DES SURFACES
  12 : C
  13 :       IMPLICIT INTEGER(I-N)
  14 :         SEGMENT IFACFO
  15 :          INTEGER   LFACT(NNBEL1)
  16 :         ENDSEGMENT
  17 :         SEGMENT LFAC
  18 :          REAL*8    FACT(NBEL2)
  19 :         ENDSEGMENT
  20 :         POINTEUR PSUR.LFAC, PLIG.LFAC
  21 :         POINTEUR MATR.IFACFO
  22 : 
  23 :         SEGACT MATR
  24 :         NBELI = MATR.LFACT(/1)
  25 :         WRITE (6,*) ' '
  26 :         WRITE (6,*) 'FACTEURS DE FORMES ou MATRICE DE RAYONNEMENT '
  27 :         WRITE (6,*) ' '
  28 :         WRITE (6,*) '     nombre d éléments : ',(NBELI-1)
  29 :         DO I=1,(NBELI-1)
  30 : 
  31 :                 WRITE (6,*) 'Elément  numéro ',I
  32 :                 PLIG = MATR.LFACT(I)
  33 :                 SEGACT PLIG 
  34 :                 N2=PLIG.FACT(/1)
  34 :                 CALL UTPRIM(PLIG.FACT,N2)
  35 :                 SEGDES PLIG
  36 : 
  37 :         ENDDO
  38 : 
  39 :         WRITE (6,*) ' '
  40 :         PSUR = MATR.LFACT(NBELI)
  41 :         SEGACT PSUR
  42 :         WRITE (6,*) 'SURFACE DES ELEMENTS'
  43 :         N2=PSUR.FACT(/1)
  43 :         CALL UTPRIM(PSUR.FACT,N2)
  44 :         WRITE (6,*) ' '
  45 : 
  46 :         SEGDES PSUR
  47 :         SEGDES MATR
  48 :         RETURN
  49 :         END
  50 :