max/iphreeqc
1
0
Fork 0
mirror of https://git.gfz-potsdam.de/naaice/iphreeqc.git synced 2025-12-16 16:44:49 +01:00
Code Issues Packages Projects Releases Wiki Activity

minor updates for manual

Browse Source 
git-svn-id: svn://136.177.114.72/svn_GW/phreeqc3/trunk@6622 1feff8c3-07ed-0310-ac33-dd36852eb9cd
This commit is contained in:
David L Parkhurst 2012-06-01 18:06:33 +00:00
parent a48d816c6f
commit 049b488cf4
4 changed files with 15 additions and 13 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
Makefile
View File

@ -17,7 +17,7 @@ endif
all: ex1.out ex2.out ex3.out ex4.out ex5.out ex6.out ex7.out ex8.out ex9.out \
ex10.out ex11.out ex12.out ex12a.out ex13a.out ex13b.out ex13c.out ex13ac.out \
ex14.out ex15.out ex16.out ex17.out ex17b.out ex18.out \
ex19.out ex19b.out ex20a.out ex20b.out ex21_radial.out
ex19.out ex19b.out ex20a.out ex20b.out ex21_radial.out ex22.out
ex1.out: $(INPUT)/ex1 $(PHREEQC) $(PHREEQCDAT)
$(PHREEQC) $(INPUT)/ex1 ex1.out $(PHREEQCDAT)

17
ex11
View File

@ -45,9 +45,9 @@ USER_GRAPH 1 Example 11
-start
10 x = (STEP_NO + 0.5) / cell_no
20 PLOT_XY x, TOT("Cl")*1000, symbol = None, color = Red
30 PLOT_XY x, TOT("Na")*1000, symbol = None
40 PLOT_XY x, TOT("K")*1000, symbol = None
50 PLOT_XY x, TOT("Ca")*1000, symbol = None
30 PLOT_XY x, TOT("Na")*1000, symbol = None, color = Blue
40 PLOT_XY x, TOT("K")*1000, symbol = None, color = Orange
50 PLOT_XY x, TOT("Ca")*1000, symbol = None, color = Magenta
60 PUT(1, 1)
-end
COPY cell 101 1-40
@ -82,10 +82,10 @@ USER_GRAPH 2 Example 11
-axis_scale y_axis 0 1.4
-plot_concentration_vs time
10 x = (STEP_NO + 0.5) / cell_no
20 PLOT_XY x, TOT("Cl")*1000, symbol = Plus, symbol_size = 2 color = Red
30 PLOT_XY x, TOT("Na")*1000, symbol = Plus, symbol_size = 2 color = Blue
40 PLOT_XY x, TOT("K") *1000, symbol = Plus, symbol_size = 2 color = Orange
50 PLOT_XY x, TOT("Ca")*1000, symbol = Plus, symbol_size = 2 color = Magenta
20 PLOT_XY x, TOT("Cl")*1000, symbol = Plus, symbol_size = 2, color = Red
30 PLOT_XY x, TOT("Na")*1000, symbol = Plus, symbol_size = 2, color = Blue
40 PLOT_XY x, TOT("K") *1000, symbol = Plus, symbol_size = 2, color = Orange
50 PLOT_XY x, TOT("Ca")*1000, symbol = Plus, symbol_size = 2, color = Magenta
# calculate Cl_analytical...
60 DATA 0.254829592, -0.284496736, 1.421413741, -1.453152027, 1.061405429, 0.3275911
70 READ a1, a2, a3, a4, a5, a6
@ -101,7 +101,8 @@ USER_GRAPH 2 Example 11
170 PLOT_XY x, y, line_width = 0, symbol = Circle, color = Red
180 d = (y - TOT("Cl")*1000)^2
190 IF EXISTS(10) THEN PUT(d + GET(10), 10) ELSE PUT(d, 10)
200 IF STEP_NO = 2 * CELL_NO THEN print 'SSQD for Cl after 2 Pore Volumes: ', GET(10), '(mmol/L)^2'
200 IF STEP_NO = 2 * CELL_NO THEN print 'SSQD for Cl after 2 Pore Volumes: ', \
GET(10), '(mmol/L)^2'
210 END
2000 REM calculate e_erfc = exp(PA) * erfc(z)...
2010 sgz = SGN(z)

6
ex12
View File

@ -51,14 +51,14 @@ TRANSPORT # Make column temperature 0C, displace Cl
-diffusion_coefficient 0.0 # No diffusion
-thermal_diffusion 1.0 # No retardation for heat
END
SOLUTION 0 Fixed temp 24C, and NaCl conc (first type boundary cond) at inlet
SOLUTION 0 Fixed temp 24C, and NaCl conc (first type bound cond) at inlet
units mol/kgw
temp 24
pH 7.0
pe 12.0 O2(g) -0.67
Na 24.e-3
Cl 24.e-3
SOLUTION 58-60 Same as soln 0 in cell 20 at closed column end (second type boundary cond)
SOLUTION 58-60 Same as soln 0 in cell 20 at closed column end (second type bound cond)
units mol/kgw
temp 24
pH 7.0
@ -79,7 +79,7 @@ TRANSPORT # Diffuse 24C, NaCl solution from column end
USER_GRAPH Example 12
-headings Na Cl Temp Analytical
-chart_title "Diffusion of Solutes and Heat"
-axis_titles "DISTANCE, IN METERS" "MILLIMOLES PER KILOGRAM WATER", "DEGREES CELSIUS"
-axis_titles "DISTANCE, IN METERS" "MILLIMOLES PER KILOGRAM WATER" "DEGREES CELSIUS"
-axis_scale x_axis 0 20
-axis_scale y_axis 0 25
-axis_scale sy_axis 0 25

3
ex9
View File

@ -124,7 +124,8 @@ SELECTED_OUTPUT
-reset false
USER_PUNCH
-headings Days Fe(2) Fe(3) pH si_goethite
10 PUNCH SIM_TIME / 3600 / 24, TOT("Fe_di")*1e6, TOT("Fe_tri")*1e6, -LA("H+"), SI("Goethite")
10 PUNCH SIM_TIME / 3600 / 24, TOT("Fe_di")*1e6, TOT("Fe_tri")*1e6, -LA("H+"), \
SI("Goethite")
USER_GRAPH Example 9
-headings _time_ Fe(2) Fe(3) pH
-chart_title "Oxidation of Ferrous Iron"