Esope source files

   1 : C CALJDC    SOURCE    CHAT      05/01/12    21:46:49     5004           
   2 :       SUBROUTINE CALJDC(FN,GR,PG,XYZ,HR,PGSQ,RPG,ND,IES,NP,NPG,
   3 :      *IAXI,AIRE)                                                                
   4 : C************************************************************************       
   5 : C                                                                               
   6 : C     CALJDC 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     ND   DIMENSION ESPACE DE DEFINITION DE L'ELEMENT (ALIAS NES)              
  17 : C     IES  DIMENSION ESPACE DE TRAVAIL >= ND                                    
  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 : C                                                                               
  28 :       REAL*8 FN(NP,NPG),GR(ND,NP,NPG),HR(IES,NP,NPG),HHR(162)
  29 :       REAL*8 PG(NPG),XYZ(IES,NP),PGSQ(NPG),RPG(NPG)                             
  30 :       REAL*8 AJ(3,3,64)
  31 : C                                                                               
  32 : C                                                                               
  33 : C     WRITE(6,*)' *** SUB CALJBC *** ND,IES',ND,IES                             
  34 :       IF(ND.EQ.IES)GOTO 20                                                      
  35 :                                                                                 
  35 :       CALL CALJCC(FN,GR,PG,XYZ,HR,PGSQ,RPG,IES,ND,NP,NPG,IAXI,AIRE,AJ,          
  36 :      &HHR)                                                                      
  37 :       RETURN                                                                    
  38 :  20   CONTINUE                                                                  
  39 :                                                                                 
  39 :       CALL CALJ22(FN,GR,PG,XYZ,HR,PGSQ,RPG,ND,NP,NPG,IAXI,AIRE,AJ)              
  40 :                                                                                 
  41 :       RETURN                                                                    
  42 :  1002 FORMAT(10(1X,1PE11.4))                                                    
  43 :  1001 FORMAT(20(1X,I5))                                                         
  44 :       END                                                                       
  45 :  
  46 :  
  47 :