…ifveR (vertical expl) and difviR (vertical impl)

…bugfix_hnode_in_Trz_and_difveR_which_tr_in_Trxy

…ix_hnode_in_Trz_and_difveR_which_tr_in_Trxy

Initialize recom variables, modify BGC tracers, call recom

Recom related tracer fields are generated by fesom.
The first step in recom is to make sure that we have correct units
and initial values. This is done in recom_init and after the ocean_setup.

… come in an separate commit

…ifveR_which_tr_in_Trxy

delete unnecessary variable

…es_diag.F90, src/oce_ale_tracer.F90, config/namelist.cvmix

Read the file names and associated variables from the namelist.recom

Atmospheric CO2, Fe and N deposition,
riverine input of alkalinity, DIC, DIN, DOC, DON, Fe
and erosion as PON and POC are read.

…order, if it is not setted in namelist.tra

…cer gradients. Its was decided to use tr^n, which means that the tr_AB can be kicked out from the DVD computation

…e first ones isnt need anymore due to the AB2/3 order update, the later one isnt need anymore due to the fixed bug from tracer synchronisation of the tracer gradient computation

…esold(1, :, :), therefor allocate new array for trstar= (trold + tr)*0.5_WP

…issolved Iron)

Dissolved iron concentration is limited by the cut-off value 0.6 (umol/m3)
in deep Eastern Equatorial Pacific.

High concentration of dissolved iron in deep ocean disperses
in all directions quickly. After several decades, mid-ocean iron
concentration reaches higher levels.
Exclude the depths shallower than 2000 m.

Read the correct file for DIC during initialisation.
Keep the initial surface ocean alkalinity in Alk_surf.

There is an option to read preindustrial or modern
Dissolved Inorganic Carbon (DIC). It is controlled
via DIC_PI switch in namelist.recom.
Alk_surf holds the initial state of surface ocean alkalinity.
Iterated alkalinity field is nudged to Alk_surf to keep
the model drift plausible.

Add REcoM and MOCSY2.0 routines for the
ecosystem and carbonate chemistry in the ocean

Basic coupled FESOM2.5-REcoM setup with 2 phytoplankton,
1 zooplankton and 1 detritus functional types. Tracers are
restructured to be operated within water column.
CO2 and O2 fluxes are calculated in recom_forcing.F90
using MOCSY2.0 routines.recom_sms.F90 has the sources
minus sink terms.

BREAKING CHANGE: FESOM2.5-REcoM setup with 2p1z1d

Use Diags in namelist.recom to compute diagnostics

Several 2D and 3D fields are saved to monitor the
model performance. 2D fields include net primary production,
respiration, chlorophyll of phytoplankton groups. 3D fields include
PAR, calcite dissolution, calcification as well as aggregation,
respiration and net primary production of small phytoplankton/diatoms
in the water column are computed.

Add second and third zooplankton classes in addition to
second detritus class to REcoM

Standard REcoM version with 2 phytoplankton (small pytoplankton
and diatoms), 3 zooplankton (micro-, meso- and macrozooplankton)
and 2 detritus (slow- and fast-sinking) functional types.
Vertical sinking for fast-sinking detritus is assumed to be constant
whereas slow-sinking detritus variable sinking has variable sinking
option.

BREAKING CHANGE: FESOM2.5-REcoM version with 2p3z2d

Add third phytoplankton (coccolithophores) and associated
output variables

Standard REcoM version with 3 phytoplankton (small pytoplankton,
diatoms and coccolithophores), 3 zooplankton (micro-, meso- and
macrozooplankton) and 2 detritus (slow- and fast-sinking) functional types.
__coccos decides the background vertical sinking speed of
intracellular calcite concentration (VCocco  in case of __coccos
otherwise VPhy) including depth dependency. The list of output files
is  extended with coccolithophores related variables.

Updates 2.5.0.4 de

CO2 dependence of photosynthesis rate, the formation and
dissolution of CaCO3 by coccolithophores and oxygen dependency
of organic matter remineralisation.

Use namelist.recom to activate ballasting, CO2lim, OmegaC_diss
and O2dep_remin. Carbonate chemistry is computed at every tracer
point.

Sinking velocity of particle (fast and slow) is scaled by particle density and/or
water viscosity (optional). Ocean carbonate chemistry is calculated
in/below the euphoric zone weekly/monthly using  MOCSY2.0
routines. The output variables are used to compute
OmegaC dependent calcite dissolution (before depth dependent) and
CO2 dependence on the photosynthesis rate of phytoplankton groups.
Organic matter remineralisation rate is not only depend on temperature
but also O2 concentration. Namelist.recom in config/bin_3p3z2d folder
is updated.

…is the improved sea ice thermodynamics for coupling with IFS/OpenIFS

Fix circular module dependency

Argument mismatch allowed on GNU Compiler on Ubuntu

This removes the most obvious traces and builds, no real test run done yet

Production de merge into refactoring

Refactoring better metadata

change some paths, remove ollie and mistral as HPCs are long gone

 stream for the forcing output bundle was defined twice

Fix FESOM compilation on Cray compilers

fix landice fw bug

fix missing/miss placed allocation

…gistry

Hey folks! I've reworked AB_tracer. Now it initializes by default, so even if you forget to set it in the namelist, it won't run wild with random values. We've identified the chaos monster!

setting default for tracers/AB_order = 2 in oce_setup_step.F90

…r_mod.F90

drop an error message for ridiculous values for AB_order in oce_tracer_mod.F90

drop an error message for ridiculous values for AB_order in oce_trace…

remove duplicate info