C MASSE     SOURCE    PV        09/03/26    21:15:30     6346           
      SUBROUTINE MASSE(ILUMP)
*
*_______________________________________________________________________
*
*            operateurs masse et lump
*
*
*           ILUMP : =1 indique si il s'agit de l'opérateur LUMP
*
*
*_______________________________________________________________________
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
*
-INC CCOPTIO
-INC SMCHAML
-INC CCHAMP
c
      CHARACTER*4  CMOT
      
C     Lump d une matrice deja constituée ( formule initiale )
      CALL LIROBJ('RIGIDITE',IRIG,0,IR1)                                        
      IF (IR1.NE.0) THEN 
      LMOT=0                                                                    
      CALL LIROBJ('LISTMOTS',LMOT,0,IR2)                                        
      IF (IERR.NE.0) RETURN       
      if (ir2.eq.0.and.ilump.eq.0) then 
       call ecrobj('RIGIDITE',IRIG)
       goto 10
      endif
                   
*                                                                               
      CALL LUMPIN(IRIG,LMOT,ILUM)                                               
      IF (IERR.NE.0) RETURN                                                     
*                                                                               
      CALL ECROBJ('RIGIDITE',ILUM)                                              
      RETURN                             
      END IF        
c
c     lecture du modele
c     
 10   IPCHE1=0
c
      CALL LIROBJ('MMODEL',IPMODL,0,IRT1)
c
      IF(IRT1.NE.0) THEN
c
c        Cas de la matrice masse sur base element finis
c        lecture du 1 ier mchaml
c
          CALL LIROBJ('MCHAML',IPCHE1,1,IRT2)
c
c         si absent erreur
c
          IF (IRT2.EQ.0) THEN 
               MOTERR(1:16) ='CARACTERISTIQUE' 
               CALL ERREUR(565)
               RETURN
          ENDIF
      call reduaf(ipche1,ipmodl,ipch,0,ir,ker)
      if (ir.ne.1) call erreur(ker)
      IF (IERR.NE.0) RETURN
      ipche1=ipch
c
c     calcul de la masse
c

          CALL MASSE1 (IPMODL,IPCHE1,IPMASS,IRET,ILUMP)
c
c     ecriture du resultat
c
          CALL ECROBJ('RIGIDITE',IPMASS)

      ELSE 
        IF (ILUMP .EQ. 0) THEN
c 
c   cas des masses additionnelles
c
          CALL LIRCHA(CMOT,0,IRT1)
          IF (IRT1.NE.0) THEN
            CALL REFUS
            CALL APPUI(2)
          ELSE 
c
c   cas de la masse sur base modale
c
            CALL MASSMO
          ENDIF
        ELSE
          CALL ERREUR(26)
        ENDIF
      ENDIF
c
c
c
      RETURN
      END