Télécharger issleq.procedur

Retour à la liste

* ISSLEQ    PROCEDUR  AF221230  17/03/28    21:15:00     9380           
DEBPROC ISSLEQ TDON1*'TABLE';
************************************************************************
*
*                            I S S L E Q
*                           ------------
*
* FONCTION:
* ---------
*
* Procedure pour la résolution des problèmes ISS et de propagation des 
*  ondes 1D avec la methode lineaire equivalente.
*
* PHRASE D'APPEL (EN GIBIANE):
* ----------------------------
*
*   ISSLEQ TDON1
*
* OPERANDES:
* ----------
*
*  TDON1               :Table des donnees (type TABLE)
*       .'SOL'             :sous-table pour definir les 
*                           donnees du profil 
*                           du sol (type TABLE)
*          .'ITERATION'       :sous table pour la partie du
*                              sol à iterer (type TABLE)
*                .I              : couche I (type TABLE)          
*                  .'MAILLAGE'        :maillage de la couche I
*                                      (type MAILLAGE)
*                  .'FRONTIERE'       :frontiere de la couche I 
*                                      (type MAILLAGE)
*                  .'MASSE_VOLUMIQUE' :masse volumique de la couche I
*                                      (type FLONNTANT)
*                  .'POISSON'         :coefficient de poisson de la 
*                                      couche I (type FLOTTANT)
*                  .'G_GAMMA'         :evolution de la courbe 
*                                      caracteristique G/G0-gamma 
*                                      (type EVOLUTION)
*                  .'H_GAMMA'         :evolution de la courbe 
*                                      caracteristique Eps-gamma 
*                                      (type EVOLUTION)
*                  .'BASE'            :frontiere inferieure de la 
*                                      derniere couche qui compose
*                                      la stratigraphie (type MAILLAGE)
*                                      a introduire une seule fois pour
*                                      la dernier couche
*          .'CONSTANTE'       :sous table pour la partie du
*                              sol à iterer (type TABLE)
*                .I              : couche I (type TABLE)    
*                  .'MAILLAGE'        :maillage de la couche I
*                                      (type MAILLAGE)
*                  .'FRONTIERE'       :frontiere de la couche I 
*                                      (type MAILLAGE)
*                  .'MASSE_VOLUMIQUE' :masse volumique de la couche I
*                                      (type FLONNTANT)
*                  .'POISSON'         :coefficient de poisson de la 
*                                      couche I (type FLOTTANT)
*                  .'AMORTISSEMENT'   :amortissement
*                                      couche I (type FLOTTANT)
*       .'STRUCTURE'       :sous-table pour definir les 
*                           donnees de la structure (type TABLE)
*          .'MAILLAGE'            :maillage de la structure 
*                                  (type MAILLAGE)
*          .'RIGIDITE'            :rigidité de la structure
*                                  (type RIGIDITE)
*          .'MASSE'               :masse de la structure
*                                  (type RIGIDITE)
*          .'AMORTISSEMENT'       :amortissement de la structure
*                                  (type RIGIDITE)
*          .'BLOCAGES_MECANIQUES' :tout type des blocages mecaniques
*                                  supplementaires (type RIGIDITE)
*       .'PARAMETRES'       :sous TABLE pour definir les parametres
*                            du calcul (type TABLE)
*          .'GAMMAO_X'      :evolution de l’acceleration selon 
*                            la direction x (type EVOLUTION)
*          .'GAMMAO_Y'      :evolution de l’acceleration selon
*                            la direction y (type EVOLUTION)
*          .'GAMMAO_Z'      :evolution de l’acceleration selon
*                            la direction z (type EVOLUTION)
*                            [seulement une des trois evolution 
*                                est necessaire]
*          .'POINT'         :point de reference (type POINT)
*          .'CRITERE'       :critere de convergence pour le 
*                            calcul lineaire equivalent (type FLOTTANT)
*                            (default 0.05)
*          .'CHI'           :coefficient pour la deformation moyenne 
*                            gamma_m (type FLOTTANT)
*          .'F1'            :premiere frequence pour modele 
*                            d'amortissement de Rayleigh 
*                            (type FLOTTANT)
*          .'F2'            :premiere frequence pour modele
*                            d'amortissement de Rayleigh 
*                            (type FLOTTANT)
*          .'FC'            :premiere de coupure du signal
*                            (type FLOTTANT)
*          .'TYPE'          :type de frontiere absorbante (LYSMER 
*                            ou WHITE) (type MOT) (defaut LYSMER)
*          .'PAR_DEC'       :sous-table à donner pour la procedure
*                            DECONV3D - voir la notice correspondante
*                            (P_GAM sous-table) (facultatif)
*          .'TYPE_CALCUL'   : type de calcul ‘ISS_COMPLET’ (type MOT)
*                                    tout le sol est à considere comme
*                                    domaine à iterer. La sous-table
*                                    ITERATION dans la table SOL est
*                                    obligatoire
*                             type de calcul ‘ISS_COMP_SIMP’ (type MOT)
*                                    une partie du sol est considere
*                                    comme domaine à iterer et une 
*                                    partie constante. La sous table
*                                    ITERATION dans la table SOL et la
*                                    sous table CONSTANTE dans la table
*                                    SOL sont obligatoires
* RESULTATS:
* ----------
*  TDON1               :Table des donnees (type TABLE)
*       .'STEP "J"'        :sous-table des resultats 
*                           pour chaque iteration (type TABLE) J=0..N
*          .'ITERATION'       :sous table pour la partie du
*                              sol à iterer (type TABLE)
*                .'T_SOL'        :sous-table avec les maillages des 
*                                 couche
*                   .I               : couche I (type MAILLAGE)
*                .'T_ESO'        :sous-table avec les mchaml pour 
*                                 les module E des couches
*                   .I               : couche I (type MCHAML)
*                .'T_MOD'        :sous-table avec les modeles pour 
*                                 les couches
*                   .I               : couche I (type MMODEL)
*                .'T_GGO'        :sous-table avec les mchaml pour 
*                                 les rapports G/G0 des couches
*                   .I               : couche I (type MCHAML)
*                .'T_MAT'        :sous-table avec les materiaux pour 
*                                 les couches
*                   .I               : couche I (type MCHAML)
*                .'T_RIG'        :sous-table avec les rigites pour 
*                                 les couches
*                   .I               : couche I (type RIGIDITE)
*                .'T_MAS'        :sous-table avec les masses pour 
*                                 les couches
*                   .I               : couche I (type RIGIDITE)
*                .'T_KSI'        :sous-table avec les amortissements
*                                 pour les couches
*                   .I               : couche I (type RIGIDITE)
*                .'T_AMO'        :sous-table avec les matrices 
*                                 d amortissement amortissements
*                                 pour les couches
*                   .I               : couche I (type MCHAML)
*                .'T_DIS'        :sous-table avec les champs des 
*                                 deformations maximales pour chaque
*                                 couche
*                   .I               : couche I (type MAILLAGE)
*                .'O_SOL'        :Assemblage des maillages pour
*                                 la partie à iterer (type MAILLAGE)
*                .'O_MOD'        :Assemblage des MMODEL pour
*                                 la partie à iterer (type MAILLAGE)
*                .'O_GGO'        :Assemblage des MCHAML G/G0 pour
*                                 la partie à iterer (type MAILLAGE)
*                .'O_MAT'        :Assemblage des materiaux pour
*                                 la partie à iterer (type MAILLAGE)
*                .'O_RIG'        :Assemblage des rigidités pour
*                                 la partie à iterer (type MAILLAGE)
*                .'O_MAS'        :Assemblage des masses pour
*                                 la partie à iterer (type MAILLAGE)
*                .'O_KSI'        :Assemblage des MCHAML Eps pour
*                                 la partie à iterer (type MAILLAGE)
*                .'O_AMO'        :Assemblage des amortissements pour
*                                 la partie à iterer (type MAILLAGE)
*                .'O_DIS'        :Assemblage des MCHAML Dmax pour
*                                 la partie à iterer (type MAILLAGE)
*          .'CONSTANTE'       :sous table pour la partie du
*                              sol constante (type TABLE)
*                .'T_SOL'        :sous-table avec les maillages des 
*                                 couche
*                   .I               : couche I (type MAILLAGE)
*                .'T_ESO'        :sous-table avec les mchaml pour 
*                                 les module E des couches
*                   .I               : couche I (type MCHAML)
*                .'T_MOD'        :sous-table avec les modeles pour 
*                                 les couches
*                   .I               : couche I (type MMODEL)
*                .'T_GGO'        :sous-table avec les mchaml pour 
*                                 les rapports G/G0 des couches
*                   .I               : couche I (type MCHAML)
*                .'T_MAT'        :sous-table avec les materiaux pour 
*                                 les couches
*                   .I               : couche I (type MCHAML)
*                .'T_RIG'        :sous-table avec les rigites pour 
*                                 les couches
*                   .I               : couche I (type RIGIDITE)
*                .'T_MAS'        :sous-table avec les masses pour 
*                                 les couches
*                   .I               : couche I (type RIGIDITE)
*                .'T_KSI'        :sous-table avec les amortissements
*                                 pour les couches
*                   .I               : couche I (type RIGIDITE)
*                .'T_AMO'        :sous-table avec les matrices 
*                                 d amortissement amortissements
*                                 pour les couches
*                   .I               : couche I (type MCHAML)
*                .'T_DIS'        :sous-table avec les champs des 
*                                 deformations maximales pour chaque
*                                 couche
*                   .I               : couche I (type MAILLAGE)
*                .'O_SOL'        :Assemblage des maillages pour
*                                 la partie à iterer (type MAILLAGE)
*                .'O_MOD'        :Assemblage des MMODEL pour
*                                 la partie à iterer (type MAILLAGE)
*                .'O_GGO'        :Assemblage des MCHAML G/G0 pour
*                                 la partie à iterer (type MAILLAGE)
*                .'O_MAT'        :Assemblage des materiaux pour
*                                 la partie à iterer (type MAILLAGE)
*                .'O_RIG'        :Assemblage des rigidités pour
*                                 la partie à iterer (type MAILLAGE)
*                .'O_MAS'        :Assemblage des masses pour
*                                 la partie à iterer (type MAILLAGE)
*                .'O_KSI'        :Assemblage des MCHAML Eps pour
*                                 la partie à iterer (type MAILLAGE)
*                .'O_AMO'        :Assemblage des amortissements pour
*                                 la partie à iterer (type MAILLAGE)
*                .'O_DIS'        :Assemblage des MCHAML Dmax pour
*                                 la partie à iterer (type MAILLAGE)
*          .'TOTALE'       :sous table pour l'assembles des deux parties
*                           CONSTANTE + ITERER (type TABLE)
*                   .I               : couche I (type MAILLAGE)
*                .'O_SOL'        :Assemblage des maillages pour
*                                 la partie à iterer (type MAILLAGE)
*                .'O_MOD'        :Assemblage des MMODEL pour
*                                 la partie à iterer (type MAILLAGE)
*                .'O_GGO'        :Assemblage des MCHAML G/G0 pour
*                                 la partie à iterer (type MAILLAGE)
*                .'O_MAT'        :Assemblage des materiaux pour
*                                 la partie à iterer (type MAILLAGE)
*                .'O_RIG'        :Assemblage des rigidités pour
*                                 la partie à iterer (type MAILLAGE)
*                .'O_MAS'        :Assemblage des masses pour
*                                 la partie à iterer (type MAILLAGE)
*                .'O_KSI'        :Assemblage des MCHAML Eps pour
*                                 la partie à iterer (type MAILLAGE)
*                .'O_AMO'        :Assemblage des amortissements pour
*                                 la partie à iterer (type MAILLAGE)
*                .'O_DIS'        :Assemblage des MCHAML Dmax pour
*                                 la partie à iterer (type MAILLAGE)
*       .'STRUCTURE'       :sous-table pour definir les 
*                           donnees de la structure (type TABLE)
*          .'MAILLAGE'            :maillage de la structure 
*                                  (type MAILLAGE)
*          .'RIGIDITE'            :rigidité de la structure
*                                  (type RIGIDITE)
*          .'MASSE'               :masse de la structure
*                                  (type RIGIDITE)
*          .'AMORTISSEMENT'       :amortissement de la structure
*                                  (type RIGIDITE)
*          .'BLOCAGES_MECANIQUES' :tout type des blocages mecaniques
*                                  supplementaires (type RIGIDITE)
*       .'RES_DYN_X'              :Table des resultats du calcul DYNAMIC
*
* COMMENTAIRES
* --------------------------------------------
* La procedure peut etre utilisée seulement pour des cas 3D
*
* AUTEUR, DATE DE CREATION ET DE MODIFICATION:
* --------------------------------------------
*
*     A. FRAU   25/11/2015 : Creation.
 
*
************************************************************************
*
*************************************************************
MESS '';
MESS '**************************************************************';
MESS '';
MESS '--------------------------------------------------------------';
MESS '---------------------- Procedure ISSLEQ ----------------------';
MESS '--   Procedure pour la realissation d un calcul ISS avec    --';
MESS '--               la methode lineaire equivalente            --';
MESS '--------------------------------------------------------------';
MESS '';
*************************************************************
*
'SI' (('VALE' 'DIME') 'NEG' 3);
  'MESS';
  'MESS' '  La procedure ISSLEQ marche seulement pour les calculs 3D';
  QUIT ISSLEQ;
'FINSI';
*
*************************************************************
***************     inizio procedura      *******************
*************************************************************
*
*
*
'SI' ('EXIST' (TDON1) 'STEP0');
*************************************************************
****************     STEP SUP A 1      **********************
*************************************************************
*
'SI' ('EXIST' (TDON1) 'COLONNE');
  'SI' ('EXIST' (TDON1) 'STRUCTURE');
     I_STEP = ('DIME' TDON1) - 5;
  'SINON';
     I_STEP = ('DIME' TDON1) - 4;
  'FINSI';
'SINON';
  'SI' ('EXIST' (TDON1) 'STRUCTURE');
     I_STEP = ('DIME' TDON1) - 4;
  'SINON';
     I_STEP = ('DIME' TDON1) - 3;
  'FINSI';
'FINSI';
*
MOT11 = 'CHAINE' 'STEP ' I_STEP;
MOT22 = 'CHAINE' 'STEP' I_STEP;
MOT33 = 'CHAINE' 'STEP' (I_STEP - 1);
*
MESS '++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++';
MESS '--                            'MOT11'                        --';
MESS '--------------------------------------------------------------';
*
TDON1.MOT22 = TABLE;
*
**************************************************************
*   extraction tables
**************************************************************
*
MESS '';
MESS '--------------------------------------------------------------';
MESS '--                 Lecture des donnees d entree             --';
MESS '--------------------------------------------------------------';
*
TABPAR = TDON1.'PARAMETRES';
TABFRO = TDON1.'FRONTIERE';
'SI' ('EXIST' (TDON1) 'STRUCTURE');
  TABSTR = TDON1.MOT33.'STRUCTURE';
'FINSI';
*
'SI' (EXIST (TDON1.'SOL') 'CONSTANTE');
  TABSOL_C = TDON1.'SOL'.'CONSTANTE';
'FINSI';
*
'SI' (EXIST (TDON1.'SOL') 'ITERATION');
  TABSOL_I = TDON1.'SOL'.'ITERATION';
'FINSI';
*
'SI' (EXIST (TDON1) 'COLONNE');
  TABCOL = TDON1.'COLONNE';
'FINSI';
*
'SI' (EXIST (TDON1.MOT33.'SOL') 'CONSTANTE');
  TABRES_C = TDON1.MOT33.'SOL'.'CONSTANTE';
*
  SOL_C = TABRES_C.T_SOL;
  ESO_C = TABRES_C.T_ESO;
  MOD_C = TABRES_C.T_MOD;
  MAT_C = TABRES_C.T_MAT;
  KSI_C = TABRES_C.T_KSI;
  RIG_C = TABRES_C.T_RIG;
  MAS_C = TABRES_C.T_MAS;
  AMO_C = TABRES_C.T_AMO;
  GGO_C = TABRES_C.T_GGO;
*
  MAI_SOLC = TABRES_C.O_SOL;
  MOD_SOLC = TABRES_C.O_MOD;
  MAT_SOLC = TABRES_C.O_MAT;
  GGO_SOLC = TABRES_C.O_GGO;
  KSI_SOLC = TABRES_C.O_KSI;
  RIG_SOLC = TABRES_C.O_RIG;
  MAS_SOLC = TABRES_C.O_MAS;
  AMO_SOLC = TABRES_C.O_AMO;
* 
  NZONE_C = DIME SOL_C;
*
'FINSI';
*
'SI' (EXIST (TDON1.MOT33.'SOL') 'ITERATION');
  TABRES_I = TDON1.MOT33.'SOL'.'ITERATION';
*
  SOL_I = TABRES_I.T_SOL;
  ESO_I = TABRES_I.T_ESO;
  MOD_I = TABRES_I.T_MOD;
  MAT_I = TABRES_I.T_MAT;
  KSI_I = TABRES_I.T_KSI;
  RIG_I = TABRES_I.T_RIG;
  MAS_I = TABRES_I.T_MAS;
  AMO_I = TABRES_I.T_AMO;
  GGO_I = TABRES_I.T_GGO;
*
  MAI_SOLI = TABRES_I.O_SOL;
  MOD_SOLI = TABRES_I.O_MOD;
  MAT_SOLI = TABRES_I.O_MAT;
  GGO_SOLI = TABRES_I.O_GGO;
  KSI_SOLI = TABRES_I.O_KSI;
  RIG_SOLI = TABRES_I.O_RIG;
  MAS_SOLI = TABRES_I.O_MAS;
  AMO_SOLI = TABRES_I.O_AMO;
*
  NZONE_I = DIME SOL_I;
*
'FINSI';
*
'SI' (EXIST (TDON1.MOT33) 'COLONNE');
  TABRESC = TDON1.MOT33.'COLONNE';
*
  SOLCOL = TABRESC.T_SOL;
  ESOCOL = TABRESC.T_ESO;
  MODCOL = TABRESC.T_MOD;
  MATCOL = TABRESC.T_MAT;
  KSICOL = TABRESC.T_KSI;
  RIGCOL = TABRESC.T_RIG;
  MASCOL = TABRESC.T_MAS;
  AMOCOL = TABRESC.T_AMO;
  GGOCOL = TABRESC.T_GGO;
*
  MAISOLC = TABRESC.O_SOL;
  MODSOLC = TABRESC.O_MOD;
  MATSOLC = TABRESC.O_MAT;
  GGOSOLC = TABRESC.O_GGO;
  KSISOLC = TABRESC.O_KSI;
  RIGSOLC = TABRESC.O_RIG;
  MASSOLC = TABRESC.O_MAS;
  AMOSOLC = TABRESC.O_AMO;
*
  NZONEC = DIME SOLCOL;
*
'FINSI';
*
F1     = TABPAR.'F1' ;
F2     = TABPAR.'F2' ;
FC     = TABPAR.'FC' ;
V_MOD = VALEUR MODE;
*
**************************************************************
*   preparation deconvolution
**************************************************************
*
MESS '';
MESS '--------------------------------------------------------------';
MESS '--                        Deconvolution                     --';
MESS '--------------------------------------------------------------';
'SI' (EGA (TYPE TABRESC) 'TABLE');
*
*  Deconvilution 3D (Cas Mode Plan et 3D) - Colonne
*  NN1 = 'DIME' (TABCOL);
  I = 1;
  'REPE' IND1 (NZONEC);
    'SI' (I 'EGA' 1);
      BOR_FRO = (TABCOL. I .'FRONTIERE');
    'SINON';
      BOR_FRO = BOR_FRO 'ET' (TABCOL. I .'FRONTIERE');
    'FINSI';
    I = I + 1;
  'FIN' IND1;
  'ELIM' 0.001 BOR_FRO;
  BOR_FON = TABCOL. NZONEC .'BASE';
  'OUBL' NN1;
  'OUBL' I;
*
  PP1 = BOR_FRO POIN PROC (10000. 10000. 10000.);
  PP2 = BOR_FRO POIN PROC (10000. 10000. -10000.);
  LL2 = BOR_FRO POIN DROIT PP1 PP2 0.001;
*
  LL2B = LL2;
  I = 1;
  NN1 = NBNO LL2;
  'REPE' IND1 (NN1);
    'SI' (I EGA NN1);
      PP1 = LL2B;
    'SINON';
      PP1 = LL2B POIN PROC (0. 0. 10000.);
    'FINSI';
    'SI' ('EGA' ('TYPE' LL2A) 'POINT');
      LL2A  = LL2A 'ET' PP1;
    'SINON';
      'SI' ('EGA' ('TYPE' LL2A) 'MAILLAGE');
        LL2A  = LL2A ET PP1;
      'SINON';
        LL2A  = PP1;
      'FINSI';
    'FINSI';
    NN2 = NN1 - (I - 1);
    J = 1;
    'SI' (NN2 'NEG' 1);
      'REPE' IND2 (NN2);
        PPX = LL2B 'POIN' J;
        VAL1 = 0.;
        VAL1 = VAL1 + (((PPX 'COOR' 1) -(PP1 'COOR' 1))**(2));
        VAL1 = VAL1 + (((PPX 'COOR' 2) -(PP1 'COOR' 2))**(2));
        VAL1 = VAL1 + (((PPX 'COOR' 3) -(PP1 'COOR' 3))**(2));
        VAL1 = ((VAL1)**(0.5));
        'SI' (VAL1 > 0.001);
          'SI' ('EGA' ('TYPE' LL2C) 'POINT');
            LL2C = LL2C 'ET' PPX;
          'SINON';
            'SI' ('EGA' ('TYPE' LL2C) 'MAILLAGE');
              LL2C = LL2C 'ET' PPX;
            'SINON';
              LL2C = PPX;
            'FINSI';
          'FINSI';
        'FINSI';
        J = J + 1;
      'FIN' IND2;
    'FINSI';
    LL2B = LL2C;
    'OUBL' LL2C;
    I = I + 1;
  'FIN' IND1;
  LL2 = LL2A;
  'OUBL' LL2A;
  'OUBL' LL2B;
  'OUBL' LL2C;
  'OUBL' NN1;
  'OUBL' NN2;
  'OUBL' VAL1;
  'OUBL' PP1;
  'OUBL' PPX;
*
  NELEM_F = ('NBNO' LL2) - 1;
  I = 1 ;
  ID2 = 'CHAINE' 'N';
  TABDV2 = 'TABLE';
  'REPE' BLOC4 NELEM_F ;
    PP1 = LL2 'POIN' I;
    PP2 = LL2 'POIN' (I + 1);
    XX1 = -10000.;
    XX2 = 10000.;
    YY1 = -10000.;
    YY2 = 10000.;
    ZZ1 = (PP1 'COOR' 3) + 0.01;
    ZZ2 = (PP2 'COOR' 3) - 0.01;
    CH1 CH2 CH3 = 'COOR'  (BOR_FRO);
    MAI2 = CH1 'POIN' 'COMPRIS' XX1 XX2;
    CH2  = 'REDU' CH2 MAI2 ;
    MAI2 = CH2 'POIN' 'COMPRIS' YY1 YY2;
    CH3  = 'REDU' CH3 MAI2 ;
    MAI2 = CH3 'POIN' 'COMPRIS' ZZ2 ZZ1;
    CH3  = 'REDU' CH3 MAI2 ;
    MAI2 = (BOR_FRO) 'ELEM' 'APPUYE' 
                     'STRICTEMENT' MAI2;
    TABDV2.I = 'TABLE' ;
    TABDV2.I.'FRONTIERE' =  MAI2;
    Z1    = (PP1 COOR 3);
    Z3    = (PP2 COOR 3);
    TABDV2.I.'EPAISSEUR' = Z1 - Z3 ;
*     ZELEM1 = MAISOLC 'ELEM' 'CONTENANT' PP1;
    ZELEM2 = MAISOLC 'ELEM' 'CONTENANT' PP2;
    PPX1 = ZELEM1 POIN PROC PP1;
    PPX2 = ZELEM1 POIN PROC PP2;
    PPY1 = ZELEM2 POIN PROC PP1;
    PPY2 = ZELEM2 POIN PROC PP2;
*     VALX1 = 0.;
    VALX2 = 0.;
    VALX1 = VALX1 + (((PPX1 COOR 1) - (PP1 COOR 1))**(2.));
    VALX1 = VALX1 + (((PPX1 COOR 2) - (PP1 COOR 2))**(2.));
    VALX1 = VALX1 + (((PPX1 COOR 3) - (PP1 COOR 3))**(2.));
    VALX1 = ((VALX1)**(0.5));
    VALX2 = VALX2 + (((PPX2 COOR 1) - (PP2 COOR 1))**(2.));
    VALX2 = VALX2 + (((PPX2 COOR 2) - (PP2 COOR 2))**(2.));
    VALX2 = VALX2 + (((PPX2 COOR 3) - (PP2 COOR 3))**(2.));
    VALX2 = ((VALX2)**(0.5));
    VALY1 = 0.;
    VALY2 = 0.;
    VALY1 = VALY1 + (((PPY1 COOR 1) - (PP1 COOR 1))**(2.));
    VALY1 = VALY1 + (((PPY1 COOR 2) - (PP1 COOR 2))**(2.));
    VALY1 = VALY1 + (((PPY1 COOR 3) - (PP1 COOR 3))**(2.));
    VALY1 = ((VALY1)**(0.5));
    VALY2 = VALY2 + (((PPY2 COOR 1) - (PP2 COOR 1))**(2.));
    VALY2 = VALY2 + (((PPY2 COOR 2) - (PP2 COOR 2))**(2.));
    VALY2 = VALY2 + (((PPY2 COOR 3) - (PP2 COOR 3))**(2.));
    VALY2 = ((VALY2)**(0.5));
    'SI' ((VALX1 + VALX2) < 0.0001);
      ZELEM = ZELEM1;
    'FINSI';
    'SI' ((VALY1 + VALY2) < 0.0001);
      ZELEM = ZELEM2;
    'FINSI';
    ZELEM = ZELEM1 'ELEM' 'CONTENANT' PP1;
    'OUBL' ZELEM1;
    'OUBL' ZELEM2;
    'OUBL' VALX1;
    'OUBL' VALX2;
    'OUBL' VALY1;
    'OUBL' VALY2;
*
    RMATE = REDU MATSOLC ZELEM ;
    RKSI  = REDU KSISOLC ZELEM ;
*
    TABDV2.I.'MASSE_VOLUMIQUE' = EXTR RMATE RHO  1 1 1 ;
    TABDV2.I.'POISSON'         = EXTR RMATE NU   1 1 1 ;
    TABDV2.I.'YOUNG'           = EXTR RMATE YOUN 1 1 1 ;
    TABDV2.I.'AMORTISSEMENT'   = EXTR RKSI  KS   1 1 1 ;
    I = I + 1 ;
  'FIN' BLOC4 ;
  'OUBL' ELEM_I;
  'OUBL' NOEU1;
  'OUBL' NOEU2;
  'OUBL' NOEU3;
  'OUBL' Z1;
  'OUBL' Z3;
  'OUBL' PP1;
  'OUBL' PP2;
  'OUBL' ID1;
*
  TDON1.MOT22.'COLONNE' = TABLE;
  'SI' ('EXIST' TABPAR 'GAMMAO_X');
    DIR    = 'UX'          ;
    GAMMAO = TABPAR.'GAMMAO_X' ;
    TYP_F = TABPAR.'TYPE' ;
    MOD_DEC = MODSOLC;
    MAT_DEC = MATSOLC;
    MAI_DEC = MAISOLC;
*
    PPX1 = 
        (TABDV2 . 1 .FRONTIERE) POIN PROC (-100000. -100000. 10000.);
    PPX2 =
        (TABDV2 . 1 .FRONTIERE) POIN PROC (100000. 100000. 10000.);
 
    PPXF = (((PPX2 COOR 1) + (PPX1 COOR 1))/(2.))
           (((PPX2 COOR 2) - (PPX1 COOR 2))/(2.))
           (((PPX2 COOR 3) - (PPX1 COOR 3))/(2.));
    P_GAM = TABLE;
    P_GAM. 1 = (TABDV2. 1 .'FRONTIERE') 'POIN' 'PROC' (10000. 10000.
                                                             10000.);
    TABRXDV2 =  DECONV3D TABDV2 (PPXF) BOR_FON MOD_DEC 
                                   MAT_DEC DIR GAMMAO F1 F2 FC P_GAM;
    TDON1.MOT22.'COLONNE'.'TAB_DEC' = TABDV2;
    TDON1.MOT22.'COLONNE'.'RES_DEC_X' = TABRXDV2;
    'OUBL' P_GAM;
    'OUBL' PPXF
    'OUBL' BOR_FRO;
    'OUBL' BOR_FON;
    'OUBL' MOD_DEC;
    'OUBL' MAT_DEC;
    'OUBL' DIR;
    'OUBL' GAMMAO;
    'OUBL' MAI_DEC;
  'FINSI';
*
  'SI' ('EXIST' TABPAR 'GAMMAO_Y');
    DIR    = 'UY'          ;
    GAMMAO = TABPAR.'GAMMAO_Y' ;
    TYP_F = TABPAR.'TYPE' ;
    MOD_DEC = MODSOLC;
    MAT_DEC = MATSOLC;
    MAI_DEC = MAISOLC;
*
    PPX1 = (TABDV2 . 1 .FRONTIERE) 
                               POIN PROC (-100000. -100000. 10000.);
    PPX2 = (TABDV2 . 1 .FRONTIERE)
                                 POIN PROC (100000. 100000. 10000.);
 
*
    PPXF = (((PPX2 COOR 1) + (PPX1 COOR 1))/(2.))
           (((PPX2 COOR 2) - (PPX1 COOR 2))/(2.))
           (((PPX2 COOR 3) - (PPX1 COOR 3))/(2.));
    P_GAM = TABLE;
    P_GAM. 1 = (TABDV2. 1 .'FRONTIERE') 'POIN' 'PROC' (10000. 10000.
                                                             10000.);
    TABRYDV2 =  DECONV3D TABDV2 (PPXF) BOR_FON MOD_DEC 
                                 MAT_DEC DIR GAMMAO F1 F2 FC P_GAM;
    TDON1.MOT22.'COLONNE'.'TAB_DEC' = TABDV2;
    TDON1.MOT22.'COLONNE'.'RES_DEC_Y' = TABRYDV2;
    'OUBL' P_GAM;
    'OUBL' PPXF
    'OUBL' BOR_FRO;
    'OUBL' BOR_FON;
    'OUBL' MOD_DEC;
    'OUBL' MAT_DEC;
    'OUBL' DIR;
    'OUBL' GAMMAO;
    'OUBL' MAI_DEC;
  'FINSI';
*
  'SI' ('EXIST' TABPAR 'GAMMAO_Z');
    DIR    = 'UZ'          ;
    GAMMAO = TABPAR.'GAMMAO_Z' ;
    TYP_F = TABPAR.'TYPE' ;
    MOD_DEC = MODSOLC;
    MAT_DEC = MATSOLC;
    MAI_DEC = MAISOLC;
*
    PPX1 = (TABDV2 . 1 .FRONTIERE) 
                               POIN PROC (-100000. -100000. 10000.);
    PPX2 = (TABDV2 . 1 .FRONTIERE) 
                                 POIN PROC (100000. 100000. 10000.);
 
*
    PPXF = (((PPX2 COOR 1) + (PPX1 COOR 1))/(2.))
           (((PPX2 COOR 2) - (PPX1 COOR 2))/(2.))
           (((PPX2 COOR 3) - (PPX1 COOR 3))/(2.));
    P_GAM = TABLE;
    P_GAM. 1 = (TABDV2. 1 .'FRONTIERE') 'POIN' 'PROC' (10000. 10000.
                                                             10000.);
    TABRZDV2 =  DECONV3D TABDV2 (PPXF) BOR_FON MOD_DEC 
                                   MAT_DEC DIR GAMMAO F1 F2 FC P_GAM;
    TDON1.MOT22.'COLONNE'.'TAB_DEC' = TABDV2;
    TDON1.MOT22.'COLONNE'.'RES_DEC_Z' = TABRZDV2;
    'OUBL' P_GAM;
    'OUBL' PPXF
    'OUBL' BOR_FRO;
    'OUBL' BOR_FON;
    'OUBL' MOD_DEC;
    'OUBL' MAT_DEC;
    'OUBL' DIR;
    'OUBL' GAMMAO;
    'OUBL' MAI_DEC;
  'FINSI';
'SINON';
*
*  Deconvilution 3D (Cas 3D)
  NN1 = 'DIME' (TABFRO.LATERALE);
  I = 1;
  'REPE' IND1 (NN1);
    'SI' (I EGA 1);
      BOR_FRO = (TABFRO.LATERALE. I .MAILLAGE);
    'SINON';
      BOR_FRO = BOR_FRO ET (TABFRO.LATERALE. I .MAILLAGE);
    'FINSI';
    I = I + 1;
  'FIN' IND1;
  'ELIM' 0.001 BOR_FRO;
  BOR_FON = TABFRO.'BASE';
  'OUBL' NN1;
  'OUBL' I;
*
  PP1 = BOR_FRO POIN PROC (10000. 10000. 10000.);
  PP2 = BOR_FRO POIN PROC (10000. 10000. -10000.);
  LL2 = BOR_FRO POIN DROIT PP1 PP2 0.001;
*
  LL2B = LL2;
  I = 1;
  NN1 = NBNO LL2;
  'REPE' IND1 (NN1);
    'SI' (I EGA NN1);
      PP1 = LL2B;
    'SINON';
      PP1 = LL2B POIN PROC (0. 0. 10000.);
    'FINSI';
    'SI' ('EGA' ('TYPE' LL2A) 'POINT');
      LL2A  = LL2A 'ET' PP1;
    'SINON';
      'SI' ('EGA' ('TYPE' LL2A) 'MAILLAGE');
        LL2A  = LL2A ET PP1;
      'SINON';
        LL2A  = PP1;
      'FINSI';
    'FINSI';
    NN2 = NN1 - (I - 1);
    J = 1;
    'SI' (NN2 'NEG' 1);
      'REPE' IND2 (NN2);
        PPX = LL2B 'POIN' J;
        VAL1 = 0.;
        VAL1 = VAL1 + (((PPX 'COOR' 1) -(PP1 'COOR' 1))**(2));
        VAL1 = VAL1 + (((PPX 'COOR' 2) -(PP1 'COOR' 2))**(2));
        VAL1 = VAL1 + (((PPX 'COOR' 3) -(PP1 'COOR' 3))**(2));
        VAL1 = ((VAL1)**(0.5));
        'SI' (VAL1 > 0.001);
          'SI' ('EGA' ('TYPE' LL2C) 'POINT');
            LL2C = LL2C 'ET' PPX;
          'SINON';
            'SI' ('EGA' ('TYPE' LL2C) 'MAILLAGE');
              LL2C = LL2C 'ET' PPX;
            'SINON';
              LL2C = PPX;
            'FINSI';
          'FINSI';
        'FINSI';
        J = J + 1;
      'FIN' IND2;
    'FINSI';
    LL2B = LL2C;
    'OUBL' LL2C;
    I = I + 1;
  'FIN' IND1;
  LL2 = LL2A;
  'OUBL' LL2A;
  'OUBL' LL2B;
  'OUBL' LL2C;
  'OUBL' NN1;
  'OUBL' NN2;
  'OUBL' VAL1;
  'OUBL' PP1;
  'OUBL' PPX;
*
  NELEM_F = ('NBNO' LL2) - 1;
  I = 1 ;
  ID2 = 'CHAINE' 'N';
  TABDV = 'TABLE';
*
  'REPE' BLOC4 NELEM_F ;
    PP1 = LL2 POIN I;
    PP2 = LL2 POIN (I + 1);
    XX1 = -10000.;
    XX2 = 10000.;
    YY1 = -10000.;
    YY2 = 10000.;
    ZZ1 = (PP1 COOR 3) + 0.01;
    ZZ2 = (PP2 COOR 3) - 0.01;
    CH1 CH2 CH3 = 'COOR'  (BOR_FRO);
    MAI2 = CH1 'POIN' 'COMPRIS' XX1 XX2;
    CH2  = 'REDU' CH2 MAI2 ;
    MAI2 = CH2 'POIN' 'COMPRIS' YY1 YY2;
    CH3  = 'REDU' CH3 MAI2 ;
    MAI2 = CH3 'POIN' 'COMPRIS' ZZ2 ZZ1;
    CH3  = 'REDU' CH3 MAI2 ;
    MAI2 = (BOR_FRO) 'ELEM' 'APPUYE' 
                     'STRICTEMENT' MAI2;
*
    TABDV.I = 'TABLE' ;
    TABDV.I.'FRONTIERE' =  MAI2;
    Z1    = (PP1 COOR 3);
    Z3    = (PP2 COOR 3);
    TABDV.I.'EPAISSEUR' = Z1 - Z3 ;
*
    'SI' ('EGA' ('TYPE' (MAI_SOLC)) 'MAILLAGE');
      PPP1 = MAI_SOLC 'POIN' 'PROC' PP1;
      VAL1 = 0.;
      VAL1 = VAL1 + (((PPP1 COOR 1) - (PP1 COOR 1))**(2.));
      VAL1 = VAL1 + (((PPP1 COOR 2) - (PP1 COOR 2))**(2.));
      VAL1 = VAL1 + (((PPP1 COOR 3) - (PP1 COOR 3))**(2.));
      VAL1 = ((VAL1)**(0.5));
    'SINON';
      VAL1 = 100000.;
    'FINSI';
    'SI' ('EGA' ('TYPE' (MAI_SOLI)) 'MAILLAGE');
      PPP2 = MAI_SOLI 'POIN' 'PROC' PP1;
      VAL2 = 0.;
      VAL2 = VAL2 + (((PPP2 COOR 1) - (PP1 COOR 1))**(2.));
      VAL2 = VAL2 + (((PPP2 COOR 2) - (PP1 COOR 2))**(2.));
      VAL2 = VAL2 + (((PPP2 COOR 3) - (PP1 COOR 3))**(2.));
      VAL2 = ((VAL2)**(0.5));
    'SINON';
      VAL2 = 100000.;
    'FINSI';
*
*
    'SI' (VAL1 < 0.0001);
      ZELEM1 = MAI_SOLC 'ELEM' 'CONTENANT' PP1;
      ZELEM2 = MAI_SOLC 'ELEM' 'CONTENANT' PP2;
      PPX1 = ZELEM1 POIN PROC PP1;
      PPX2 = ZELEM1 POIN PROC PP2;
      PPY1 = ZELEM2 POIN PROC PP1;
      PPY2 = ZELEM2 POIN PROC PP2;
*
      VALX1 = 0.;
      VALX2 = 0.;
      VALX1 = VALX1 + (((PPX1 COOR 1) - (PP1 COOR 1))**(2.));
      VALX1 = VALX1 + (((PPX1 COOR 2) - (PP1 COOR 2))**(2.));
      VALX1 = VALX1 + (((PPX1 COOR 3) - (PP1 COOR 3))**(2.));
      VALX1 = ((VALX1)**(0.5));
      VALX2 = VALX2 + (((PPX2 COOR 1) - (PP2 COOR 1))**(2.));
      VALX2 = VALX2 + (((PPX2 COOR 2) - (PP2 COOR 2))**(2.));
      VALX2 = VALX2 + (((PPX2 COOR 3) - (PP2 COOR 3))**(2.));
      VALX2 = ((VALX2)**(0.5));
      VALY1 = 0.;
      VALY2 = 0.;
      VALY1 = VALY1 + (((PPY1 COOR 1) - (PP1 COOR 1))**(2.));
      VALY1 = VALY1 + (((PPY1 COOR 2) - (PP1 COOR 2))**(2.));
      VALY1 = VALY1 + (((PPY1 COOR 3) - (PP1 COOR 3))**(2.));
      VALY1 = ((VALY1)**(0.5));
      VALY2 = VALY2 + (((PPY2 COOR 1) - (PP2 COOR 1))**(2.));
      VALY2 = VALY2 + (((PPY2 COOR 2) - (PP2 COOR 2))**(2.));
      VALY2 = VALY2 + (((PPY2 COOR 3) - (PP2 COOR 3))**(2.));
      VALY2 = ((VALY2)**(0.5));
      'SI' ((VALX1 + VALX2) < 0.0001);
        ZELEM = ZELEM1;
      'FINSI';
      'SI' ((VALY1 + VALY2) < 0.0001);
        ZELEM = ZELEM2;
      'FINSI';
*
      'OUBL' ZELEM1;
      'OUBL' ZELEM2;
      'OUBL' VALX1;
      'OUBL' VALX2;
      'OUBL' VALY1;
      'OUBL' VALY2;
      ID1 = 'CHAINE' 'C';
      'SI' ('EGA' ID2 'C');
        ID2 = 'CHAINE' 'C';
      'SINON';
        'SI' ('EGA' ID2 'I');
          ID2 = 'M';
        'SINON';
          ID2 = 'C';
        'FINSI';
      'FINSI';
    'SINON';
      'SI' (VAL2 < 0.0001);
        ZELEM1 = MAI_SOLI 'ELEM' 'CONTENANT' PP1;
        ZELEM2 = MAI_SOLI 'ELEM' 'CONTENANT' PP2;
        PPX1 = ZELEM1 POIN PROC PP1;
        PPX2 = ZELEM1 POIN PROC PP2;
        PPY1 = ZELEM2 POIN PROC PP1;
        PPY2 = ZELEM2 POIN PROC PP2;
*
        VALX1 = 0.;
        VALX2 = 0.;
        VALX1 = VALX1 + (((PPX1 COOR 1) - (PP1 COOR 1))**(2.));
        VALX1 = VALX1 + (((PPX1 COOR 2) - (PP1 COOR 2))**(2.));
        VALX1 = VALX1 + (((PPX1 COOR 3) - (PP1 COOR 3))**(2.));
        VALX1 = ((VALX1)**(0.5));
        VALX2 = VALX2 + (((PPX2 COOR 1) - (PP2 COOR 1))**(2.));
        VALX2 = VALX2 + (((PPX2 COOR 2) - (PP2 COOR 2))**(2.));
        VALX2 = VALX2 + (((PPX2 COOR 3) - (PP2 COOR 3))**(2.));
        VALX2 = ((VALX2)**(0.5));
        VALY1 = 0.;
        VALY2 = 0.;
        VALY1 = VALY1 + (((PPY1 COOR 1) - (PP1 COOR 1))**(2.));
        VALY1 = VALY1 + (((PPY1 COOR 2) - (PP1 COOR 2))**(2.));
        VALY1 = VALY1 + (((PPY1 COOR 3) - (PP1 COOR 3))**(2.));
        VALY1 = ((VALY1)**(0.5));
        VALY2 = VALY2 + (((PPY2 COOR 1) - (PP2 COOR 1))**(2.));
        VALY2 = VALY2 + (((PPY2 COOR 2) - (PP2 COOR 2))**(2.));
        VALY2 = VALY2 + (((PPY2 COOR 3) - (PP2 COOR 3))**(2.));
        VALY2 = ((VALY2)**(0.5));
        'SI' ((VALX1 + VALX2) < 0.0001);
          ZELEM = ZELEM1;
        'FINSI';
        'SI' ((VALY1 + VALY2) < 0.0001);
          ZELEM = ZELEM2;
        'FINSI';
*
        ID1 = 'CHAINE' 'I';
        'SI' ('EGA' ID2 'I');
          ID2 = 'CHAINE' 'I';
        'SINON';
          'SI' ('EGA' ID2 'C');
            ID2 = 'M';
          'SINON';
            ID2 = 'I';
          'FINSI';
        'FINSI';
      'SINON';
        'MESS' 'Error contruction table DECONV';
        'QUIT' ISSLEQ;
      'FINSI';
    'FINSI';
    'OUBL' ZELEM1;
    'OUBL' ZELEM2;
*
    'SI' ('EGA' ID1 'C');
      RMATE = REDU MAT_SOLC ZELEM ;
      RKSI  = REDU KSI_SOLC ZELEM ;
    'FINSI';
    'SI' ('EGA' ID1 'I');
      RMATE = REDU MAT_SOLI ZELEM ;
      RKSI  = REDU KSI_SOLI ZELEM ;
    'FINSI';
*
    TABDV.I.'MASSE_VOLUMIQUE' = EXTR RMATE RHO  1 1 1 ;
    TABDV.I.'POISSON'         = EXTR RMATE NU   1 1 1 ;
    TABDV.I.'YOUNG'           = EXTR RMATE YOUN 1 1 1 ;
    TABDV.I.'AMORTISSEMENT'   = EXTR RKSI  KS   1 1 1 ;
    I = I + 1 ;
  'FIN' BLOC4 ;
  'OUBL' ELEM_I;
  'OUBL' NOEU1;
  'OUBL' NOEU2;
  'OUBL' NOEU3;
  'OUBL' Z1;
  'OUBL' Z3;
  'OUBL' PP1;
  'OUBL' PP2;
  'OUBL' ID1;
*
  TDON1.MOT22.'SOL' = TABLE;
  'SI' ('EXIST' TABPAR 'GAMMAO_X');
    DIR    = 'UX'          ;
    GAMMAO = TABPAR.'GAMMAO_X' ;
    TYP_F = TABPAR.'TYPE' ;
    'SI' ('EGA' ID2 'C');
      MOD_DEC = MOD_SOLC;
      MAT_DEC = MAT_SOLC;
      MAI_DEC = MAI_SOLC;
    'FINSI';
*
    'SI' ('EGA' ID2 'I');
      MOD_DEC = MOD_SOLI;
      MAT_DEC = MAT_SOLI;
      MAI_DEC = MAI_SOLI;
    'FINSI';
*
    'SI' ('EGA' ID2 'M');
      MOD_DEC = MOD_SOLC ET MOD_SOLI;
      MAT_DEC = MAT_SOLC ET MAT_SOLI;
      MAI_DEC = MAI_SOLC ET MAI_SOLI;
    'FINSI';
*
    P_GAM = TABPAR.'PAR_DEC';
    TABRXDV =  DECONV3D TABDV (TABPAR.POINT) BOR_FON MOD_DEC 
                      MAT_DEC DIR GAMMAO F1 F2 FC P_GAM;
*
    TDON1.MOT22.'SOL'.'TAB_DEC' = TABDV;
    TDON1.MOT22.'SOL'.'RES_DEC_X' = TABRXDV;
  'FINSI';
*
  'SI' ('EXIST' TABPAR 'GAMMAO_Y');
    DIR    = 'UY'          ;
    GAMMAO = TABPAR.'GAMMAO_Y' ;
    TYP_F = TABPAR.'TYPE' ;
    'SI' ('EGA' ID2 'C');
      MOD_DEC = MOD_SOLC;
      MAT_DEC = MAT_SOLC;
      MAI_DEC = MAI_SOLC;
    'FINSI';
*
    'SI' ('EGA' ID2 'I');
      MOD_DEC = MOD_SOLI;
      MAT_DEC = MAT_SOLI;
      MAI_DEC = MAI_SOLI;
    'FINSI';
*
    'SI' ('EGA' ID2 'M');
      MOD_DEC = MOD_SOLC ET MOD_SOLI;
      MAT_DEC = MAT_SOLC ET MAT_SOLI;
      MAI_DEC = MAI_SOLC ET MAI_SOLI;
    'FINSI';
*
    P_GAM = TABLE;
    P_GAM. 1 = (TABDV. 1 .'FRONTIERE') 'POIN' 'PROC' (10000. 10000.
                                                              10000.);
    TABRYDV =  DECONV3D TABDV (TABPAR.POINT) BOR_FON MOD_DEC 
                                MAT_DEC DIR GAMMAO F1 F2 FC P_GAM;
*
    TDON1.MOT22.'SOL'.'TAB_DEC' = TABDV;
    TDON1.MOT22.'SOL'.'RES_DEC_Y' = TABRYDV;
  'FINSI';
*
  'SI' ('EXIST' TABPAR 'GAMMAO_Z');
    DIR    = 'UZ'          ;
    GAMMAO = TABPAR.'GAMMAO_Z' ;
    TYP_F = TABPAR.'TYPE' ;
    'SI' ('EGA' ID2 'C');
      MOD_DEC = MOD_SOLC;
      MAT_DEC = MAT_SOLC;
      MAI_DEC = MAI_SOLC;
    'FINSI';
*
    'SI' ('EGA' ID2 'I');
      MOD_DEC = MOD_SOLI;
      MAT_DEC = MAT_SOLI;
      MAI_DEC = MAI_SOLI;
    'FINSI';
*
    'SI' ('EGA' ID2 'M');
      MOD_DEC = MOD_SOLC ET MOD_SOLI;
      MAT_DEC = MAT_SOLC ET MAT_SOLI;
      MAI_DEC = MAI_SOLC ET MAI_SOLI;
    'FINSI';
*
    P_GAM = TABLE;
    P_GAM. 1 = (TABDV. 1 .'FRONTIERE') 'POIN' 'PROC' (10000. 10000.
                                                               10000.);
    TABRZDV =  DECONV3D TABDV (TABPAR.POINT) BOR_FON MOD_DEC 
                                   MAT_DEC DIR GAMMAO F1 F2 FC P_GAM;
*
    TDON1.MOT22.'SOL'.'TAB_DEC' = TABDV;
    TDON1.MOT22.'SOL'.'RES_DEC_Z' = TABRZDV;
  'FINSI';
  'OUBL' P_GAM;
  'OUBL' BOR_FRO;
  'OUBL' BOR_FON;
  'OUBL' MOD_DEC;
  'OUBL' MAT_DEC;
  'OUBL' DIR;
  'OUBL' GAMMAO;
  'OUBL' MAI_DEC;
'FINSI';
*
 
'SI' ('EXIST' (TDON1) 'STRUCTURE');
  TDON1.MOT22.'STRUCTURE' = 'COPI' (TDON1.'STRUCTURE');
'FINSI';
*
**************************************************************
*   resolution dynamique
**************************************************************
*
MESS '';
MESS '--------------------------------------------------------------';
MESS '--                     Calcul Dynamic                       --';
MESS '--------------------------------------------------------------';
*
*  Detection plage de temps
A1 = 0;
'SI' ('EXIST' TABPAR 'GAMMAO_X');
  DT1 = EXTR (EXTR (TABPAR.'GAMMAO_X') ABSC) 2;
  AA1 AA2 A3 = 'MAXI' (TABPAR.'GAMMAO_X');
  'SI' (AA1 > A1);
    A1 = AA1;
  'FINSI'; 
  'OUBL' AA1;
  'OUBL' AA2;
  'OUBL' AA3;
'FINSI';
'SI' ('EXIST' TABPAR 'GAMMAO_Y');
  DT1 = EXTR (EXTR (TABPAR.'GAMMAO_Y') ABSC) 2;
  AA1 AA2 A3 = 'MAXI' (TABPAR.'GAMMAO_Y');
  'SI' (AA1 > A1);
    A1 = AA1;
  'FINSI'; 
  'OUBL' AA1;
  'OUBL' AA2;
  'OUBL' AA3;
'FINSI';
'SI' ('EXIST' TABPAR 'GAMMAO_Z');
  DT1 = EXTR (EXTR (TABPAR.'GAMMAO_Z') ABSC) 2;
  AA1 AA2 A3 = 'MAXI' (TABPAR.'GAMMAO_Z');
  'SI' (AA1 > A1);
    A1 = AA1;
  'FINSI'; 
  'OUBL' AA1;
  'OUBL' AA2;
  'OUBL' AA3;
'FINSI';
A2 = (ENTI((TABPAR.'DELTA_T')/DT1))/2;
NB_PAS = A1 + A2;
NN1 = A1 + A2;
NN2 = A1 - A2;
'SI' (NN2 < 0);
  NN2 = 0;
'FINSI';
'OUBL' A2;
'OUBL' A1;
'OUBL' DT1;
NN1 = ((NN1 - 1)*(4)) + 1;
NN2 = ((NN2 - 1)*(4)) + 1;
NB_PAS = NN1 + 4;
*
'SI' ('EGA' ('TYPE' TABCOL) 'TABLE');
*
*   resolution dynamique sur la colonne
  'SI' ('EGA' V_MOD 'TRID' ) ; 
    'SI' (('DIME' ('EXTR' RIGSOLC 'COMP')) 'EGA' 3);
      CH_DEPI = 'MANU' 'CHPO' 3 MAISOLC UX 0.0 UY 0.0 UZ 0.0 ;
      CH_VITI = 'MANU' 'CHPO' 3 MAISOLC UX 0.0 UY 0.0 UZ 0.0 ;
    'SINON';
      CH_DEPI = 'MANU' 'CHPO' 6 MAISOLC UX 0.0 UY 0.0 UZ 0.0 
                                       RX 0.0 RY 0.0 RZ 0.0;
      CH_VITI = 'MANU' 'CHPO' 6 MAISOLC UX 0.0 UY 0.0 UZ 0.0
                                       RX 0.0 RY 0.0 RZ 0.0;
    'FINSI';
  'FINSI';
  'SI' ('EXIST' TABPAR 'GAMMAO_X');
    TAB_DYXC = TABLE ;
    TAB_DYXC.'DEPL' = CH_DEPI ;
    TAB_DYXC.'VITE' = CH_VITI ;
    TAB_DYXC.'RIGI' = RIGSOLC ;
    TAB_DYXC.'MASS' = MASSOLC;
    TAB_DYXC.'AMOR' = AMOSOLC 'ET' (TABRXDV2.'AMOR');
    TAB_DYXC.'CHAR' = TABRXDV2.'CHAR' ;
    TAB_DYXC.'FREQ' = TABRXDV2.'FCDYN' ;
    TAB_DYXC.'DEBU' = 0.0 ;
    TAB_DYXC.'INST' = 'PROG' 0.0 PAS ((1)*(TABRXDV2.'PAS'))
                                  NPAS ( NB_PAS - 1 ) ;
    TRESUXC = DYNAMIC TAB_DYXC ;
    'TEMPS' ;
    'MENAGE' ;
    TDON1.MOT22.'COLONNE'.'RES_DYN_X' = 'COPI' TRESUXC;
  'FINSI';
  'SI' ('EXIST' TABPAR 'GAMMAO_Y');
    TAB_DYYC = TABLE ;
    TAB_DYYC.'DEPL' = CH_DEPI ;
    TAB_DYYC.'VITE' = CH_VITI ;
    TAB_DYYC.'RIGI' = RIGSOLC ;
    TAB_DYYC.'MASS' = MASSOLC;
    TAB_DYYC.'AMOR' = AMOSOLC 'ET' (TABRYDV2.'AMOR');
    TAB_DYYC.'CHAR' = TABRYDV2.'CHAR' ;
    TAB_DYYC.'FREQ' = TABRYDV2.'FCDYN' ;
    TAB_DYYC.'DEBU' = 0.0 ;
    TAB_DYYC.'INST' = 'PROG' 0.0 PAS ((1)*(TABRYDV2.'PAS'))
                                  NPAS ( NB_PAS - 1 ) ;
    TRESUYC = DYNAMIC TAB_DYYC ;
    'TEMPS' ;
    'MENAGE' ;
    TDON1.MOT22.'COLONNE'.'RES_DYN_Y' = 'COPI' TRESUYC;
  'FINSI';
  'SI' ('EXIST' TABPAR 'GAMMAO_Z');
    TAB_DYZC = TABLE ;
    TAB_DYZC.'DEPL' = CH_DEPI ;
    TAB_DYZC.'VITE' = CH_VITI ;
    TAB_DYZC.'RIGI' = RIGSOLC ;
    TAB_DYZC.'MASS' = MASSOLC;
    TAB_DYZC.'AMOR' = AMOSOLC 'ET' (TABRZDV2.'AMOR');
    TAB_DYZC.'CHAR' = TABRZDV2.'CHAR' ;
    TAB_DYZC.'FREQ' = TABRZDV2.'FCDYN' ;
    TAB_DYZC.'DEBU' = 0.0 ;
    TAB_DYZC.'INST' = 'PROG' 0.0 PAS ((1)*(TABRZDV2.'PAS'))
                                  NPAS ( NB_PAS - 1 ) ;
    TRESUZC = DYNAMIC TAB_DYZC ;
    'TEMPS' ;
    'MENAGE' ;
    TDON1.MOT22.'COLONNE'.'RES_DYN_Z' = 'COPI' TRESUZC;
  'FINSI';
'SINON';
*
*   resolution dynamique du systeme sol+structure
  'SI' ('EGA' ('TYPE' MAI_SOLC) 'MAILLAGE' ) ;
    MAITOT = MAI_SOLC;
    MODTOT = MOD_SOLC;
    RIGTOT = RIG_SOLC;
    MASTOT = MAS_SOLC;
    AMOTOT = AMO_SOLC;  
    'SI' ('EGA' ('TYPE' MAI_SOLI) 'MAILLAGE' ) ;
      MAITOT = MAITOT 'ET' MAI_SOLI;
      MODTOT = MODTOT 'ET' MOD_SOLI;
      RIGTOT = RIGTOT 'ET' RIG_SOLI;
      MASTOT = MASTOT 'ET' MAS_SOLI;
      AMOTOT = AMOTOT 'ET' AMO_SOLI;  
    'FINSI';
  'SINON';
    'SI' ('EGA' ('TYPE' MAI_SOLI) 'MAILLAGE' );
      MAITOT = MAI_SOLI;
      MODTOT = MOD_SOLI;
      RIGTOT = RIG_SOLI;
      MASTOT = MAS_SOLI;
      AMOTOT = AMO_SOLI;  
    'FINSI';
  'FINSI';
  'SI' ('EGA' ('TYPE' TABSTR) 'TABLE' );
      MAITOT = MAITOT 'ET' (TABSTR.'MAILLAGE');
      RIGTOT = RIGTOT 'ET' (TABSTR.'RIGIDITE');
      MASTOT = MASTOT 'ET' (TABSTR.'MASSE');
      AMOTOT = AMOTOT 'ET' (TABSTR.'AMORTISSEMENT');  
    'SI' ('EXIST' (TABSTR) 'BLOCAGES_MECANIQUES' );
       RIGTOT = RIGTOT 'ET' (TABSTR.'BLOCAGES_MECANIQUES');
    'FINSI';
  'FINSI';
*  'ELIM' 0.001 MAITOT;
*
  'SI' ('EGA' V_MOD 'TRID' ) ; 
    'SI' (('DIME' ('EXTR' RIGTOT 'COMP')) 'EGA' 3);
      CH_DEPI = 'MANU' 'CHPO' 3 MAITOT UX 0.0 UY 0.0 UZ 0.0 ;
      CH_VITI = 'MANU' 'CHPO' 3 MAITOT UX 0.0 UY 0.0 UZ 0.0 ;
    'SINON';
      CH_DEPI = 'MANU' 'CHPO' 6 MAITOT UX 0.0 UY 0.0 UZ 0.0 
                                       RX 0.0 RY 0.0 RZ 0.0;
      CH_VITI = 'MANU' 'CHPO' 6 MAITOT UX 0.0 UY 0.0 UZ 0.0
                                       RX 0.0 RY 0.0 RZ 0.0;
    'FINSI';
  'FINSI';
*
  'SI' ('EXIST' TABPAR 'GAMMAO_X');
    TAB_DYX = TABLE ;
    TAB_DYX.'DEPL' = CH_DEPI ;
    TAB_DYX.'VITE' = CH_VITI ;
    TAB_DYX.'RIGI' = RIGTOT ;
    TAB_DYX.'MASS' = MASTOT;
    TAB_DYX.'AMOR' = AMOTOT 'ET' (TABRXDV.'AMOR');
    TAB_DYX.'CHAR' = TABRXDV.'CHAR' ;
    TAB_DYX.'FREQ' = TABRXDV.'FCDYN' ;
    TAB_DYX.'DEBU' = 0.0 ;
    TAB_DYX.'INST' = 'PROG' 0.0 PAS ((1)*(TABRXDV.'PAS'))
                                NPAS ( NB_PAS - 1 ) ;
    TRESUX = DYNAMIC TAB_DYX ;
    'TEMPS' ;
    'MENAGE' ;
    TDON1.MOT22.'RES_DYN_X' = 'COPI' TRESUX;
  'FINSI';
  'SI' ('EXIST' TABPAR 'GAMMAO_Y');
    TAB_DYY = TABLE ;
    TAB_DYY.'DEPL' = CH_DEPI ;
    TAB_DYY.'VITE' = CH_VITI ;
    TAB_DYY.'RIGI' = RIGTOT ;
    TAB_DYY.'MASS' = MASTOT;
    TAB_DYY.'AMOR' = AMOTOT 'ET' (TABRYDV.'AMOR');
    TAB_DYY.'CHAR' = TABRYDV.'CHAR' ;
    TAB_DYY.'FREQ' = TABRYDV.'FCDYN' ;
    TAB_DYY.'DEBU' = 0.0 ;
    TAB_DYY.'INST' = 'PROG' 0.0 PAS ((1)*(TABRYDV.'PAS'))
                                NPAS ( NB_PAS - 1 ) ;
    TRESUY = DYNAMIC TAB_DYY ;
    'TEMPS' ;
    'MENAGE' ;
    TDON1.MOT22.'RES_DYN_Y' = 'COPI' TRESUY;
  'FINSI';
  'SI' ('EXIST' TABPAR 'GAMMAO_Z');
    TAB_DYZ = TABLE ;
    TAB_DYZ.'DEPL' = CH_DEPI ;
    TAB_DYZ.'VITE' = CH_VITI ;
    TAB_DYZ.'RIGI' = RIGTOT ;
    TAB_DYZ.'MASS' = MASTOT;
    TAB_DYZ.'AMOR' = AMOTOT 'ET' (TABRZDV.'AMOR');
    TAB_DYZ.'CHAR' = TABRZDV.'CHAR' ;
    TAB_DYZ.'FREQ' = TABRZDV.'FCDYN' ;
    TAB_DYZ.'DEBU' = 0.0 ;
    TAB_DYZ.'INST' = 'PROG' 0.0 PAS ((1)*(TABRZDV.'PAS'))
                                NPAS ( NB_PAS - 1 ) ;
    TRESUZ = DYNAMIC TAB_DYZ ;
    'TEMPS' ;
    'MENAGE' ;
    TDON1.MOT22.'RES_DYN_Z' = 'COPI' TRESUZ;
  'FINSI';
'FINSI';
*
**************************************************************
*   deductions reponse
**************************************************************
*
*
*
**************************************************************
*   deduction du champ de deformation maximale
**************************************************************
*
MESS '';
MESS '--------------------------------------------------------------';
MESS '--           Calcul Deformation Maximale                    --';
MESS '--------------------------------------------------------------';
'SI' ('EGA' ('TYPE' TABSOL_C) 'TABLE');
  DEF_C = 'TABLE';
  'SI' ('EGA' ('TYPE' TABCOL) 'TABLE');
     DEFC = 'TABLE';
     II1 = 1;
     'REPE' IND1 (NZONE_C);
       DEF_C. II1 = 'TABLE';
       JJ1 = 1;
       'REPE' IND2 ('DIME' (SOL_C. II1));
         DEF_C. II1 . JJ1 = 'MANU' 'CHML' (MOD_C. II1 . JJ1) 'DMAX' 0.0;
         JJ1 = JJ1 + 1;
       'FIN' IND2;
       II1 = II1 + 1;
     'FIN' IND1;
     II1 = 1;
     'REPE' IND1 (NZONEC);
       DEFC. II1 = 'TABLE';
       JJ1 = 1;
       'REPE' IND2 ('DIME' (SOLCOL. II1));
         DEFC. II1 . JJ1 = 'MANU' 'CHML' (MODCOL. II1 . JJ1) 'DMAX' 0.0;
         JJ1 = JJ1 + 1;
       'FIN' IND2;
       II1 = II1 + 1;
     'FIN' IND1;
  'SINON';
    II1 = 1;
    'REPE' IND1 (NZONE_C);
      DEF_C. II1  = 'MANU' 'CHML' (MOD_C. II1) 'DMAX' 0.0;
      II1 = II1 + 1;
    'FIN' IND1;
  'FINSI';
'FINSI';
*
'SI' ('EGA' ('TYPE' TABSOL_I) 'TABLE');
  DEF_I = 'TABLE';
  II1 = 1;
  'REPE' IND1 ('DIME' SOL_I);
*    DEF_I. II1 = 'TABLE';
*    JJ1 = 1;
*    'REPE' IND2 ('DIME' (SOL_I. II1));
*      DEF_I. II1 . JJ1 = 'MANU' 'CHML' (MOD_I. II1 . JJ1) 'DMAX' 0.0;
*      JJ1 = JJ1 + 1;
*    'FIN' IND2;
    DEF_I. II1 = 'MANU' 'CHML' (MOD_I. II1) 'DMAX' 0.0;
    II1 = II1 + 1;
  'FIN' IND1;
'FINSI';
*
'SI' ('EGA' ('TYPE' TABCOL) 'TABLE');
  I = 1;
  'REPE' IND1 (NN1);
    'SI' (I '>' (NN2));
      'SI' ('EGA' ('TYPE' CH1) 'CHPOINT');
        'SI' ('EGA' ('TYPE' TRESUXC) 'TABLE');
          CH1 = CH1 + (TRESUXC. I .'DEPL');
        'FINSI';
        'SI' ('EGA' ('TYPE' TRESUYC) 'TABLE');
          CH1 = CH1 + (TRESUYC. I .'DEPL');
        'FINSI';
        'SI' ('EGA' ('TYPE' TRESUZC) 'TABLE');
          CH1 = CH1 + (TRESUZC. I .'DEPL');
        'FINSI';
      'SINON';
        'SI' ('EGA' ('TYPE' TRESUXC) 'TABLE');
          CH1 = (TRESUXC. I .'DEPL');
        'FINSI';
        'SI' ('EGA' ('TYPE' TRESUYC) 'TABLE');
          CH1 = (TRESUYC. I .'DEPL');
        'FINSI';
        'SI' ('EGA' ('TYPE' TRESUZC) 'TABLE');
          CH1 = (TRESUZC. I .'DEPL');
        'FINSI';
      'FINSI';
 
*      DEPSOL = 'REDU' (TRESUC. I .'DEPL') ('EXTR' MODSOLC 'MAIL');
      DEPSOL = 'REDU' (CH1) ('EXTR' MODSOLC 'MAIL');
      'OUBL' CH1;
 
      EPSSOL = 'EPSI' DEPSOL MODSOLC;
      EPSSOL   = 'CHAN' 'GRAVITE' MODSOLC EPSSOL;
*
*  Cas 3D - COLONNE
      'SI' ('EGA' V_MOD 'TRID' ) ; 
         II1 = 1;
         'REPE' IND3 (NZONEC);
           JJ1 = 1;
           'REPE' IND4 ('DIME' (SOLCOL. II1));
             MOD1 = 'REDU' MODSOLC (SOLCOL. II1 . JJ1);
             EPS1 = 'REDU' EPSSOL (SOLCOL. II1 . JJ1);
             E11  =   EXCO 'EPXX' EPS1 'SMXX' ;
             E22  =   EXCO 'EPYY' EPS1 'SMYY' ;
             E33  =   EXCO 'EPZZ' EPS1 'SMZZ' ;
             G12  =   EXCO 'GAXY' EPS1 'SMXY' ;
             G13  =   EXCO 'GAXZ' EPS1 'SMXZ' ;
             G23  =   EXCO 'GAYZ' EPS1 'SMYZ' ;
             E12  =   G12 / 2.0 ;
             E13  =   G13 / 2.0 ;
             E23  =   G23 / 2.0 ;
             EPS2 = E11;
             EPS2 = EPS2 + (E22);
             EPS2 = EPS2 + (E33);
             EPS2 = EPS2 + (E12);
             EPS2 = EPS2 + (E13);
             EPS2 = EPS2 + (E23);
             EPS2 = CHAN 'TYPE' EPS2 'CONTRAINTES' ;
             EPS2 = 'PRIN' EPS2 MOD1;
             'DETR' E11;
             'DETR' E22;
             'DETR' E33;
             'DETR' G12;
             'DETR' G13;
             'DETR' G23;
             'DETR' E12;
             'DETR' E13;
             'DETR' E23;
             S11  = 'EXCO' SI11 EPS2 'DMAX' ;
             S33  = 'EXCO' SI33 EPS2 'DMAX' ;
             CHDEFM = (S11 - S33)*(0.5);
             VAL1 = 'EXTR' CHDEFM 'DMAX' 1 1 1;
             VAL2 = 'EXTR' (DEFC. II1 . JJ1) 'DMAX' 1 1 1;
             'SI' (VAL1 > VAL2);
                DEFC. II1 . JJ1= CHDEFM;
             'FINSI';
             JJ1 = JJ1 + 1;
           'FIN' IND4;
           II1 = II1 + 1;
         'FIN' IND3;
       'FINSI';
    'FINSI';
    I = I + 1;
  'FIN' IND1;
'SINON';
*
  CHDEFM_I = 'TABLE';
  CHDEFM_C = 'TABLE';
  I = 1;
  'REPE' IND1 (NN1);
    'SI' (I '>' (NN2));
      'SI' ('EGA' ('TYPE' CH1) 'CHPOINT');
        'SI' ('EGA' ('TYPE' TRESUX) 'TABLE');
          CH1 = CH1 + (TRESUX. I .'DEPL');
        'FINSI';
        'SI' ('EGA' ('TYPE' TRESUY) 'TABLE');
          CH1 = CH1 + (TRESUY. I .'DEPL');
        'FINSI';
        'SI' ('EGA' ('TYPE' TRESUZ) 'TABLE');
          CH1 = CH1 + (TRESUZ. I .'DEPL');
        'FINSI';
      'SINON';
        'SI' ('EGA' ('TYPE' TRESUX) 'TABLE');
          CH1 = (TRESUX. I .'DEPL');
        'FINSI';
        'SI' ('EGA' ('TYPE' TRESUY) 'TABLE');
          CH1 = (TRESUY. I .'DEPL');
        'FINSI';
        'SI' ('EGA' ('TYPE' TRESUZ) 'TABLE');
          CH1 = (TRESUZ. I .'DEPL');
        'FINSI';
      'FINSI';     
*      DEPSOL = 'REDU' (TRESU. I .'DEPL') ('EXTR' MODTOT 'MAIL');
      DEPSOL = 'REDU' (CH1) ('EXTR' MODTOT 'MAIL');
      'OUBL' CH1;
      EPSSOL = 'EPSI' DEPSOL MODTOT;
      EPSSOL   = 'CHAN' 'GRAVITE' MODTOT EPSSOL;
*
*  Cas 3D
      'SI' ('EGA' V_MOD 'TRID' ) ; 
*
********  VECCHIA VERSIONE
*         'SI' ('EGA' ('TYPE' TABSOL_C) 'TABLE');
*           II1 = 1;
*           'REPE' IND3 (NZONE_C);
*             MOD1 = 'REDU' MODTOT (SOL_C. II1);
*             EPS1 = 'REDU' EPSSOL (SOL_C. II1);
*             E11  =   EXCO 'EPXX' EPS1 'SMXX' ;
*             E22  =   EXCO 'EPYY' EPS1 'SMYY' ;
*             E33  =   EXCO 'EPZZ' EPS1 'SMZZ' ;
*             G12  =   EXCO 'GAXY' EPS1 'SMXY' ;
*             G13  =   EXCO 'GAXZ' EPS1 'SMXZ' ;
*             G23  =   EXCO 'GAYZ' EPS1 'SMYZ' ;
*             E12  =   G12 / 2.0 ;
*             E13  =   G13 / 2.0 ;
*             E23  =   G23 / 2.0 ;
*             EPS2 = E11;
*             EPS2 = EPS2 + (E22);
*             EPS2 = EPS2 + (E33);
*             EPS2 = EPS2 + (E12);
*             EPS2 = EPS2 + (E13);
*             EPS2 = EPS2 + (E23);
*             EPS2 = CHAN 'TYPE' EPS2 'CONTRAINTES' ;
*             EPS2 = 'PRIN' EPS2 MOD1;
*             'DETR' E11;
*             'DETR' E22;
*             'DETR' E33;
*             'DETR' G12;
*             'DETR' G13;
*             'DETR' G23;
*             'DETR' E12;
*             'DETR' E13;
*             'DETR' E23;
*             S11  = 'EXCO' SI11 EPS2 'DMAX' ;
*             S33  = 'EXCO' SI33 EPS2 'DMAX' ;
*             CHDEFM = (S11 - S33)*(0.5);
*             VAL1 = 'EXTR' CHDEFM 'DMAX' 1 1 1;
*             VAL2 = 'EXTR' (DEF_C. II1) 'DMAX' 1 1 1;
*             'SI' (VAL1 > VAL2);
*                DEF_C. II1 = CHDEFM;
*             'FINSI';
*           II1 = II1 + 1;
*           'FIN' IND3;
*         'FINSI';
*         'SI' ('EGA' ('TYPE' TABSOL_I) 'TABLE');
*           II1 = 1;
*           'REPE' IND3 (NZONE_I);
*             JJ1 = 1;
*             'REPE' IND4 ('DIME' SOL_I. II1);
*               MOD1 = 'REDU' MODTOT (SOL_I. II1 . JJ1);
*               EPS1 = 'REDU' EPSSOL (SOL_I. II1 . JJ1);
*               E11  =   'EXCO' 'EPXX' EPS1 'SMXX' ;
*               E22  =   'EXCO' 'EPYY' EPS1 'SMYY' ;
*               E33  =   'EXCO' 'EPZZ' EPS1 'SMZZ' ;
*               G12  =   'EXCO' 'GAXY' EPS1 'SMXY' ;
*               G13  =   'EXCO' 'GAXZ' EPS1 'SMXZ' ;
*               G23  =   'EXCO' 'GAYZ' EPS1 'SMYZ' ;
*               E12  =   G12 / 2.0 ;
*               E13  =   G13 / 2.0 ;
*               E23  =   G23 / 2.0 ;
*               EPS2 = E11;
*               EPS2 = EPS2 + (E22);
*               EPS2 = EPS2 + (E33);
*               EPS2 = EPS2 + (E12);
*               EPS2 = EPS2 + (E13);
*               EPS2 = EPS2 + (E23);
*               EPS2 = 'CHAN' 'TYPE' EPS2 'CONTRAINTES' ;
*               EPS2 = 'PRIN' EPS2 MOD1;
*               'DETR' E11;
*               'DETR' E22;
*               'DETR' E33;
*               'DETR' G12;
*               'DETR' G13;
*               'DETR' G23;
*               'DETR' E12;
*               'DETR' E13;
*               'DETR' E23;
*               EPS2 = 'PRIN' EPS1 MOD1;
*               S11  = 'EXCO' SI11 EPS2 'DMAX' ;
*               S33  = 'EXCO' SI33 EPS2 'DMAX' ;
*               CHDEFM = (S11 - S33)*(0.5);
*               VAL1 = 'EXTR' CHDEFM 'DMAX' 1 1 1;
*               VAL2 = 'EXTR' (DEF_I. II1 . JJ1) 'DMAX' 1 1 1;
*               'SI' (VAL1 > VAL2);
*                  DEF_I. II1 . JJ1 = CHDEFM;
*               'FINSI';
*               JJ1 = JJ1 + 1;
*             'FIN' IND4;
*             II1 = II1 + 1;
*           'FIN' IND3;
*         'FINSI';
*
********  NUOVA VERSIONE
         'SI' ('EGA' ('TYPE' TABSOL_C) 'TABLE');
*           CHDEFM_C = 'TABLE';
           II1 = 1;
           'REPE' IND3 (NZONE_C);
             'SI' (I EGA (NN2 + 1));
             MOD1 = 'REDU' MODTOT (SOL_C. II1);
             EPS1 = 'REDU' EPSSOL (SOL_C. II1);
             E11  =   EXCO 'EPXX' EPS1 'SMXX' ;
             E22  =   EXCO 'EPYY' EPS1 'SMYY' ;
             E33  =   EXCO 'EPZZ' EPS1 'SMZZ' ;
             G12  =   EXCO 'GAXY' EPS1 'SMXY' ;
             G13  =   EXCO 'GAXZ' EPS1 'SMXZ' ;
             G23  =   EXCO 'GAYZ' EPS1 'SMYZ' ;
             E12  =   G12 / 2.0 ;
             E13  =   G13 / 2.0 ;
             E23  =   G23 / 2.0 ;
             EPS2 = E11;
             EPS2 = EPS2 + (E22);
             EPS2 = EPS2 + (E33);
             EPS2 = EPS2 + (E12);
             EPS2 = EPS2 + (E13);
             EPS2 = EPS2 + (E23);
             EPS2 = CHAN 'TYPE' EPS2 'CONTRAINTES' ;
             EPS2 = 'PRIN' EPS2 MOD1;
             'DETR' E11;
             'DETR' E22;
             'DETR' E33;
             'DETR' G12;
             'DETR' G13;
             'DETR' G23;
             'DETR' E12;
             'DETR' E13;
             'DETR' E23;
             S11  = 'EXCO' SI11 EPS2 'DMAX' ;
             S33  = 'EXCO' SI33 EPS2 'DMAX' ;
*             CHDEFM_C. II1 = (S11 - S33)*(0.5);
             CHDEFM_C. II1 = (S11 - S33)*(1.0);
             'SINON';
             MOD1 = 'REDU' MODTOT (SOL_C. II1);
             EPS1 = 'REDU' EPSSOL (SOL_C. II1);
             E11  =   EXCO 'EPXX' EPS1 'SMXX' ;
             E22  =   EXCO 'EPYY' EPS1 'SMYY' ;
             E33  =   EXCO 'EPZZ' EPS1 'SMZZ' ;
             G12  =   EXCO 'GAXY' EPS1 'SMXY' ;
             G13  =   EXCO 'GAXZ' EPS1 'SMXZ' ;
             G23  =   EXCO 'GAYZ' EPS1 'SMYZ' ;
             E12  =   G12 / 2.0 ;
             E13  =   G13 / 2.0 ;
             E23  =   G23 / 2.0 ;
             EPS2 = E11;
             EPS2 = EPS2 + (E22);
             EPS2 = EPS2 + (E33);
             EPS2 = EPS2 + (E12);
             EPS2 = EPS2 + (E13);
             EPS2 = EPS2 + (E23);
             EPS2 = CHAN 'TYPE' EPS2 'CONTRAINTES' ;
             EPS2 = 'PRIN' EPS2 MOD1;
             'DETR' E11;
             'DETR' E22;
             'DETR' E33;
             'DETR' G12;
             'DETR' G13;
             'DETR' G23;
             'DETR' E12;
             'DETR' E13;
             'DETR' E23;
             S11  = 'EXCO' SI11 EPS2 'DMAX' ;
             S33  = 'EXCO' SI33 EPS2 'DMAX' ;
*             CHDEFM_C. II1 = MAXI (CHDEFM_C. II1) ((S11 - S33)*(0.5));
             CH_111 = (CHDEFM_C. II1);
*             CH_222 = ((S11 - S33)*(0.5));
             CH_222 = ((S11 - S33)*(1.0));
             DIF_11 = CH_111 - CH_222;
             MASQ1 = MASQ DIF_11 'INFEURIEUR' 0.D0;
             MASQ2 = MASQ DIF_11 'EGSUPE' 0.D0;
             CHDEFM_C. II1 = (CH_111 '*' MASQ2) + (CH_222 '*' MASQ1);
             'FINSI';
             DEF_C. II1 = CHDEFM_C. II1;
             II1 = II1 + 1;
           'FIN' IND3;
         'FINSI';
*
         'SI' ('EGA' ('TYPE' TABSOL_I) 'TABLE');
*           CHDEFM_I = 'TABLE';
           II1 = 1;
           'REPE' IND3 (NZONE_I);
             'SI' (I EGA (NN2 + 1));
             MOD1 = 'REDU' MODTOT (SOL_I. II1);
             EPS1 = 'REDU' EPSSOL (SOL_I. II1);
             E11  =   'EXCO' 'EPXX' EPS1 'SMXX' ;
             E22  =   'EXCO' 'EPYY' EPS1 'SMYY' ;
             E33  =   'EXCO' 'EPZZ' EPS1 'SMZZ' ;
             G12  =   'EXCO' 'GAXY' EPS1 'SMXY' ;
             G13  =   'EXCO' 'GAXZ' EPS1 'SMXZ' ;
             G23  =   'EXCO' 'GAYZ' EPS1 'SMYZ' ;
             E12  =   G12 / 2.0 ;
             E13  =   G13 / 2.0 ;
             E23  =   G23 / 2.0 ;
             EPS2 = E11;
             EPS2 = EPS2 + (E22);
             EPS2 = EPS2 + (E33);
             EPS2 = EPS2 + (E12);
             EPS2 = EPS2 + (E13);
             EPS2 = EPS2 + (E23);
             EPS2 = 'CHAN' 'TYPE' EPS2 'CONTRAINTES' ;
             EPS2 = 'PRIN' EPS2 MOD1;
             'DETR' E11;
             'DETR' E22;
             'DETR' E33;
             'DETR' G12;
             'DETR' G13;
             'DETR' G23;
             'DETR' E12;
             'DETR' E13;
             'DETR' E23;
*             EPS2 = 'PRIN' EPS1 MOD1;
*             EPS2 = 'PRIN' EPS2 MOD1;
             S11  = 'EXCO' SI11 EPS2 'DMAX' ;
             S33  = 'EXCO' SI33 EPS2 'DMAX' ;
*             CHDEFM_I. II1 = (S11 - S33)*(0.5);
             CHDEFM_I. II1 = (S11 - S33)*(1.0);
             'SINON';
             MOD1 = 'REDU' MODTOT (SOL_I. II1);
             EPS1 = 'REDU' EPSSOL (SOL_I. II1);
             E11  =   'EXCO' 'EPXX' EPS1 'SMXX' ;
             E22  =   'EXCO' 'EPYY' EPS1 'SMYY' ;
             E33  =   'EXCO' 'EPZZ' EPS1 'SMZZ' ;
             G12  =   'EXCO' 'GAXY' EPS1 'SMXY' ;
             G13  =   'EXCO' 'GAXZ' EPS1 'SMXZ' ;
             G23  =   'EXCO' 'GAYZ' EPS1 'SMYZ' ;
             E12  =   G12 / 2.0 ;
             E13  =   G13 / 2.0 ;
             E23  =   G23 / 2.0 ;
             EPS2 = E11;
             EPS2 = EPS2 + (E22);
             EPS2 = EPS2 + (E33);
             EPS2 = EPS2 + (E12);
             EPS2 = EPS2 + (E13);
             EPS2 = EPS2 + (E23);
             EPS2 = 'CHAN' 'TYPE' EPS2 'CONTRAINTES' ;
             EPS2 = 'PRIN' EPS2 MOD1;
             'DETR' E11;
             'DETR' E22;
             'DETR' E33;
             'DETR' G12;
             'DETR' G13;
             'DETR' G23;
             'DETR' E12;
             'DETR' E13;
             'DETR' E23;
*             EPS2 = 'PRIN' EPS1 MOD1;
*             EPS2 = 'PRIN' EPS2 MOD1;
             S11  = 'EXCO' SI11 EPS2 'DMAX' ;
             S33  = 'EXCO' SI33 EPS2 'DMAX' ;
*             CHDEFM_I. II1 = MAXI (CHDEFM_I. II1) ((S11 - S33)*(0.5));
             CH_111 = (CHDEFM_I. II1);
*             CH_222 = ((S11 - S33)*(0.5));
             CH_222 = ((S11 - S33)*(1.0));
             DIF_11 = CH_111 - CH_222;
             MASQ1 = MASQ DIF_11 'INFEURIEUR' 0.D0;
             MASQ2 = MASQ DIF_11 'EGSUPE' 0.D0;
             CHDEFM_I. II1 = (CH_111 '*' MASQ2) + (CH_222 '*' MASQ1);
             'FINSI';
             DEF_I. II1 = CHDEFM_I. II1;
             II1 = II1 + 1;
           'FIN' IND3;
         'FINSI';
      'FINSI';
    'FINSI';
    I = I + 1;
  'FIN' IND1;
'FINSI';
'OUBL' NN1;
'OUBL' NN2;
'OUBL' DEPSOL;
'OUBL' EPSSOL;
'OUBL' MOD1;
'OUBL' EPS1;
'OUBL' EPS2;
'OUBL' S11;
'OUBL' S33;
'OUBL' CHDEFM;
'OUBL' VAL1;
'OUBL' VAL2;
'OUBL' II1;
'OUBL' JJ1;
*
'SI' ('EGA' ('TYPE' DEF_I) 'TABLE');
  DEF_SOLI = 'VIDE' 'MCHAML';
*  NN1 = 'DIME' DEF_I;
  I = 1;
  'REPE' IND1 (NZONE_I);
*    J = 1;
*    NN2 = 'DIME' (DEF_I. I);
*    'REPE' IND2 (NN2);
*      DEF_SOLI = DEF_SOLI 'ET' (DEF_I. I . J);
*      J = J + 1;
*    'FIN' IND2;
    DEF_SOLI = DEF_SOLI 'ET' (DEF_I. I);
    I = I + 1;
  'FIN' IND1;
'FINSI';
*
'SI' ('EGA' ('TYPE' DEF_C) 'TABLE');
  DEF_SOLC = 'VIDE' 'MCHAML';
*  NN1 = 'DIME' DEF_C;
  I = 1;
  'REPE' IND1 (NZONE_C);
    DEF_SOLC = DEF_SOLC 'ET' (DEF_C. I);
    I = I + 1;
  'FIN' IND1;
'FINSI';
*
'SI' ('EGA' ('TYPE' DEFC) 'TABLE');
  DEF_SOLI = 'VIDE' 'MCHAML';
*  NN1 = 'DIME' DEFC;
  I = 1;
  'REPE' IND1 (NZONEC);
    J = 1;
    NN2 = 'DIME' (DEFC. I);
    'REPE' IND2 (NN2);
      DEFCOL = DEFCOL 'ET' (DEFC. I . J);
      J = J + 1;
    'FIN' IND2;
    I = I + 1;
  'FIN' IND1;
'FINSI';
*
**************************************************************
*   Modification des matrices
**************************************************************
*
MESS '';
MESS '--------------------------------------------------------------';
MESS '--           Modification termes K et C                     --';
MESS '--------------------------------------------------------------';
  ALPHA  = 4.0 * PI * F1 * F2 / ( F1 + F2 ) ;
  BETA   = 1.0 / ( PI * ( F1 + F2 ) )       ;
 
'SI' (EXIST (TDON1.MOT33.'SOL') 'CONSTANTE');
*
  ESO_C_N = TABLE;
  MAT_C_N = TABLE;
  KSI_C_N = TABLE;
  RIG_C_N = TABLE;
  AMO_C_N = TABLE;
  GGO_C_N = TABLE;
*
'FINSI';
*
'SI' (EXIST (TDON1.MOT33.'SOL') 'ITERATION');
  TABRES_I = TDON1.MOT33.'SOL'.'ITERATION';
*
  ESO_I_N = TABLE;
  MAT_I_N = TABLE;
  KSI_I_N = TABLE;
  RIG_I_N = TABLE;
  AMO_I_N = TABLE;
  GGO_I_N = TABLE;
*
'FINSI';
*
'SI' (EXIST (TDON1.MOT33) 'COLONNE');
  TABRESC = TDON1.MOT33.'COLONNE';
*
  ESOCOL_N = TABLE;
  MATCOL_N = TABLE;
  KSICOL_N = TABLE;
  RIGCOL_N = TABLE;
  AMOCOL_N = TABLE;
  GGOCOL_N = TABLE;
*
'FINSI';
*
'SI' ('EXIST' (TDON1.MOT33.'SOL') 'ITERATION');
*  NN1 = 'DIME' DEF_I;
  I = 1;
  'REPE' IND1 (NZONE_I);
*
*   VECCHIA VERSIONE
*    ESO_I_N. I = TABLE;
*    MAT_I_N. I = TABLE;
*    KSI_I_N. I = TABLE;
*    RIG_I_N. I = TABLE;
*    AMO_I_N. I = TABLE;
*    GGO_I_N. I = TABLE;
*    J = 1;
*    'REPE' IND2 ('DIME' (DEF_I. I));
*      LLGAM1 = 'EXTR' (TABSOL_I. I .'G_GAMMA') ABSC;
*      LLGAM2 = 'EXTR' (TABSOL_I. I .'H_GAMMA') ABSC;
*      LLGG01 = 'EXTR' (TABSOL_I. I .'G_GAMMA') ORDO; 
*      LLEPS1 = 'EXTR' (TABSOL_I. I .'H_GAMMA') ORDO;
*      VDEF1 = (TABPAR.'CHI')*('EXTR' (DEF_I. I .J) 
*                                      'DMAX' 1 1 1);
*      'SI' (VDEF1 < ('EXTR' LLGAM1 1));
*        VGG01 = 'EXTR' LLGG01 1;
*      'SINON';
*        'SI' (VDEF1 > ('EXTR' LLGAM1 ('DIME' LLGAM1)));
*          VGG01 = 'EXTR' LLGG01 ('DIME' LLGAM1);
*        'SINON';
*          VGG01 = 'EXTR' (IPOL (LOG ('PROG' VDEF1))
*                               (LOG LLGAM1) LLGG01) 1 ;
*
*        'FINSI';
*      'FINSI';
*      'SI' (VDEF1 < ('EXTR' LLGAM2 1));
*        VEPS1 = 'EXTR' LLEPS1 1;
*      'SINON';
*        'SI' (VDEF1 > ('EXTR' LLGAM2 ('DIME' LLGAM2)));
*          VEPS1 = 'EXTR' LLEPS1 ('DIME' LLGAM2);
*        'SINON';
*          VEPS1 = EXTR (IPOL (LOG (PROG VDEF1))
*                               (LOG LLGAM2) LLEPS1) 1 ;
*        'FINSI';
*      'FINSI';
*
*      ESO_I_N. I . J = VGG01*(TABSOL_I. I .'MODULE_E');
*      MAT_I_N. I . J = 'MATE' (MOD_I. I . J ) 
*                 'YOUN' (ESO_I_N. I . J) 'NU' (TABSOL_I. I .'POISSON')
*                             'RHO' (TABSOL_I. I .'MASSE_VOLUMIQUE');
*      GGO_I_N. I . J = 'MANU' 'CHML' MOD_I.I.J 'G_GO' VGG01;
*      KSI_I_N. I . J = 'MANU' 'CHML' MOD_I.I.J 'KS' VEPS1;
*      RIG_I_N. I . J = 'RIGI' (MOD_I.I.J) (MAT_I_N.I.J) ;
*      AMO_I_N. I . J = (VEPS1) * (((ALPHA) * (MAS_I. I . J))
*                          'ET' ((BETA) * (RIG_I_N. I . J)));
*
*      J = J + 1;
*    'FIN' IND2;
*
*   NUOVA VERSIONE
    LLGAM1 = 'EXTR' (TABSOL_I. I .'G_GAMMA') ABSC;
    LLGAM2 = 'EXTR' (TABSOL_I. I .'H_GAMMA') ABSC;
    LLGAM1 = LOG(LLGAM1);
    LLGAM2 = LOG(LLGAM2);
    LLGG01 = 'EXTR' (TABSOL_I. I .'G_GAMMA') ORDO; 
    LLEPS1 = 'EXTR' (TABSOL_I. I .'H_GAMMA') ORDO;
 
    EV_CO1 = EVOL MANU 'DMAX' LLGAM1 'G/G0' LLGG01;
    EV_CO2 = EVOL MANU 'DMAX' LLGAM2 'EPS' LLEPS1;
 
 
    GGO_I_N. I = 'EXCO' (VARI (MOD_I. I) 
                 (LOG((TABPAR.'CHI')*((DEF_I. I)))) EV_CO1 'GRAVITE')
                  'SCAL' 'G_GO';
 
    KSI_I_N. I = 'EXCO' (VARI (MOD_I. I) 
                 (LOG((TABPAR.'CHI')*((DEF_I. I)))) EV_CO2 'GRAVITE')
                  'SCAL' 'KS';
 
    ESO_I_N. I = (GGO_I_N. I)*(TABSOL_I. I .'MODULE_E');
    MAT_I_N. I = 'MATE' (MOD_I. I) 
               'YOUN' (ESO_I_N. I) 'NU' (TABSOL_I. I .'POISSON')
                           'RHO' (TABSOL_I. I .'MASSE_VOLUMIQUE');
    RIG_I_N. I = 'RIGI' (MOD_I.I) (MAT_I_N.I) ;
 
    MAT_PRO = 'MATE' (MOD_I. I) 'YOUN' 
                 (((ESO_I_N. I)*(KSI_I_N. I))*(BETA))
                  'NU' (TABSOL_I. I .'POISSON')
             'RHO' 
             (((TABSOL_I. I .'MASSE_VOLUMIQUE')*(KSI_I_N. I))*(ALPHA));
    AMO_I_N. I  = (MASS (MOD_I.I) MAT_PRO) ET (RIGI (MOD_I.I) MAT_PRO);
    I = I + 1;
  'FIN' IND1;
'FINSI';
*
'SI' (EXIST (TDON1.MOT33.'SOL') 'CONSTANTE');
*  NN1 = 'DIME' DEF_C;
  I = 1;
  'REPE' IND1 (NZONE_C);
    VGG01 = 'EXTR' (GGO_C. I) 'G_GO' 1 1 1;
    VEPS1 = 'EXTR' (KSI_C. I) 'KS' 1 1 1;
*
    ESO_C_N. I  = VGG01*(TABSOL_C. I .'MODULE_E');
    MAT_C_N. I  = 'MATE' (MOD_C. I) 'YOUN' (ESO_C_N. I)
                    'NU' (TABSOL_C. I .'POISSON')
                    'RHO' (TABSOL_C. I .'MASSE_VOLUMIQUE');
    GGO_C_N. I = 'MANU' 'CHML' MOD_C.I 'G_GO' VGG01;
    KSI_C_N. I = 'MANU' 'CHML' MOD_C.I 'KS' VEPS1;
    RIG_C_N. I = 'RIGI' (MOD_C.I) (MAT_C_N.I) ;
    AMO_C_N. I = (VEPS1) * (((ALPHA) * (MAS_C. I))
      'ET' ((BETA) * (RIG_C_N. I)));
   I = I + 1;
 'FIN' IND1;
'FINSI';
*
**************************************************************
*   Controle convergence
**************************************************************
*
MESS '';
MESS '--------------------------------------------------------------';
MESS '--           Critere de convergence                         --';
MESS '--------------------------------------------------------------';
*'SI' ( 'NON' ( 'EXIS' TOL_EPS ) ) ;
*  TOL_EPS = 0.05 ;
*'FINSI' ;
TOL_EPS = TABPAR.'CRITERE';
*
K = 0;
'SI' ('EXIST' (TDON1.MOT33.'SOL') 'ITERATION') ;
*
  I = 1;
  'REPE' IND1 (NZONE_I);
     L_ELE1 = (TABSOL_I.I.MAILLAGE) 'ELEM' 'TYPE';
     J = 1;
     'REPE' IND2 (DIME L_ELE1);
       JJ1 = 1;
       MAI1 = ((TABSOL_I.I.'MAILLAGE') 'ELEM' ('EXTR'L_ELE1 J));
       'REPE' IND3 (NBEL MAI1);
         ELE1 = MAI1 ELEM JJ1;
         VAL0 = EXTR  (REDU (GGO_I. I) ELE1) 'G_GO' 1 1 1 ;
         VAL1 = EXTR  (REDU (GGO_I_N. I) ELE1) 'G_GO' 1 1 1 ;
         DELTA = ( 'ABS' ( VAL1 - VAL0 ) ) / VAL1 ;
         'SI' ( DELTA >EG 0.05 ) ;
           K = K + 1 ;
         'FINSI' ;
         JJ1 = JJ1 + 1;
       'FIN' IND3;
       J = J + 1;
     'FIN' IND2;
*    J = 1;
*    'REPE' IND2 ('DIME' (DEF_I. I ));
*      VAL0 = EXTR  (GGO_I. I . J) 'G_GO' 1 1 1 ;
*      VAL1 = EXTR  (GGO_I_N. I . J) 'G_GO' 1 1 1 ;
*      DELTA = ( 'ABS' ( VAL1 - VAL0 ) ) / VAL1 ;
*      'SI' ( DELTA >EG 0.05 ) ;
*        K = K + 1 ;
*      'FINSI' ;
*      J = J + 1;
*    'FIN' IND2;
    I = I + 1;
  'FIN' IND1;
*
  'MESS' '';
  'MESS' 
       '--------------------------------------------------------------';
  TIT11 = 'CHAINE' '--           Non convergence sur '
                 K ' elements                --';
  'MESS' TIT11;
  'MESS' 
       '--------------------------------------------------------------';
'FINSI';
*
**************************************************************
*   Assemblage
**************************************************************
*
'MESS' '';
'MESS' '--------------------------------------------------------------';
'MESS' '--                         Assemblage                       --';
'MESS' '--------------------------------------------------------------';
I = 1;
'SI' ('EGA' ('TYPE' TABSOL_I) 'TABLE');
  'REPE' BLOC3 NZONE_I;
*    NELEM_I = (DIME SOL_I.I);
*
* VECCHIA VERSIONE
*    J = 1;
*    'REPE' BLOC31 NELEM_I ;
*      'SI' ( ( I 'EGA' 1 ) 'ET' ( J 'EGA' 1 ) ) ;
*        'SI' ('EXIST' (SOL_I. I) J);
*          GGO_SI_N = GGO_I_N. I . J ;
*          MAT_SI_N = MAT_I_N. I . J ;
*          RIG_SI_N = RIG_I_N. I . J ;
*          KSI_SI_N = KSI_I_N. I . J ;
*          AMO_SI_N = AMO_I_N. I . J ;
*        'FINSI';
*      'SINON';
*        'SI' ('EXIST' (SOL_I. I) J);
*          GGO_SI_N = GGO_SI_N 'ET' (GGO_I_N. I . J);
*          MAT_SI_N = MAT_SI_N 'ET' (MAT_I_N. I . J);
*          RIG_SI_N = RIG_SI_N 'ET' (RIG_I_N. I . J);
*          KSI_SI_N = KSI_SI_N 'ET' (KSI_I_N. I . J);
*          AMO_SI_N = AMO_SI_N 'ET' (AMO_I_N. I . J);
*        'FINSI';
*      'FINSI';
*      J = J + 1;
*    'FIN' BLOC31;
*
* NUOVA VERSIONE
    'SI' ( ( I 'EGA' 1 )) ;
      GGO_SI_N = GGO_I_N. I;
      MAT_SI_N = MAT_I_N. I;
      RIG_SI_N = RIG_I_N. I;
      KSI_SI_N = KSI_I_N. I;
      AMO_SI_N = AMO_I_N. I;
    'SINON';
      GGO_SI_N = GGO_SI_N 'ET' (GGO_I_N. I);
      MAT_SI_N = MAT_SI_N 'ET' (MAT_I_N. I);
      RIG_SI_N = RIG_SI_N 'ET' (RIG_I_N. I);
      KSI_SI_N = KSI_SI_N 'ET' (KSI_I_N. I);
      AMO_SI_N = AMO_SI_N 'ET' (AMO_I_N. I);
    'FINSI';
    I = I + 1;
  'FIN' BLOC3;
'FINSI';
*
'SI' ('EGA' ('TYPE' TABSOL_C) 'TABLE');
'SI' ('EGA' ('TYPE' TABCOL) 'TABLE');
  I = 1;
  'REPE' BLOC3 NZONEC;
    NELEMC = (DIME SOLCOL.I);
    J = 1;
    'REPE' BLOC31 NELEMC ;
      'SI' ( ( I 'EGA' 1 ) 'ET' ( J 'EGA' 1 ) ) ;
        'SI' ('EXIST' (SOLCOL. I) J);
          GGOSC_N = GGOCOL_N. I . J ;
          MATSC_N = MATCOL_N. I . J ;
          RIGSC_N = RIGCOL_N. I . J ;
          KSISC_N = KSICOL_N. I . J ;
          AMOSC_N = AMOCOL_N. I . J ;
        'FINSI';
      'SINON';
        'SI' ('EXIST' (SOLCOL. I) J);
          GGOSC_N = GGOSC_N 'ET' (GGOCOL_N. I . J);
          MATSC_N = MATSC_N 'ET' (MATCOL_N. I . J);
          RIGSC_N = RIGSC_N 'ET' (RIGCOL_N. I . J);
          KSISC_N = KSISC_N 'ET' (KSICOL_N. I . J);
          AMOSC_N = AMOSC_N 'ET' (AMOCOL_N. I . J);
        'FINSI';
      'FINSI';
      J = J + 1;
    'FIN' BLOC31;
    I = I + 1;
  'FIN' BLOC3;
*
  I = 1;
  'REPE' BLOC3 NZONE_C;
    NELEM_C = (DIME SOL_C.I);
    J = 1;
    'REPE' BLOC31 NELEM_C ;
      'SI' ( ( I 'EGA' 1 ) 'ET' ( J 'EGA' 1 ) ) ;
        'SI' ('EXIST' (SOL_C. I) J);
          GGO_SC_N = GGO_C_N. I . J ;
          MAT_SC_N = MAT_C_N. I . J ;
          RIG_SC_N = RIG_C_N. I . J ;
          KSI_SC_N = KSI_C_N. I . J ;
          AMO_SC_N = AMO_C_N. I . J ;
        'FINSI';
      'SINON';
        'SI' ('EXIST' (SOL_C. I) J);
          GGO_SC_N = GGO_SC_N 'ET' (GGO_C_N. I . J);
          MAT_SC_N = MAT_SC_N 'ET' (MAT_C_N. I . J);
          RIG_SC_N = RIG_SC_N 'ET' (RIG_C_N. I . J);
          KSI_SC_N = KSI_SC_N 'ET' (KSI_C_N. I . J);
          AMO_SC_N = AMO_SC_N 'ET' (AMO_C_N. I . J);
        'FINSI';
      'FINSI';
      J = J + 1;
    'FIN' BLOC31;
    I = I + 1;
  'FIN' BLOC3;
'SINON';
  I = 1;
  'REPE' BLOC3 NZONE_C;
      'SI' ( ( I 'EGA' 1 )) ;
          GGO_SC_N = GGO_C_N. I ;
          MAT_SC_N = MAT_C_N. I ;
          RIG_SC_N = RIG_C_N. I ;
          KSI_SC_N = KSI_C_N. I ;
          AMO_SC_N = AMO_C_N. I ;
      'SINON';
          GGO_SC_N = GGO_SC_N 'ET' (GGO_C_N. I );
          MAT_SC_N = MAT_SC_N 'ET' (MAT_C_N. I );
          RIG_SC_N = RIG_SC_N 'ET' (RIG_C_N. I );
          KSI_SC_N = KSI_SC_N 'ET' (KSI_C_N. I );
          AMO_SC_N = AMO_SC_N 'ET' (AMO_C_N. I );
      'FINSI';
    I = I + 1;
  'FIN' BLOC3;
'FINSI';
'FINSI';
*
'SI' ('EGA' ('TYPE' TABSOL_I) 'TABLE');
  'ELIM' MAI_SOLI 0.001 ;
  DEF_SOLT = DEF_SOLI;
  MAI_SOLT = MAI_SOLI;
  MOD_SOLT = MOD_SOLI;
  GGO_ST_N = GGO_SI_N;
  MAT_ST_N = MAT_SI_N;
  RIG_ST_N = RIG_SI_N;
  MAS_SOLT = MAS_SOLI;
  KSI_ST_N = KSI_SI_N;
  AMO_ST_N = AMO_SI_N;
  'SI' ('EGA' ('TYPE' TABSOL_C) 'TABLE');
    'SI' ('EGA' ('TYPE' TABCOL) 'TABLE');
      DEF_SOLT = DEF_SOLC 'ET' DEF_SOLI;
      MAI_SOLT = MAI_SOLI 'ET' MAI_SOLC;
      MOD_SOLT = MOD_SOLI 'ET' MOD_SOLC;
      MAS_SOLT = MAS_SOLI 'ET' MAS_SOLC;
      GGO_ST_N = GGO_ST_N 'ET' GGO_SOLC;
      MAT_ST_N = MAT_ST_N 'ET' MAT_SOLC;
      RIG_ST_N = RIG_ST_N 'ET' RIG_SOLC;
      KSI_ST_N = KSI_ST_N 'ET' KSI_SOLC;
      AMO_ST_N = AMO_ST_N 'ET' AMO_SOLC;
*
      GGOST_C = GGOSC_N;
      MATST_C = MATSC_N;
      RIGST_C = RIGSC_N;
      KSIST_C = KSISC_N;
      AMOST_C = AMOSC_N;
    'SINON';
      DEF_SOLT = DEF_SOLC 'ET' DEF_SOLI;
      MAI_SOLT = MAI_SOLI 'ET' MAI_SOLC;
      MOD_SOLT = MOD_SOLI 'ET' MOD_SOLC;
      GGO_ST_N = GGO_ST_N 'ET' GGO_SOLC;
      MAT_ST_N = MAT_ST_N 'ET' MAT_SOLC;
      RIG_ST_N = RIG_ST_N 'ET' RIG_SOLC;
      MAS_SOLT = MAS_SOLI 'ET' MAS_SOLC;
      KSI_ST_N = KSI_ST_N 'ET' KSI_SOLC;
      AMO_ST_N = AMO_ST_N 'ET' AMO_SOLC;
    'FINSI';
  'FINSI';
'SINON';
  'SI' ('EGA' ('TYPE' TABSOL_C) 'TABLE');
    'SI' ('EGA' ('TYPE' TABCOL) 'TABLE');
      DEF_SOLT = DEF_SOLC;
      MAI_SOLT = MAI_SOLC;
      MOD_SOLT = MOD_SOLC;
      MAS_SOLT = MAS_SOLC;
      GGO_ST_N = GGO_SOLC;
      MAT_ST_N = MAT_SOLC;
      RIG_ST_N = RIG_SOLC;
      KSI_ST_N = KSI_SOLC;
      AMO_ST_N = AMO_SOLC;
*
      GGOST_C = GGOSC_N;
      MATST_C = MATSC_N;
      RIGST_C = RIGSC_N;
      KSIST_C = KSISC_N;
      AMOST_C = AMOSC_N;
    'SINON';
      DEF_SOLT = DEF_SOLC;
      GGO_ST_N = GGO_SOLC;
      MAT_ST_N = MAT_SOLC;
      RIG_ST_N = RIG_SOLC;
      KSI_ST_N = KSI_SOLC;
      AMO_ST_N = AMO_SOLC;
    'FINSI';
  'FINSI';
'FINSI';
*
ELIM 0.001 MAI_SOLT;
*
**************************************************************
*   Sauvegarde
**************************************************************
*
MESS '';
MESS '--------------------------------------------------------------';
MESS '--                         Stockage                         --';
MESS '--------------------------------------------------------------';
*
'SI' ('EGA' ('TYPE' TABSOL_C) 'TABLE');
  'SI' ('EGA' ('TYPE' TABCOL) 'TABLE');
    TABRC_N = TABLE ;
    TABRC_N.'T_SOL' = SOLCOL ;
    TABRC_N.'T_ESO' = ESOCOL ;
    TABRC_N.'T_MOD' = MODCOL ;
    TABRC_N.'T_GGO' = GGOCOL_N ;
    TABRC_N.'T_MAT' = MATCOL_N;
    TABRC_N.'T_RIG' = RIGCOL_N ;
    TABRC_N.'T_MAS' = MASCOL ;
    TABRC_N.'T_KSI' = KSICOL_N ;
    TABRC_N.'T_AMO' = AMOCOL_N ;
    TABRC_N.'T_DIS' = DEFC ;
    TABRC_N.'O_SOL' = MAISOLC ;
    TABRC_N.'O_MOD' = MODSOLC ;
    TABRC_N.'O_GGO' = GGOST_C ;
    TABRC_N.'O_MAT' = MATST_C ;
    TABRC_N.'O_RIG' = RIGST_C ;
    TABRC_N.'O_MAS' = MASSOLC ;
    TABRC_N.'O_KSI' = KSIST_C ;
    TABRC_N.'O_AMO' = AMOST_C ;
    TABRC_N.'O_DIS' = DEFCOL ;
*
    TABR_C_N = TABLE ;
    TABR_C_N.'T_SOL' = SOL_C ;
    TABR_C_N.'T_ESO' = ESO_C_N ;
    TABR_C_N.'T_MOD' = MOD_C ;
    TABR_C_N.'T_GGO' = GGO_C_N ;
    TABR_C_N.'T_MAT' = MAT_C_N ;
    TABR_C_N.'T_RIG' = RIG_C_N ;
    TABR_C_N.'T_MAS' = MAS_C ;
    TABR_C_N.'T_KSI' = KSI_C_N ;
    TABR_C_N.'T_AMO' = AMO_C_N ;
    TABR_C_N.'T_DIS' = DEF_C ;
    TABR_C_N.'O_SOL' = MAI_SOLC ;
    TABR_C_N.'O_MOD' = MOD_SOLC ;
    TABR_C_N.'O_GGO' = GGO_SC_N ;
    TABR_C_N.'O_MAT' = MAT_SC_N ;
    TABR_C_N.'O_RIG' = RIG_SC_N ;
    TABR_C_N.'O_MAS' = MAS_SOLC ;
    TABR_C_N.'O_KSI' = KSI_SC_N ;
    TABR_C_N.'O_AMO' = AMO_SC_N ;
    TABR_C_N.'O_DIS' = DEF_SOLC ;
  'SINON';
    TABR_C_N = TABLE ;
    TABR_C_N.'T_SOL' = SOL_C ;
    TABR_C_N.'T_ESO' = ESO_C_N ;
    TABR_C_N.'T_MOD' = MOD_C ;
    TABR_C_N.'T_GGO' = GGO_C_N ;
    TABR_C_N.'T_MAT' = MAT_C_N ;
    TABR_C_N.'T_RIG' = RIG_C_N ;
    TABR_C_N.'T_MAS' = MAS_C ;
    TABR_C_N.'T_KSI' = KSI_C_N ;
    TABR_C_N.'T_AMO' = AMO_C_N ;
    TABR_C_N.'T_DIS' = DEF_C ;
    TABR_C_N.'O_SOL' = MAI_SOLC ;
    TABR_C_N.'O_MOD' = MOD_SOLC ;
    TABR_C_N.'O_GGO' = GGO_SC_N ;
    TABR_C_N.'O_MAT' = MAT_SC_N ;
    TABR_C_N.'O_RIG' = RIG_SC_N ;
    TABR_C_N.'O_MAS' = MAS_SOLC ;
    TABR_C_N.'O_KSI' = KSI_SC_N ;
    TABR_C_N.'O_AMO' = AMO_SC_N ;
    TABR_C_N.'O_DIS' = DEF_SOLC ;
  'FINSI';
'FINSI';
'SI' ('EGA' ('TYPE' TABSOL_I) 'TABLE');
  TABR_I_N = TABLE ;
  TABR_I_N.'T_SOL' = SOL_I ;
  TABR_I_N.'T_ESO' = ESO_I_N ;
  TABR_I_N.'T_MOD' = MOD_I ;
  TABR_I_N.'T_GGO' = GGO_I_N ;
  TABR_I_N.'T_MAT' = MAT_I_N ;
  TABR_I_N.'T_RIG' = RIG_I_N ;
  TABR_I_N.'T_MAS' = MAS_I ;
  TABR_I_N.'T_KSI' = KSI_I_N ;
  TABR_I_N.'T_AMO' = AMO_I_N ;
  TABR_I_N.'T_DIS' = DEF_I ;
  TABR_I_N.'O_SOL' = MAI_SOLI ;
  TABR_I_N.'O_MOD' = MOD_SOLI ;
  TABR_I_N.'O_GGO' = GGO_SI_N ;
  TABR_I_N.'O_MAT' = MAT_SI_N ;
  TABR_I_N.'O_RIG' = RIG_SI_N ;
  TABR_I_N.'O_MAS' = MAS_SOLI ;
  TABR_I_N.'O_KSI' = KSI_SI_N ;
  TABR_I_N.'O_AMO' = AMO_SI_N ;
  TABR_I_N.'O_DIS' = DEF_SOLI ;
'FINSI';
*
TABR_T_N = TABLE;
TABR_T_N.'O_SOL' = MAI_SOLT ;
TABR_T_N.'O_MOD' = MOD_SOLT ;
TABR_T_N.'O_GGO' = GGO_ST_N ;
TABR_T_N.'O_MAT' = MAT_ST_N ;
TABR_T_N.'O_RIG' = RIG_ST_N ;
TABR_T_N.'O_MAS' = MAS_SOLT ;
TABR_T_N.'O_KSI' = KSI_ST_N ;
TABR_T_N.'O_AMO' = AMO_ST_N ;
TABR_T_N.'O_DIS' = DEF_SOLT ;
*
TDON1.MOT22 .'SOL' = TABLE;
'SI' ('EGA' ('TYPE' TABSOL_I) 'TABLE');
  TDON1.MOT22 .'SOL'.'ITERATION' = TABR_I_N;
'FINSI';
*
'SI' ('EGA' ('TYPE' TABSOL_C) 'TABLE');
  'SI' ('EGA' ('TYPE' TABCOL) 'TABLE');
    TDON1.MOT22.'COLONNE' = TABRC_N;
    TDON1.MOT22.'SOL'.'CONSTANTE' = TABR_C_N;
  'SINON';
    TDON1.MOT22.'SOL'.'CONSTANTE' = TABR_C_N;
  'FINSI';
'FINSI';
*
TDON1.MOT22.'SOL'.'TOTALE' = TABR_T_N;
'SI' ('EXIST' (TDON1) 'STRUCTURE');
  TDON1.MOT22.'STRUCTURE' = TABSTR;
'FINSI';
*
'TRAC' (TABR_T_N.'O_GGO') (TABR_T_N.'O_MOD')
  'TITR' ('CHAINE' 'Non Convergence sur ' K ' elements') 'NCLK';
'TRAC' (TABR_T_N.'O_KSI') (TABR_T_N.'O_MOD')
  'TITR' ('CHAINE' 'Non Convergence sur ' K ' elements') 'NCLK';
*
'OUBL' I;
'OUBL' J;
'OUBL' ZZ;
'OUBL' MOD_C;
'OUBL' MODC;
'OUBL' MAT_C;
'OUBL' MATC;
'OUBL' RIG_C;
'OUBL' RIGC;
'OUBL' MAS_C;
'OUBL' MASC;
'OUBL' AMO_C;
'OUBL' AMOC;
'OUBL' KSI_C;
'OUBL' KSIC;
'OUBL' SOL_C;
'OUBL' SOLC;
'OUBL' ESO_C;
'OUBL' ESOC;
'OUBL' GGO_C;
'OUBL' GGOC;
'OUBL' MOD_I;
'OUBL' MAT_I;
'OUBL' RIG_I;
'OUBL' MAS_I;
'OUBL' AMO_I;
'OUBL' KSI_I;
'OUBL' SOL_I;
'OUBL' ESO_I;
'OUBL' GGO_I;
'OUBL' MAI_SOLC;
'OUBL' MAT_SOLC;
'OUBL' GGO_SOLC;
'OUBL' KSI_SOLC;
'OUBL' MOD_SOLC;
'OUBL' RIG_SOLC;
'OUBL' MAS_SOLC;
'OUBL' AMO_SOLC;
'OUBL' MAI_SOLI;
'OUBL' MAT_SOLI;
'OUBL' GGO_SOLI;
'OUBL' KSI_SOLI;
'OUBL' MOD_SOLI;
'OUBL' RIG_SOLI;
'OUBL' MAS_SOLI;
'OUBL' AMO_SOLI;
'OUBL' MAISOLC;
'OUBL' MATSOLC;
'OUBL' GGOSOLC;
'OUBL' KSISOLC;
'OUBL' MODSOLC;
'OUBL' RIGSOLC;
'OUBL' MASSOLC;
'OUBL' AMOSOLC;
'OUBL' MAI_SOLT;
'OUBL' MAT_SOLT;
'OUBL' GGO_SOLT;
'OUBL' KSI_SOLT;
'OUBL' MOD_SOLT;
'OUBL' RIG_SOLT;
'OUBL' MAS_SOLT;
'OUBL' AMO_SOLT;
'OUBL' MAISOLT;
'OUBL' MATSOLT;
'OUBL' GGOSOLT;
'OUBL' KSISOLT;
'OUBL' MODSOLT;
'OUBL' RIGSOLT;
'OUBL' MASSOLT;
'OUBL' AMOSOLT;
'OUBL' TABSOL_C;
'OUBL' TABCOL;
'OUBL' TABSOL_I;
'OUBL' TABPAR;
'OUBL' TABSOL;
'OUBL' TABSTR;
'OUBL' TABFRO;
'OUBL' H_SAND;
'OUBL' G_SAND;
'OUBL' H_CLAY;
'OUBL' G_CLAY;
'OUBL' H_SAND;
'OUBL' EH_SAND;
'OUBL' EG_SAND;
'OUBL' EH_CLAY;
'OUBL' EG_CLAY;
'OUBL' EH_SAND;
'OUBL' MAIL1C;
'OUBL' GAM_EFF;
'OUBL' TABRES_T;
'OUBL' TABRESC;
'OUBL' TABRES_C;
'OUBL' TABRES_I;
'OUBL' ZONE_C;
'OUBL' ZONE_I;
'OUBL' ZONEC;
'OUBL' NN1;
'OUBL' NELEM_C;
'OUBL' NELEMC;
'OUBL' NELEM1C;
'OUBL' NELEM1_I;
'OUBL' NELEM_I;
'OUBL' NZONE_C;
'OUBL' NZONE_I;
'OUBL' NZONEC;
'OUBL' F1;
'OUBL' F2;
'OUBL' FC;
'OUBL' ALPHA;
'OUBL' BETA;
'OUBL' MAIL1_I;
'OUBL' GGO_SI_N;
'OUBL' KSI_SI_N;
'OUBL' MAT_SI_N;
'OUBL' KSI_SOLI;
'OUBL' GGO_SOLI;
'OUBL' TABR_T_N;
'OUBL' TABR_I_N;
'OUBL' GGO_I_N;
'OUBL' AMO_I_N;
'OUBL' RIG_I_N;
'OUBL' KSI_I_N;
'OUBL' MAT_I_N;
'OUBL' ESO_I_N;
'OUBL' DEF_I;
'OUBL' TRESUX;
'OUBL' TAB_DYX;
'OUBL' TABRXDV;
'OUBL' TABDV;
'OUBL' LLEPS1;
'OUBL' LLGG01;
'OUBL' LLGAM2;
'OUBL' LLGAM1;
'OUBL' L_GG03;
'OUBL' L_GG02;
'OUBL' L_GG01;
*
*************************************************************
***************     FIN  STEP SUP A 1    ********************
*************************************************************
*
'SINON';
*************************************************************
********************     STEP0      *************************
*************************************************************
*
MESS '++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++';
MESS '--                            STEP 0                        --';
MESS '--------------------------------------------------------------';
MESS '';
MESS '--------------------------------------------------------------';
MESS '--                Controle des donnees d entree             --';
MESS '--------------------------------------------------------------';
MESS '';
*
*************  controle general *****************************
'SI' ('EXIST' (TDON1) 'SOL');
  'SI' ('NEG' (TYPE (TDON1.'SOL')) 'TABLE');
    'MESS' ;
    'MESS' '   Erruer....l indice de la table des donnees';
    'MESS' '   SOL doit etre du type TABLE';
    'QUIT' ISSLEQ;
  'SINON';
    TABSOL = TDON1.'SOL';
  'FINSI';
'SINON';
  'MESS' ;
  'MESS' '   Erruer....Il manque la table du sol';
  'MESS' '   dans la table des donn�es';
  'QUIT' ISSLEQ;
'FINSI';
'SI' ('EXIST' (TDON1) 'PARAMETRES');
  'SI' ('NEG' (TYPE (TDON1.'PARAMETRES')) 'TABLE');
    'MESS' ;
    'MESS' '   Erruer....l indice de la table des donnees';
    'MESS' '   PARAMETRE doit etre du type TABLE';
    'QUIT' ISSLEQ;
  'SINON';
    TABPAR = TDON1.'PARAMETRES';
  'FINSI';
'SINON';
  'MESS' ;
  'MESS' '   Erruer....Il manque la table des parametres';
  'MESS' '   dans la table des donn�es';
  'QUIT' ISSLEQ;
'FINSI';
V_STR = FAUX;
'SI' ('EXIST' (TDON1) 'STRUCTURE');
  'SI' ('NEG' (TYPE (TDON1.'STRUCTURE')) 'TABLE');
    'MESS' ;
    'MESS' '   Erruer....l indice de la table des donnees';
    'MESS' '   STRUCTURE doit etre du type TABLE';
    'QUIT' ISSLEQ;
  'SINON';
    TABSTR = TDON1.'STRUCTURE';
    V_STR = VRAI;
  'FINSI';
'SINON';
  'MESS' ;
  'MESS' '   Attention...calcul sans la structure';
'FINSI';
*
*************  controle parametre *****************************
'SI' ('EXIST' (TABPAR) 'F1');
  'SI' ('NEG' ('TYPE' (TABPAR.'F1')) 'FLOTTANT');
    'MESS' ;
    'MESS' '   Erruer....F1 n est pas un flottant'
    'MESS' '   dans la table des parametres';
    'QUIT' ISSLEQ;
  FINSI;
'SINON';
  'MESS' ;
  'MESS' '   Erruer....il manque F1 dans la table des parametres';
  'QUIT' ISSLEQ;
'FINSI';
'SI' ('EXIST' (TABPAR) 'F2');
  'SI' ('NEG' ('TYPE' (TABPAR.'F2')) 'FLOTTANT');
    'MESS' ;
    'MESS' '   Erruer....F2 n est pas un flottant'
    'MESS' '   dans la table des parametres';
    'QUIT' ISSLEQ;
  FINSI;
'SINON';
  'MESS' ;
  'MESS' '   Erruer....il manque F2 dans la table des parametres';
  'QUIT' ISSLEQ;
'FINSI';
'SI' ('EXIST' (TABPAR) 'FC');
  'SI' ('NEG' ('TYPE' (TABPAR.'FC')) 'FLOTTANT');
    'MESS' ;
    'MESS' '   Erruer....FC n est pas un flottant'
    'MESS' '   dans la table des parametres';
    'QUIT' ISSLEQ;
  FINSI;
'SINON';
  'MESS' ;
  'MESS' '   Erruer....il manque FC dans la table des parametres';
  'QUIT' ISSLEQ;
'FINSI';
'SI' ('EXIST' (TABPAR) 'TYPE');
  'SI' (('NEG' TABPAR.'TYPE' 'LYSMER') ET 
         ('NEG' TABPAR.'TYPE' 'WHITE'));
     'MESS' ;
     'MESS' '   Erruer....type de frontiere mauvaise';
     'QUIT' ISSLEQ;
  'FINSI';
'SINON';
  TABPAR.'TYPE' = CHAINE 'LYSMER';
'FINSI';
*'SI' ('EXIST' (TABPAR) 'REPONSE');
*  I = 1;
*  'REPE' IND1 ('DIME' (TABPAR.'REPONSE'));
*    'SI' ('NEG' (TYPE (TABPAR.'REPONSE'. I .'NOM')) 'MOT');
*      'MESS' ;
*      'MESS' '   Erruer dans la tables des parametres';
*      'MESS' '   dans la table des reponse a l indice ' I;
*      'MESS' '   le sous-indice NOM doit etre un mot';
*      'QUIT' ISSLEQ; 
*    'FINSI';
*    'SI' ('NEG' (TYPE (TABPAR.'REPONSE'. I .'NOEUD')) 'POINT');
*      'MESS' ;
*      'MESS' '   Erruer dans la tables des parametres';
*      'MESS' '   dans la table des reponse a l indice ' I;
*      'MESS' '   le sous-indice NOEUD doit etre un point';
*      'QUIT' ISSLEQ; 
*    'FINSI';
*    'SI' ('NEG' (TYPE (TABPAR.'REPONSE'. I .'COMP')) 'MOT');
*      'MESS' ;
*      'MESS' '   Erruer dans la tables des parametres';
*      'MESS' '   dans la table des reponse a l indice ' I;
*      'MESS' '   le sous-indice COMP doit etre un mot';
*      'QUIT' ISSLEQ; 
*    'FINSI';
*    I = I + 1;
*  'FIN' IND1;
*'SINON';
*  'MESS' ;
*  'MESS' '   Erruer....il manque la table des point de reponse';
*  'QUIT' ISSLEQ;
*'FINSI';
'SI' ('EXIST' (TABPAR) 'TYPE_CALCUL');
  'SI' (('NEG' (TABPAR.'TYPE_CALCUL') (CHAINE 'ISS_COMPLET')) ET
        ('NEG' (TABPAR.'TYPE_CALCUL') 'ISS_CHAMPLOIN') ET
        ('NEG' (TABPAR.'TYPE_CALCUL') (CHAINE 'ISS_COMP_SIMP')));
    'MESS' ;
    'MESS' '   Erruer....type de calcul non prevu';
    'QUIT' ISSLEQ;
  'FINSI';
'SINON';
  TABPAR.'TYPE_CALCUL' = CHAINE 'ISS_COMPLET';
'FINSI';
T_CAL = TABPAR.'TYPE_CALCUL';
'SI' ('EXIST' (TABPAR) 'CRITERE');
  'SI' ('NEG' ('TYPE' (TABPAR.'CRITERE')) 'FLOTTANT');
    'MESS';
    'MESS' '   Erruer....CRITERE n est pas un flottant'
    'MESS' '   dans la table des parametres';
    'QUIT' ISSLEQ;
  'FINSI';
'SINON';
  'MESS';
  'MESS' '   CRITERE fixe à 0.05'
  'MESS' '   dans la table des parametres';
  TABPAR.'CRITERE' = 0.05;
'FINSI';
*
*************  controle sol **********************************
'SI' (('VALEUR' 'DIME') 'EGA' 2);
'OPTI' 'ELEM' SEG3;
  F_BASE = 'VIDE' MAILLAGE;
  F_LATE = 'VIDE' MAILLAGE;
'OPTI' 'ELEM' QUA8;
'SINON';
'OPTI' 'ELEM' QUA4;
  F_BASE = 'VIDE' MAILLAGE;
  F_LATE = 'VIDE' MAILLAGE;
'OPTI' 'ELEM' CUB8;
'FINSI';
*
* serie de checks pour le calcul complet
'SI' ('EGA' T_CAL 'ISS_COMPLET');
  'SI' ('EXIST' TABSOL 'CONSTANTE');
    'MESS' ;
    'MESS' '   Erreur....pour le calcul COMPLET';
    'MESS' '   il n y a pas besoin de la partie du sol constante';
    'QUIT' ISSLEQ;
  'FINSI';
  'SI' ('EXIST' TABSOL 'ITERATION');
    'SI' ('NEG'