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