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