noaa-ocs-modeling / EnsemblePerturbation

perturbation of coupled model input over a space of input variables
https://ensembleperturbation.readthedocs.io
Creative Commons Zero v1.0 Universal
7 stars 3 forks source link

tidal perturbations #25

Open krober10nd opened 3 years ago

krober10nd commented 3 years ago

What's the strategy you guys use to perturb the tides to capture arbitrary storm-tide timing?

ghost commented 3 years ago

oh whoops, this function call is outdated and will not work. Good catch @krober10nd , I need to update the examples

the function generate_adcirc_configuration is now in CoupledModelDriver.

The function works by reading JSON files from the specified directory and extracting model configuration from them. The JSON files are written by creating an AdcircRunConfiguration or NEMSADCIRCRunConfiguration object and calling .write_directory(directory)

Perturbations can be designated by editing configure_modeldriver.json (which is wrapped by a ModelDriverJSON object) and adding any number of cases in which any number of variables can be perturbed to any value.

krober10nd commented 3 years ago

Thanks Zach for the clarification.

I'll try these out. I guess I was looking for an existing methodology that made sense to perturb start times to "adequately" sample the storm-tide uncertainty.

This is of course a science question that constantly arises in our field.

ghost commented 3 years ago

@krober10nd here is the current example to perturb ATCF tracks by 4 variables (2 scalar, 2 spatial) and build configurations per-value: https://github.com/noaa-ocs-modeling/EnsemblePerturbation/blob/6e4b873f70589b2412f8d2e8a9fa8c8fa6dcedef/examples/adcirc/perturb_tracks.py#L88-L108

I will look into how we might perturb tidal forcing. I'm not an oceanographer but I imagine we could perturb phase + timings of tides after they are read from the TPXO dataset