MARS3D / AGRIF model configuration for the Bay of Brest

Date 2022-02-18
Author(s) Petton SebastienORCID1, Dumas Franck2
Affiliation(s) 1 : Ifremer, Univ Brest, CNRS, IRD, LEMAR, F‐29280, Plouzané, France
2 : SHOM, Service Hydrographique et Océanographique de la Marine, 13 rue de Chatellier, CS592803, 29228 Brest CEDEX2, France
DOI 10.17882/86400
Publisher SEANOE
Keyword(s) Bay of Brest, hydrodynamic model, coastal dynamics

This repository provides the source code written in Fortran 90 language based on MARS3D model code (Lazure and Dumas, 2008) implemented with the AGRIF library (Adaptive Grid Refinement In Fortran from Debreu et al. (2008)). It also contains the bathymetric grids and the input namelist files.

The MARS3D/AGRIF model is setup over the Iroise sea (geographic limits 47.74°N - 48.82°N and 4.08°W - 5.55°W) with a horizontal grid resolution of 250m. A zoom over the Bay of Brest is introduced with a resolution of 50m (geographic limits 48.20°N - 48.44°N and 4.09°W - 4.72°W). The time and space refinement factors are both equal to 5. The vertical discretization is performed with 20 equidistant σ-layers in both grids. The timestep for mother grid is set to 30s which consequently leads to 6s for the child grid.

The bathymetries are interpolated from a combination of different digital terrain models (SHOM, Ifremer, IGN). The Iroise model is forced by harmonic components from the SHOM CST-France model (Le Roy and Simon, 2003). The 3D open boundary conditions for temperature and salinity are imposed at hourly frequency from a regional configuration based on a hindcast (Caillaud et al., 2016). Freshwater inputs for the four main rivers in the Bay were taken from the French HYDRO database and corrected with corresponding watershed surface rates

The horizontal turbulent closure consists of a Laplacian operator with a constant turbulent viscosity coefficient. This coefficient was set to 0.5 m2 s−1. The vertical turbulent closure is performed using the generalised length-scale parametrisation with a two-equation k-ε model.

For the child grid, the surface drag coefficient was set uniformly to obtain the optimum validation of the currents in the Bay of Brest. However, the bottom drag coefficient relied on the Prandtl theory of the bottom boundary layer and its logarithmical shape. It includes in its formulation a classical roughness scale, which was spatialised according to a recent sedimentological map (Gregoire et al. 2016).

Licence CC-BY-NC-ND
File Size Format Processing Access
Source code 13 MB Fortran 90 Open access
Top of the page

How to cite 

Petton Sebastien, Dumas Franck (2022). MARS3D / AGRIF model configuration for the Bay of Brest. SEANOE.