NWM 1km LDAS kerchunked dataset

recipes/NWM-1km/meta.yaml

@@ -0,0 +1,22 @@
+title: 'NWM-2.1-grid1km-LDAS'
+description: LDAS variables at 3-hourly intervals from the US National Water Model Version 2.1 Retrospective (1980-). There are 21 LDAS variables including soil moisture, rain, snow, evapotranspiration
+pangeo_forge_version: '0.9.0'
+pangeo_notebook_version: '2022.07.13'
+recipes:
+  - id: NWM-2.1-grid1km-LDAS
+    object: 'recipe:recipe'
+provenance:
+  providers:
+    - name: 'NOAA Office of Water Prediction'
+      description: 'US NOAA OWP'
+      roles:
+        - provider
+        - licensor
+      url: https://water.noaa.gov/about/nwm
+  license: 'CC0'
+maintainers:
+  - name: 'Rich Signell'
+    orcid: '0000-0003-0682-9613'
+    github: rsignell
+bakery:
+  id: 'pangeo-ldeo-nsf-earthcube'

recipes/NWM-1km/recipe.py

@@ -0,0 +1,19 @@
+import fsspec
+import kerchunk
+
+from pangeo_forge_recipes.patterns import pattern_from_file_sequence
+from pangeo_forge_recipes.recipes.reference_hdf_zarr import HDFReferenceRecipe
+
+fs = fsspec.filesystem('s3', anon=True)
+
+all_paths = sorted(fs.glob('noaa-nwm-retrospective-2-1-pds/model_output/*/*LDAS*'))
+
+pattern = pattern_from_file_sequence(['s3://' + path for path in all_paths], 'time')
+
+
+recipe = HDFReferenceRecipe(
+    pattern,
+    netcdf_storage_options={'anon': True},
+    identical_dims=['x', 'y'],
+    preprocess=kerchunk.combine.drop('reference_time'),
+)