C POLA      SOURCE    CHAT      06/08/24    21:58:20     5529           
      SUBROUTINE POLA
*---------------------------------------------------------------------          
*                                                                               
*            DECOMPOSITION POLAIRE D'UN CHAMP DE GRADIENT
*                          F = R U
*                                                                               
*     CHR CHU  = POLA  MOD CHEL (CAR) ( MOT );
*                                                                               
*           MOD   MODELE DE CALCUL, TYPE MMODEL                                 
*           CHEL  OBJET DE TYPE MCHAML DE SOUS TYPE GRADIENT
*           CHR   OBJET DE TYPE MCHAML DE SOUS TYPE
*           CHU   OBJET DE TYPE MCHAML DE SOUS TYPE
*                                                                               
*--------------------------------------------------------------------           
*                                                                               
      IMPLICIT INTEGER(I-N)
-INC CCOPTIO                                                                    
      CHARACTER*4 MOCLE(2)
      DATA MOCLE/'GEOM','DEPL'/
*                                                                               
      CALL LIROBJ('MMODEL',IPMODL,1,IRTM)
      IF(IERR.NE.0) RETURN
*                                                                               
      CALL LIROBJ('MCHAML',IPCHE1,1,IRT1)
      IF(IERR.NE.0) RETURN
*                                                                               
*                                                                               
*     LECTURE D'UN MOT SI NECESSAIRE                                            
*                                                                               
      IMIL=0                                                                    
      CALL LIRMOT(MOCLE,2,IMOT,0)
      IF(IMOT.EQ.2) IMIL=1
*                                                                               
      CALL POLA1(IPMODL,IPCHE1,IPCHE2,IPCHE3,IMIL)
      IF(IERR.NE.0) RETURN
*                                                                               
      CALL ECROBJ('MCHAML',IPCHE3)
      CALL ECROBJ('MCHAML',IPCHE2)
      RETURN                                                                    
      END