Esope source files

   1 : C GRAD      SOURCE    AM        08/02/14    21:23:43     6050           
   2 : 
   3 : C=======================================================================
   4 : C=                              G R A D                                =
   5 : C=                              -------                                =
   6 : C=                                                                     =
   7 : C=  OPERATEUR CAST3M "GRAD" DE CALCUL DU GRADIENT D'UN CHAMP :         =
   8 : C=  ----------------------------------------------------------         =
   9 : C=   GRA1  =  'GRAD'  MODL1  CHP1 ou CHE1  ( CAR1 )  ;                 =
  10 : C=                                                                     =
  11 : C=  Cet operateur calcule le gradient d'un champ de temperatures ou    =
  12 : C=  de deplacements.                                                   =
  13 : C=                                                                     =
  14 : C=  ARGUMENTS :                                                        =
  15 : C=  -----------                                                        =
  16 : C=   MODL1   (MMODEL)   Modele associe a la structure etudiee          =
  17 : C=   CHP1    (CHPOINT)  | Champ de deplacements           (NOEUDS)     =
  18 : C=   CHE1    (MCHAML)   | ou de temperatures              (ELEMENT)    =
  19 : C=   CAR1    (MCHAML)   Caracteristiques geometriques     (facultatif) =
  20 : C=                      Sous-type 'CARACTERISTIQUES'                   =
  21 : C=                                                                     =
  22 : C=  RESULTAT :                                                         =
  23 : C=  ----------                                                         =
  24 : C=   GRA1    (MCHAML)   Champ de gradients sur la structure            =
  25 : C=                                                                     =
  26 : C=  SUO X.Z., le 16 novemebre 1986.                                    =
  27 : C=  Christian LE BRETON et Denis ROBERT-MOUGIN, le 31 juillet 1989.    =
  28 : C=  Modifications aux nouvelles normes par I. MONNIER, le 28 mai 1990. =
  29 : C=======================================================================
  30 : 
  31 :       SUBROUTINE GRAD
  32 : 
  33 :       IMPLICIT INTEGER(I-N)
  34 :       IMPLICIT REAL*8 (A-H,O-Z)
  35 : 
  36 : -INC CCOPTIO
  37 : -INC CCHAMP
  38 : 
  39 : C  1 - LECTURE DES ARGUMENTS DE L'OPERATEUR
  40 : C ==========================================
  41 : C  1.1 - Lecture OBLIGATOIRE du modele (IPMODL)
  42 : C =====
  42 :       CALL LIROBJ('MMODEL',IPMODL,1,IRET)
  43 :       IF (IERR.NE.0) RETURN
  44 : C =====
  45 : C  1.2 - Syntaxe 1 : lecture d'un CHPOINT de depl. ou T (IPCHP1)
  46 : C        Transformation du CHPOINT en MCHAML aux noeuds (IPCHE2)
  47 : C =====
  48 :       IRET1=0
  48 :       CALL LIROBJ('CHPOINT',IPCHP1,0,IRET1)
  49 :       IF (IERR.NE.0) RETURN
  50 :       IRET2=0
  51 :       IF (IRET1.NE.0) THEN
  52 :         IPCHL1=IPCHP1
  52 :         CALL CHAME1(0,IPMODL,IPCHP1,' ',IPCHE2,1)
  53 :         IF (IERR.NE.0) RETURN
  54 : C =====
  55 : C  1.3 - Syntaxe 3 : lecture d'un MCHAML de depl. ou T (IPCHE2)
  56 : C        ERREUR si le mode de calcul est "generalise"
  57 : C =====
  58 :       ELSE
  59 :         IF (IFOUR.EQ.-3.OR.
  60 :      .      (IFOUR.GE.7.AND.IFOUR.LE.11).OR.IFOUR.EQ.14) THEN
  61 :           MOTERR(1:8)='CHPOINT '
  61 :           CALL ERREUR(37)
  62 :           RETURN
  63 :         ENDIF
  64 :         IPCHL1=0
  64 :         CALL LIROBJ('MCHAML',IPCHE2,1,IRET2)
  65 :         IF (IERR.NE.0) RETURN
  66 :       ENDIF
  67 : C =====
  68 : C  1.4 - Lecture FACULTATIVE des caracteristiques geometriques (IPCHE1)
  69 : C =====
  70 :       IRET3=0
  70 :       CALL LIROBJ('MCHAML',IPCHE1,0,IRET3)
  71 :       IF (IERR.NE.0) RETURN
  72 :       IF (IRET3.NE.0.AND.IRET2.NE.0) THEN
  72 :         CALL RNGCHA(IPCHE1,IPCHE2,'CARACTERISTIQUES',
  73 :      .              'DEPLACEMENTS',IPCHA1,IPCHA2)
  74 :         IF (IPCHA1.EQ.0.OR.IPCHA2.EQ.0) THEN
  74 :           CALL ERREUR(805)
  75 :           RETURN
  76 :         ENDIF
  77 :         IPCHE1=IPCHA1
  78 :         IPCHE2=IPCHA2
  79 :       ENDIF
  80 : C =====
  81 : C  1.5 - Lecture FACULTATIVE du mot 'CONS'
  82 : C =====
  83 :       ISCA=0
  83 :       CALL LIRMOT('CONS',1,ISCA,0)
  84 : 
  85 :       IF (ISCA.EQ.0) THEN
  85 :         CALL GRAD1(IPMODL,IPCHE2,IPCHE1,IPCHL1,IRET)
  86 :       ELSE
  86 :         CALL GRAD2(IPMODL,IPCHE2,IPCHE1,IPCHL1,IRET)
  87 :       ENDIF
  88 : 
  89 :       IF (IRET.EQ.1) THEN
  89 :         IF (IRET1.NE.0) CALL DTCHAM(IPCHE2)
  89 :         CALL ECROBJ('MCHAML',IPCHL1)
  90 :       ENDIF
  91 : 
  92 :       RETURN
  93 :       END
  94 : 
  95 :  
  96 :  
  97 :  
  98 :  
  99 :