Esope source files

   1 : C RELRIG    SOURCE    GOUNAND   12/07/23    21:15:09     7441           
   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 : -INC CCOPTIO
  33 : -INC SMCHPOI
  34 : -INC SMRIGID
  35 :       POINTEUR MLIN.MRIGID
  36 :       POINTEUR MCON.MRIGID
  37 : * Includes persos
  38 : CBEGININCLUDE SMMATASS
  39 :       SEGMENT MATASS
  40 :       POINTEUR KJPOPA.MLENTI
  41 :       POINTEUR LINCPA.MLMOTS
  42 :       POINTEUR MINCPA.MINC
  43 :       POINTEUR KJPODA.MLENTI
  44 :       POINTEUR LINCDA.MLMOTS
  45 :       POINTEUR MINCDA.MINC
  46 :       POINTEUR PROFMA.PMORS
  47 :       POINTEUR VALMA.IZA
  48 :       ENDSEGMENT
  49 : CENDINCLUDE SMMATASS
  50 : *
  51 :       INTEGER IMPR,IRET
  52 : *     
  53 : * Executable statements
  54 : *
  55 :       IMPR=0
  56 :       IF (IMPR.GT.2) WRITE(IOIMP,*) 'Entrée dans relrig'
  57 : *
  58 : * Lecture de la matrice des contraintes
  59 : *
  59 :       CALL LIROBJ('RIGIDITE',MLIN,0,IRET)
  60 : * Sinon, on va essayer de lire un MATRIK
  61 :       IF (IRET.EQ.0) THEN
  61 :          CALL RELRI2
  62 :       ELSE
  63 : * Lecture d'un chpoint éventuel précisant les blocages
  63 :          CALL LIROBJ('CHPOINT',MCHPOI,0,IRET)
  64 : *
  65 : * Assemblage de la matrice
  66 : *
  66 :          CALL RELR10(MLIN,MATASS,IMPR,IRET)
  67 :          IF (IRET.NE.0) GOTO 9999
  68 : C*
  69 : C* Création de la matrice des contraintes avec deux multiplicateurs
  70 : C*
  70 : C         CALL RELR20(MATASS,MCON,IMPR,IRET)
  71 : C         IF (IRET.NE.0) GOTO 9999
  72 : *
  73 : * Création de la matrice des contraintes avec un multiplicateur
  74 : *
  74 :          CALL RELR30(MATASS,MCHPOI,MCON,MCHPO2,IMPR,IRET)
  75 :          IF (IRET.NE.0) GOTO 9999
  76 : *
  77 : * Ecriture de la matrice des contraintes
  78 : *
  79 :          IF (MCHPO2.NE.0) THEN
  79 :             CALL ECROBJ('CHPOINT',MCHPO2)
  80 :          ENDIF
  80 :          CALL ECROBJ('RIGIDITE',MCON)
  81 :       ENDIF
  82 : *
  83 : * Normal termination
  84 : *
  85 :       IRET=0
  86 :       RETURN
  87 : *
  88 : * Format handling
  89 : *
  90 : *
  91 : * Error handling
  92 : *
  93 :  9999 CONTINUE
  94 :       IRET=1
  95 :       WRITE(IOIMP,*) 'An error was detected in subroutine relrig'
  96 :       RETURN
  97 : *
  98 : * End of subroutine RELRIG
  99 : *
 100 :       END
 101 : 
 102 :  
 103 :