Esope source files

   1 : C KSOMT     SOURCE    GOUNAND   06/04/27    21:15:06     5430           
   2 :       SUBROUTINE KSOMT
   3 :       IMPLICIT INTEGER(I-N)
   4 :       IMPLICIT REAL*8 (A-H,O-Z)
   5 : C     
   6 : C      CALCUL DE LA SOMME DES TERMES D UN TABLEAU
   7 : C      COMPOSANTE PAR COMPOSANTE
   8 : C     
   9 : -INC CCOPTIO
  10 : -INC SMCHPOI
  11 :       CHARACTER*8 TYPE
  12 : 
  12 :       CALL LIROBJ('CHPOINT ',MCHPOI,1,IRETOU)
  13 :       IF(IRETOU.EQ.0) THEN
  14 : C   37 2
  15 : C On ne trouve pas d'objet de type %m1:8
  16 :          MOTERR(1:8)='CHPOINT '
  16 :          CALL ERREUR(37)
  17 :          RETURN
  18 :       ENDIF
  18 :       CALL LICHT(MCHPOI,MPOVAL,TYPE,IGEOM)
  19 :       IF (MPOVAL.EQ.0) THEN
  20 :          U=0.D0 
  20 :          CALL ECRREE(U)
  21 :          RETURN
  22 :       ENDIF
  23 :       NC=VPOCHA(/2)
  24 :       IF(NC.NE.1)THEN
  25 : C  180 2
  26 : C Il faut spécifier un champ par point avec une seule composante
  26 :          CALL ERREUR(180)
  27 :          RETURN
  28 :       ENDIF
  29 :       NPT=VPOCHA(/1)
  30 :       U=0.D0
  31 :       DO 1 I=1,NPT
  32 :          U=U+VPOCHA(I,1)
  33 :  1    CONTINUE
  34 :       SEGDES MCHPOI,MPOVAL
  34 :       CALL ECRREE(U)
  35 :       RETURN
  36 :       END
  37 :  
  38 :