Tags: NOAA-OWP/t-route
Tags
Ten diffusive domains w schism boundary (#573) * f90 files * step2 * run all ten diffusive domains in serial with schism depth boundary data * code update 6/29/2022 * run all ten diffusive domains in serial with schism boundary * more clean up * comment out unused maxCourant variable * removed private file paths * further cleanup * change dt_schisma and dt_db_g and cleanup * remove coastal_domain_testing.yaml * remove old files in LowererColorado_TX folder * remove 2 more old files in LowerColorado_TX folder * updated hydrofabric & natural xsec files * remove private path in yaml * change default value of diffusivity limits * remove list() * list() removed * remove two natural xsec nc files * updated mask file only for Lower Colorado River * remove data of HurricaneLaura Co-authored-by: Dong Kim <[email protected]>
Expedite time slice interpolation routine (#536) * lake-gage crosswalk creation * use cross walk dataframes and parallelize interpolation * name index * return hybrid lake-gage crosswalks * save link_gage_df for preprocess * pass cross walk dataframe to get_obs_from_timeslices * simplify crosswalk with dataframes * fix link_lake_df creation * change waterbody types to 1 if USGS or USACE hybrid obs are unavailable * save usgs and usace crosswalk dfs to preprocess dict * add reservoir crosswalk dfs to list of things to unpack from preprocess and return * add crosswalk to list of returns from unpack_nwm_preprocess Co-authored-by: awlostowski-noaa <[email protected]>
Add Nudging decay (including reading LastObs files) (#310) * added last obs df including discharge and ids * lastobs nc file folder * formatting and cleanup * added both discharge and model discharge * added prediction delta timesteps * da decay with exp decay incorporated, unblackened for readability * added check between obs file results and our fvd output * removed pdbs and set to run on lastobs * added verbose statements * restructuring of lastobs df to simplify process * generalized last timestep index call to automatically determine from various input sizes * github requested fixes to if statements and cleanup * updated variable * working ncar da decay prototype need to bring in real lastobs data inside equation * getting da values in order through classic da assim technique * pushing ids and values to mc reach * fixed gage id matching, cython is broken cant compile correctly, need to print values to identify if da is properly working * saving changes, trying to fix old DA function * restructing da timeslice file read to use datetime, not generalized before * added last obs df including discharge and ids * lastobs nc file folder * formatting and cleanup * added both discharge and model discharge * added prediction delta timesteps * da decay with exp decay incorporated, unblackened for readability * added check between obs file results and our fvd output * removed pdbs and set to run on lastobs * added verbose statements * restructuring of lastobs df to simplify process * generalized last timestep index call to automatically determine from various input sizes * github requested fixes to if statements and cleanup * updated variable * generalized da naming conventional and date timeframe with improved interpolation * removed extra comments * remove dependence on not-yet-created flowveldepth * name "last_obs_file" * include data_assimilation_parameters to yaml * include empty dict for data_assimilation_parameters in yaml * black * added paths to shared drive locations and blackened files * quick merge changes * files working * sync for merge * sync for merge * sync merge to upstream * move last obs function next to usgs_da function * add TODOs * add lastobs to other parallel modes * move last_obs above usgs_df * cimports and cdefs for last_obs * fixed broken usgs_df names were changed to folder in some areas and filter in others * moved da functions into single wrapper in nhd_network_utilities also harmonized inputs a bit. * add da yaml file (DATA NOT YET ADDED) * drop non-5-minute entries from time_slices Also harmonizing inputs for merge. * add function for finding the tailwater for a given segment * add masks * use pandas date_range * cleanup * add comment showing possible handling of extended DA * Revert "add comment showing possible handling of extended DA" This reverts commit 0455466. * temporarily disable last_obs * Update example yaml with inputs that work * temporarily disable last_obs * update comment * adjust DA for perfect match * removed filter list * use efficient shape call for usgs_positions_list length * add gage_maxtime and pseudocode for lastobs * Identified major hard-code issue in structured- and structured-obj * update lastobs comments and pseudocode * update da test yaml file with additional gage options * use "reindex" to fill/eliminate columns for usgs_df * functions in place for decay but last obs file is behaving incorrectly and changing parity check even when all uses are turned off * added decay timestep count * Use new fields in flowveldepth to simplify initial condition handling To facilitate this, added a constant qvd_ts_w (flowveldepth timestep width) to define the standard column width. * add two additional segments for parity checking * reconfigure reach splitting to consider gages * update diffusive call signature * yaml updates for test * black Co-authored-by: James Halgren <[email protected]>