Télécharger clorite.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : clorite.dgibi
  2. SAUT PAGE ;
  3. *********************************************************************
  4. * Utilisation des modules CHI1 et CHI2
  5. * Jeu de données pour tester la chlorite d'Aspe
  6. *********************************************************************
  7. *option echo 1 ;
  8. OPTION ELEM QUA4 ;
  9. n1 = 1 ;
  10. n2 = 1 ;
  11. *
  12. * POINT SERVANT A DEFINIR LE CONTOUR
  13. *
  14. A = 0.0 0.0 ;
  15. B = 1. 0.0 ;
  16. C = 1. 2. ;
  17. D = 0. 2. ;
  18. *
  19. * DEFINITION DE SEGMENTS
  20. *
  21. *
  22. AB = A DROIT N1 B ;
  23. BC = B DROIT N2 C ;
  24. CD = C DROIT N1 D ;
  25. DA = D DROIT N2 A ;
  26. *
  27. * DEFINITION DU MAILLAGE
  28. *
  29. *
  30. GP = AB BC CD DA DALL 'PLAN' ;
  31. elim 0.001 gp ;
  32. ab= chang poi1 ab ;
  33. g1= chang poi1 gp ;
  34. g0= diff g1 ab ;
  35. *********************************************************************
  36. *
  37. * La reaction de formation de la chlorite de formule structurale
  38. *
  39. * Si Al O Al Fe3+ Fe2+ Mg OH
  40. * 3,06 0,94 10 1,10 0,32 1,43 2,89 8
  41. *
  42. * s'ecrit
  43. *
  44. * chlo +3,192 H20 +7.56 H+ +0.32 e- <-> 2,89 Mg +3,06 Si +1,75 Fe2+
  45. * +2,04Al(OH)4-
  46. * Table des données de CHI1
  47. * on choisi les solides devant intervenir dans la reaction CHXMX
  48. ***********************************************************************
  49. TABDON=TABLE ;
  50. TABDON.BDD= 'STRASBG' ;
  51. TABDON.IDEN= LECT 1 2 6 7 50 60 61 99 101 165 ;
  52. TABDON.CHXMX= LECT 2157 2166 2141 2147 2150 2152 2154 2148 2139 ;
  53. TABESP1= TABLE ;
  54. TABESP1.IDEN= 2175 ;
  55. TABESP1.LOGK= 9.78 ;
  56. TABESP1.ITYP= 5 ;
  57. TABESP1.COMP= LECT 2 61 165 7 50 99 ;
  58. TABESP1.STOECH= PROG 2.89 3.06 2.04 1.75 -7.56 -0.32 ;
  59. TABDON.NVESP= TABLE ;
  60. TABDON.NVESP.1= TABESP1 ;
  61. TABCLIM=TABLE ;
  62. TABCLIM.TYP3= LECT 50 99 ;
  63. TABDON.CLIM=TABCLIM ;
  64. *
  65. TB1=CHI1 TABDON
  66. COMP '/u/castem/divers/COMPSM'
  67. LOGK '/u/castem/divers/COMPSM' ;
  68. * ----------------------------
  69. * table des paramètres de CHI2
  70. * ----------------------------
  71. TBPAR2= TABLE ;
  72. TBPAR2.'SOUSTYPE'='DONNEES_CHIMIQUES' ;
  73. TBPAR2.FIONI= MANU CHPO G1 1 SCAL 0.0 ;
  74. TBPAR2.LOGC= MANU CHPO G1 10 X001 -2.7 X002 -3.8 X006 -15.
  75. X007 -8. X050 -6.8 X060 -6.9 X061 -4. X099 3. X101 -5.4 X165 -8.5 ;
  76. TOTCA1= MANU CHPO G1 1 X001 1.31D-2 ;
  77. TOTCA2= MANU CHPO AB 1 X001 (4.82D-4 - 1.31D-2) ;
  78. TOTCA= TOTCA1+ TOTCA2 ;
  79. TOTMG1= MANU CHPO G1 1 X002 1.99D-2 ;
  80. TOTMG2= MANU CHPO AB 1 X002 (1.50D-4 - 1.99D-2) ;
  81. TOTMG= TOTMG1+ TOTMG2 ;
  82. TOFE31= MANU CHPO G1 1 X006 2.D-2 ;
  83. TOFE32= MANU CHPO AB 1 X006 (3.66116D-12 - 2.D-2 ) ;
  84. TOTFE3= TOFE31+ TOFE32 ;
  85. TOTFE2 =MANU CHPO G1 1 X007 0.D0 ;
  86. TOTOH= MANU CHPO G1 1 X060 0.D0 ;
  87. TOTHSI1= MANU CHPO G1 1 X061 2.D-1 ;
  88. TOTHSI2= MANU CHPO AB 1 X061 (1.0022474D-4 - 2.D-1) ;
  89. TOTHSI= TOTHSI1 + TOTHSI2 ;
  90. TOTE= MANU CHPO G1 1 X099 0.D0 ;
  91. TOTCO1= MANU CHPO G1 1 X101 1.7D-2 ;
  92. TOTCO2= MANU CHPO AB 1 X101 (2.796D-2 - 1.7D-2) ;
  93. TOTCO= TOTCO1+ TOTCO2 ;
  94. TOTAL1= MANU CHPO G1 1 X165 5.D-2 ;
  95. TOTAL2= MANU CHPO AB 1 X165 (2.28366D-7 - 5.D-2) ;
  96. TOTAL= TOTAL1 + TOTAL2 ;
  97. TOTH= MANU CHPO G1 1 X050 0. ;
  98. TBPAR2.TOT=
  99. TOTCA+TOTMG+TOTFE3+TOTFE2+TOTOH+TOTHSI+TOTE+TOTCO+TOTAL+TOTH ;
  100. CLIMPE1= MANU CHPO G1 1 W026 -3. ;
  101. CLIMPE2= MANU CHPO AB 1 W026 (14. + 3. ) ;
  102. CLIMPE = CLIMPE1 + CLIMPE2 ;
  103. CLIMPH1= MANU CHPO G1 1 W027 7. ;
  104. CLIMPH2= MANU CHPO AB 1 W027 ( 5. - 7. ) ;
  105. CLIMPH = CLIMPH1 + CLIMPH2 ;
  106. TBPAR2.CLIM= CLIMPH+CLIMPE ;
  107. TBPARM= TABLE ;
  108. TBPARM.ITMAX = 80;
  109. TBPARM.EPS= 1.D-6 ;
  110. TBPARM.NFI= 6 ;
  111. TBPARM.IMPRIM = LECT 1 3 ;
  112. TBPARM.SORTIE= MOTS 'PREC' 'TYP5' 'FION' ;
  113. *
  114. TB3= CHI2 TB1 TBPARM TBPAR2 ;
  115. *
  116. * controle des résultats
  117. *
  118. FIONTE1=MANU CHPO G1 1 SCAL 8.07816D-3 'NATURE' DISCRET ;
  119. FIONTE2=MANU CHPO AB 1 SCAL (1.87502D-3 - 8.07816D-3 )
  120. 'NATURE' DISCRET ;
  121. FIONTES= FIONTE1 + FIONTE2 ;
  122. VERR1= ( ABS ( FIONTES - TB3.FION )) MASQUE SUPERIEUR 5.D-8 SOMME ;
  123. *LIST VERR1 ;
  124. PRECTE1= MANU CHPO G1 10 W030 0. W031 0. W032 0. W033 1.02796E-02
  125. W034 0. W035 0. W036 2.64965E-03 W037 1.42876E-01 W038 1.79765E-02
  126. W039 6.88576E-03 'NATURE' DISCRET ;
  127. PRECTE2= MANU CHPO AB 10 W030 0. W031 0. W032 0. W033 -1.02796E-02
  128. W034 0. W035 0. W036 -2.64965E-03 W037 -1.42876E-01
  129. W038 -1.79765E-02 W039 -6.88576E-03 'NATURE' DISCRET ;
  130. PRECTES= PRECTE1+PRECTE2 ;
  131. *LIST PRECTES ;
  132. *LIST TB3.PREC ;
  133. PRECD=PRECTES / 50. ;
  134. VERR2= ( ABS ( PRECTES - TB3.PREC )) MASQUE SUPERIEUR PRECD SOMME ;
  135. *LIST VERR2 ;
  136. *LIST TB3.TYP5 ;
  137. TY5TE1= MANU CHPO G1 10 W030 2.63183E-03 W031 1.70579E-05
  138. W032 2.60854E-07 W033 0. W034 2.28431E-02 W035 1.00674E-04 W036 0.
  139. W037 0. W038 0. W039 0. 'NATURE' DISCRET ;
  140. TY5TE2= MANU CHPO AB 10 W030 ( 3.80867E-11 - 2.63183E-03)
  141. W031 ( 9.52143E-01 - 1.70579E-05) W032 (6.45597E-05 - 2.60854E-07)
  142. W033 5.22263E-04 W034 (1.43073E-13 - 2.28431E-02)
  143. W035 ( 1.63439E-07 - 1.00674E-04) W036 1.74267E-08 W037 9.99926E-01
  144. W038 9.55357E-01 W039 1.50703E-20 'NATURE' DISCRET ;
  145. TY5TES= TY5TE1+TY5TE2 ;
  146. *LIST TY5TES ;
  147. TY5D= TY5TES/ 50. ;
  148. VERR3= ( ABS ( TY5TES - TB3.TYP5 )) MASQUE SUPERIEUR TY5D SOMME ;
  149. *LIST VERR3 ;
  150. ********************************************************************
  151. * On traite le meme cas que precedemment mais au lieu de
  152. * preciser les solides qui interviennent on precise ceux
  153. * qui n'interviennent pas ( CLIM.TYP6 )
  154. ********************************************************************
  155. TABDO2=TABLE ;
  156. TABDO2.BDD= 'STRASBG' ;
  157. TABDO2.IDEN= LECT 1 2 6 7 50 60 61 99 101 165 ;
  158. TABDO2.NVESP= TABLE ;
  159. TABDO2.NVESP.1= TABESP1 ;
  160. TABCLIM.TYP6= LECT 2130 2140 2143 2144 2149 2185 2231
  161. 2153 2160 2161 2167 2168 2169 2176 2177 2178 2179 2180 2181 2182
  162. 2183 2186 2188 2192 2197 2198 2200 2201 2203 2204 2205 2208 2216
  163. 2217 2224 2225 2234 2235 2236 2237 2249 2252 2254 2255 2272 2273
  164. 2224 2225 2234 2235 2236 2237 2249 2252 2254 2255 2272 2273 2274
  165. 2277 2278 2279 2281 2282 2283 2286 2287 2288 2300 2311 2313
  166. 2316 2317 2321 2322 2349 2359 2360
  167. 2371 2372 2374 2375 2377 2378 2380 2381 ;
  168. TABDO2.CLIM=TABCLIM ;
  169. *OPTION DONN 5 ;
  170. TB2=CHI1 TABDO2
  171. COMP '/u/castem/divers/COMPSM' ;
  172. *OPTION DONN 5 ;
  173. * ----------------------------
  174. * la table des paramètres de CHI2 est celle du cas précedent
  175. * mais on reinitialise la force ionique avec la valeur calculée
  176. * ----------------------------
  177. TBPAR2.FIONI=TB3.FION ;
  178. TB4= CHI2 TB2 TBPARM TBPAR2 ;
  179. *
  180. * controle des résultats
  181. *
  182. FIONTE1=MANU CHPO G1 1 SCAL 8.08143D-3 'NATURE' DISCRET ;
  183. FIONTE2=MANU CHPO AB 1 SCAL (1.87502D-3 - 8.08143D-3 )
  184. 'NATURE' DISCRET ;
  185. FIONTES= FIONTE1 + FIONTE2 ;
  186. VERR4= ( ABS ( FIONTES - TB4.FION )) MASQUE SUPERIEUR 5.D-8 SOMME ;
  187. *LIST VERR4 ;
  188. PRECTE1= MANU CHPO G1 10 W030 6.88576E-03 W031 0. W032 0. W033 0.
  189. W034 1.79765E-02 W035 2.64965E-03 W036 0. W037 1.02795E-02
  190. W038 1.42876E-01 W039 0. 'NATURE' DISCRET ;
  191. PRECTE2= MANU CHPO AB 10 W030 -6.88576E-03 W031 0. W032 0. W033 0.
  192. W034 -1.79765E-02 W035 -2.64965E-03 W036 0. W037 -1.02795E-02
  193. W038 -1.42876E-01 W039 0. 'NATURE' DISCRET ;
  194. PRECTES= PRECTE1+PRECTE2 ;
  195. PRECD=PRECTES / 50. ;
  196. VERR5= ( ABS ( PRECTES - TB4.PREC )) MASQUE SUPERIEUR PRECD SOMME ;
  197. TY5TE1= MANU CHPO G1 10 W030 0. W031 2.63201E-03 W032 1.00675E-04
  198. W033 1.70591E-05 W034 0. W035 0. W036 2.60859E-07 W037 0.
  199. W038 0. W039 2.28433E-02 'NATURE' DISCRET ;
  200. TY5TE2= MANU CHPO AB 10 W030 1.50703E-20
  201. W031 (3.80867E-11 - 2.63201E-03) W032 ( 1.63439E-07 - 1.00675E-04)
  202. W033 (9.52143E-01 - 1.70591E-05) W034 9.55356E-01 W035 1.74267E-08
  203. W036 ( 6.45597E-05 - 2.60859E-07) W037 5.22263E-04
  204. W038 9.99926E-01 W039 ( 1.43073E-13 - 2.28433E-02) 'NATURE' DISCRET ;
  205. TY5TES= TY5TE1+TY5TE2 ;
  206. TY5D= TY5TES/ 50. ;
  207. VERR6= ( ABS ( TY5TES - TB4.TYP5 )) MASQUE SUPERIEUR TY5D SOMME ;
  208. ********************************************************************
  209. * CAS No 3
  210. * Presence de redox
  211. ********************************************************************
  212. TABDO3=TABLE ;
  213. TABDO3.BDD= 'STRASBG' ;
  214. TABDO3.IDEN= LECT 1 2 6 7 50 60 61 99 101 165 ;
  215. TABDO3.NVESP= TABLE ;
  216. TABDO3.NVESP.1= TABESP1 ;
  217. TABDO3.CLIM= TABLE ;
  218. TABDO3.CLIM.TYP6= LECT 99 ;
  219. TABDO3.CHXMX= LECT 2157 2166 2141 2147 2150 2152 2154 2148 2139 ;
  220. TB5=CHI1 TABDO3
  221. COMP '/u/castem/divers/COMPSM' ;
  222. * ----------------------------
  223. * table des paramètres de CHI2
  224. * ----------------------------
  225. TBPAR3= TABLE ;
  226. TBPAR3.ITMAX = 80;
  227. TBPAR3.EPS= 1.D-6 ;
  228. TBPAR3.NFI= 6 ;
  229. *TBPAR3.ITERSOLI=35 ;
  230. TBPAR3.IMPRIM= LECT 1 3 ;
  231. TBPAR4= TABLE ;
  232. TBPAR4.'SOUSTYPE'='DONNEES_CHIMIQUES' ;
  233. TBPAR4.FIONI= TB4.FION ;
  234. TBPAR4.LOGC= MANU CHPO G1 10 X001 -2.7 X002 -3.8 X006 -15.
  235. X007 -8. X050 -6.8 X060 -6.9 X061 -4. X099 3. X101 -5.4 X165 -8.5 ;
  236. TOTE1= MANU CHPO G0 1 X099 1.2497250381D-2 ;
  237. TOTE2= MANU CHPO AB 1 X099 9.57287D-16 ;
  238. TOTE= TOTE1 + TOTE2 ;
  239. TOTH1= MANU CHPO G0 1 X050 -2.948552539D-2 ;
  240. TOTH2= MANU CHPO AB 1 X050 5.4659D-2 ;
  241. TOTH= TOTH1 + TOTH2 ;
  242. TBPAR4.TOT=
  243. TOTCA+TOTMG+TOTFE3+TOTFE2+TOTOH+TOTHSI+TOTE+TOTCO+TOTAL+TOTH ;
  244. TBPAR3.SORTIE= MOTS 'PREC' 'TYP5' 'FION' ;
  245. TB6= CHI2 TB5 TBPAR4 TBPAR3 ;
  246. *
  247. * controle des résultats
  248. *
  249. FIONTE1=MANU CHPO G1 1 SCAL 8.08125D-3 'NATURE' DISCRET ;
  250. FIONTE2=MANU CHPO AB 1 SCAL (1.88536D-3 - 8.08125D-3 )
  251. 'NATURE' DISCRET ;
  252. FIONTES= FIONTE1 + FIONTE2 ;
  253. VERR7= ( ABS ( FIONTES - TB6.FION )) MASQUE SUPERIEUR 5.D-8 SOMME ;
  254. *LIST VERR7 ;
  255. PRECTE1= MANU CHPO G1 10 W029 6.88576E-03 W030 0. W031 0. W032 0.
  256. W033 1.02796E-02 W034 0. W035 0. W036 2.64965E-03
  257. W037 1.42876E-01 W038 1.79765E-02 'NATURE' DISCRET ;
  258. PRECTE2= MANU CHPO AB 10 W029 -6.88576E-03 W030 0. W031 0. W032 0.
  259. W033 -1.02796E-02 W034 0. W035 0. W036 -2.64965E-03
  260. W037 -1.42876E-01 W038 -1.79765E-02 'NATURE' DISCRET ;
  261. PRECTES= PRECTE1+PRECTE2 ;
  262. PRECD=PRECTES / 50. ;
  263. VERR8= ( ABS ( PRECTES - TB6.PREC )) MASQUE SUPERIEUR PRECD SOMME ;
  264. TY5TE1= MANU CHPO G1 10 W029 0. W030 2.63200E-03 W031 1.70595E-05
  265. W032 2.60862E-07 W033 0. W034 2.28415E-02
  266. W035 1.00671E-04 W036 0. W037 0. W038 0. 'NATURE' DISCRET ;
  267. TY5TE2= MANU CHPO AB 10 W029 1.67592E-20
  268. W030 ( 3.99616E-11 - 2.63200E-03) W031 (9.85139E-01 - 1.70595E-05 )
  269. W032 ( 6.56605E-05 - 2.60862E-07) W033 5.30806E-04
  270. W034 (1.45031E-13 - 2.28415E-02) W035 (1.68843E-07 - 1.00671E-04)
  271. W036 1.82753E-08 W037 9.99926E-01 W038 9.86350E-01 'NATURE' DISCRET ;
  272. TY5TES= TY5TE1+TY5TE2 ;
  273. TY5D= TY5TES/ 50. ;
  274. VERR9= ( ABS ( TY5TES - TB6.TYP5 )) MASQUE SUPERIEUR TY5D SOMME ;
  275. *
  276. VERR= VERR1+VERR2+VERR3+VERR4+VERR5+VERR6+VERR7+VERR8+VERR9 ;
  277. SI (VERR EGA 0 ) ;
  278. ERRE 0 ;
  279. SINO ;
  280. ERRE 5 ;
  281. FINSI ;
  282. FIN ;
  283.  
  284.  
  285.  
  286.  
  287.  
  288.  
  289.  
  290.  
  291.  
  292.  
  293.  

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