Calculate model fitness and stratification statistics

Calculate model fitness and stratification statistics for each model output.

Usage

calc_stats(
  obs,
  model,
  depths,
  folder = ".",
  NH,
  flake_nml,
  out_time,
  par_file,
  start,
  out_hour = 0
)

obs: dataframe; observation data in the LER standardised format
model: character; Model for which output will be analysed. Options include c('GOTM', 'GLM', 'Simstrat', 'FLake')
depths: vector; Depths which are to be analysed in model output.
folder: filepath; to folder which contains the model folders generated by export_config()
NH: boolean; northern hemisphere? TRUE or FALSE. Defaults to true
flake_nml: filepath; To FLake nml file. Only used if model = 'FLake'
out_time: vector; of output time values to subset data by passed to `read_flake_out`.
par_file: filepath; To Simstrat par file. Only used if model = 'Simstrat'
start: character; of start to set as origin for Simstrat output
out_hour: numeric; hour of output time values to subset data to be passed to `read_flake_out`. Only used for FLake if model time step is 86400s.