Procedures | Cast3M

   1 : * PAS_ETAT  PROCEDUR  CHAT      13/03/21    21:15:07     7735           
   2 : 'DEBPROC' PAS_ETAT PRECED*'TABLE' TI*'FLOTTANT'   ;
   3 : *
   4 : * cette procedure cree un MCHAML qui contient tout ce que
   5 : * l'on connait pour un instant donne (TI)
   6 : * Les informations sont tirees du chargement et, en plus,
   7 : * on  interpole entre le dernier pas (dans la table estimation)
   8 : * et les dernieres valeurs convergees (dans continuation).
   9 : *
  10 : WTAB=PRECED.'WTABLE';
  11 : *
  12 : IE=FAUX;
  13 : 'SI' ('EXIS' WTAB 'CHARGEMENT'); 
  14 :    IE=VRAI;
  15 :    MCH= 'ADET' 'NOUV' WTAB.'MO_TOT' WTAB.'CHARGEMENT' TI ;
  16 : 'FINS';
  17 : 
  18 : CONTI=PRECED.'CONTINUATION';ESTIM=PRECED.'ESTIMATION';
  19 : NCONT=DIME CONTI;
  20 : 
  21 : 'SI' (NCONT NEG 0);
  22 :    IND= 'INDEX' CONTI;
  23 :    TCONV= CONTI . 'TEMPS';
  24 :    TCONT= ESTIM.'TEMPS';
  25 :    IAD=FAUX; IAC=FAUX;
  26 :    'SI' ( 'NEG' TCONT  TCONV);
  27 :       XPRO= (TI - TCONV)/ (TCONT - TCONV); 
  28 :       UNXPRO= 1.d0 - XPRO;
  29 :       IAD = ('ABS' UNXPRO) '>' 1.E-10 ;
  30 :       IAC = ('ABS' XPRO) '>' 1.E-10 ;
  31 :    'FINSI';   
  32 :    'REPE' BOU NCONT;
  33 :       MO= 'MOT' IND . &BOU;
  34 : *      mess 'tctc mio' mo;
  35 :       II='TEXT' ' ';
  36 :       SI ( 'EGA'   MO 'TEMPS' ) ; II= 'TEXT' MO;'FINSI';
  37 : * pas decide si grave d'avoir les reactions !!!! ni les temperatures
  38 : * ni les vitesses ni...
  39 :       'SI' ( 'EGA' MO 'GEODEB' ) ; 'ITER' BOU; 'FINSI';
  40 :       'SI' ( 'EGA' MO 'GEOFIN' ) ; 'ITER' BOU; 'FINSI';      
  41 :       'SI' ( 'EGA' MO 'REACTIONS' ) ; 'ITER' BOU; 'FINSI';
  42 :       'SI' ( 'EGA' MO 'DEFORMATIONS_INELASTIQUES');'ITER' BOU;'FINSI';
  43 :       'SI' ( 'EGA' MO 'VITESSES');'ITER' BOU;'FINSI'; 
  44 : *      'SI' ( 'EGA' MO 'TEMPERATURES');'ITER' BOU;'FINSI';
  45 :       'SI' ( 'EGA' MO 'DEFORMATIONS');'ITER' BOU;'FINSI';
  46 :       'SI' ( 'EGA' MO 'CARA_ACTU');'ITER' BOU;'FINSI'; 
  47 :       'SI' IAD;
  48 :          VAD = CONTI. MO;
  49 :          'SI' IAC ;
  50 :             VAD = (VAD * UNXPRO) + (ESTIM.MO * XPRO);
  51 :          'FINSI';
  52 :       'SINON' ;
  53 :          VAD = ESTIM. MO;
  54 :       'FINSI';
  55 :       'SI' IE; 
  56 :         'SI' ( 'EGA' MO 'TEMPS' );
  57 :            MCH= MCH 'ADET' WTAB.'MO_TOT'  MO VAD ;
  58 :          'SINON';
  59 :            MCH= MCH 'ADET' WTAB.'MO_TOT'   VAD ;
  60 :          'FINSI';
  61 :       'SINON';
  62 :          IE=VRAI;
  63 :          'SI' ( 'EGA' MO 'TEMPS' );
  64 :            MCH= 'ADET' 'NOUV'   WTAB.'MO_TOT'  MO VAD ;
  65 :          'SINON';
  66 :            MCH=  'ADET' 'NOUV' WTAB.'MO_TOT'   VAD ;
  67 :          'FINSI';
  68 :       'FINSI';
  69 :    'FIN' BOU;
  70 : 'FINSI';
  71 : 
  72 : 'FINP' MCH;
  73 :  
  74 :  
  75 :  
  76 :  
  77 :