This directory contains runscripts applied for optimisations ECCO*
(directory OptECCO) and optimisations UHigh*, U20* and U17.5* (directory
OptUVic).

Iris Kriest, 14 January 2020

In each subdirectory is a runscript for the optimiser CMAES
("serial.job"), which is the same for all optimisations, and for 10
parallel model simulations ("parallel.job), which differ.

Differences in "parallel.job" are related to 

(1) Timestepping:

      ECCO* uses 720 ocean timesteps per year (for circulation) resulting 
      in an ocean timestep length of 1/720 (years) = 0.0013888888888889 = deltat_clock.
      For biogeochemistry we give the timestep length in seconds: 
      1/720 (y) * 360 (days/y) * 86400 (seconds/day) = 43200 = -biogeochem_deltat
      Length of simulation is 3000 (years) * 720 (timesteps/year) = 2160000 = max_steps.
      A snapshot of the simulation will be written at the end: -write_steps 2160000
      Time averaging for annual average tracers and diagnostic output is over a year 
      (-avg_time_steps = -diag_timesteps = 720), 
      starting in year 2999 (-avg_start_time_step, -diag_start_time_step).
      Burial is intergrated over one year (-burial_sum_steps 720).
      MOPS in ECCO applies 8 biogeochemical time steps per ocean time step 
      (-num_biogeochem_steps_per_ocean_step 8), i.e. 5760 BGC time steps in a year or 16 time steps in a day.
       
      UHigh* (U20*, U17.5*) uses 1095 ocean timesteps per year (for circulation) resulting 
      in an ocean timestep length of 1/1095 (years) = 0.0009132420091324 = deltat_clock.
      For biogeochemistry we give the timestep length in seconds: 
      1/1095 (y) * 365 (days/y) * 86400 (seconds/day) = 28800 = -biogeochem_deltat
      Length of simulation is 3000 (years) * 1095 (timesteps/year) = 3285000 = max_steps.
      A snapshot of the simulation will be written at the end: -write_steps 3285000
      Time averaging for annual average tracers and diagnostic output is over a year 
      (-avg_time_steps = -diag_timesteps = 1095), 
      starting in year 2999 (-avg_start_time_step, -diag_start_time_step).
      Burial is intergrated over one year (-burial_sum_steps 1095).
      MOPS in UVOK applies 5 biogeochemical time steps per ocean time step 
      (-num_biogeochem_steps_per_ocean_step 5), i.e. 5475 BGC time steps in a year, or 15 time steps in a day.
      
      The new driver code used for UHigh, U20, U17.5 
            - allows to prescribe the year length (-days_per_year 365.0)
            - instead of using -matrix_cycle_step 0.0833333333333333 it applies -matrix_num_per_period 12
            - instead of -periodic_biogeochem_cycle_step 0.08333333333333 is applies -periodic_biogeochem_num_per_period 12 
            - uses a  different nomenclature for the misfit time stepper: 
              instead of -cost_start_time_step and -cost_time_steps is applies 
              -misfit_start_time_step and -misfit_time_steps 
            - does not take the option -average_cost any more

(2) Vertical domains: MOPS in 

      ECCO* considers the upper 6 vertical layers as euphotic zone (-nzeuph 6) (= upper 100m)                        
      UHigh* (U20*, U17.5*) considers the upper 2 vertical layers as euphotic zone (-nzeuph 2) (= upper 130m)                        
            

Other changes relate to the naming of the input files, executable and computational ressources.