2_observations.yml

target_default: 2_observations

packages:
  - yaml
  - whisker
  - dplyr
  - tidyr
  - meddle
  - readr
  - rgdal
  - stringr
  - glmtools
  - sf


sources:
  - src/file_utils.R

targets:
  2_observations:
    depends:
      - out_data/lake_surface_temp_obs_preds.csv

  # add all_outputs_and_obs_error_estimate.csv, obs and pred for all three models

  lake_surface_temp_preds:
    command: convert_preds_tibble('in_data/all_outputs_and_obs_081021_wBachmann.csv')

  obs_metadata:
    command: subset_tbl(metadata_tbl, I('site_id'),
      I('lake_lon_deg'), I('lake_lat_deg'))

 # ERA5 data from https://cds.climate.copernicus.eu/cdsapp#!/dataset/10.24381/cds.e2161bac %>%
 # UTC 5PM to be an approx noon for US central location
  out_data/obs_pred.csv:
    command: match_era5_grid2obs(target_name, obs_pred = lake_surface_temp_preds,
      nc_fl = 'in_data/adaptor.mars.internal-1639759493.2063537-10416-6-f31d02b2-53d8-467c-844f-f62fae559a98.nc',
      centroids_sf = modeled_centroids_sf,
      cell_res = I(0.1))

# also add de-biased ERA5* per L&O:L review
  out_data/lake_surface_temp_obs_preds.csv:
    command: add_source_info_obs(target_name,
      obs_pred_fl = 'out_data/obs_pred.csv',
      source_fl = 'out_data/lake_surface_temp_obs.csv')

  out_data/lake_surface_temp_obs.csv:
    command: convert_feather_file(target_name, feather_fl = 'in_data/surface_lake_temp_daily.feather')