C MOTAN     SOURCE    JK148537  24/11/05    21:15:07     12067          
      SUBROUTINE MOTAN
*_______________________________________________________________________
*
*           Operateur permettant de calculer le module tangent
*
*   SYNTAXE :
*   ---------
*
*    ETAN=MOTAN MOD1 SI1 VA1 ( CA1 ) (PRECIS) ;
*
*        PRECIS  Flottant (par defaut 1.D-3)
*        MOD1    Modele de calcul, type MMODEL
*        SI1     Champ par element de contraintes,type MCHAML
*        VA1     Champ par element de variables internes,type MCHAML
*        CA1     Champ par element de caracteristiques (materielle et/ou
*                geometriques),type MCHAML
*        ETAN    Champ par element de type scalaire contenant les
*                modules tangent aux points des elements,type MCHAML
*
*    CODE L.EBERSOLT NOVEMBRE 86
*
*    passage aux nouveaux CHAMELEMs par JM CAMPENON LE 06/91
*_______________________________________________________________________
*
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
*

-INC PPARAM
-INC CCOPTIO
*
*     Lecture d'un modele
*
      CALL LIROBJ('MMODEL  ',IPMODL,1,IRETOU)
      CALL ACTOBJ('MMODEL  ',IPMODL,1)
      IF (IERR.NE.0) RETURN
*
*     Lecture du MCHAML de CONTRAINTES
*
      CALL LIROBJ('MCHAML  ',IPIN,1,IRT1)
      CALL ACTOBJ('MCHAML  ',IPIN,1)
      IF (IERR.NE.0) RETURN
      CALL REDUAF(IPIN,IPMODL,IPCHE1,0,IR,KER)
      IF(IR   .NE. 1) CALL ERREUR(KER)
      IF(IERR .NE. 0) RETURN
*
*     Lecture du MCHAML de VARIABLES INTERNES
*
      CALL LIROBJ('MCHAML',IPIN,1,IRT2)
      CALL ACTOBJ('MCHAML',IPIN,1)
      IF (IERR.NE.0) RETURN
      CALL REDUAF(IPIN,IPMODL,IPCHE2,0,IR,KER)
      IF(IR   .NE. 1) CALL ERREUR(KER)
      IF(IERR .NE. 0) RETURN
*
*     Lecture du MCHAML de CARACTERISTIQUES
*
      CALL LIROBJ('MCHAML',IPIN,1,IRT3)
      CALL ACTOBJ('MCHAML',IPIN,1)
      IF (IERR.NE.0) RETURN
      CALL REDUAF(IPIN,IPMODL,IPCHE3,0,IR,KER)
      IF(IR   .NE. 1) CALL ERREUR(KER)
      IF(IERR .NE. 0) RETURN
*
*     Lecture eventuelle de la precision des iteration internes
*
      CALL LIRREE(PRECIS,0,IRETPR)
      IF (IERR.NE.0) RETURN
      IF (IRETPR.EQ.0) PRECIS = 0.9
*
      CALL MOTANA(IPMODL,IPCHE1,IPCHE2,IPCHE3,PRECIS,IPSCAL,IRET)
*
*     Ecriture de la matrice de raideur tangente
*
      IF (IRET.EQ.1) THEN
        CALL ACTOBJ('MCHAML  ',IPSCAL,1)
        CALL ECROBJ('MCHAML  ',IPSCAL)
      ENDIF

      END