C DECO      SOURCE    AM        08/02/14    21:20:05     6050           
      SUBROUTINE DECO                                                           
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)                                                  
*--------------------------------------------------------------------*          
*                                                                    *          
*                     Opérateur DEnsite de COurant                   * 
*                     ____________________________                   *          
*                                                                    *          
*     Calcul de la densite de courant de Foucault                    *          
*                                                                    *          
* Syntaxe n›1:  DNS1 = DECO MOD1 DEP1 ( CAR1 ) ;                    *          
*                                                                    *          
*     Entr{es:                                                       *          
*                                                                    *          
*        MOD1      objet MMODEL                                      *          
*        DEP1      CHPOINT (FONCTION DE COURANT)                     *          
*        CAR1      MCHAML de sous-type CARACTERISTIQUES              *          
*                                                                    *          
*     Sortie:                                                        *          
*                                                                    *          
*        DNS1      MCHAML de sous-type GRADIENT                      *          
*                                                                    *          
*                                                                    *          
*      Auteurs, date de cr{ation:                                    *          
*                                                                    *          
*                                                                    *          
*        Y. STEPHAN le 22/09/97 (copie de GRAD)                      *
*--------------------------------------------------------------------*          
-INC CCOPTIO                                                                    
-INC CCHAMP                                                                     
*
*     NINFOS est le nombre d'informations contenues dans INFELE                 
*     ( VOIR LE S-P ELQUOI )                                                    
*                                                                               
      PARAMETER ( NINFOS=15 )                                                   
*                                                                                                                                               
* MC 19/01/98 : Fonction indisponible en défo planes généralisées
      IF(IFOUR.EQ.-3) THEN
         CALL ERREUR(710)
         RETURN
      ENDIF
*
      IPCHL1=0
      IPCHE1=0
      IRET1=0 
      IRET2=0
      IRET3=0                                                                     
      ICONV=0                                                                   
*                                                                               
*                                                                               
*     LECTURE D'UN MMODEL                                                       
*                                                                               
      CALL LIROBJ('MMODEL',IPMODL,1,IRET)                                       
      IF (IERR.NE.0) RETURN                                                     
*                                                                               
      IF (IRET.NE.0) THEN                                                      
*                                                                               
*       1-}RE SYNTAXE                                                           
*       _____________                                                           
*                                                                               
*                                                                               
*       LECTURE D'UN CHPOINT                                                    
*                                                                               
        CALL LIROBJ('CHPOINT',IPCHP1,0,IRET1)                                    
        IF(IERR.NE.0) RETURN
        IF (IRET1.NE.0) THEN
*                                                                               
*    CHPOINT  ---> CHAMELEM AUX NOEUDS                               
*                                                                                
           CALL CHAME1(0,IPMODL,IPCHP1,' ',IPCHE2,1)
           IF (IERR.NE.0) RETURN
        ENDIF
*                                                                               
*     LECTURE D'UN  MCHAML (CARACTERISTIQUE)                                    
*                                                                               
        CALL LIROBJ('MCHAML',IPCHE1,0,IRET3)                                
        IF (IERR.NE.0) RETURN
*                                                                               
      ENDIF                                                                     
*                                                                               
*     GESTION DES OP{RANDES                                                     
*     IPCHL1 = 0 (MCHAML calcule)
*                                
      CALL DECO1(IPMODL,IPCHE2,IPCHE1,IPCHL1,IRET)                              
*                                                                               
      IF (IRET.EQ.1) THEN                                                       
        IF (ICONV.EQ.0) THEN
          IF (IRET1.NE.0) CALL DTCHAM(IPCHE2)
          CALL ECROBJ('MCHAML',IPCHL1)                                          
        ELSE                                                                    
          CALL POPEL3(IPAFF1,1,IPCHL1,IPCHEL2)                                  
          IF (IPCHE1.NE.0) CALL DTCHAM(IPCHE1)                                  
          CALL DTCHAM(IPCHL1)                                                   
          CALL DTMODL(IPMODL)                                                   
          CALL ECROBJ('CHAMELEM',IPCHEL2)                                       
        ENDIF                                                                   
      ENDIF                                                                     
      RETURN                                                                    
      END