Missing nodata value in DataCatalog for worldclim
The function precip() from forcing.py gave an error when using the worldclim data because the nodata value was not set. This could be added by modifying the argument for no data in the current yml file of the data catalog
---- code piece that didn't work ---- #Testing with clim argument cat = DataCatalog() p_precip = cat.get_rasterdataset("era5", variables="precip") #era horuly
#We create a more refined grid
p_transform = p_precip.raster.transform
p_shape = (21,18)
grid = full_from_transform(transform = p_transform, shape = p_shape, crs = 4326)
p_clim = cat.get_rasterdataset("worldclim") p_clim.raster.set_nodata(-999.0) #without this line, ti does not work
pout_clim = precip(p_precip, grid, clim=p_clim)
@couasnonanais For clarity, this is with the "artifact_data" data catalog right?
Both the deltares_data and artifact_data catalogs need updating