Assumption initial conditions are on z-levels
The initial conditions for runs are assumed to be on z-levels based on gdept_1d
. This is not immediately apparent from reading the manual. This is a particular issue when running with ln_sco=.true.
as the user would have to update gdept_1d
in the domain_cfg.nc file to match that of the initial condition vertical levels.
A few options would be to either:
- modify code such that depth information is required in the initial condition file. This would allow an initial condition file, comprising a different number of levels and/or grid type/or depth levels, to be read in and interpolated onto
gdept_0
. - remove the interpolation loop completely (bearing in mind the on-going simplification process) and impose that the initial conditions read in have to be on
gdept_0
, requiring a pre-processing step.
In either case, explicit reference to the method employed needs to be made in section 3.3 of the reference manual.