Télécharger uferdx.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : uferdx.dgibi
  2. SAUT PAGE ;
  3. *********************************************************************
  4. * Utilisation des modules CHI1 et CHI2
  5. * avec presence de redox
  6. *********************************************************************
  7. *
  8. *
  9. * DEFINITION DU MAILLAGE
  10. *
  11. n1 = 1 ;
  12. n2 = 1 ;
  13. a = 0.0 0.0 ;
  14. b = 0.0 0.1 ;
  15. c = 2. 0.1 ;
  16. d = 2. 0. ;
  17. *
  18. *
  19. option elem qua4 ;
  20. *
  21. ab = a droit n1 b ;
  22. bc = b droit n2 c ;
  23. cd = c droit n1 d ;
  24. da = d droit n2 a ;
  25. *
  26. *
  27. GP = AB BC CD DA DALL 'PLAN' ;
  28. ELIM 0.001 GP ;
  29. AB= CHANG POI1 AB ;
  30. G1= CHANG POI1 GP ;
  31. G0= DIFF G1 AB ;
  32. *
  33. * TABLE ENTREE DE CHI1
  34. TABDON=TABLE ;
  35. TABDON.IDEN= LECT 5 4 103 99 6 7 32 39 40 41 50 ;
  36. TABDON.CLIM= TABLE ;
  37. TABDON.CLIM . TYP6 = LECT 99 ;
  38. TABDON.CHXMX = LECT 40007 20310 40501 32530 32520 21460
  39. 32009 32011 ;
  40. * on rajoute 3 espèces ne figurant pas dans
  41. * la base de données
  42. TABDON.NVESP= TABLE ;
  43. TABDON.NVESP.1 =TABLE ;
  44. TABDON.NVESP.1 .IDEN= 40501 ;
  45. TABDON.NVESP.1 .LOGK= -92.81 ;
  46. TABDON.NVESP.1 .ITYP= 5 ;
  47. TABDON.NVESP.1 .COMP= LECT 40 99 ;
  48. TABDON.NVESP.1 .STOECH= PROG 1. 4. ;
  49. TABDON.NVESP.2 =TABLE ;
  50. TABDON.NVESP.2 .IDEN= 32530 ;
  51. TABDON.NVESP.2 .LOGK= 13.03 ;
  52. TABDON.NVESP.2 .ITYP= 5 ;
  53. TABDON.NVESP.2 .COMP= LECT 32 40 50 ;
  54. TABDON.NVESP.2 .STOECH= PROG 1. 3. -14. ;
  55. TABDON.NVESP.3 =TABLE ;
  56. TABDON.NVESP.3 .IDEN= 32520 ;
  57. TABDON.NVESP.3 .LOGK= -2.46 ;
  58. TABDON.NVESP.3 .ITYP= 5 ;
  59. TABDON.NVESP.3 .COMP= LECT 32 40 50 ;
  60. TABDON.NVESP.3 .STOECH= PROG 2. 1. -8. ;
  61. *OPTION DONN 5 ;
  62. *
  63. TB1=CHI1 TABDON
  64. COMP '/u/castem/divers/COMPOM' ;
  65. *
  66. * TABLE DES PARAMETRES DE CHI2
  67. *OPTION DONN 5 ;
  68. TBPAR2= TABLE ;
  69. TBPAR2.'SOUSTYPE'='DONNEES_CHIMIQUES' ;
  70. TBPAR2.FIONI= MANU CHPO G1 1 SCAL 0.001 ;
  71. TBPAR2.LOGC= MANU CHPO G1 11 X005 -3. X004 -11. X103 -3. X099 -5.
  72. X006 -10. X007 -10. X032 -10. X039 -10. X040 -10. X041 -10.
  73. X050 -7. ;
  74. *
  75. * on forme le CHPOIN des concentrations totales
  76. *
  77. TOTNA1= MANU CHPO G1 1 X005 1.0008D-3 ;
  78. TOTNA2= MANU CHPO AB 1 X005 (1.D-2 - 1.0008D-3) ;
  79. TOTNA= TOTNA1 + TOTNA2 ;
  80. TOTK1=MANU CHPO G1 1 X004 1.D-11 ;
  81. TOTK2=MANU CHPO AB 1 X004 (1.D-5 - 1.D-11) ;
  82. TOTK= TOTK1 + TOTK2 ;
  83. TOTCL1=MANU CHPO G1 1 X103 1.0915D-3 ;
  84. TOTCL2=MANU CHPO AB 1 X103 (1.0022057D-2 - 1.0915D-3) ;
  85. TOTCL=TOTCL1 + TOTCL2 ;
  86. TOTE1= MANU CHPO G1 1 X099 4.5619624684D-5 ;
  87. TOTE2= MANU CHPO AB 1 X099 (1.49708535D-12 - 4.5619624684D-5) ;
  88. TOTE= TOTE1 + TOTE2 ;
  89. TOTFE3= MANU CHPO G1 1 X006 4.559E-5 ;
  90. TOTFE5= MANU CHPO AB 1 X006 (1.D-8 - 4.559E-5) ;
  91. TOTFE3= TOTFE3 + TOTFE5 ;
  92. TOTFE2= MANU CHPO G1 1 X007 0.D0 ;
  93. TOTU6 = MANU CHPO G1 1 X032 2.D-8 ;
  94. TOTU = MANU CHPO AB 1 X032 ( 1.D-6 - 2.D-8) ;
  95. TOTU6 = TOTU6 + TOTU ;
  96. TOTU3 = MANU CHPO G1 1 X039 0.D0 ;
  97. TOTU4 = MANU CHPO G1 1 X040 0.D0 ;
  98. TOTU5 = MANU CHPO G1 1 X041 0.D0 ;
  99. TOTH1= MANU CHPO G1 1 X050 -4.8578064D-7 ;
  100. TOTH2= MANU CHPO AB 1 X050 (1.0028126D-5 + 4.8578064D-7) ;
  101. TOTH = TOTH1 + TOTH2 ;
  102. TBPAR2.TOT= TOTNA + TOTK + TOTCL + TOTE + TOTFE3 + TOTFE2 + TOTU6 +
  103. TOTU5 + TOTU4 + TOTU3 + TOTH ;
  104. TBPARM= TABLE ;
  105. TBPARM.ITMAX = 80;
  106. TBPARM.EPS= 1.D-6 ;
  107. TBPARM.NFI= 4 ;
  108. *TBPARM.ITERSOLI=15 ;
  109. *
  110. * on précise les elements de la table CHIMI2
  111. * que l'on souhaite conserver
  112. *
  113. TBPARM.SORTIE= MOTS 'TYP5' 'PREC' 'FION' ;
  114. *
  115. TB3= CHI2 TB1 TBPARM TBPAR2 ;
  116. * CONTROLE DES RESULTATS
  117.  
  118. FIONTE1=MANU CHPO G1 1 SCAL 1.13636E-03 'NATURE' DISCRET ;
  119. FIONTE2=MANU CHPO AB 1 SCAL (1.00229E-02 - 1.13636E-03 )
  120. 'NATURE' DISCRET ;
  121. FIONTES= FIONTE1 + FIONTE2 ;
  122. VERR1= ( ABS ( FIONTES - TB3.FION )) MASQUE SUPERIEUR 5.D-8 SOMME ;
  123. PRECTE1= MANU CHPO G1 8 W050 0. W051 0. W052 5.18766E-09 W053 0.
  124. W054 0. W055 1.99605E-08 W056 0. W057 0. 'NATURE' DISCRET ;
  125. PRECTE2= MANU CHPO AB 8 W050 0. W051 0. W052 -5.18766E-09 W053 0.
  126. W054 0. W055 -1.99605E-08 W056 0. W057 0. 'NATURE' DISCRET ;
  127. PRECTES= PRECTE1+PRECTE2 ;
  128. PRECD=PRECTES / 50. ;
  129. VERR2= ( ABS ( PRECTES - TB3.PREC )) MASQUE SUPERIEUR PRECD SOMME ;
  130. TY5TE1= MANU CHPO AB 8 W050 8.28223E-16 W051 4.85086E-02
  131. W052 2.07294E-03 W053 1.31814E-04 W054 1.65945E-02 W055 2.17882E-16
  132. W056 2.31600E-180 W057 1.16045E-44 'NATURE' DISCRET ;
  133. TY5TE2= MANU CHPO G0 8 W050 1.15633E-36 W051 8.45130E-07 W052 0.
  134. W053 7.27012E-23 W054 9.15253E-21 W055 0. W056 7.36067E-97
  135. W057 6.18787E-16 'NATURE' DISCRET ;
  136. TY5TES= TY5TE1+TY5TE2 ;
  137. TY5D= TY5TES/ 50. ;
  138. * DANS LE CAS DES CONCENTRATIONS TRES FAIBLES
  139. * ON ADMETTRA UNE ERREUR PLUS IMPORTANTE SUR LE RESULTAT
  140. TT1= TY5TES MASQUE DIFFERENT 0. ;
  141. TTT= TT1 MASQUE INFERIEUR 1.E-40 ;
  142. TY5D= TY5D + ( TTT / 2. ) ;
  143. VERR3= ( ABS ( TY5TES - TB3.TYP5 )) MASQUE SUPERIEUR TY5D SOMME ;
  144. *OPTION DONN 5 ;
  145. * TABLE ENTREE DE CHI1
  146. TABLIM2= TABLE ;
  147. TABLIM2. TYP3 = LECT 50 99 ;
  148. TABDON.CLIM = TABLIM2 ;
  149. *OPTION DONN 5 ;
  150. TB4=CHI1 TABDON
  151. COMP '/u/castem/divers/COMPOM'
  152. LOGK '/u/castem/divers/COMPOM' ;
  153. *
  154. * TABLE ENTREE DE CHI2
  155. *
  156. TBPAR3= TABLE ;
  157. TBPAR3.'SOUSTYPE'='DONNEES_CHIMIQUES' ;
  158. TBPAR2.FIONI= MANU CHPO G1 1 SCAL 0.001 ;
  159. TBPAR3.LOGC= MANU CHPO G1 11 X005 -3. X004 -11. X103 -3. X099 -5.
  160. X006 -10. X007 -10. X032 -10. X039 -10. X040 -10. X041 -10.
  161. X050 -7. ;
  162. TBPAR3.FIONI= TB3.FION ;
  163. PH1=MANU CHPO G1 1 W046 7.D0 ;
  164. PH2=MANU CHPO AB 1 W046 (5.D0 - 7.D0) ;
  165. PH=PH1+PH2 ;
  166. PE1=MANU CHPO G1 1 W045 -7.5D0 ;
  167. PE2=MANU CHPO AB 1 W045 (12.D0 + 7.5D0) ;
  168. PE=PE1+PE2 ;
  169. TBPAR3.CLIM= PH + PE ;
  170. TOTE=MANU CHPO G1 1 X099 0.D0 ;
  171. TOTH=MANU CHPO G1 1 X050 0.D0 ;
  172. TBPAR3.TOT= TOTNA + TOTK + TOTCL + TOTE + TOTFE3 + TOTFE2 + TOTU6 +
  173. TOTU5 + TOTU4 + TOTU3 + TOTH ;
  174. TB5= CHI2 TB4 TBPARM TBPAR3 ;
  175. * CONTROLE DES RESULTATS
  176.  
  177. FIONTE1=MANU CHPO G1 1 SCAL 1.13652E-03 'NATURE' DISCRET ;
  178. FIONTE2=MANU CHPO AB 1 SCAL (1.00178E-02 - 1.13652E-03 )
  179. 'NATURE' DISCRET ;
  180. FIONTES= FIONTE1 + FIONTE2 ;
  181. FIOND= FIONTES/ 2.D5 ;
  182. VERR4= ( ABS ( FIONTES - TB5.FION )) MASQUE SUPERIEUR FIOND SOMME ;
  183. PRECTE1= MANU CHPO G1 8 W051 0. W052 0. W053 0. W054 0.
  184. W055 1.99815E-08 W056 0. W057 0. W058 0. 'NATURE' DISCRET ;
  185. PRECTE2= MANU CHPO AB 8 W051 0. W052 0. W053 0. W054 0.
  186. W055 -1.99815E-08 W056 0. W057 0. W058 0. 'NATURE' DISCRET ;
  187. PRECTES= PRECTE1+PRECTE2 ;
  188. PRECD=PRECTES / 50. ;
  189. VERR5= ( ABS ( PRECTES - TB5.PREC )) MASQUE SUPERIEUR PRECD SOMME ;
  190. TY5TE1= MANU CHPO AB 8 W051 4.45529E-02 W052 7.02472E-04
  191. W053 1.13068E-04 W054 1.42344E-02 W055 3.57552E-17
  192. W056 1.39104E-182 W057 4.39903E-47 W058 1.00004E-16
  193. 'NATURE' DISCRET ;
  194. TY5TE2= MANU CHPO G0 8 W051 1.22641E-07 W052 4.63350E-03
  195. W053 3.16228E-23 W054 3.98107E-21 W055 0. W056 3.89045E-96
  196. W057 2.69153E-16 W058 2.18776E-37 'NATURE' DISCRET ;
  197. TY5TES= TY5TE1+TY5TE2 ;
  198. TY5D= TY5TES/ 50. ;
  199. * DANS LE CAS DES CONCENTRATIONS TRES FAIBLES
  200. * ON ADMETTRA UNE ERREUR PLUS IMPORTANTE SUR LE RESULTAT
  201. TT1= TY5TES MASQUE DIFFERENT 0. ;
  202. TTT= TT1 MASQUE INFERIEUR 1.E-40 ;
  203. TY5D= TY5D + ( TTT / 2. ) ;
  204. VERR6= ( ABS ( TY5TES - TB5.TYP5 )) MASQUE SUPERIEUR TY5D SOMME ;
  205. *
  206. VERR= VERR1+VERR2+VERR3+VERR4+VERR5+VERR6 ;
  207. SI (VERR EGA 0 ) ;
  208. ERRE 0 ;
  209. SINO ;
  210. ERRE 5 ;
  211. FINSI ;
  212. FIN ;
  213.  
  214.  
  215.  
  216.  
  217.  
  218.  
  219.  
  220.  
  221.  

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