CROCO – training 2019
Introduction to CROCO model CROCO
Coastal and Regional Ocean COmmunity model
CROCO – training 2019 - Brest
CROCO history CROCO
Coastal and Regional Ocean COmmunity model
SCRUM
ROMS UCLA / RUTGERS
ROMS Regional Oceanic Modeling System AGRIF
ROMS AGRIF NH-NBQ of Symphonie-NH
Mutli-resolution
MUSTANG sediment module of MARS
CROCO
Coastal and Regional Ocean COmmunity model
CROCO philosophy CROCO
Coastal and Regional Ocean COmmunity model
Community development High-level numerical schemes
Continuous development Realistic complex modelling
Stable releases: every 1 / Circulation forced by waves
1.5 year
Coupling with atmosphere and
Help/support through a waves
forum Tools for an easy
built of regional and Coupling with biogeochemistry
coastal and ecology
configurations
LES / DNS
Adapted to IRD
CROCO examples CROCO
Coastal and Regional Ocean COmmunity model
For starting, here are a few examples of use of CROCO
Gula et al.
High resolution in
the Gulf Stream
NBQ mode
CROCO examples CROCO
Coastal and Regional Ocean COmmunity model
For starting, here are a few examples of use of CROCO
Ponte & Klein, 2015,, internal Marchesiello et al. 2015, Rip current
tides and eddies
Runup min
0.2
0.1
Sea level
0
-0.1
CROCO
OBS
-0.2
0 10 20 30 40 50 60 70 80 90
Setup Hs
CROCO wave-to-wave (Marchesiello, Benshila)
0.35 Runup max
0.2
0.03 CROCO
OBS 0.3
0.25 0.1
0.02
0.2
Sea level
0.01 0
0.15
0.1
0 -0.1
CROCO
0.05 CROCO
OBS
OBS
-0.01 0
10 20 30 40 50 60 70 80 90 10 20 30 40 50 60 70 80 90 -0.2
0 10 20 30 40 50 60 70 80 90
CROCO numerics CROCO
Coastal and Regional Ocean COmmunity model
• Solves the Primitive Equations in an Earth-centered rotating environment:
* momentum conservation
* continuity
* tracer conservation
* equation of state
• Boussinesq hystrostatic mode, and non-hydrostatic, non-Boussinesq mode (NBQ) available
• Split-explicit time-stepping: see dedicated course
* short time steps are used to advance the surface elevation and barotropic momentum
* much larger time step used for temperature, salinity, and baroclinic momentum
* for NBQ mode: barotropic mode solver is replaced by a fully 3D fast mode solver,
resolving all waves down to acoustic waves (with sound speed that can be decreased to
the maximum wave velocity one wants to solve)
CROCO numerics CROCO
Coastal and Regional Ocean COmmunity model
CROCO grid is discretized in coastline- and terrain-following curvilinear coordinates with free-
surface, on an Arakawa-C grid
free-surface vertical stretching coord. Water column thickness
(changed in sediment applications with erosion/deposition)
z(x, y, σ, t) = ζ.(1+ σ) + hc .σ + (h – hc).C(σ)
parameter controlling stretching Stretching function C=f( 𝝷s 𝝷b)
between surface and bottom levels
CROCO numerics CROCO
Coastal and Regional Ocean COmmunity model
High-order numerics
High-order numerical schemes: 3rd and 5th-order advection schemes
Rotated tensors to reduce diapycnal mixing
see dedicated
KPP and GLS mixing parameterizations courses
Optimization
Parallelization by two-dimensional subdomain partitioning
OPEN-MP and MPI
CROCO configurations CROCO
Coastal and Regional Ocean COmmunity model
Idealized conditions
• Several idealized test cases are provided
• Analytical forcing and boundary conditions
can be set
Regional configurations
• Open boundaries: active, implicit,
upstream-biased radiation conditions
• Climatological or interannual surface and
boundary conditions
• Bulk formulations for atmospheric forcing
• Rivers, and tidal forcing available
CROCO nesting CROCO
Coastal and Regional Ocean COmmunity model
Nesting with AGRIF
• Grid refinement with the AGRIF library
(developed at Inria)
• 1-way (coarse grid force the finer grid) and 2-
way (feedback of the finer grid to the coarse
grid) nesting capabilities
Towards multi-grid / multi-resolution (in dev.)
• Exchanges between grids of the same level
• Refinement criteria
• Good CPU load balance
• Management of numerous grid outputs
CROCO modules and coupling CROCO
Coastal and Regional Ocean COmmunity model
Wave-current interactions
• Based on the vortex force formalism (Uchiyama et
al. 2010):
• Impact of evolving water level on waves
OCEAN CIRCULATION
• Impact of current on waves evolution (refraction, etc)
• Wave-induced circulation (stokes drift and transport, Wave-averaged
acceleration by breaking) WAVES U,V, Tracer equations
(WKB, WW3)
• Enhanced mixing due to wave breaking
• Surface and bottom streaming (wave-induced thin Hs, T, Dir
viscous boundary layer) BOTTOM LAYER
• Mass flux due to wave rollers Bottom stress due to waves
• Wave-induced pressure effects and currents
• Wave-induced additional diffusivity DATA
• Wave-induced setup WAVES
Boundary forcing and
• WKB module breaking parameters
• Coupling with a wave model through OASIS3-MCT
library (developed at CERFACS)
CROCO modules and coupling CROCO
Coastal and Regional Ocean COmmunity model
Sediment modules
• USGS Sediment model (Blaas et al. (2007); Warner
et al. (2008))
• Wave input (specified, WKB, or WW3)
• Wave-current combined bottom stress (Soulsby, 1995)
OCEAN CIRCULATION
• Erosion (armoring), deposition, suspended transport
• Bedload transport and flux divergence Wave-averaged
• Bed model (sand, mud, or mixed) WAVES U,V, Tracer equations
(WKB, WW3)
• Morphological evolution (with acceleration factor)
• Wetting and drying Hs, T, Dir
BOTTOM LAYER
• Positive-definite advection schemes (WENO,TVD)
• Sediment influence on density Bottom stress due to waves
and currents
• MUSTANG (Mud and Sand Transport Modeling, Le DATA
Hir et al., 2011, in dev. by Ifremer/DHYSED) WAVES
SEDIMENT model
Boundary forcing and
breaking parameters
Morphodynamics Deposition, erosion, suspended and
SEDIMENTS bedload transport, bed evolution
Currents-sediment coupling (Warner et al. 2008): Sediment and bed
• Volume and constancy preserving scheme properties
• Speed-up equilibration: morpho. factor (Roelvink, 2006)
CROCO modules and coupling CROCO
Coastal and Regional Ocean COmmunity model
Ocean-wave-atmosphere coupling
• Current feedback (CFB) option available
ATSMOPHERE (WRF)
• Coupling with an atmospheric model through
OASIS3-MCT library (developed at CERFACS) OCEAN CIRCULATION
Wave-averaged
WAVES U,V, Tracer equations
(WKB, WW3)
⇒ Need to download Hs, T, Dir
BOTTOM LAYER
and compile OASIS
Bottom stress due to waves
and coupled models and currents
DATA
WAVES
SEDIMENT model
Boundary forcing and
breaking parameters
Deposition, erosion, suspended and
SEDIMENTS bedload transport, bed evolution
Sediment and bed
properties
CROCO modules and coupling CROCO
Coastal and Regional Ocean COmmunity model
Biogeochemistry
• PISCES module (Aumont and Bopp, 2006)
• BioEBUS module (Gutknecht et al., 2013)
• NPZD
Coupling with lagrangian and ecosystemic models
• ARIANE
• ICHTYOP
• OSMOSE
• APECOSM
CROCO tools and facilities CROCO
Coastal and Regional Ocean COmmunity model
Matlab CROCO_TOOLS Python CROCO_TOOLS
• Climatological pre-processing • Pre-processing: in dev.
• Interannual pre-processing • Visualization
• Visualization
Online diagnostics XIOS (dev. at ISPL)
• Online temperature / vorticity / energy balance • Outputs facilities
• Diagnostics facilities
⇒ Need to download and compile XIOS
CROCO help CROCO
Coastal and Regional Ocean COmmunity model
CROCO team CROCO
Coastal and Regional Ocean COmmunity model
CROCO training program CROCO
Coastal and Regional Ocean COmmunity model
CROCO training program CROCO
Coastal and Regional Ocean COmmunity model
