C CERBSG SOURCE CHAT 05/01/12 21:52:17 5004 SUBROUTINE CERBSG(FORCE,COORD,EFFORT,KERRE) C----------------------------------------------------------------------- C CALCUL DES FORCES INTERNES B*SIGMA C POUR LA CERCE C C ENTREE C EFFORT L'EFFORT DANS LA CERCE C COORD COORDONNEES DU NOEUD C C SORTIE C FORCE FORCES INTERNES C KERRE INDICE D'ERREUR ( 0 SI PAS DE PB ) C C----------------------------------------------------------------------- IMPLICIT INTEGER(I-N) IMPLICIT REAL*8(A-H,O-Z) -INC CCREEL -INC CCOPTIO DIMENSION COORD(3,*),FORCE(*) * KERRE=0 IF(IFOUR.EQ.0.OR.(IFOUR.EQ.1.AND.NIFOUR.EQ.0)) THEN FAC=2*XPI*EFFORT ELSE FAC=XPI*EFFORT ENDIF IF(IFOUR.EQ.0) THEN FORCE(1)=FAC FORCE(2)=0.D0 ELSE IF(IFOUR.EQ.1) THEN AN=DBLE(NIFOUR) FORCE(1)=FAC FORCE(2)=0.D0 FORCE(3)=FAC*AN ELSE KERRE=2 RETURN ENDIF RETURN END