Esope source files

   1 : C DEFFIN    SOURCE    CHAT      05/01/12    22:40:22     5004           
   2 :       SUBROUTINE DEFFIN(INPLAS,TSOM,NSOM,NCOMP,NINV,NINCMA,
   3 :      .           TCAR,DTOPTI,IPOTAB,KERRE)  
   4 : *  
   5 : ************************************************************
   6 : *  ENTREES  
   7 : ************************************************************  
   8 : *  
   9 : *  INPLAS  : numéro du matériau inélastique 
  10 : *  TSOM    :  
  11 : *  NSOM    :
  12 : *  NCOMP   :  
  13 : *  NINV    : 
  14 : *  NINCMA  :  
  15 : *  TCAR    : 
  16 : *  IPOTAB  : pointeur sur segment table 
  17 : *   
  18 : ************************************************************ 
  19 : *  SORTIES  
  20 : ************************************************************  
  21 : * 
  22 : *  DTOPTI  : pas de temps optimal ( modèles visqueux) 
  23 : *  KERRE   : indice d'erreur  
  24 : *
  25 : ************************************************************ 
  26 : *  
  27 :       IMPLICIT INTEGER(I-N)
  28 :       IMPLICIT REAL*8(A-H,O-Z)
  29 : *
  30 : -INC CCOPTIO
  31 : *  
  32 :       LOGICAL LOGIN,LOGRE 
  33 :       CHARACTER*72 CHARRE
  34 : *
  35 :       IF ( INPLAS .EQ. 17 .OR.
  36 :      2   ( INPLAS .GE. 19 .AND. INPLAS .LE. 25) .OR.
  37 :      4     INPLAS .EQ. 61 .OR.  
  38 :      1     INPLAS .EQ. 53 .OR.  
  39 :      1     INPLAS .EQ. 65 .OR.        
  40 :      1     INPLAS .EQ. 29 .OR.   
  41 :      2     INPLAS .EQ. 44 .OR.
  42 :      9     INPLAS .EQ. 45 .OR. 
  43 :      9     INPLAS .EQ. 76 .OR.  INPLAS .EQ. 77  .OR.
  44 :      9     INPLAS .EQ. 70 ) THEN
  45 : C
  46 : *    9     INPLAS .EQ. 70 .OR.
  47 : *    9     INPLAS .EQ. 74     ) THEN
  48 :          TMOY = TSOM/NCOMP
  49 :          NMOY = NSOM/NCOMP
  50 :          IF (IIMPI.GE.1) THEN
  51 :          WRITE(IOIMP,*)' NBRE DE SS PAS : MOYENNE ',NMOY,
  52 :      1                    ' MAX ',NINCMA
  53 :             IF (NINV.NE.0) THEN
  54 :             WRITE(IOIMP,*)' NBRE D"INVERSION ',NINV,'/',NCOMP
  55 :             END IF
  56 :          ENDIF
  57 :          TECAR = SQRT(ABS(TCAR/NCOMP - TMOY*TMOY))
  58 :          DTOPTI = MAX(TMOY-TECAR/2.D0,DTOPTI)
  58 :          CALL ECCTAB(IPOTAB,'MOT     ',IVALIN,XVALIN,'DTOPTI',
  59 :      1              LOGIN,IOBIN,'FLOTTANT',IVALRE,DTOPTI,CHARRE,
  60 :      2                                               LOGRE,IOBRE)
  61 :          IF(IERR.NE.0) THEN
  62 :             KERRE=999
  63 :          END IF
  64 :       ENDIF 
  65 :       RETURN 
  66 :       END  
  67 :  
  68 :  
  69 :  
  70 :  
  71 :  
  72 :  
  73 :  
  74 :