C MOCUR     SOURCE    PV090527  24/09/03    21:15:02     11998          
      SUBROUTINE MOCUR
**********************************************************************
*
*     OPERATEUR MOCU  (MOment/CUrvature)
*
**********************************************************************
*     Pierre Pegon (ISPRA) Juillet/Aout 1993
**********************************************************************
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
C
      DIMENSION DEPSI(6),SIGMA(6),CRIGI(12),CMASS(12)
      logical zveri
C

-INC PPARAM
-INC CCOPTIO
-INC SMLREEL
-INC CCGEOME
-INC CCREEL
-INC SMMODEL
-INC SMCHAML
      POINTEUR IPCURY.MLREEL,IPCURZ.MLREEL,IPAXIA.MLREEL
      POINTEUR IPMOMY.MLREEL,IPMOMZ.MLREEL,IPEPAX.MLREEL
      CHARACTER*5 MOVERI(1)
      DATA MOVERI/'VERIF'/
C
      NMAXIT=50
      NSTRS2 = 6
      zveri=.false.
C
      ICOUL1=IDCOUL
C
C
C------- READ THE INPUT --
      IF ((IFOUR.EQ.-2).OR.(IFOUR.EQ.-1)) THEN
C
        CALL LIROBJ ('LISTREEL',IPCURZ,1,IRET)
        IF(IRET.EQ.0) RETURN
C
        CALL LIROBJ ('LISTREEL',IPAXIA,1,IRET0)
        IF(IRET.EQ.0) RETURN
C
        CALL LIROBJ ('MMODEL  ',IPMODL,1,IRET)
        IF(IRET.EQ.0) RETURN
        CALL ACTOBJ('MMODEL',IPMODL,1)
C
        CALL LIROBJ ('MCHAML  ',IPCAR,1,IRET)
        IF(IRET.EQ.0) RETURN
        CALL ACTOBJ('MCHAML',IPCAR ,1)
C
        CALL LIRREE (TOL,1,IRET)
        IF(IRET.EQ.0) RETURN
C
        CALL LIRMOT(MOVERI,1,iret,0)
        if (iret.eq.1) zveri=.true.
C
        CALL MOCU2D(IPCURZ,IPAXIA,IPMODL,IPCAR,IPEPAX,
     $       IPMOMZ,TOL,mtable,zveri)
C
        if (zveri) Call ecrobj('TABLE', mtable)
        CALL ECROBJ('LISTREEL',IPEPAX)
        CALL ECROBJ('LISTREEL',IPMOMZ)

      ELSE
C
        CALL LIROBJ ('LISTREEL',IPCURY,1,IRET)
        IF(IRET.EQ.0) RETURN
C
        CALL LIROBJ ('LISTREEL',IPCURZ,1,IRET)
        IF(IRET.EQ.0) RETURN
C
        CALL LIROBJ ('LISTREEL',IPAXIA,1,IRET0)
        IF(IRET0.EQ.0) RETURN
C
        CALL LIROBJ ('MMODEL  ',IPMODL,1,IRET)
        IF(IRET.EQ.0) RETURN
        CALL ACTOBJ('MMODEL',IPMODL,1)
C
        CALL LIROBJ ('MCHAML  ',IPCAR,1,IRET)
        IF(IRET.EQ.0) RETURN
        CALL ACTOBJ('MCHAML',IPCAR ,1)
C
        CALL LIRREE (TOL,1,IRET)
        IF(IRET.EQ.0) RETURN
C
        CALL LIRMOT(MOVERI,1,iret,0)
        if (iret.eq.1) zveri=.true.

C
         CALL MOCU3D(IPCURY,IPCURZ,IPAXIA,IPMODL,IPCAR,IPEPAX,
     $ IPMOMY,IPMOMZ,TOL,mtable,zveri)
C
        if (zveri) Call ecrobj('TABLE', mtable)
        CALL ECROBJ('LISTREEL',IPEPAX)
        CALL ECROBJ('LISTREEL',IPMOMZ)
        CALL ECROBJ('LISTREEL',IPMOMY)

C
      ENDIF
C
      RETURN
      END