Esope source files

   1 : C CERBSG    SOURCE    CHAT      05/01/12    21:52:17     5004           
   2 :       SUBROUTINE CERBSG(FORCE,COORD,EFFORT,KERRE)
   3 : C-----------------------------------------------------------------------
   4 : C             CALCUL DES FORCES INTERNES B*SIGMA
   5 : C                    POUR LA CERCE
   6 : C     
   7 : C  ENTREE
   8 : C          EFFORT  L'EFFORT DANS LA CERCE
   9 : C          COORD   COORDONNEES DU NOEUD
  10 : C     
  11 : C  SORTIE
  12 : C          FORCE   FORCES INTERNES
  13 : C          KERRE   INDICE D'ERREUR ( 0 SI PAS DE PB )
  14 : C     
  15 : C-----------------------------------------------------------------------
  16 :       IMPLICIT INTEGER(I-N)
  17 :       IMPLICIT REAL*8(A-H,O-Z)
  18 : -INC CCREEL
  19 : -INC CCOPTIO
  20 :       DIMENSION COORD(3,*),FORCE(*)
  21 : *     
  22 :       KERRE=0
  23 :       IF(IFOUR.EQ.0.OR.(IFOUR.EQ.1.AND.NIFOUR.EQ.0)) THEN
  24 :          FAC=2*XPI*EFFORT
  25 :       ELSE
  26 :          FAC=XPI*EFFORT
  27 :       ENDIF
  28 :       IF(IFOUR.EQ.0) THEN
  29 :          FORCE(1)=FAC
  30 :          FORCE(2)=0.D0
  31 :       ELSE IF(IFOUR.EQ.1) THEN
  32 :          AN=DBLE(NIFOUR)
  33 :          FORCE(1)=FAC
  34 :          FORCE(2)=0.D0
  35 :          FORCE(3)=FAC*AN
  36 :       ELSE
  37 :          KERRE=2
  38 :          RETURN
  39 :       ENDIF
  40 :       RETURN
  41 :       END
  42 : 
  43 : 
  44 : 
  45 :  
  46 :  
  47 :