Esope source files

   1 : C DFOURI    SOURCE    FANDEUR   10/12/14    21:15:52     6812           
   2 :       SUBROUTINE DFOURI                                                         
   3 : C=======================================================================        
   4 : C                                                                               
   5 : C          OPERATEUR DFOURI                                                     
   6 : C                                                                               
   7 : C     SYNTAXE   I CHPT I  =  DFOURI  I CHPT I  ANGL                             
   8 : C               I CHEL I             I CHEL I                                   
   9 : C                                                                               
  10 : C     CHPT    : CHPOINT DE TYPE FORCES OU DEPLACEMENTS                          
  11 : C     CHEL    : CHAMP/ELEMENT DE CONTRAINTES OU DEFORMATIONS (MCHAML)           
  12 : C     ANGL    : VALEUR DE L ANGLE                                               
  13 : C                                                                               
  14 : C   DANS LE CAS D UNE ANALYSE EN SERIE DE FOURIER , CALCULE LES VALEURS         
  15 : C   DU CHPOINT OU DU MCHALM POUR L ANGLE DONNE                                  
  16 : C                                                                               
  17 : C     J BROCHARD MARS 87                                                        
  18 : C     PASSAGE AUX NOUVEAUX MCHAMLS  P DOWLATYARI SEPT. 91                       
  19 : C=======================================================================        
  20 :       IMPLICIT INTEGER(I-N)
  21 :       IMPLICIT REAL*8(A-H,O-Z)                                                  
  22 : -INC CCOPTIO  
  23 : -INC SMCHPOI                                                                  
  24 : C                                                                               
  25 : C ********  ON LIT L ANGLE                                                      
  26 : C                                                                               
  26 :       CALL LIRREE(ANGL,1,IRETOU)                                                
  27 :       IF (IERR.NE.0) RETURN                                                     
  28 : C                                                                               
  28 :       CALL LIROBJ('CHPOINT',IPCHP1,0,IRT1)                                      
  29 :       IF(IERR.NE.0)RETURN  
  30 :                                                      
  31 :       IF(IRT1.NE.0) THEN                                                        
  32 : C                                                                               
  33 : C    C'EST BIEN UN CHPOINT QU'IL FAUT LIRE                                      
  34 : C                                                                               
  34 :         CALL DFOUR1(IPCHP1,ANGL,IPCHP2)                                          
  35 :         IF (IERR.NE.0) RETURN
  36 : *        ATTRIBUTION D'UNE NATURE IDENTIQUE AU CHPO ARGUMENT
  37 :         MCHPOI = IPCHP1
  38 :         MCHPO1 = IPCHP2
  39 :         SEGACT ,MCHPOI,MCHPO1
  40 :         NAT = MAX ( 1 , MCHPO1.JATTRI(/1) )
  41 :         NSOUPO=MCHPO1.IPCHP(/1)
  42 :         SEGADJ MCHPO1
  43 :         IF ( JATTRI(/1) .GE. 1 )  THEN
  44 :              MCHPO1.JATTRI(1) = JATTRI(1)
  45 :         ELSE
  46 :              JATTRI(1) = 0
  47 :         ENDIF
  48 :         SEGDES ,MCHPOI,MCHPO1
  48 :         CALL ECROBJ('CHPOINT',IPCHP2)                               
  49 :       ELSE                                                                      
  50 : C                                                                               
  51 : C     LE CHAMP EST UN CHAMP/ELEMENT                                             
  52 : C                                                                               
  52 :          CALL LIROBJ('MCHAML',IPCHE1,1,IRT3)                                    
  53 :          IF(IERR.NE.0)RETURN                                                    
  53 :       CALL DFOUR2(IPCHE1,ANGL,IPCHE2)                                           
  53 :          IF (IERR.EQ.0) CALL ECROBJ('MCHAML',IPCHE2)                            
  54 :       ENDIF                                                                     
  55 :       RETURN                                                                    
  56 :       END                                                                       
  57 :  
  58 :