diff --git a/reports/manuscript_semic.tex b/reports/manuscript_semic.tex
index 7bf8e63a3ec9c3774b55731c14cf81351e70ff20..a723c7e8aba3e61341cd3fe46bc0f501d7f08d4a 100644
--- a/reports/manuscript_semic.tex
+++ b/reports/manuscript_semic.tex
@@ -2,6 +2,8 @@
% \documentclass[tc]{copernicus}
\usepackage[utf8]{inputenc}
+\usepackage{authblk}
+
\usepackage{graphicx}
\usepackage{amsmath}
@@ -14,7 +16,7 @@
\linenumbers
% notes and comments
-\usepackage[nosilent,draft]{fixme}
+\usepackage[nosilent,final]{fixme}
\FXRegisterAuthor{mk}{amk}{Mario}
\fxusetheme{colorsig}
@@ -32,7 +34,13 @@
%opening
\title{\textit{SEMIC}: An efficient surface energy and mass balance model applied to the Greenland ice sheet}
-\author{Mario Krapp, Alexander Robinson, and Andrey Ganopolski}
+\author[1,2]{Mario Krapp\thanks{mario.krapp@pik-potsdam.de}}
+\author[3,1]{Alexander Robinson}
+\author[1]{Andrey Ganopolski}
+\affil[1]{Potsdam Institute for Climate Impact Research}
+\affil[2]{Dept. of Zoology, Univ. of Cambridge}
+\affil[3]{Dpto. AstrofÃsica y CC de la Atmósfera, Universidad Complutense de Madrid}
+\date{}
% FIX for lineno/amsmath line numbering issues:
\newcommand*\patchAmsMathEnvironmentForLineno[1]{%
@@ -134,13 +142,13 @@ The parameter $c_\text{eff}$ denotes the effective heat capacity of the snowpack
In a strict sense of the term ``energy balance'' the left-hand-side of Eq.~(\ref{eq:temp}) should be zero.
Here, we assume that surface temperature and the energy are not in equilibrium because the snowpack or surface exerts some thermal inertia.
-Over snow and ice, the temperature cannot exceed 0\textdegree~C.
+Temperatures of snow- and ice-covered surfaces cannot exceed 0\textdegree~C.
However, for computational purposes, we initially assume that $T_s$ represents the potential temperature, which would be observed in the absence of phase transitions, i.e., melting or refreezing.
Once, melting and refreezing has been computed (see Sect.~\ref{subsec:melt_refr}), the residual heat flux $Q_\text{M/R}$ in Eq.~(\ref{eq:temp}) keeps track of any heat flux surplus or deficit and is added back to the energy balance.
This way, $T_s$ never exceeds 0\textdegree~C.
-For coupling to an ice sheet model, the surface mass balance for ice ($SMB_i$) must be known.
-SEMIC separates the total surface mass balance into the surface mass balance for snow and for ice:
+For coupling to an ice sheet model, the surface mass balance for ice ($SMB_i$) is computed by SEMIC.
+It separates the total surface mass balance into the surface mass balance for snow and for ice:
\begin{align}
SMB &= SMB_s + SMB_i = P_s - SU - M + R,\\
SMB_s &= P_s - SU - M_\text{snow} - C_{si},\\
@@ -212,7 +220,7 @@ Fortunately, an analytical solution to this problem exists.
We calculate the roots of the cosine function and then integrate between the roots to solve for average above- and below-freezing mean surface temperatures $T_s^+$ and $T_s^-$. %\mknote{$f_\text{avg} = \frac{1}{b-a}\int_a^b f(x)dx$}
The roots are
\begin{equation}
-t_1 = \frac{24}{2\pi}\cos^{-1}(\frac{T_s}{A}), \quad
+t_1 = \frac{24}{2\pi}\arccos(\frac{T_s}{A}), \quad
t_2 = 24 - t_1. \notag
\end{equation}
Thus, the time span for temperatures above and below freezing is
@@ -224,9 +232,9 @@ These are the integrals of the cosine function
\begin{subequations}
\begin{align}
T_s^+ &= \frac{1}{\Delta t_+}\int_{t_1}^{t_2}T(t)dt \\
- &= \frac{24}{\pi \Delta t_+} \left[-T_s \cos^{-1}(\frac{T_s}{A}) + A\sqrt{1-\frac{T_s^2}{A^2}} + \pi T_s\right] \notag \\
+ &= \frac{24}{\pi \Delta t_+} \left[-T_s \arccos(\frac{T_s}{A}) + A\sqrt{1-\frac{T_s^2}{A^2}} + \pi T_s\right] \notag \\
T_s^- &= \frac{1}{\Delta t_-}\int_{0}^{t_1}T(t) dt + \int_{t_2}^{24} T(t) dt \\
- &= \frac{24}{\pi \Delta t_-} \left[T_s \cos^{-1}(\frac{T_s}{A}) - A\sqrt{1-\frac{T_s^2}{A^2}}\right]. \notag
+ &= \frac{24}{\pi \Delta t_-} \left[T_s \arccos(\frac{T_s}{A}) - A\sqrt{1-\frac{T_s^2}{A^2}}\right]. \notag
\end{align}
\label{eq:above_below}
\end{subequations}
@@ -329,8 +337,8 @@ The actual surface albedo $\alpha$ is then the average of snow albedo $\alpha_s$
\begin{align}
\alpha = \alpha_{s} - f_a (\alpha_s - \alpha_{bg}) \quad &\text{where} \quad \alpha_{bg} =
\begin{cases}
- \alpha_i \quad \text{over ice-covered or} \\
- \alpha_l \quad \text{over ice-free land}
+ \alpha_i \quad \text{for ice-covered or} \\
+ \alpha_l \quad \text{for ice-free land}
\label{eq:alb}
\end{cases} \\ &\text{and} \quad f_a = \exp(-h_s/h_\text{crit}).\notag
\end{align}
@@ -545,7 +553,7 @@ We have developed SEMIC as a coupler between interactive ice sheet models and EM
SEMIC realistically represents the energy transfer between atmosphere and surface as radiation and turbulent mixing of heat and water vapour, thus providing a general solution to the surface energy balance that is applicable for different climates and time scales.
Ice-free land and ice-covered land are treated differently in SEMIC because of the different physical processes involved.
-For example, the surface temperature over ice- and snow-free land has no upper limit as is the case for surface temperatures over ice, which is always lower than or equal to the freezing point.
+For example, the surface temperature of ice- and snow-free land has no upper limit as is the case for surface temperatures of ice, which is always lower than or equal to the freezing point.
Generally, land albedo is much more variable than as described by the single bare land albedo used in SEMIC.
Different land and vegetation types have different effects on the radiation budget.
Consequently, net shortwave radiation errors in SEMIC are larger over ice-free land than over the ice sheet (Fig.~\ref{fig:comp}).
@@ -604,12 +612,16 @@ M.K. was funded by the Deutsche Forschungsgemeinschaft (DFG) Project "Modeling t
We hereby acknowledge, support, and encourage research that follows standards with respect to scientific reproducibility, transparency, and data availability.
Any model source code and the authors' manuscript source (typeset in \LaTeX) is freely available and accessible online.
-The project infrastructure covering individuals step starting from data download and preparation, model source code compilation, running the optimisation, running the calibrated model, running the model with historical and RCP8.5 scenario data, as well as the source code of this manuscript with its figures can be downloaded from the repository website \url{https://gitlab.pik-potsdam.de/krapp/semic-project} or cloned using \texttt{git}\mknote{Latest version will be tagged for submission to journal}:
+The project infrastructure covering individuals step starting from data download and preparation, model source code compilation, running the optimisation, running the calibrated model, running the model with historical and RCP8.5 scenario data, as well as the source code of this manuscript with its figures can be downloaded from the repository website \url{https://gitlab.pik-potsdam.de/krapp/semic-project}.
+See the project website's \texttt{README.md} for details.
+The project can also be cloned using \texttt{git}:
\begin{verbatim}
-git clone git@gitlab.pik-potsdam.de:krapp/semic-project.git
+git clone -b v1.0 git@gitlab.pik-potsdam.de:krapp/semic-project.git
\end{verbatim}
+
+
%%% Tables
\begin{table*}