Esope source files

   1 : C TAILLE    SOURCE    FANDEUR   10/08/30    21:15:48     6729           
   2 : c     mr millard      93    20:00:16     354
   3 : 
   4 :       SUBROUTINE TAILLE
   5 : 
   6 : c=====================================================================
   7 : c    operateur calculant les parametres de taille par element
   8 : c        aux points de gauss (tenseur de taille)
   9 : c
  10 : c    ancienne syntaxe :
  11 : c    ------------------
  12 : c
  13 : c    chel = tail modl
  14 : c            ob1 : modele de calcul, type mmodel
  15 : c            chel: chamelem de sous type scalaires contenant
  16 : c                    les parametres de taille aux points de gauss
  17 : c
  18 : c=====================================================================
  19 : c
  20 :       IMPLICIT INTEGER(I-N)
  21 :       IMPLICIT REAL*8(A-H,O-Z)
  22 : 
  23 : -INC CCOPTIO
  24 : *
  25 :       CHARACTER*4 MUNIF(3)
  26 :       DATA MUNIF/'DIAM','DIRE','UNIF'/
  27 : c
  28 : c           ipmodl          pointeur sur un modele mmodel
  29 : c           ipche           pointeur sur un chamelem  de
  30 : c                           sous type scalaires contenant
  31 : c                           les parametres de taille
  32 : c
  33 :       IPMODL = 0
  34 :       IUNIF  = 0
  35 : c
  36 : c Lecture du mot cle
  37 : c
  37 :       CALL LIRMOT(MUNIF,3,IUNIF,0)
  38 :       IF (IERR.NE.0) RETURN
  39 : c
  40 : c Lecture du modele
  41 : c
  41 :       CALL LIROBJ('MMODEL',IPMODL,1,IRET)
  42 :       IF (IERR.NE.0) RETURN
  43 : C
  44 : C Option : 'DIAMETRE_MIN'
  45 :       IF (IUNIF.EQ.1) THEN
  46 :         IPCHA1 = 0
  47 :         IPCHA2 = 0
  48 :         IPCHA3 = 0
  49 :         IPCHA4 = 0
  49 :         CALL LIROBJ('MCHAML',IPCHA1,0,IRET)
  50 :         IF (IERR.NE.0) RETURN
  51 :         ICAS = 5
  51 :         CALL CFL1(IPMODL,IPCHA1,IPCHA2,IPCHA3,IPCHA4,ICAS)
  51 :         IF (IPCHA4 .NE. 0) CALL ECROBJ('MCHAML',IPCHA4)
  52 : 
  53 : C Option : 'DIRECTION' & Cas ou le mot cle UNIF est seul
  54 :       ELSE
  55 :         IF (IUNIF.EQ.2) THEN
  56 : C Lecture eventuelle du mot cle 'UNIF' apres le mot cle 'DIRECTION'
  57 :           IUNIF = 0
  57 :           CALL LIRMOT(MUNIF(3),1,IUNIF,0)
  58 :           IF (IERR.NE.0) RETURN
  59 :         ENDIF
  60 : c Calcul des parametres de taille par element a chaque point de Gauss
  61 : c Si IUNIF non nul 0, toutes les composantes sont nulles.
  62 :         IPCHE = 0
  63 :         IRET  = 0
  63 :         CALL TAILPO(IPMODL,IPCHE,IUNIF,IRET)
  63 :         IF (IRET.NE.0) CALL ECROBJ('MCHAML',IPCHE)
  64 :       ENDIF
  65 : 
  66 :       RETURN
  67 :       END
  68 : 
  69 :