Esope source files

   1 : C NCOQ4     SOURCE    CHAT      06/03/29    21:27:51     5360           
   2 : C                                                                               
   3 :       SUBROUTINE NCOQ4(IGAU,XEL,SHPTOT,SHP,BGENE,DJAC,iarr)                     
   4 : C                                                                               
   5 : C ====================================================================          
   6 : C    CARICANEMTO MATRICE N PER COQ4 AD USO OPERATORE MASSE                      
   7 : C      GIANNUZZI 15/12/86                                                       
   8 : C ====================================================================          
   9 :       IMPLICIT INTEGER(I-N)
  10 :       IMPLICIT REAL*8(A-H,O-Z)                                                  
  11 :       PARAMETER (LRE=24,NST=6,NBNO=4)                                           
  12 :       DIMENSION XEL(3,*),BGENE(NST,*),SHP(6,*),SHPTOT(6,NBNO,*)                 
  13 :       DATA  ZER/0.D0/                                                           
  14 : C                                                                               
  15 :       iarr=0                                                                    
  15 :       CALL GINT2                                                                
  15 :       CALL ZERO(BGENE,NST,LRE)                                                  
  16 :       DO 101 NP=1,NBNO                                                          
  17 :       SHP(1,NP)=SHPTOT(1,NP,IGAU)                                               
  18 :       SHP(2,NP)=SHPTOT(2,NP,IGAU)                                               
  19 :       SHP(3,NP)=SHPTOT(3,NP,IGAU)                                               
  20 : 101   CONTINUE                                                                  
  20 :       CALL JACOBI(XEL,SHP,2,NBNO,DJAC)                                          
  21 :       IF(DJAC.LE.ZER) iarr=1                                                    
  22 :       IF(DJAC.LE.ZER) RETURN                                                    
  23 : C                                                                               
  24 : C     CARICAMENTO MATRICE  - N -                                                
  25 : C                                                                               
  26 :       K=0                                                                       
  27 :       DO 102 NP=1,NBNO                                                          
  28 :       DO 103 I =1,NST                                                           
  29 :       BGENE(I,K+I)=SHP(1,NP)                                                    
  30 : 103   CONTINUE                                                                  
  31 :       K=K+6                                                                     
  32 : 102   CONTINUE                                                                  
  33 :       RETURN                                                                    
  34 :       END                                                                       
  35 :  
  36 :