-
Notifications
You must be signed in to change notification settings - Fork 63
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
utility script for making animations of forcing inputs #695
base: master
Are you sure you want to change the base?
Conversation
crs = ccrs.epsg(epsg) | ||
crs = ccrs.PlateCarree() |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Which one is crs
meant to be?
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
The pyproj 5070 epsg code enforces a projection boundary that clips off the southern part of the CONUS domain, so for now using the PlateCarree. See comment on line 160. I'll clean this up a bit, just wanted to get it published in some form.
vmin = forcing['RAINRATE'].min().item() | ||
vmin = 0.0009 |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Another duplicate assignment
|
||
forcing = forcing.set_index({"catchment ids":"catchment_ids"}).to_array().squeeze().drop_vars('variable') | ||
#Get the catchment geopandas index, strip to int for xarray index | ||
cat_ids = catchments['divide_id'].apply(lambda x: int(x.split("-")[1])) |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Presumably true at the time of writing, but does the hydrofabric work guarantee unique numeric IDs of catchments across the various kinds that potentially prefix them?
Actually, are there currently or plans for catchment IDs prefixed other than cat-NNNN
, like we have for nexuses?
Script used to generate forcing animations for AGU 2023 fall meeting.