# This file was created with only commonly used settings, plus those that differ from the default specified by GOTM.
# You can generate a configuration with every possible setting with: gotm --write_yaml <OUTFILE> --detail full
# To see only the settings that differ from the default, use: gotm --write_yaml <OUTFILE> --detail minimal
version: 7                                # version of configuration file [default=7]
title: OWS Papa                           # simulation title used in output [default=GOTM simulation]
location:
   name: Ocean Weather Station Papa       # station name used in output [default=GOTM site]
   latitude: 50.0                         # latitude [degrees North; min=-90.0; max=90.0; default=0.0]
   longitude: -145.0                      # longitude [degrees East; min=-360.0; max=360.0; default=0.0]
   depth: 150.0                           # water depth [m; min=0.0; default=100.0]
time:
   start: 1961-03-25 00:00:00             # start date and time [yyyy-mm-dd HH:MM:SS; default=2017-01-01 00:00:00]
   stop: 1962-03-25 00:00:00              # stop date and time [yyyy-mm-dd HH:MM:SS; default=2018-01-01 00:00:00]
   dt: 3600.0                             # time step for integration [s; min=0.0; default=3600.0]
grid:
   nlev: 250                              # number of layers [min=1; default=100]
   method: analytical                     # layer thickness specification [analytical=equal by default with optional zooming, file_sigma=prescribed relative fractions, file_h=prescribed thicknesses; default=analytical]
   ddu: 0.0                               # surface zooming [dimensionless; min=0.0; default=0.0]
   ddl: 0.0                               # bottom zooming [dimensionless; min=0.0; default=0.0]
   file:                                  # path to file with layer thicknesses [default=]
temperature:                              # temperature profile used for initialization and optionally relaxation [Celsius]
   method: file                           # method [off, file=from file, constant, two_layer=two layers with linear gradient in between, buoyancy=from salinity and buoyancy frequency; default=off]
   constant_value: 20.0                   # value to use throughout the simulation [Celsius; default=0.0]
   file: t_prof_legacy.dat                # path to file with series of profiles [default=]
   column: 1                              # index of column to read from [default=1]
   type: in-situ                          # temperature measure [in-situ, potential, conservative; default=in-situ]
   two_layer:
      z_s: 30.0                           # depth where upper layer ends [m; min=0.0; default=0.0]
      t_s: 20.0                           # upper layer temperature [Celsius; min=-2.0; max=40.0; default=0.0]
      z_b: 40.0                           # depth where lower layer begins [m; min=0.0; default=0.0]
      t_b: 15.0                           # lower layer temperature [Celsius; min=-2.0; max=40.0; default=0.0]
   NN: 2.56000000E-04                     # square of buoyancy frequency [s^-2; min=0.0; default=0.0]
   relax:                                 # relax model temperature to observed/prescribed value
      tau: 1.00000000E+15                 # time scale for interior layer [s; min=0.0; default=1.00000000E+15]
salinity:                                 # salinity profile used for initialization and optionally relaxation [psu]
   method: file                           # method [off, file=from file, constant, two_layer=two layers with linear gradient in between, buoyancy=from temperature and buoyancy frequency; default=off]
   constant_value: 20.0                   # value to use throughout the simulation [psu; min=0.0; default=0.0]
   file: s_prof_legacy.dat                # path to file with series of profiles [default=]
   column: 1                              # index of column to read from [default=1]
   type: practical                        # salinity measure [practical, absolute; default=practical]
   two_layer:
      z_s: 30.0                           # depth where upper layer ends [m; min=0.0; default=0.0]
      s_s: 20.0                           # upper layer salinity [g/kg; min=0.0; max=40.0; default=0.0]
      z_b: 40.0                           # depth where lower layer begins [m; min=0.0; default=0.0]
      s_b: 15.0                           # lower layer salinity [g/kg; min=0.0; max=40.0; default=0.0]
   NN: 2.56000000E-04                     # square of buoyancy frequency [s^-2; min=0.0; default=0.0]
   relax:                                 # relax model salinity to observed/prescribed value
      tau: 86400.0                        # time scale for interior layer [s; min=0.0; default=1.00000000E+15]
      tau_s: 86400.0                      # time scale for surface layer [s; min=0.0; default=1.00000000E+15]
      tau_b: 86400.0                      # time scale for bottom layer [s; min=0.0; default=1.00000000E+15]
surface:
   fluxes:                                # heat and momentum fluxes
      method: off                         # method to calculate fluxes from meteorological conditions [off=use prescribed fluxes, kondo=Kondo (1975), fairall=Fairall et al. (1996); default=off]
      heat:                               # prescribed total heat flux (sensible, latent and net back-radiation) [W/m^2]
         method: file                     # method [constant, file=from file; default=constant]
         constant_value: 0.0              # value to use throughout the simulation [W/m^2; default=0.0]
         file: heatflux_legacy.dat        # path to file with time series [default=]
         column: 1                        # index of column to read from [default=1]
      tx:                                 # prescribed momentum flux in West-East direction [Pa]
         method: file                     # method [constant, file=from file; default=constant]
         constant_value: 0.0              # value to use throughout the simulation [Pa; default=0.0]
         file: momentumflux_legacy.dat    # path to file with time series [default=]
         column: 1                        # index of column to read from [default=1]
      ty:                                 # prescribed momentum flux in South-North direction [Pa]
         method: file                     # method [constant, file=from file; default=constant]
         constant_value: 0.0              # value to use throughout the simulation [Pa; default=0.0]
         file: momentumflux_legacy.dat    # path to file with time series [default=]
         column: 2                        # index of column to read from [default=1]
   u10:                                   # wind speed in West-East direction @ 10 m [m/s]
      method: file                        # method [constant, file=from file; default=constant]
      constant_value: 0.0                 # value to use throughout the simulation [m/s; default=0.0]
      file:                               # path to file with time series [default=]
      column: 1                           # index of column to read from [default=1]
   v10:                                   # wind speed in South-North direction @ 10 m [m/s]
      method: file                        # method [constant, file=from file; default=constant]
      constant_value: 0.0                 # value to use throughout the simulation [m/s; default=0.0]
      file:                               # path to file with time series [default=]
      column: 2                           # index of column to read from [default=1]
   ssuv_method: absolute                  # wind treatment [absolute=use absolute wind speed, relative=use wind speed relative to current velocity; default=relative]
   airp:                                  # air pressure [Pa]
      method: file                        # method [constant, file=from file; default=constant]
      constant_value: 0.0                 # value to use throughout the simulation [Pa; default=0.0]
      file:                               # path to file with time series [default=]
      column: 3                           # index of column to read from [default=1]
      scale_factor: 100.0                 # scale factor to be applied to values read from file [default=1.0]
   airt:                                  # air temperature @ 2 m [Celsius or K]
      method: file                        # method [constant, file=from file; default=constant]
      constant_value: 0.0                 # value to use throughout the simulation [Celsius or K; default=0.0]
      file:                               # path to file with time series [default=]
      column: 4                           # index of column to read from [default=1]
   hum:                                   # humidity @ 2 m
      method: file                        # method [constant, file=from file; default=constant]
      constant_value: 0.0                 # value to use throughout the simulation [default=0.0]
      file:                               # path to file with time series [default=]
      column: 5                           # index of column to read from [default=1]
      type: relative                      # humidity metric [relative=relative humidity (%), wet_bulb=wet-bulb temperature, dew_point=dew point temperature, specific=specific humidity (kg/kg); default=relative]
   cloud:                                 # cloud cover [1]
      method: file                        # method [constant, file=from file; default=constant]
      constant_value: 0.0                 # value to use throughout the simulation [fraction; min=0.0; max=1.0; default=0.0]
      file:                               # path to file with time series [default=]
      column: 6                           # index of column to read from [default=1]
   precip:                                # precipitation [m/s]
      method: constant                    # method [constant, file=from file; default=constant]
      constant_value: 0.0                 # value to use throughout the simulation [m/s; default=0.0]
      file:                               # path to file with time series [default=]
      column: 1                           # index of column to read from [default=1]
      flux_impact: false                  # include effect on fluxes of sensible heat and momentum [default=false]
   calc_evaporation: false                # calculate evaporation from meteorological conditions [default=false]
   swr:                                   # shortwave radiation [W/m^2]
      method: file                        # method [constant, file=from file, calculate=from time, location and cloud cover; default=constant]
      constant_value: 0.0                 # value to use throughout the simulation [W/m^2; min=0.0; default=0.0]
      file: swr_legacy.dat                # path to file with time series [default=]
      column: 1                           # index of column to read from [default=1]
   longwave_radiation:                    # net longwave radiation [W/m^2]
      method: clark                       # method [file=from file, clark=Clark et al. (1974), hastenrath_lamb=Hastenrath and Lamb (1978), bignami=Bignami et al. (1995), berliand_berliand=Berliand and Berliand (1952), josey1=Josey et al. (2003) - 1, josey2=Josey et al. (2003) - 2; default=clark]
      file:                               # path to file with time series [default=]
      column: 1                           # index of column to read from [default=1]
   albedo:
      method: constant                    # method to compute albedo [constant, payne=Payne (1972), cogley=Cogley (1979); default=payne]
      constant_value: 0.0                 # constant value to use throughout the simulation [fraction; min=0.0; max=1.0; default=0.0]
   sst:                                   # observed surface temperature [Celsius]
      method: file                        # method [constant, file=from file; default=constant]
      file: sst_legacy.dat                # path to file with time series [default=]
   ice:
      model: no_ice                       # model [no_ice=no ice model, simple=freezing temperature, basal_melt=Basal Melt , lebedev=Lebedev (1938), mylake=Saloranta and Andersen (2007), winton=Michael Winton (2000); default=simple]
      H: 0.0                              # initial ice thickness [m; default=0.0]
   roughness:
      charnock: false                     # use Charnock (1955) roughness adaptation [default=false]
      charnock_val: 1400.0                # empirical constant for roughness adaptation [dimensionless; min=0.0; default=1400.0]
      z0s_min: 0.02                       # hydrodynamic roughness (minimum value if Charnock adaptation is used) [m; min=0.0; default=0.02]
bottom:
   h0b: 0.05                              # physical bottom roughness [m; min=0.0; default=0.05]
light_extinction:
   method: jerlov-i                       # water type [jerlov-i=Jerlov type I, jerlov-1-50m=Jerlov type 1 (upper 50 m), jerlov-ia=Jerlov type IA, jerlov-ib=Jerlov type IB, jerlov-ii=Jerlov type II, jerlov-iii=Jerlov type III, custom; default=jerlov-i]
   A:                                     # non-visible fraction of shortwave radiation [1]
      method: constant                    # method [constant, file=from file; default=constant]
      constant_value: 0.7                 # value to use throughout the simulation [fraction; min=0.0; max=1.0; default=0.7]
      file:                               # path to file with time series [default=]
      column: 1                           # index of column to read from [default=1]
   g1:                                    # e-folding depth of non-visible shortwave radiation [m]
      method: constant                    # method [constant, file=from file; default=constant]
      constant_value: 0.4                 # value to use throughout the simulation [m; min=0.0; default=0.4]
      file:                               # path to file with time series [default=]
      column: 1                           # index of column to read from [default=1]
   g2:                                    # e-folding depth of visible shortwave radiation [m]
      method: constant                    # method [constant, file=from file; default=constant]
      constant_value: 8.0                 # value to use throughout the simulation [m; min=0.0; default=8.0]
      file:                               # path to file with time series [default=]
      column: 1                           # index of column to read from [default=1]
turbulence:
   turb_method: second_order              # turbulence closure [no_model=constant nuh and num, first_order=first-order, second_order=second-order, cvmix; default=second_order]
   tke_method: tke                        # turbulent kinetic energy equation [local_eq=algebraic length scale equation, tke=differential equation for tke (k-eps or k-w style), mellor_yamada=differential equation for q^2/2 (Mellor-Yamada style); default=tke]
   len_scale_method: dissipation          # dissipative length scale [parabolic, triangular, xing_davies=Xing and Davies (1995), robert_ouellet=Robert and Ouellet (1987), blackadar=Blackadar (two boundaries) (1962), bougeault_andre=Bougeault and Andre (1986), dissipation=dynamic dissipation rate equation, mellor_yamada=dynamic Mellor-Yamada q^2 l-equation, gls=generic length scale (GLS), omega=dynamic omega equation; default=dissipation]
   stab_method: constant                  # stability functions [constant, munk_anderson=Munk and Anderson (1954), schumann_gerz=Schumann and Gerz (1995); default=constant]
   bc:                                    # boundary conditions
      k_ubc: neumann                      # upper boundary condition for k-equation [dirichlet, neumann; default=neumann]
      k_lbc: neumann                      # lower boundary condition for k-equation [dirichlet, neumann; default=neumann]
      psi_ubc: neumann                    # upper boundary condition for length-scale equation [dirichlet, neumann; default=neumann]
      psi_lbc: neumann                    # lower boundary condition for length-scale equation [dirichlet, neumann; default=neumann]
      ubc_type: logarithmic               # upper boundary layer [logarithmic=logarithmic law of the wall, tke_injection=tke-injection (breaking waves); default=logarithmic]
      lbc_type: logarithmic               # lower boundary layer [logarithmic=logarithmic law of the wall; default=logarithmic]
   turb_param:
      cm0_fix: 5.47700000E-01             # value of the stability function in the log-law [dimensionless; min=0.0; default=5.47700000E-01]
      Prandtl0_fix: 0.74                  # turbulent Prandtl-number [dimensionless; min=0.0; default=0.74]
      cw: 100.0                           # constant of the wave-breaking model [dimensionless; min=0.0; default=100.0]
      compute_kappa: true                 # compute von Karman constant from model parameters [default=true]
      kappa: 0.4                          # von Karman constant [dimensionless; min=0.0; default=0.4]
      compute_c3: true                    # compute c3 (E3 for Mellor-Yamada) from steady-state Richardson number [default=true]
      Ri_st: 0.25                         # desired steady-state Richardson number [dimensionless; min=0.0; default=0.25]
      length_lim: true                    # apply Galperin et al. (1988) length scale limitation [default=true]
      galp: 0.27                          # coefficient for length scale limitation [dimensionless; min=0.0; default=0.27]
      const_num: 5.00000000E-04           # constant eddy diffusivity [m^2/s; min=0.0; default=5.00000000E-04]
      const_nuh: 5.00000000E-04           # constant heat diffusivity [m^2/s; min=0.0; default=5.00000000E-04]
      k_min: 1.00000000E-06               # minimum turbulent kinetic energy [m^2/s^2; min=0.0; default=1.00000000E-10]
      eps_min: 1.00000000E-12             # minimum dissipation rate [m^2/s^3; min=0.0; default=1.00000000E-12]
      kb_min: 1.00000000E-10              # minimum buoyancy variance [m^2/s^4; min=0.0; default=1.00000000E-10]
      epsb_min: 1.00000000E-14            # minimum buoyancy variance destruction rate [m^2/s^5; min=0.0; default=1.00000000E-14]
   generic:                               # generic length scale (GLS) model
      compute_param: false                # compute the model parameters [default=false]
      gen_m: 1.5                          # exponent for k [dimensionless; default=1.5]
      gen_n: -1.0                         # exponent for l [dimensionless; default=-1.0]
      gen_p: 3.0                          # exponent for cm0 [dimensionless; default=3.0]
      cpsi1: 1.44                         # empirical coefficient cpsi1 in psi equation [dimensionless; default=1.44]
      cpsi2: 1.92                         # empirical coefficient cpsi2 in psi equation [dimensionless; default=1.92]
      cpsi3minus: 0.0                     # cpsi3 for stable stratification [dimensionless; default=0.0]
      cpsi3plus: 1.0                      # cpsi3 for unstable stratification [dimensionless; default=1.0]
      cpsi4: 0.0                          # empirical coefficient cpsi4 in psi equation [dimensionless; default=0.0]
      sig_kpsi: 1.0                       # Schmidt number for TKE diffusivity [dimensionless; default=1.0]
      sig_psi: 1.3                        # Schmidt number for psi diffusivity [dimensionless; default=1.3]
      gen_d: -1.2                         # temporal decay rate d in homogeneous turbulence [dimensionless; default=-1.2]
      gen_alpha: -2.0                     # decay exponent alpha [dimensionless; default=-2.0]
      gen_l: 0.2                          # slope L of the length scale in shear-free turbulence [dimensionless; default=0.2]
   keps:                                  # k-epsilon model
      ce1: 1.44                           # empirical coefficient ce1 in dissipation equation [dimensionless; default=1.44]
      ce2: 1.92                           # empirical coefficient ce2 in dissipation equation [dimensionless; default=1.92]
      ce3minus: -0.4                      # ce3 for stable stratification [dimensionless; default=0.0]
      ce3plus: 1.5                        # ce3 for unstable stratification [dimensionless; default=1.5]
      ce4: 0.0                            # empirical coefficient ce4 in dissipation equation [dimensionless; default=0.0]
      sig_k: 1.0                          # Schmidt number for TKE diffusivity [dimensionless; default=1.0]
      sig_e: 1.3                          # Schmidt number for dissipation diffusivity [dimensionless; default=1.3]
      sig_peps: false                     # use Burchard (2001) wave breaking parameterisation [default=false]
   kw:                                    # k-omega model
      cw1: 0.555                          # empirical coefficient cw1 in omega equation [dimensionless; default=0.555]
      cw2: 0.833                          # empirical coefficient cw2 in omega equation [dimensionless; default=0.833]
      cw3minus: 0.0                       # cw3 for stable stratification [dimensionless; default=0.0]
      cw3plus: 0.5                        # cw3 for unstable stratification [dimensionless; default=0.5]
      cw4: 0.15                           # empirical coefficient cw4 in omega equation [dimensionless; default=0.15]
      sig_kw: 2.0                         # Schmidt number for TKE diffusivity [dimensionless; default=2.0]
      sig_w: 2.0                          # Schmidt number for omega diffusivity [dimensionless; default=2.0]
   my:                                    # Mellor-Yamada model
      e1: 1.8                             # coefficient e1 in q^2 l equation [dimensionless; default=1.8]
      e2: 1.33                            # coefficient e2 in q^2 l equation [dimensionless; default=1.33]
      e3: 1.8                             # coefficient e3 in q^2 l equation [dimensionless; default=1.8]
      e6: 4.0                             # coefficient e6 in q^2 l equation [dimensionless; default=4.0]
      sq: 0.2                             # turbulent diffusivities of q^2 (= 2k) [dimensionless; default=0.2]
      sl: 0.2                             # turbulent diffusivities of q^2 l [dimensionless; default=0.2]
      length: linear                      # barotropic length scale in q^2 l equation [parabolic, triangular, linear=linear from surface; default=parabolic]
   scnd:                                  # second-order model
      method: quasi_eq                    # method [quasi_eq=quasi-equilibrium, weak_eq_kb_eq=weak equilibrium with algebraic buoyancy variance, quasi_eq_h15=quasi-equilibrium with Langmuir (Harcourt, 2015); default=weak_eq_kb_eq]
      scnd_coeff: canuto-a                # coefficients of second-order model [custom, gibson_launder=Gibson and Launder (1978), mellor_yamada=Mellor and Yamada (1982), kantha_clayson=Kantha and Clayson (1994), luyten=Luyten et al. (1996), canuto-a=Canuto et al. (2001) (version A), canuto-b=Canuto et al. (2001) (version B), cheng=Cheng et al. (2002); default=canuto-a]
   iw:                                    # internal wave mixing
      method: large                       # method [off, mellor=Mellor (1989), large=Large et al. (1994); default=off]
      alpha: 0.7                          # coefficient for Mellor internal wave model [dimensionless; default=0.0]
mimic_3d:                                 # effects of horizontal gradients
   ext_pressure:                          # external pressure
      type: elevation                     # pressure metric [elevation=horizontal gradient in surface elevation, velocity=horizontal velocities at given height above bed, average_velocity=vertically averaged horizontal velocities; default=elevation]
      dpdx:                               # pressure in West-East direction
         method: constant                 # method [constant, tidal=from tidal constituents, file=from file; default=constant]
         constant_value: 0.0              # value to use throughout the simulation [default=0.0]
         file:                            # path to file with time series [default=]
         column: 1                        # index of column to read from [default=1]
         tidal:                           # tidal constituents
            amp_1: 0.0                    # amplitude of 1st harmonic [default=0.0]
            phase_1: 0.0                  # phase of 1st harmonic [s; default=0.0]
            amp_2: 0.0                    # amplitude of 2nd harmonic [default=0.0]
            phase_2: 0.0                  # phase of 2nd harmonic [s; default=0.0]
      dpdy:                               # pressure in South-North direction
         method: constant                 # method [constant, tidal=from tidal constituents, file=from file; default=constant]
         constant_value: 0.0              # value to use throughout the simulation [default=0.0]
         file:                            # path to file with time series [default=]
         column: 1                        # index of column to read from [default=1]
         tidal:                           # tidal constituents
            amp_1: 0.0                    # amplitude of 1st harmonic [default=0.0]
            phase_1: 0.0                  # phase of 1st harmonic [s; default=0.0]
            amp_2: 0.0                    # amplitude of 2nd harmonic [default=0.0]
            phase_2: 0.0                  # phase of 2nd harmonic [s; default=0.0]
      h:                                  # height above bed [m]
         method: constant                 # method [constant, file=from file; default=constant]
         constant_value: 0.0              # value to use throughout the simulation [m; min=0.0; default=0.0]
         file:                            # path to file with time series [default=]
         column: 1                        # index of column to read from [default=1]
      period_1: 44714.0                   # period of 1st tidal harmonic (eg. M2-tide) [s; default=44714.0]
      period_2: 43200.0                   # period of 2nd tidal harmonic (eg. S2-tide) [s; default=43200.0]
   int_pressure:                          # internal pressure
      type: none                          # method [none, gradients=prescribed horiztonal gradients of T and S, plume=surface or bottom plume; default=none]
      gradients:                          # horizontal salinity and temperature gradients
         dtdx:                            # temperature gradient in West-East direction [Celsius/m]
            method: off                   # method [off, constant, file=from file; default=off]
            constant_value: 0.0           # value to use throughout the simulation [Celsius/m; default=0.0]
            file:                         # path to file with series of profiles [default=]
            column: 1                     # index of column to read from [default=1]
         dtdy:                            # temperature gradient in South-North direction [Celsius/m]
            method: off                   # method [off, constant, file=from file; default=off]
            constant_value: 0.0           # value to use throughout the simulation [Celsius/m; default=0.0]
            file:                         # path to file with series of profiles [default=]
            column: 1                     # index of column to read from [default=1]
         dsdx:                            # salinity gradient in West-East direction [psu/m]
            method: off                   # method [off, constant, file=from file; default=off]
            constant_value: 0.0           # value to use throughout the simulation [psu/m; default=0.0]
            file:                         # path to file with series of profiles [default=]
            column: 1                     # index of column to read from [default=1]
         dsdy:                            # salinity gradient in South-North direction [psu/m]
            method: off                   # method [off, constant, file=from file; default=off]
            constant_value: 0.0           # value to use throughout the simulation [psu/m; default=0.0]
            file:                         # path to file with series of profiles [default=]
            column: 1                     # index of column to read from [default=1]
      plume:                              # surface or bottom plume
         type: bottom                     # plume type [surface=buoyant surface-attached, bottom=dense bottom-attached; default=bottom]
         x_slope: 0.0                     # plume slope in West-East direction [dimensionless; default=0.0]
         y_slope: 0.0                     # plume slope in South-North direction [dimensionless; default=0.0]
      t_adv: false                        # horizontally advect temperature [default=false]
      s_adv: false                        # horizontally advect salinity [default=false]
   zeta:                                  # surface elevation [m]
      method: constant                    # method [constant, tidal=from tidal constituents, file=from file; default=constant]
      constant_value: 0.0                 # value to use throughout the simulation [m; default=0.0]
      file:                               # path to file with time series [default=]
      column: 1                           # index of column to read from [default=1]
      tidal:                              # tidal constituents
         period_1: 44714.0                # period of 1st harmonic (eg. M2-tide) [s; default=44714.0]
         amp_1: 1.0                       # amplitude of 1st harmonic [m; default=0.0]
         phase_1: 0.0                     # phase of 1st harmonic [s; default=0.0]
         period_2: 43200.0                # period of 2nd harmonic (eg. S2-tide) [s; default=43200.0]
         amp_2: 0.5                       # amplitude of 2nd harmonic [m; default=0.0]
         phase_2: 0.0                     # phase of 2nd harmonic [s; default=0.0]
   w:                                     # vertical velocity
      max:                                # maximum velocity [m/s]
         method: off                      # method [off, constant, file=from file; default=off]
         constant_value: 0.0              # value to use throughout the simulation [m/s; default=0.0]
         file:                            # path to file with time series [default=]
         column: 1                        # index of column to read from [default=1]
      height:                             # height of maximum velocity [m]
         method: constant                 # method [constant, file=from file; default=constant]
         constant_value: 0.0              # value to use throughout the simulation [m; default=0.0]
         file:                            # path to file with time series [default=]
         column: 1                        # index of column to read from [default=1]
      adv_discr: p2_pdm                   # vertical advection scheme [upstream=first-order upstream, p2=third-order upstream-biased polynomial, superbee=third-order TVD with Superbee limiter, muscl=third-order TVD with MUSCL limiter, p2_pdm=third-order TVD with ULTIMATE QUICKEST limiter; default=p2_pdm]
o2:                                       # oxygen
   scale_factor: 31.25                    # scale factor to be applied to values read from file [default=1.0]
cvmix:
   surface_layer:                         # surface layer mixing
      use: true                           # compute surface layer mixing coefficients [default=true]
      langmuir_method: none               # method of Langmuir turbulence pararmeterization [none, lwf16=Li et al. (2016), lf17=Li and Fox-Kemper (2017), rwh16=Reichl et al. (2016); default=none]
      surface_layer_extent: 0.1           # extent of surface layer in fraction of the boundary layer [dimensionless; min=0.0; max=1.0; default=0.1]
      Ri_c: 0.3                           # critical Richardson number [dimensionless; min=0.0; default=0.3]
      check_Ekman_length: false           # limit the OBL by the Ekman depth [default=false]
      check_MonOb_length: false           # limit the OBL by the Monin-Obukhov depth [default=false]
      use_enhanced_diff: true             # enhance diffusivity at OBL [default=true]
      use_noDGat1: true                   # zero gradient of the shape function at OBL [default=true]
      match_technique: simple             # matching technique of shape functions with the ocean interior [simple=simple shapes, gradient=gradient term, both=both gradient and nonlocal terms, parabolic=parabolic nonlocal term; default=simple]
      bulk_Ri_interp_type: quadratic      # interpolation type for the bulk Richardson number [linear, quadratic, cubic; default=quadratic]
      OBL_interp_type: lmd94              # interpolation type for diffusivity and viscosity at OBL [linear, quadratic, cubic, lmd94=Large et al. (1994); default=lmd94]
   bottom_layer:                          # bottom layer mixing
      use: false                          # compute bottom layer mixing coefficients [default=false]
      surface_layer_extent: 0.1           # extent of surface layer in fraction of the boundary layer [dimensionless; min=0.0; max=1.0; default=0.1]
      Ri_c: 0.3                           # critical Richardson number [dimensionless; min=0.0; default=0.3]
      check_Ekman_length: false           # limit the OBL by the Ekman depth [default=false]
      check_MonOb_length: false           # limit the OBL by the Monin-Obukhov depth [default=false]
      use_noDGat1: true                   # zero gradient of the shape function at OBL [default=true]
      match_technique: simple             # matching technique of shape functions with the ocean interior [simple=simple shapes, gradient=gradient term, both=both gradient and nonlocal terms, parabolic=parabolic nonlocal term; default=simple]
      bulk_Ri_interp_type: quadratic      # interpolation type for the bulk Richardson number [linear, quadratic, cubic; default=quadratic]
      OBL_interp_type: lmd94              # interpolation type for diffusivity and viscosity at OBL [linear, quadratic, cubic, lmd94=Large et al. (1994); default=lmd94]
   interior:                              # interior mixing
      use: false                          # compute interior mixing coefficients [default=false]
      background:
         use: false                       # use interior background mixing coefficients [default=false]
         diffusivity: 1.00000000E-05      # background diffusivity [m^2/s; default=1.00000000E-05]
         viscosity: 1.00000000E-04        # background viscosity [m^2/s; default=1.00000000E-04]
      shear:
         use: false                       # compute interior shear mixing coefficients [default=false]
         num_smooth_Ri: 1                 # number of iterations to smooth the gradient Richardson number [default=1]
         mix_scheme: kpp                  # shear mixing scheme [pp=Pacanowski and Philander (1981), kpp=Large et al. (1994); default=kpp]
         PP_nu_zero: 0.005                # numerator in viscosity term in PP [m^2/s; default=0.005]
         PP_alpha: 5.0                    # coefficient of Ri in denominator of visc / diff terms [dimensionless; default=5.0]
         PP_exp: 2.0                      # exponent of denominator in viscosity term [dimensionless; default=2.0]
         KPP_nu_zero: 0.005               # leading coefficient of the KPP shear mixing formula [m^2/s; default=0.005]
         KPP_Ri_zero: 0.7                 # critical Richardson number for KPP shear mixing [dimensionless; default=0.7]
         KPP_exp: 3.0                     # exponent of unitless factor of diffusivity [dimensionless; default=3.0]
      convection:
         use: false                       # compute interior convective mixing coefficients [default=false]
         diffusivity: 1.0                 # convective diffusivity [m^2/s; default=1.0]
         viscosity: 1.0                   # convective viscosity [m^2/s; default=1.0]
         basedOnBVF: true                 # triger convection based on the squared Brunt-Vaisala frequency [default=true]
         triggerBVF: 0.0                  # threshold of squared Brunt-Vaisala frequency [1/s^2; default=0.0]
      tidal_mixing:
         use: false                       # compute interior tidal mixing coefficients [default=false]
      double_diffusion:
         use: false                       # compute interior double diffusion mixing coefficients [default=false]
fabm:                                     # Framework for Aquatic Biogeochemical Models
   use: false                             # enable FABM [default=false]
   freshwater_impact: true                # enable dilution/concentration by precipitation/evaporation [default=true]
   feedbacks:                             # feedbacks to physics
      shade: false                        # interior light absorption [default=false]
      albedo: false                       # surface albedo [default=false]
      surface_drag: false                 # surface drag [default=false]
   repair_state: false                    # clip state to minimum/maximum boundaries [default=false]
   input:
equation_of_state:                        # equation of state
   method: full_teos-10                   # density formulation [full_teos-10=TEOS-10, linear_teos-10=linearized from user-specified T0, S0, p0 (alpha, beta, cp calculated), linear_custom=linearized from user-specified T0, S0, rho0, alpha, beta, cp; default=full_teos-10]
   rho0: 1027.0                           # reference density [kg/m3; default=1027.0]
   linear:
      T0: 10.0                            # reference temperature [Celsius; min=-2.0; default=15.0]
      S0: 35.0                            # reference salinity [g/kg; min=0.0; default=35.0]
      p0: 0.0                             # reference pressure [dbar; min=0.0; default=0.0]
      alpha: 2.00000000E-04               # thermal expansion coefficient [1/K; default=2.00000000E-04]
      beta: 7.50000000E-04                # saline contraction coefficient [kg/g; default=7.50000000E-04]
      cp: 3.99186796E+03                  # specific heat capacity [J/(kg/K); min=0.0; default=3.99186796E+03]
seagrass:                                 # calculate seagrass effect on turbulence
   method: 0                              #  [0=off, 1=from file; default=0]
   file: seagrass.dat                     # path to file with grass specifications [default=seagrass.dat]
   alpha: 0.0                             # efficiency of leafes turbulence production [default=0.0]
restart:
   load: false                            # initialize simulation with state stored in restart.nc [default=false]
output:
   ows_papa_keps_legacy:                  # path of output file, excluding extension
      time_unit: dt                       # time unit [second, hour, day, month, year, dt=model time step; default=day]
      time_step: 24                       # number of time units between output [min=1; default=1]
      time_method: point                  # treatment of time dimension [point=instantaneous, mean, integrated; default=point]
      variables:
      - source: /*                        # variable name in model