Télécharger la plaquette Cast3M


fr en

Télécharger clorite.dgibi

Retour à la liste

Numérotation des lignes : 

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

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