3_ncdf.yml

target_default: 3_ncdf

packages:
  - tidyverse
  - scipiper
  - sf
  - ncdfgeom

sources:
  - src/ncdf_utils.R

targets:
  3_ncdf:
    depends:
      - out_data/netcdf_filehashes.yaml


  out_data/netcdf_filehashes.yaml:
    command: combine_hash_files(target_name,
      'out_data/weather_data_filehashes.yaml', 'out_data/prediction_data_filehashes.yml')


  predict_id_groups:
    command: subset_tbl(metadata_tbl, I('group_id'),
      I('site_id'), I('lake_lon_deg'), I('lake_lat_deg'), I('elevation_m'), I('group_bbox'))

  weather_metadata:
    command: tibble(
      name = I(c('dlwrfsfc','dswrfsfc','tmp2m','ugrd10m','vgrd10m')),
      units = I(c('W/m^2','W/m^2','K','m/s','m/s')),
      long_name = I(c(
        'lw radiation flux downwards (surface) [w/m^2]',
        'sw radiation flux downwards (surface) [w/m^2]',
        '2-m above ground temperature [k]',
        '10-m above ground zonal wind speed [m/s]',
        '10-m above ground meridional wind speed [m/s]')),
      precision = I(c('.1','.1','.2','.3','.3')))

  predict_metadata:
    command: tibble(
      name = I('surftemp'),
      units = I('°C'),
      long_name = I('Surface water temperature [°C]'),
      precision = I('.2'))

  weather_id_groups:
    command: subset_tbl(metadata_tbl, I('group_id'), I('weather_id'), I('x'), I('y'),
      I('weather_lon_deg'), I('weather_lat_deg'), I('group_bbox'))

  out_data/weather_data_filehashes.yaml:
    command: build_driver_nc(target_name, 'in_data/NLDAS_nc_dailies.yaml',
      export_range = I(c('1980-01-01','2020-12-31')),
      out_pattern = I('tmp/%s_weather_%s.nc'),
      weather_id_groups,
      weather_metadata)

  out_data/prediction_data_filehashes.yml:
    command: build_prediction_nc(target_name,
      I('/Volumes/ThunderBlade/SWT_results'),
      export_range = I(c('1980-01-01','2020-12-31')),
      out_pattern = I('tmp/%s_predicted_temp_%s.nc'),
      dummy = I('2021-07-15'),
      predict_id_groups,
      predict_metadata)