Procedures | Cast3M

   1 : * ELECNEUT  PROCEDUR  AURI      98/12/29    21:16:17     3419           
   2 : DEBPROC ELECNEUT TB3*TABLE TBC2*OBJET TBPAR2*OBJET LIMODI/LISTENTI ;
   3 : *********************************************************************
   4 : *
   5 : *      recherche de la neutralité électrique
   6 : *
   7 : *********************************************************************
   8 : MSINEUT=VRAI ;
   9 : **
  10 : 'MESS''Voulez afficher les speciations ?' ;
  11 : 'OBTE' ISPECI ;
  12 : IIMPI= VALE 'IMPI' ;
  13 : IIMP1= VALE 'IMPI' ;
  14 : SI ('EGA' ISPECI 'OUI' ) ;
  15 : IIMP1= 1 ;
  16 : 'FINSI' ;
  17 : OPTION IMPI IIMP1 ;
  18 :  TB4 = CHI2 TB3  TBC2 TBPAR2 ;
  19 : OPTION IMPI IIMPI ;
  20 : **
  21 : EPS= 1.D-3 ;
  22 : CHACHA= 'NEUT' TB3 TB4.'SOLU' ;
  23 :       'MESS' 'Voulez vous afficher le bilan électrique' ;
  24 :       'OBTE' ISUITE ;
  25 :       'SI' ('EGA' ISUITE 'OUI') ;
  26 :        NTAB CHACHA 'TITR' 'BILAN ELECTRIQUE' ;
  27 :       'FINSI' ;
  28 : CHAPLUS= 'EXCO' CHACHA 'CATI' 'SCAL' ;
  29 : CHAMOIN= 'EXCO' CHACHA 'ANIO' 'SCAL' ;
  30 : CHANEUT= (CHAPLUS +CHAMOIN)*2. / ( CHAPLUS -CHAMOIN);
  31 : TEST= (ABS CHANEUT) MASQ 'SUPERIEUR' 'SOMME' EPS ;
  32 : SI ( TEST 'EGA' 0) ;
  33 :              MESS ' Neutralité électrique vérifiée' ;
  34 :              MSINEUT=FAUX ;
  35 : SINON ;
  36 : 'SI' ('NON'( 'EXIS' LIMODI)) ;
  37 :         MESS ' Entrez la liste des composants dont on peut modifier '
  38 :         'les concentrations totales ' ;
  39 :          OBTE LIMODI*LISTENTI ;
  40 : 'FINSI' ;
  41 :  DDZZ=  CHAPLUS+ CHAMOIN ;
  42 :  NBIDX=DIME TB3.DESCHI.IDX ;
  43 :  TOTO= TBC2%ACCES GTOT ;
  44 :  TOT1= COPI TOTO   ;
  45 :  NBMODI= DIME LIMODI ;
  46 :  TABCHA= TABLE ;
  47 :  'REPE' BOUNEUT 20 ;
  48 :  CHATOT= 0 ;
  49 :  'REPE' BCHAR NBMODI ;
  50 :  NN1= EXTR LIMODI &BCHAR ;
  51 :  REPE BOUDX NBIDX ;
  52 :        NUII= EXTR (TB3.DESCHI.IDX) &BOUDX ;
  53 :        SI( 'EGA' NUII NN1) ;
  54 :           CHACDX= EXTR(TB3.DESCHI.CHARGE) &BOUDX ;
  55 :           TABCHA.NN1= CHACDX ;
  56 :        FINSI ;
  57 :  FIN BOUDX ;
  58 :  CHATOT= CHATOT+(ABS  CHACDX)  ;
  59 :  FIN BCHAR ;
  60 :  'REPE' BNTOT NBMODI ;
  61 :  NN1= EXTR LIMODI &BNTOT ;
  62 :  CHACDX= TABCHA.NN1 ;
  63 :  MAPLU= MASQ ( DDZZ*CHACDX) SUPERIEUR 0. ;
  64 :  MANEG= MASQ  (DDZZ*CHACDX) INFERIEUR 0. ;
  65 :  TOTN= ((ABS DDZZ)/CHATOT ) * MANEG ;
  66 :  TOTP= ((ABS DDZZ)/CHATOT ) * MAPLU ;
  67 : *MESS ' TOTI ' ; LIST TOTI ;
  68 :  MO4 MO3 NU2 = NOCOMCHI TB3 'NUMCOMP' NN1 ;
  69 :  TOTN= NOMC MO3 TOTN ;
  70 :  TOTP= NOMC MO3 TOTP ;
  71 :  TOTO=TOTO+TOTN - TOTP ;
  72 :  SI (NEG NN1 50) ;
  73 :  TOT2= 'EXCO' TOTO MO3  MO3 ;
  74 :  TES2= TOT2 MASQ 'INFERIEUR' 'SOMME' 0 ;
  75 :  SI (TES2 > 0) ;
  76 :  MSINEUT= VRAI ;
  77 :  QUITTER BOUNEUT ;
  78 :  FINSI ;
  79 :  FINSI ;
  80 :  FIN BNTOT ;
  81 :  TBC2%GTOT  TOTO  ;
  82 : *TB4 = SPECHI2 TB3  TBC2 TBPAR2 ;
  83 : **
  84 : 'MESS''Voulez afficher les speciations ?' ;
  85 : 'OBTE' ISPECI ;
  86 : IIMPI= VALE 'IMPI' ;
  87 : IIMP1= VALE 'IMPI' ;
  88 : SI ('EGA' ISPECI 'OUI' ) ;
  89 : IIMP1= 1 ;
  90 : 'FINSI' ;
  91 : OPTION IMPI IIMP1 ;
  92 :  TB4 = CHI2 TB3  TBC2 TBPAR2 ;
  93 : OPTION IMPI IIMPI ;
  94 : **
  95 : CHACHA= 'NEUT' TB3 TB4.'SOLU' ;
  96 :       'MESS' 'Voulez vous afficher le bilan électrique' ;
  97 :       'OBTE' ISUITE ;
  98 :       'SI' ('EGA' ISUITE 'OUI') ;
  99 :        NTAB CHACHA 'TITR' 'BILAN ELECTRIQUE' ;
 100 :       'FINSI' ;
 101 : CHAPLUS= 'EXCO' CHACHA 'CATI' 'SCAL' ;
 102 : CHAMOIN= 'EXCO' CHACHA 'ANIO' 'SCAL' ;
 103 : CHANEUT= (CHAPLUS +CHAMOIN)*2. / ( CHAPLUS -CHAMOIN);
 104 : TEST= (ABS CHANEUT) MASQ 'SUPERIEUR' 'SOMME' EPS ;
 105 : SI ( TEST 'EGA' 0) ;
 106 :        MSINEUT=FAUX ;
 107 :        MESS ' Neutralité électrique vérifiée' ;
 108 :        OUBLIE TOT1 ;
 109 :         QUITTER BOUNEUT ;
 110 : FINSI ;
 111 :    DDZZ=  CHAPLUS+ CHAMOIN ;
 112 :  'FIN' BOUNEUT ;
 113 :  SI MSINEUT ;
 114 :  TBC2%GTOT TOT1 ;
 115 :  'MESS' ' On ne peut atteindre l équilibre electrique en '
 116 :   'modifiant les concentrations de ces composants' ;
 117 :   OUBLIE TOTO ;
 118 :  FINSI ;
 119 :  FINSI ;
 120 : 'FINP' MSINEUT ;
 121 :