Esope source files

   1 : C JACO      SOURCE    PV        09/03/12    21:25:41     6325           
   2 :       SUBROUTINE JACO                                                           
   3 : C=====================================================================          
   4 : C    OPERATEUR POUR CALCULER LE JACOBIEN AUX POINTS DE GAUSS                    
   5 : C                                                                               
   6 : C    CHAM = JACO  MODL                                                          
   7 : C              MODL : MODELE de calcul , type MMODEL                            
   8 : C              CHAM : CHAMELEM de sous type SCALAIRE contenant                  
   9 : C                       les JACOBIENS aux POINTS DE GAUSS                       
  10 : C                                                                               
  11 : C=====================================================================          
  12 : C                                                                               
  13 :       IMPLICIT INTEGER(I-N)
  14 :       IMPLICIT REAL*8(A-H,O-Z)                                                  
  15 :       PARAMETER(XZER=0.D0,UNDEMI=.5D0,UN=1.D0)                                  
  16 : -INC CCOPTIO                                                                    
  17 : -INC SMCHAML                                                                    
  18 : C                                                                               
  19 : C                                                                               
  20 : C                                                                               
  21 : C                                                                               
  22 : C           IPMODL          MODELE MMODEL                                       
  23 : C           IPCHE           CHAMELEM  DE SOUS TYPE SCALAIRE                     
  24 : C                           CONTENANT LE JACOBIEN                               
  25 : C                                                                               
  26 : C                                                                               
  27 :       IPMODL=0                                                                  
  28 :       IPCHE=0                                                                   
  29 :       IRET=0                                                                    
  30 : C                                                                               
  31 : C     LECTURE DU MODELE                                                         
  32 : C                                                                               
  32 :       CALL LIROBJ('MMODEL',IPMODL,1,IRTM)
  33 :       IF (IERR.NE.0) RETURN                                                     
  34 : C                                                                               
  35 : C     CALCUL DU JACOBIEN
  36 : C                                                                               
  36 :       CALL JACOPO(IPMODL,IPCHE,IRET)                                            
  37 : C
  37 :       IF(IRET.NE.0) CALL ECROBJ('MCHAML',IPCHE)
  38 :       RETURN                                                                    
  39 :       END                                                                       
  40 :  
  41 :  
  42 :  
  43 :