Esope source files

   1 : C TRJCOU    SOURCE    CHAT      05/01/13    03:49:24     5004           
   2 :       SUBROUTINE TRJCOU(UE,IZCOU,IZFACE,IEL,JEL,NUCENT)
   3 : C                                                                               
   4 :       IMPLICIT INTEGER(I-N)
   5 :       IMPLICIT REAL*8 (A-H,O-Z)                                                 
   6 : C                                                                               
   7 : CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC        
   8 : C                                                                               
   9 : C     CALCUL DE DTCO(JEL) POUR AVOIR UN NOMBRE DE COURANT DE 1 DANS IEL
  10 : C     UE  VITESSES MAX POUR L'ELEMENT JEL
  11 : C     IZCOU POINTEUR DU SEGMENT CONTENANT LE NOMBRE DE COURANT
  12 : C     IZFACE POINTEUR DU MAILLAGE CONTENANT LES CENTRES DES FACES
  13 : C     IEL NUMERO DE L ELEMENT CONSIDERE DANS LE SOUS MAILLAGE
  14 : C     JEL NUMERO DE L ELEMENT CONSIDERE DANS LE MAILLAGE COMPLET
  15 : C     NUCENT NUMERO DU NOEUD CENTRE DE L ELEMENT
  16 : C                                                                               
  17 : C                                                                               
  18 : CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC        
  19 : C                                                                               
  20 : C                                                                               
  21 : -INC CCOPTIO
  22 : -INC SMCOORD
  23 : -INC SMELEME
  24 :       POINTEUR IZFACE.MELEME
  25 :       SEGMENT IZUMAX
  26 :            REAL*8 UMAX(NBREL)
  27 :       ENDSEGMENT
  28 : C                                                                               
  29 :       SEGMENT  IZCOU 
  30 :            REAL*8 DTCO(NEL),COU
  31 :       ENDSEGMENT
  32 : C                                                                               
  33 : C                                                                               
  34 :       NF=IZFACE.NUM(/1)
  35 : C        ON VA DETERMINER DIAM : LA PLUS PETITE DISTANCE ENTRE LE
  36 : C        POINT CENTRE DE L ELEMENT ET LE MILIEU DES FACES
  37 :       IPCENT=(NUCENT-1)*(IDIM+1)
  38 :       NF1=(IZFACE.NUM(1,IEL)-1)*(IDIM+1)
  39 :       VINT=0.D0
  40 :       DO 5 I=1,IDIM
  41 :         VINT=VINT+(XCOOR(IPCENT+I)-XCOOR(NF1+I))**2
  42 :     5 CONTINUE
  43 :       DIAM=VINT
  44 :       DO 20 J=2,NF
  45 :            NF1=(IZFACE.NUM(J,IEL)-1)*(IDIM+1)
  46 :            VINT=0.D0
  47 :            DO 15 I=1,IDIM
  48 :              VINT=VINT+(XCOOR(IPCENT+I)-XCOOR(NF1+I))**2
  49 :    15      CONTINUE
  50 :            IF(VINT.LT.DIAM)DIAM=VINT
  51 :    20 CONTINUE
  52 :       DIAM=SQRT(DIAM)*2.D0
  53 : C
  54 : C
  55 : C     UE=UMAX(JEL)
  56 : C     write(6,*)' DIAM UE ',DIAM,UE,JEL
  57 :       DTCO(JEL)=DIAM/UE
  58 :       RETURN                                                                    
  59 :       END                                                                       
  60 :  
  61 :