Open durack1 opened 1 year ago
@lauradomar, I'll await these new location deets before prepping the next round of data
@durack1 The maximum extent limits for the seagrass groups are (format: max lat, max lon, min lat, min lon):
-9, 160, -46, 113
We could just add a boundary around that to add a few degrees either side as previously I think?
Let me know of any questions and thanks! Laura
@lauradomar here is a first attempt at infilling
Native output
Zonally infilled output
I am unfamiliar with the subregions that are the focus, so it would be great to know where these seagrass beds lie, so we can ascertain that the infilling is providing a useable and representative field (e.g. the top of Australia and Gulf of Carpentaria is not usable for this infill method)
@lauradomar just pinging you on this, so we can get these fields finalized and prepped for model runs
@lauradomar just pinging you here, I am yet to provide you data, so we better get on with this!
Thanks! Sorry, had deadlines and long holiday weekend. Here the info: The maximum extent limits for the seagrass groups are (format: max lat, max lon, min lat, min lon):
-9, 160, -46, 113
The different individual extents for the genera are:
-9, 160, -44, 113 (Halophila)
-9, 160, -46, 114 (Zostera)
-25, 153, -44, 113 (Posidonia)
-9, 149, -22, 114 (Thalassia)
-9, 147, -20, 130 (Enhalus)
We could just drop the ones that are in the North if the data would not be usable? The output otherwise looks good, I would just need to have the land based on your map removed like before so only cells in the ocean would have values that are put into the model.
@lauradomar wow, these regions are far more expansive than I would have assumed, I thought we'd be constraining to tight coastal boxes in representative zones, rather than a box that includes the whole of Australia for e.g. - red box. red - full extent blue - Halophila black - Zostera white - Posidonia magenta - Thalassia yellow - Enhalus
By eyeball, if we constrained the data usage South of 15deg S, then I think we could get away with it. My intention would be to provide the complete infilled array, as I had presumed you'd be cookie cutting out small boxes where the seagrass beds are, not Aust-wide regions. This doesn't seem to be the use case that you have in mind.
There is a need to create higher (~0.1deg) resolution data for the seagrass, near-coastal locations. These will need to be defined, so if they can be included below, ideal:
See https://github.com/durack1/Roethigetal21NatClimChg/blob/main/readNcWriteAsc.py#L200-L212
ping @lauradomar