Sources Esope | Cast3M

ncoq4
C NCOQ4     SOURCE    CHAT      06/03/29    21:27:51     5360
C
      SUBROUTINE NCOQ4(IGAU,XEL,SHPTOT,SHP,BGENE,DJAC,iarr)
C
C ====================================================================
C    CARICANEMTO MATRICE N PER COQ4 AD USO OPERATORE MASSE
C      GIANNUZZI 15/12/86
C ====================================================================
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
      PARAMETER (LRE=24,NST=6,NBNO=4)
      DIMENSION XEL(3,*),BGENE(NST,*),SHP(6,*),SHPTOT(6,NBNO,*)
      DATA  ZER/0.D0/
C
      iarr=0
      CALL GINT2
      CALL ZERO(BGENE,NST,LRE)
      DO 101 NP=1,NBNO
      SHP(1,NP)=SHPTOT(1,NP,IGAU)
      SHP(2,NP)=SHPTOT(2,NP,IGAU)
      SHP(3,NP)=SHPTOT(3,NP,IGAU)
101   CONTINUE
      CALL JACOBI(XEL,SHP,2,NBNO,DJAC)
      IF(DJAC.LE.ZER) iarr=1
      IF(DJAC.LE.ZER) RETURN
C
C     CARICAMENTO MATRICE  - N -
C
      K=0
      DO 102 NP=1,NBNO
      DO 103 I =1,NST
      BGENE(I,K+I)=SHP(1,NP)
103   CONTINUE
      K=K+6
102   CONTINUE
      RETURN
      END