Télécharger la plaquette Cast3M


fr en

Télécharger eauacti.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  eauacti.dgibi
  2. * repertoire des fichiers "divers"
  3. DIVERS = VENV 'CASTEM_DIVERS';
  4. *
  5.  SAUT PAGE ;
  6. ********************************************************************
  7. *  Test de bon fonctionnement des operateurs LOGK COAC FION et NEUT
  8. *               eau + temperature  80.   + pression CO2
  9. *********************************************************************
  10. OPTION DIME 2 ;                                                         
  11. A= 0. 0. ;                                                              
  12. B= 5. 0. ;
  13. OPTION  ELEM QUA4 ;                                                     
  14.  AB= A DROIT 1 B ;                                                      
  15.  MT = CHANGE AB POI1  ;
  16.   XX YY = COOR MT ;
  17. *
  18. *                table de données pour CHI1
  19. *
  20. TABDON=TABLE ;                                                          
  21. TABDON.BDD= 'STRASBG' ;
  22. TABDON.IDEN= LECT 1 50 60 61 101 165 ;
  23.  TABDON.CHXMX= LECT 2144 2148 2157 2166 2168 2192 2198 2200
  24.     2224 2231 2249 2278 2281 2372 ;
  25.  TABCLIM=TABLE ;
  26.  TABCLIM.TYP6= LECT 2002 ;
  27.  TABCLIM.TYP3= LECT 2400 ;
  28.  TABCLIM.COMP3= LECT 101 ;
  29.  TABDON.CLIM=TABCLIM ;
  30.  TABDON.TEMPERATURE = 'OUI' ;
  31. *
  32.    TB1=CHI1 TABDON    
  33.   LOGK ('CHAINE' DIVERS '/COMPSM')
  34.   ENTH  ('CHAINE' DIVERS '/COMPSM')
  35.   COMP  ('CHAINE' DIVERS '/COMPSM') ;
  36. *
  37. *
  38. *                   table de données pour CHI2
  39. *
  40.  TBPARM= TABLE ;
  41.  TBPAR2= TABLE ;
  42.  TBPAR2.'SOUSTYPE'= 'DONNEES_CHIMIQUES' ;
  43.  TBPARM.ITMAX = 95;
  44.  TBPARM.ITERSOLI  = 10 ;
  45.  TBPARM.NFI  = 8;
  46.  TBPARM.EPS= 1.D-4 ;
  47.  TBPAR2.LOGC= MANU CHPO MT 6 X001 -3.6 X050 -6.0 X060 -5.0
  48.               X061 -3. X101 -6.5  X165 -6.5  ;
  49.   TOT001= MANU CHPO MT 1 X001  6.D-4 ;
  50.   TOT050= MANU CHPO MT 1 X050  -1.15D-3  ;
  51.   TOT060= MANU CHPO MT 1 X060  0.D0 ;
  52.   TOT061= MANU CHPO MT 1 X061  2.8D-2 ;
  53.   TOT101= MANU CHPO MT 1 X101  0.D0 ;
  54.   TOT165= MANU CHPO MT 1 X165  5.5D-5 ;
  55.  TBPAR2.TOT= TOT001+ TOT050 + TOT060 + TOT061 + TOT101 + TOT165 ;
  56. *
  57. *     définition du champ de temperature
  58. *
  59.   TMPE= MANU CHPO MT 1 SCAL 80. ;
  60.   TBPAR2.TEMPE= TMPE ;
  61. *
  62. *
  63. *  KCO2= MANU CHPO MT 1 W016 20.45 ;
  64.    KCO2= MANU CHPO MT 1 W016 2.32 ;
  65.   TBPAR2.CLIM= KCO2 ;
  66.  TBPARM.SORTIE= MOTS 'PREC' 'FION' 'TYP5' 'SOLU' ;
  67.  TBPARM.IMPRIM= LECT 1   ;
  68. *
  69. * option impi 1 ;
  70.  TB3= CHI2 TB1 TBPARM TBPAR2 ;
  71.   option impi 0 ;
  72. *
  73. *                  Opérateur LOGK
  74. *
  75.   ONI= TB3.FION ;
  76.  LLK= LOGK TB1  FORCEIONI ONI TEMPERATURE TMPE ;
  77. *  LIST LLK ;
  78. *
  79. *                 Verification du résultat
  80. *    Les valeurs LKREF sont celles figurant dans la colonne LOGK
  81. *       dans le tableau sorti par CHIMI2 avec l'option  impi 1
  82. *
  83.  LKREF = MANU CHPO MT 32  W001 0.  W002 0. W003 0.  W004 0.  W005  0.
  84.  W006 0. W007  14.799 W008 1.41924 W009  3.5127 W010 12.063
  85.  W011 -8.982 W012 12.103 W013 6.6868 W014 10.049 W015 16.374
  86.  W016 18.266 W017 12.564  W018  12.793 W019 8.871 W020  3.2849
  87.  W021  33.481 W022 38.633 W023  27.960  W024  32.431 W025 41.452
  88.  W026 29.842 W027 2.4878 W028 12.715 W029  42.9798 W030 47.0985
  89.  W031  26.582 W032 18.491   'NATURE' DISCRET ;
  90.  VERLK=( ABS( LLK-LKREF)) MASQ SUPERIEUR  0.01 SOMME  ;
  91.  LIST VERLK ;
  92.  SI ( VERLK EGA 0) ;
  93.      ERRE 0 ;
  94.  SINO ;
  95.      ERRE 5 ;
  96.  FINSI ;
  97. *
  98. *                   Opérateur NEUT
  99. *
  100.  CHA= NEUT TB1 ( TB3.SOLU ) ;
  101. *LIST CHA ;
  102. *
  103. *                  Vérification du résultat
  104. *      Les valeurs de référence pour NCAT et NANI sont les valeurs
  105. *    de cations et anions dans le tableau imprimé par CHIMI2 avec
  106. *    l'option impi 1
  107. *
  108.  NCAT= EXTR CHA  'CATI' ((EXTR  CHA MAIL) POINT 1) ;
  109.  NANI= EXTR CHA  'ANIO' ((EXTR  CHA MAIL) POINT 1) ;
  110.  LIST NANI ;
  111.  SI (ABS(NANI + 0.001136 ) > 0.000001) ;
  112.      ERRE 5 ;
  113.  FINSI ;
  114.  LIST NCAT ;
  115.  SI (ABS(NCAT - 0.001131 ) > 0.000001) ;
  116.      ERRE 5 ;
  117.  FINSI ;
  118. *
  119. *                Opérateur FION
  120. *
  121.  FIO = FION TB1 ( TB3.SOLU ) ;
  122. *LIST FIO ;
  123. *
  124. *              vérification du résultat
  125. *
  126.  FFF= RESU( ABS(TB3.FION - FIO)) ;
  127.  VFF=  EXTR FFF 'SCAL' ((EXTR FFF MAIL) POINT 1) ;
  128.  LIST VFF ;
  129.  SI (VFF > 1.D-10) ;
  130.       ERRE 5 ;
  131.  FINSI ;
  132. *
  133. *                Opérateur COAC
  134. *
  135.  ACT = COAC TB1 FORCEIONI ONI TEMPERAT TMPE ;
  136. *LIST ACT ;
  137.  VAC=  EXTR ACT 'SCAL' ((EXTR ACT MAIL) POINT 1) ;
  138.  LIST VAC ;
  139.  SI (ABS(VAC + 0.02229 ) > 0.0001) ;
  140.      ERRE 5 ;
  141.  FINSI ;
  142. *
  143. *
  144. *
  145. fin ;
  146.  
  147.  
  148.  
  149.  
  150.  
  151.  
  152.  
  153.  
  154.  

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