4_model_forecast.yml

target_default: 4_model_forecast

packages:
  - yaml
  - dplyr
  - tidyr
  - scipiper

sources:
  - 4_model/src/EnKF_functions.R
  - 4_model/src/get_sntemp_values.R
  - 4_model/src/run_sntemp.R
  - 4_model/src/update_sntemp.R
  - 4_model/src/set_sntemp_output.R
  - 4_model_forecast/src/forecast.R


targets:
  4_model_forecast:
    depends:
      - forecast_horizon

  forecast_settings:
    command: read_yaml('4_model_forecast/cfg/forecast_settings.yml')

  forecast_horizon:
    command: forecast_settings[[I('forecast_horizon')]]

  forecast_start:
    command: forecast_settings[[I('start')]]

  forecast_stop:
    command: forecast_settings[[I('stop')]]

  forecast_param_groups:
    command: forecast_settings[[I('param_groups')]]

  forecast_n_en:
    command: forecast_settings[[I('n_en')]]

  # applies to a set of forecasts - equivalent to a DOI and shouldn't change unless algorithms change
  forecast_project_id:
    command: forecast_settings[[I('forecast_project_id')]]

  2_3_model_parameters/out/forecast_params.nc:
    command: nc_create_cal_params(
      n_en = forecast_n_en,
      forecast_project_id = forecast_project_id,
      vars = forecast_param_groups,
      nc_name_out = target_name,
      model_run_loc = I('4_model_calibrate/tmp'),
      param_default_file = I('control/delaware.control.par_name'))

  # model_spinup:
  #   command: model_spinup(
  #     n_en = n_en,
  #     start = start,
  #     stop = stop,
  #     time_step = I('days'),
  #     model_run_loc = I('4_model_forecast/tmp'),
  #     spinup_days = I(730))
  #
  # 4_model/out/model_out_gwsum_sssum_subbasin_4182_inf_factor.rds.ind:
  #   command: EnKF(
  #     ind_file = target_name,
  #     n_en = n_en,
  #     start = start,
  #     stop = stop,
  #     time_step = I('days'),
  #     subbasin_file = '4_model_calibrate/out/drb_subbasins.rds',
  #     subbasin_outlet_file = '4_model_calibrate/cfg/subbasin_outlets.yml',
  #     subbasin_outlet_id = I('4182'),
  #     model_fabric_file = '20191002_Delaware_streamtemp/GIS/Segments_subset.shp',
  #     obs_file = '3_observations/in/obs_temp_full.rds',
  #     init_param_file = '2_3_model_parameters/out/init_params.rds',
  #     model_run_loc = I('4_model/tmp'),
  #     orig_model_loc = I('20200207_Delaware_streamtemp_state_adj'),
  #     state_names = state_names,
  #     n_params_obs = I(0),
  #     obs_cv = I(0.05),
  #     param_cv = I(0.2),
  #     init_cond_cv = I(0.1),
  #     covar_inf_factor = TRUE)