Télécharger relrig.eso

Retour à la liste

Numérotation des lignes :

relrig
  1. C RELRIG SOURCE CB215821 20/11/25 13:38:58 10792
  2. SUBROUTINE RELRIG()
  3. IMPLICIT REAL*8 (A-H,O-Z)
  4. IMPLICIT INTEGER (I-N)
  5. C***********************************************************************
  6. C NOM : RELRIG
  7. C DESCRIPTION : Lecture des données et écriture des résultats
  8. C pour la création des contraintes associées à un
  9. C opérateur linéaire
  10. C
  11. C LANGAGE : ESOPE
  12. C AUTEUR : Stéphane GOUNAND (CEA/DRN/DMT/SEMT/LTMF)
  13. C mél : gounand@semt2.smts.cea.fr
  14. C***********************************************************************
  15. C APPELES :
  16. C APPELES (E/S) :
  17. C APPELE PAR :
  18. C***********************************************************************
  19. C SYNTAXE GIBIANE :
  20. C ENTREES :
  21. C SORTIES :
  22. C***********************************************************************
  23. C VERSION : v1, 26/06/03, version initiale
  24. C HISTORIQUE : v1, 26/06/03, création
  25. C HISTORIQUE :
  26. C HISTORIQUE :
  27. C***********************************************************************
  28. C Prière de PRENDRE LE TEMPS de compléter les commentaires
  29. C en cas de modification de ce sous-programme afin de faciliter
  30. C la maintenance !
  31. C***********************************************************************
  32.  
  33. -INC PPARAM
  34. -INC CCOPTIO
  35. -INC SMCHPOI
  36. -INC SMRIGID
  37. POINTEUR MLIN.MRIGID
  38. POINTEUR MCON.MRIGID
  39. * Includes persos
  40. CBEGININCLUDE SMMATASS
  41. SEGMENT MATASS
  42. POINTEUR KJPOPA.MLENTI
  43. POINTEUR LINCPA.MLMOTS
  44. POINTEUR MINCPA.MINC
  45. POINTEUR KJPODA.MLENTI
  46. POINTEUR LINCDA.MLMOTS
  47. POINTEUR MINCDA.MINC
  48. POINTEUR PROFMA.PMORS
  49. POINTEUR VALMA.IZA
  50. ENDSEGMENT
  51. CENDINCLUDE SMMATASS
  52. *
  53. INTEGER IMPR,IRET
  54. *
  55. * Executable statements
  56. *
  57. IMPR=0
  58. IF (IMPR.GT.2) WRITE(IOIMP,*) 'Entrée dans relrig'
  59. *
  60. * Lecture de la matrice des contraintes
  61. *
  62. CALL LIROBJ('RIGIDITE',MLIN,0,IRET)
  63. * Sinon, on va essayer de lire un MATRIK
  64. IF (IRET.EQ.0) THEN
  65. CALL RELRI2
  66. ELSE
  67. * Lecture d'un chpoint éventuel précisant les blocages
  68. CALL LIROBJ('CHPOINT',MCHPOI,0,IRET)
  69. *
  70. * Assemblage de la matrice
  71. *
  72. CALL RELR10(MLIN,MATASS,IMPR,IRET)
  73. IF (IRET.NE.0) GOTO 9999
  74. C*
  75. C* Création de la matrice des contraintes avec deux multiplicateurs
  76. C*
  77. C CALL RELR20(MATASS,MCON,IMPR,IRET)
  78. C IF (IRET.NE.0) GOTO 9999
  79. *
  80. * Création de la matrice des contraintes avec un multiplicateur
  81. *
  82. CALL RELR30(MATASS,MCHPOI,MCON,MCHPO2,IMPR,IRET)
  83. IF (IRET.NE.0) GOTO 9999
  84. *
  85. * Ecriture de la matrice des contraintes
  86. *
  87. IF (MCHPO2.NE.0) THEN
  88. CALL ECROBJ('CHPOINT',MCHPO2)
  89. ENDIF
  90. CALL ECROBJ('RIGIDITE',MCON)
  91. ENDIF
  92. *
  93. * Normal termination
  94. *
  95. IRET=0
  96. RETURN
  97. *
  98. * Format handling
  99. *
  100. *
  101. * Error handling
  102. *
  103. 9999 CONTINUE
  104. IRET=1
  105. WRITE(IOIMP,*) 'An error was detected in subroutine relrig'
  106. RETURN
  107. *
  108. * End of subroutine RELRIG
  109. *
  110. END
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  

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