Télécharger la plaquette Cast3M


fr en

Télécharger kcacu2.eso

Retour à la liste

Numérotation des lignes : 

kcacu2

  1. C KCACU2    SOURCE    CB215821  16/04/21    21:17:29     8920
  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 = SQRT(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.  
  86.  
  87.  
  88.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales