C DECO      SOURCE    PV        20/09/12    21:15:07     10711          
      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 PPARAM
-INC CCOPTIO
-INC CCHAMP
-INC SMCOORD
*
*     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
      segact mcoord
*
      IPCHL1=0
      IPCHE1=0
      IRET1 =0
      IRET3 =0
*
*
*     LECTURE D'UN MMODEL
*
      CALL LIROBJ('MMODEL  ',IPMODL,1,IRET)
      CALL ACTOBJ('MMODEL  ',IPMODL,1)
      IF (IERR.NE.0) RETURN
*
      IF (IRET.NE.0) THEN
*
*       1-ERE SYNTAXE
*       _____________
*
*
*       LECTURE D'UN CHPOINT
*
        CALL LIROBJ('CHPOINT ',IPCHP1,0,IRET1)
        IF(IERR.NE.0) RETURN
        IF (IRET1.NE.0) THEN
        CALL ACTOBJ('CHPOINT ',IPCHP1,1)
*
*    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  ',IPIN,0,IRET3)
        IF (IERR.NE.0) RETURN
        IPCHE1 = 0
        IF (IRET3 .EQ. 1) THEN
          CALL ACTOBJ('MCHAML  ',IPIN,1)
          CALL REDUAF(IPIN,IPMODL,IPCHE1,0,IR,KER)
          IF(IR   .NE. 1) CALL ERREUR(KER)
          IF(IERR .NE. 0) RETURN
        ENDIF
*
      ENDIF
*
*     GESTION DES OP{RANDES
*     IPCHL1 = 0 (MCHAML calcule)
*
      CALL DECO1(IPMODL,IPCHE2,IPCHE1,IPCHL1,IRET)
*
      IF (IRET.EQ.1) THEN
        CALL ACTOBJ('MCHAML  ',IPCHL1,1)
        CALL ECROBJ('MCHAML  ',IPCHL1)
      ENDIF
      END