Esope source files

   1 : C CALJBC    SOURCE    CHAT      05/01/12    21:46:40     5004           
   2 :       SUBROUTINE CALJBC(FN,GR,PG,XYZ,HR,PGSQ,RPG,NES,IES,NP,NPG,
   3 :      *IAXI,AIRE)                                                                
   4 : C************************************************************************       
   5 : C                                                                               
   6 : C     CALJBC APPEL CALJ22                                                       
   7 : C                                                                               
   8 : C            LA DIFFERENCE AVEC CALJAC EST LE RANGEMENT DE HR                   
   9 : C                        HR(IES,NP)                                             
  10 : C                                                                               
  11 : C                                                                               
  12 : C     CALCUL DE L'INVERSE DU JACOBIEN  AJ=1/J                                   
  13 : C     CALCUL DE L'ELEMENT D'AIRE       SQ=DET(J)                                
  14 : C     DANS LES CAS 1D 2D ET 3D                                                  
  15 : C                                                                               
  16 : C     NES  DIMENSION ESPACE DE DEFINITION DE L'ELEMENT
  17 : C     IES  DIMENSION ESPACE
  18 : C     NP   NOMBRE DE NOEUDS DE L'ELEMENT                                        
  19 : C     NPG  NOMBRE DE POINTS D'INTEGRATION                                       
  20 : C                                                                               
  21 : C     XYZ  COORDONNEES                                                          
  22 : C     GR   GRADIENT                                                             
  23 : C                                                                               
  24 : C************************************************************************       
  25 :       IMPLICIT INTEGER(I-N)
  26 :       IMPLICIT REAL*8 (A-H,O-Z)                                                 
  27 : -INC CCREEL
  28 : C                                                                               
  29 :       REAL*8 FN(NP,NPG),GR(NES,NP,NPG),HR(IES,NP,NPG)
  30 :       REAL*8 PG(NPG),XYZ(IES,NP),PGSQ(NPG),RPG(NPG)                             
  31 :       REAL*8 AJ(3,3,64)
  32 : C                                                                               
  32 :       CALL CALJBR
  33 :      &(FN,GR,PG,XYZ,HR,PGSQ,RPG,NES,IES,NP,NPG,IAXI,AIRE,AJ,SGN)
  34 :       IF(IAXI.NE.0)THEN
  35 :       DEUPI=2.D0*XPI
  36 :       DO 1 L=1,NPG
  37 :       PGSQ(L)=PGSQ(L)*DEUPI*RPG(L)
  38 :  1    CONTINUE
  39 :       ENDIF
  40 :                                                                                 
  41 :       RETURN                                                                    
  42 :  1002 FORMAT(10(1X,1PE11.4))                                                    
  43 :  1001 FORMAT(20(1X,I5))                                                         
  44 :       END                                                                       
  45 :  
  46 :  
  47 :  
  48 :  
  49 :