Esope source files

   1 : C POLA      SOURCE    CHAT      06/08/24    21:58:20     5529           
   2 :       SUBROUTINE POLA
   3 : *---------------------------------------------------------------------          
   4 : *                                                                               
   5 : *            DECOMPOSITION POLAIRE D'UN CHAMP DE GRADIENT
   6 : *                          F = R U
   7 : *                                                                               
   8 : *     CHR CHU  = POLA  MOD CHEL (CAR) ( MOT );
   9 : *                                                                               
  10 : *           MOD   MODELE DE CALCUL, TYPE MMODEL                                 
  11 : *           CHEL  OBJET DE TYPE MCHAML DE SOUS TYPE GRADIENT
  12 : *           CHR   OBJET DE TYPE MCHAML DE SOUS TYPE
  13 : *           CHU   OBJET DE TYPE MCHAML DE SOUS TYPE
  14 : *                                                                               
  15 : *--------------------------------------------------------------------           
  16 : *                                                                               
  17 :       IMPLICIT INTEGER(I-N)
  18 : -INC CCOPTIO                                                                    
  19 :       CHARACTER*4 MOCLE(2)
  20 :       DATA MOCLE/'GEOM','DEPL'/
  21 : *                                                                               
  21 :       CALL LIROBJ('MMODEL',IPMODL,1,IRTM)
  22 :       IF(IERR.NE.0) RETURN
  23 : *                                                                               
  23 :       CALL LIROBJ('MCHAML',IPCHE1,1,IRT1)
  24 :       IF(IERR.NE.0) RETURN
  25 : *                                                                               
  26 : *                                                                               
  27 : *     LECTURE D'UN MOT SI NECESSAIRE                                            
  28 : *                                                                               
  29 :       IMIL=0                                                                    
  29 :       CALL LIRMOT(MOCLE,2,IMOT,0)
  30 :       IF(IMOT.EQ.2) IMIL=1
  31 : *                                                                               
  31 :       CALL POLA1(IPMODL,IPCHE1,IPCHE2,IPCHE3,IMIL)
  32 :       IF(IERR.NE.0) RETURN
  33 : *                                                                               
  33 :       CALL ECROBJ('MCHAML',IPCHE3)
  33 :       CALL ECROBJ('MCHAML',IPCHE2)
  34 :       RETURN                                                                    
  35 :       END                                                                       
  36 :  
  37 :  
  38 :