This zip-file is the Online Supplement to: "Version 1 of a sea ice module for the physics based, detailed, multi-layer SNOWPACK model" by Wever et al. and contains all the necessary input files to run the simulations shown in the mansucript.


SOURCE CODE:
============
The source code and compile instructions can be obtained from https://models.slf.ch. The version required to execute the simulations corresponds to revision 2382 of MeteoIO (https://models.slf.ch/svn/meteoio/trunk) and revision 1751 of SNOWPACK (https://models.slf.ch/svn/snowpack/branches/dev). The LAPACK library is required for reproducibility, and needs to be selected when compiling SNOWPACK using the ENABLE_LAPACK option.


SIMULATION FILES:
=================
The following folder structure will be generated when unpacking the zip-file:
cfgfiles: Contains the model configuration files
input: Contains the model initial state files
smet: Contains the files with meteorological driving data for the model
output: After executing the simulations, the output files will appear here


SETTINGS:
=========
Below follows explanations of sea ice specific keys in the configuration files (*.ini in the cfgfiles folder):

VARIANT	in [SNOWPACKADVANCED] section:			SEAICE: specifies the SNOWPACK model to run the sea ice module.
USEREFERENCELAYER in [OUTPUT] section:			TRUE:   output is referenced with respect to the marked reference level in the *.sno file.
							FALSE:  output is referenced to sea level.
							In the *pro output files, the reference z-value is 4 m.
ENFORCE_MEASURED_SNOW_HEIGHTS in [SNOWPACK] section:	TRUE:   use measured snow depths to derive snowfall events. This rquires the marking of a
							        reference level in the *sno file (mk > 9000)
							FALSE:  use precipitation sum (PSUM) to derive precipitation events
GEO_HEAT in [SNOWPACK] section:				Can be used to prescribe the ocean heat flux. A time series of ocean heat flux may be provided
							by the variable GEO_HEAT in the *smet input file.


To use the explicit brine dynamics, the following settings are mandatory:
WATERTRANSPORTMODEL_SNOW	= RICHARDSEQUATION		(in [SNOWPACKADVANCED] section)
WATERTRANSPORTMODEL_SOIL	= RICHARDSEQUATION		(in [SNOWPACKADVANCED] section)
LB_COND_WATERFLUX		= SEAICE			(in [SNOWPACKADVANCED] section)
SALINITYPROFILE			= NONE				(in [SNOWPACKSEAICE] section)
SALINITYTRANSPORT_SOLVER	= IMPLICIT			(in [SNOWPACKSEAICE] section)
Note that it is possible to set other combinations of settings and start the simulation without receiving warnings or errors, but the simulation setup may not make sense. For example, other options for LB_COND_WATERFLUX are available, but only useful in other contexts, such as alpine snow covers.


To use fixed brine dynamics, the following settings are mandatory:
WATERTRANSPORTMODEL_SNOW	= BUCKET			(in [SNOWPACKADVANCED] section)
WATERTRANSPORTMODEL_SOIL	= BUCKET			(in [SNOWPACKADVANCED] section)
SALINITYPROFILE			= SINUSSAL			(in [SNOWPACKSEAICE] section)




SIMULATIONS:
============
Below follows a list of simulations which were used to produce the figures.


Figure 2:
---------
Execute: snowpack -c cfgfiles/dynamicsal_psum_drytowet.ini -e 2014-01-19T00:32


Figure 3a and 4:
----------------
Execute: snowpack -c cfgfiles/dynamicsal_hs_S31_marked.ini -e NOW


Figure 3b:
----------
Execute: snowpack -c cfgfiles/dynamicsal_psum_S31_marked.ini -e NOW


Figure 5:
---------
Execute: snowpack -c cfgfiles/fixedsal_hs_S31_marked.ini -e NOW


Figure 6 and 7:
---------------
Execute: snowpack -c cfgfiles/dynamicsal_hs_S12_marked.ini -e NOW


Figure 8 and 9:
---------------
Execute: snowpack -c cfgfiles/dynamicsal_psum_warming_marked.ini -e NOW


Figure 10:
----------
Execute: snowpack -c cfgfiles/dynamicsal_psum_cooling_marked.ini -e NOW


Figure 11:
----------
Change L29 in snowpack/snowpackCore/SalinityTransport.cc to:
  >>>  static const bool ZeroFluxUpperBoundary_advection = true;
and recompile SNOWPACK.
Execute: snowpack -c cfgfiles/dynamicsal_thinice.ini -e 2014-05-19T01:00
