Télécharger la plaquette Cast3M


fr en

Télécharger pas_resu.procedur

Retour à la liste

Numérotation des lignes : 

  1. * PAS_RESU  PROCEDUR  JK148537  24/04/19    21:15:03     11906          
  2. 'DEBP' PAS_RESU TEMP0*'FLOTTANT' IAUTO*'ENTIER' PRECED*'TABLE' ;
  3.  
  4. WTAB  = PRECED.'WTABLE';
  5. *
  6. * Recuperation de la table 'ESTIMATION' des champs en fin de pas :
  7. estim = PRECED.'ESTIMATION' ;
  8. *
  9. * Remplissage de la table de 'CONTINUATION' a partir de 'ESTIMATION'
  10. conti= PRECED.'CONTINUATION';
  11. *
  12. 'SI' ('NEG' estim.'TEMPS' TEMP0 ('VALE' 'PREC'));
  13.   'MESS' 'Temps differents dans PAS_RESU ?' ;
  14.   estim.'TEMPS' = TEMP0 ;
  15. 'FINS';
  16.  
  17. * Transfert des indices de la table "estim" dans "conti"
  18. DIM1   = 'DIME' estim ;
  19. INDEX1 = 'INDE' estim ;
  20.  
  21. 'REPE' SURi DIM1 ;
  22.   ii     = &SURi;
  23.   TYP1   ='TYPE' INDEX1.ii;
  24.   'SI' ('NEG' TYP1 'MOT');
  25.     MOT_tmp = INDEX1.ii ;
  26.   'SINO';
  27.     MOT_tmp ='MOT' INDEX1.ii;
  28.   'FINS';
  29.   conti.MOT_tmp = estim.MOT_tmp ;
  30. 'FIN' SURi;
  31.  
  32. ************************************************************************
  33. *************** Sauvegarde eventuelle des resultats ********************
  34. ************************************************************************
  35. *La sauvegarde n'est faite que si on a eu convergence des calculs :
  36.    ISAUV = FAUX;
  37.   'SI' WTAB.'CONV' ;
  38.      ISAUV ='DANS' WTAB.'IPRESU' TEMP0 ;
  39.  
  40.     'SI' ('MULT' IAUTO WTAB.'NAURES' ) ;
  41.        ISAUV = VRAI ;
  42.     'FINSI' ;
  43.   'FINSI' ;
  44.  
  45.   'SI' ISAUV ;
  46.      NSOR        = WTAB.'NSOR' + 1 ;
  47.      WTAB.'NSOR' = NSOR            ;
  48.  
  49. *    Cas particuliers de 'DEFIN' et 'DEFTO'
  50.      TABSAUV = WTAB.'TABSAUV' ;
  51.     'SI' ('NON' TABSAUV.'DEFTO');
  52.       'SI' ('EXIS' estim 'DEFORMATIONS');
  53.          estim ='ENLE' estim 'DEFORMATIONS';
  54.       'FINS';
  55.     'FINSI' ;
  56.  
  57.     'SI' ('NON' TABSAUV.'DEFIN') ;
  58.       'SI' ('EXIS' estim 'DEFORMATIONS_INELASTIQUES');
  59.          estim ='ENLE' estim 'DEFORMATIONS_INELASTIQUES';
  60.       'FINS';
  61.     'FINSI' ;
  62.  
  63.     'SI' ('EXIS' estim 'VITESSES_FLUIDE_0') ;
  64.        estim ='ENLE' estim 'VITESSES_FLUIDE_0';
  65.     'FINSI' ;
  66.  
  67. *    Transfert des indices de la table "estim" dans PRECED
  68.      DIM1  ='DIME' estim;
  69.      INDEX1='INDE' estim;
  70.  
  71.     'REPE' SURi DIM1 ;
  72.        ii     = &SURi;
  73.        TYP1   ='TYPE' INDEX1.ii;
  74.       'SI' ('NEG' TYP1 'MOT');
  75.          MOT_tmp = INDEX1.ii     ;
  76.          OBJi    = estim.MOT_tmp ;
  77.  
  78.       'SINO';
  79.          MOT_tmp ='MOT' INDEX1.ii;
  80.          OBJi    = estim.MOT_tmp ;
  81.         'SI' ('EGA' TYP1 'TABL');
  82.           'SI' ('EXIS' estim.MOT_tmp 'SOUSTYPE');
  83.             'SI' ('EGA' estim.MOT_tmp.'SOUSTYPE' 'ESCLAVE');
  84.                ERRE 21;
  85.                OBJi ='ETG' estim.MOT_tmp ;
  86.             'FINS';
  87.           'FINS';
  88.         'FINS';
  89.       'FINS';
  90.       'SI'('NON' ('EXIS' PRECED MOT_tmp));
  91.          PRECED.MOT_tmp ='TABL' ;
  92.       'FINS';
  93.        PRECED.MOT_tmp.NSOR = estim.MOT_tmp ;
  94.     'FIN' SURi;
  95.  
  96. *-----------------------------------------------------------------------
  97. *---  1- Sauvegarde incrementale sous certaines conditions
  98. *---  2- Fantomisation (option 'ECONOMIQUE'=VRAI de PASAPAS) pour liberer la memoire
  99. *-----------------------------------------------------------------------
  100.     'SI' WTAB.'ISAUVEGL';
  101.        ltpssauv = WTAB.'ISAUVEGA' ;
  102.       'SI' ('DANS' ltpssauv TEMP0) ;
  103.          EPOQ1 ='DATE' 'EPOCH' ;
  104. *        Sauvegarde : Dernier pas OU toutes les 300 secondes OU si DATE ne marche pas
  105. *        Fantomisation si on sauvegarde effectivement !
  106.         'SI' (     ('EGA' TEMP0 ('EXTR' ltpssauv ('DIME' ltpssauv)))
  107.               'OU' ((EPOQ1 '-' WTAB.'EPOQUE') '>EG' 300.D0)
  108.               'OU' (EPOQ1 'EGA' 0)) ;
  109.  
  110.           'SAUV' 'MUET' PRECED ;
  111.            WTAB.'EPOQUE' =EPOQ1 ;
  112.  
  113.           'SI' WTAB.'LECONOMI';
  114.              NFFSOR = NSOR '-' WTAB.'NFSOR' '-' 1 ;
  115.             'SI' ( NFFSOR '>' 0) ;
  116.  
  117. *              Boucle sur les instants SAUVES depuis la derniere fantomisation
  118.               'REPE' DETRES  NFFSOR ;
  119.                  NSORM2 = WTAB.'NFSOR' - 1 + &DETRES ;
  120.                 'SI' (NSORM2 '<' 0) ;
  121.                   'ITER' DETRES ;
  122.                 'FINSI' ;
  123.  
  124. *                Boucle sur les indices a fantomiser
  125.                 'REPE' SURi DIM1 ;
  126.                    ii     = &SURi;
  127.                    TYP1   ='TYPE' INDEX1.ii;
  128.                   'SI' ('NEG' TYP1 'MOT');
  129.                      MOT_tmp = INDEX1.ii     ;
  130.                   'SINO';
  131.                      MOT_tmp ='MOT' INDEX1.ii;
  132.                   'FINS';
  133.                   'FANTOM' PRECED.MOT_tmp NSORM2 ;
  134.                 'FIN' SURi;
  135.               'FIN' DETRES ;
  136.                WTAB.'NFSOR'= WTAB.'NFSOR' + NFFSOR ;
  137.             'FINSI' ;
  138.           'FINS' ;
  139.         'FINSI';
  140.       'FINSI';
  141.     'FINSI';
  142.   'FINSI' ;
  143. 'FINP' ISAUV ;
  144.  
  145.  
  146.  
  147.  
  148.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales