Esope source files

   1 : C KCACU2    SOURCE    CHAT      05/01/13    00:52:11     5004           
   2 :       SUBROUTINE  KCACU2(NES,NR,G,VP,NFACE,KRO,KSI,IR,KA,IM)
   3 :       IMPLICIT INTEGER(I-N)
   4 :       IMPLICIT REAL*8 (A-H,O-Z)                                                 
   5 :       DIMENSION VP(NES,NR),G(NR),IR(1),KA(NFACE)
   6 :       DIMENSION KRO(NFACE,NES),KSI(NFACE,NES),IM(NFACE,NFACE)
   7 : C
   8 : C     DIM 2
   9 : C     CALCUL DES FF ET DES DIRECTIONS POUR TOUTE LA SURFACE DE REFERENCE
  10 : C        ( PLAN Z = +1)
  11 : C
  12 : C        NF = 1    X = +1        NF = 2     X = -1
  13 : C        NF = 3    Y = +1        NF = 4     Y = -1
  14 : C
  15 : C     UTILISATION :
  16 : C
  17 : C            D(NF,K,I) = KSI(NF,K) * V(KRO(NF,K),I)
  18 : C
  19 : C           NF : NUMERO DE FACE
  20 : C           K  : DIRECTION DE ESPACE
  21 : C           I  : INDICES DE DISCRETISATION LOCAL A CHAQUE FACE
  22 : C
  23 : C    REM: VP ORIENTE CU CENTRE VERS LA CELLULE
  24 : C    ***
  25 : C
  26 : C
  27 :       NR2 = NR/2
  28 : 
  29 :       DO 30 I = 1,NR2
  30 :            IR(I)      = NR2-I+1
  31 :            IR(I+ NR2) = I
  32 :    30     CONTINUE
  33 : 
  34 :       DO 1  I = 1, NR
  35 :             VP(1,I) = (IR(I)-0.5)/NR
  36 :             RP = VP(1,I)*VP(1,I) + 0.5*0.5
  37 :             RP = DSQRT(RP)
  38 :             IS = 1
  39 :             IF (I.LE.NR2)   IS = -1
  40 :             VP(1,I) = IS * VP(1,I)/RP
  41 :             VP(2,I) =  0.5/RP
  42 :             G(I) = 0.5/(RP*RP)/NR/2
  43 :    1  CONTINUE
  44 : C
  45 : C     TABLEAU DE CORRESPONDANCE
  46 : C                                                                               
  47 :       KRO(1,1) = 2
  48 :       KRO(1,2) = 1
  49 :       KRO(2,1) = 2
  50 :       KRO(2,2) = 1
  51 :       KRO(3,1) = 1
  52 :       KRO(3,2) = 2
  53 :       KRO(4,1) = 1
  54 :       KRO(4,2) = 2
  55 : C
  56 :       DO 4 NF = 1,NFACE
  57 :       DO 5 K = 1,NES
  58 :          KSI(NF,K) = 1
  59 :   5   CONTINUE
  60 :   4   CONTINUE
  61 : C
  62 :       KSI(2,1) = -1
  63 :       KSI(4,2) = -1
  64 : C
  65 :       DO 50 K1=1,NFACE                                                          
  66 :       DO 51 K2=1,NFACE                                                          
  67 :       IM(K1,K2) = 0                                                             
  68 :  51   CONTINUE                                                                  
  69 :  50   CONTINUE                                                                  
  70 :                                                                                 
  71 :       IM(1,4) = 1                                                               
  72 :       IM(1,3) = NR                                                              
  73 :       IM(2,4) = 1                                                               
  74 :       IM(2,3) = NR                                                              
  75 :       IM(3,2) = 1                                                               
  76 :       IM(3,1) = NR                                                              
  77 :       IM(4,2) = 1                                                               
  78 :       IM(4,1) = NR                                                              
  79 :       KA(1) = 2                                                                 
  80 :       KA(2) = 2                                                                 
  81 :       KA(3) = 1                                                                 
  82 :       KA(4) = 1                                                                 
  83 :       RETURN
  84 :       END
  85 :