Télécharger askstb.eso

Retour à la liste

Numérotation des lignes :

askstb
  1. C ASKSTB SOURCE CHAT 07/10/22 21:15:10 5921
  2. SUBROUTINE ASKSTB(XKSTAB,XK11,XK22,XK33,XK12,XK21,XK13,XK23,
  3. ! XK31,XK32)
  4. IMPLICIT REAL*8 (A-H,O-Z)
  5. IMPLICIT INTEGER (I-N)
  6.  
  7.  
  8. * SUBROUTINE ASSE_KSTAB(KSTAB,K11,K22,K33,K12,K21,K13,K23,
  9. * ! K31,K32)
  10. C ANTOINE LEGAY LE 31/01/01
  11. * REAL *8 KSTAB(24,24),K11(8,8),K22(8,8),K33(8,8),K12(8,8)
  12. * REAL *8 K21(8,8),K13(8,8),K23(8,8),K31(8,8),K32(8,8)
  13. * INTEGER I,J
  14. CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
  15. C
  16. C ASSEMBLAGE DE LA MATRICE DE STABILISATION
  17. C AVEC LES 3 MATRICES 8*8 XK11 XK22 ET XK33
  18. C
  19.  
  20.  
  21. DIMENSION XKSTAB(24,24),XK11(8,8),XK22(8,8),XK33(8,8),XK12(8,8)
  22. DIMENSION XK21(8,8),XK13(8,8),XK23(8,8),XK31(8,8),XK32(8,8)
  23. CALL ZDANUL(XKSTAB,576)
  24. DO I=1,8
  25. DO J=1,8
  26. XKSTAB(I,J)=XK11(I,J)
  27. END DO
  28. END DO
  29. DO I=9,16
  30. DO J=9,16
  31. XKSTAB(I,J)=XK22(I-8,J-8)
  32. END DO
  33. END DO
  34. DO I=17,24
  35. DO J=17,24
  36. XKSTAB(I,J)=XK33(I-16,J-16)
  37. END DO
  38. END DO
  39. DO I=1,8
  40. DO J=9,16
  41. XKSTAB(I,J)=XK12(I,J-8)
  42. END DO
  43. END DO
  44. DO I=9,16
  45. DO J=1,8
  46. XKSTAB(I,J)=XK21(I-8,J)
  47. END DO
  48. END DO
  49. DO I=1,8
  50. DO J=17,24
  51. XKSTAB(I,J)=XK13(I,J-16)
  52. END DO
  53. END DO
  54. DO I=9,16
  55. DO J=17,24
  56. XKSTAB(I,J)=XK23(I-8,J-16)
  57. END DO
  58. END DO
  59. DO I=17,24
  60. DO J=1,8
  61. XKSTAB(I,J)=XK31(I-16,J)
  62. END DO
  63. END DO
  64. DO I=17,24
  65. DO J=9,16
  66. XKSTAB(I,J)=XK32(I-16,J-8)
  67. END DO
  68. END DO
  69. RETURN
  70. END
  71.  
  72.  
  73.  

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