Skip to content
Snippets Groups Projects
Commit 20034cd3 authored by Guillaume Samson's avatar Guillaume Samson :snowman2:
Browse files

Merge branch '16-update-time-dependant-reference-pressure-option-in-sbc_apr' into 'main'

Resolve "update time-dependant reference pressure option in sbc_apr"

Closes #16

See merge request !26
parents b0f8339e 37917368
No related branches found
No related tags found
1 merge request!26Resolve "update time-dependant reference pressure option in sbc_apr"
......@@ -1066,9 +1066,12 @@ an equivalent inverse barometer sea surface height, $\eta_{ib}$, using:
\eta_{ib} = - \frac{1}{g\,\rho_o} \left( P_{atm} - P_o \right)
\]
where $P_{atm}$ is the atmospheric pressure and $P_o$ a reference atmospheric pressure.
A value of $101,000~N/m^2$ is used unless \np{ln_ref_apr}{ln\_ref\_apr} is set to true.
In this case, $P_o$ is set to the value of $P_{atm}$ averaged over the ocean domain,
\ie\ the mean value of $\eta_{ib}$ is kept to zero at all time steps.
3 different options are available to define this reference atmospheric pressure using \np{nn_ref_apr}{nn\_ref\_apr}:
\begin{description}
\item $\bullet$ 0 sets $P_o$ to a constant value of $101,000~N/m^2$
\item $\bullet$ 1 sets $P_o$ to the value of $P_{atm}$ averaged over the ocean domain (the mean value of $\eta_{ib}$ is kept to zero at all time steps)
\item $\bullet$ 2 reads a time-dependant $P_o$ value from a 1D file defined via \np{sn_apref}{sn\_apref} structure (usefull for regional configurations)
\end{description}
The gradient of $\eta_{ib}$ is added to the RHS of the ocean momentum equation (see \mdl{dynspg} for the ocean).
For sea-ice, the sea surface height, $\eta_m$, which is provided to the sea ice model is set to $\eta - \eta_{ib}$
......
Subproject commit 0e8152c2f6ad129e5bf1bd1a2c85173a2916bb78
Subproject commit 43a9e90e73133b4d705f0576a474c5c0b04c2204
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment