Télécharger eauacti.dgibi

Retour à la liste

Numérotation des lignes :

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

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