input('test-data', name: 'data') as data
->
# conceptually this join could be used to convert rotated pole grid coords
# to lat/lon.
join(on: data.LAT = site.LAT and data.LON = site.LON) as data_with_site

input('test-site', name: 'site') as sites
->
data_with_site.rhs

data_with_site
->
select({
  data.*,
  site.SITE as SITE
})
->
save('output', format: 'csv')
