Generally reasonable input parameters for running the model include:



	[Hg]atm.- 1-20 ng/m^3

	[Hg]aq.	-2-200 ng/m^3 (this is input using up to a 5th order polynomial; 

		for constant concentrations ,enter a value for the a0 coefficient

		and set the remaining coefficients to 0)

	Tatm.	- circa 20 deg. C

	Tatm. range - (Tmax - Tmin; circa 10 deg. C for diel simulations)



	Taq.	- circa 20 deg. C

	Taq. range  - (Tmax - Tmin; circa 2 deg. C for diel simulations)



If you wish, you can enter the polynomial coeficients in the manuscript to generate 

the figures in the text. So far as I can tell, everything is operating properly. 

One limitation with this version is that you can not enter values of 0 for either 

the atmospheric or aqueous elemental mercury concentrations (e.g., you presumably 

could enter a value of 1e-6 ng/m^3). A future version also will set the flux to 

zero when the aqueous temperature drops to zero (or perhaps, below the freezing 

point of sea water).