Procedures | Cast3M

   1 : * CALACTIV  PROCEDUR  AURI      98/12/29    21:16:01     3419           
   2 : DEBPROC CALACTIV  TB3*TABLE TB4*TABLE  LILI*LISTENTI ;
   3 : **********************************************************************
   4 : *     CALCUL DE L'ACTIVITE DES ESPECES DE LA LISTE LILI
   5 : *      TB3 table issue de CHI1
   6 : *      TB4 table issue de CHI2
   7 : **********************************************************************
   8 :  ACTIV='COAC' TB3 'FORCEION' TB4.FION ;
   9 :  NIDX= DIME(TB3.DESCHI.IDX) ;
  10 :  NIDY= DIME(TB3.DESCHI.IDY) ;
  11 :  NLILI= DIME LILI ;
  12 : *                        initialisation du résultat
  13 :  NUESP= 'EXTR' LILI 1 ;
  14 :  MO1 NO1 = NOESPCHI TB3  NUESP ;
  15 :  CCACT= 'EXCO' MO1 TB4.'SOLU' MO1 ;
  16 :  CCACT= CCACT- CCACT ;
  17 : *
  18 :  'REPE' BOU1 NLILI ;
  19 :  NUESP= 'EXTR' LILI &BOU1 ;
  20 :  MO1 NO1 = NOESPCHI TB3  NUESP ;
  21 :  'SI' ('EGA' ('DIME' ( LILI 'SAUF'('LECT' NUESP)))  NLILI ) ;
  22 :           'MESS' NUESP ' n est pas une espèce en solution'
  23 :            'QUIT' BOU1 ;
  24 :  'FINSI' ;
  25 :  REPE BDY NIDY ;
  26 :  NUNU1='EXTR' (TB3.DESCHI.IDY) &BDY ;
  27 :  'SI' ('EGA' NUNU1 NUESP) ;
  28 :  NUDY=&BDY ;
  29 :  'FINS' ;
  30 :  'FIN' BDY ;
  31 :  CHESP=0. ;
  32 :  'REPE' BDX NIDX ;
  33 :        NUAA= NIDY*(&BDX-1)+NUDY ;
  34 :        AAA= 'EXTR'(TB3.DESCHI.MATRICEA) NUAA ;
  35 :        CHXX= 'EXTR' (TB3.DESCHI.CHARGE) &BDX ;
  36 :        CHESP=CHESP+(AAA*CHXX) ;
  37 :  'FIN' BDX ;
  38 :  CONC= 'EXCO' MO1 TB4.'SOLU' 'SCAL' ;
  39 :  COESP=(('ABS' CHESP)*ACTIV)+(('LOG' CONC)/('LOG' 10.))  ;
  40 :  COESP= NOMC MO1 COESP ;
  41 :  CCACT= CCACT+COESP ;
  42 : FIN BOU1 ;
  43 : FINP CCACT ;
  44 :