espin-2020 / coastalteam Goto Github PK

Simulate passive particle transport on tidal flows using dorado

Ok so when I call and use tidal_flow_calculator.py, it only creates 3 fields:

['topographic__elevation', 'mean_water__depth'] *at cells
['roughness', 'flood_tide_flow__velocity', 'ebb_tide_flow__velocity'] *at links

1. mean_water__depth
2. flood_tide_flow__velocity
3. ebb_tide_flow_velocity

sadly the rest of them are created by us initially. the point I'm making is that internally there are a number of fields calculated (like half-tide depth) that we would want to use in subsequent scripts/functions. But these aren't exported...so not sure what we want to do exactly here. It might make sense to change this code? Or create a new version that does create more fields for grids.

Variable names are in () for python, and italics for matlab. * means it can be an array (vs. single value)

Inputs already coded in python

• input grids (grid) A

1. bed elevation (z) z
2. Bed Roughness (roughness) MANN +optional, can be single value
   TO ADD GRIDS
3. hydroperiod (hydroperiod) fTide

• tidal range (tidal_range) Trange
• tidal period (tidal_period) Ttide
• bed roughness (roughness)*
• number of rows (nrows)
• number of columns (ncols)
• spacing of grid cells (grid_spacing) dx
• depth of grid cells (mean_depth)

unsure of status

• minimum water depth (min_water_depth) kro
• coefficient of tidal disperson (_diffusion_coef_at_links) DiffS
• how many total time steps () tmax
• length of one time step () dt0
• how many time steps to plot () tINT
• time array () time
• mean sea level (mean_sea_level) msl
• hydroperiod () fTide
• water depth at mean high water () dHW

inputs to be coded in python

• tidal range of vegetation (tidal_range_veg?) TrangeVEG
• Suspended sediment at open boundary, assumes sea (ssc_bc) co2
• rate of relative SLR (RSLR) RSLR
• mud settling velocity (mud_settle_veloc) ws2
• dry bulk density of mud (mud_dry_bulk_density) rbulk2
• mud erodability (mud_erod) me
• linear critical shear stress increase below mean low water (tau_cr_grad) tcrgradeint
• critical shear stress of sediment (tau_cr) taucr
• unvegetated creep coefficient () crMUD ^^no idea what this correlates to in paper...something to do with divergence of bed evolution, see bedcreep.m
• vegetated creep coefficient () crMARSH ^^see above notes
• Lower limit of vegetation growth () dBlo
• Upper limit of vegetation growth () dBup
• Manning coefficient roughness vegetated areas (roughness_deep) Cv ^^see above for grid inputs
• Manning coefficient roughness for unveg. areas (roughness_shallow) Cb ^^see above
  -settling velocity in vegetated area (veg_settle_veloc) wsB
• Critical shear stress of sediment in vegetated areas (tau_cr_veg) taucrVEG
• organic sediment production () Korg

has been coded in python:

• flow velocity () U
• flow velocity in x-dir Ux
• flow velocity in y-dir Uy
• flood tide flow velocity (flood_tide_vel) Ux/Uy ***look at def map_velocity_components_to_nodes
• flood tide ebb velocity (ebb_tide_vel) Ux/Uy ***look at def map_velocity_components_to_nodes

does Not need to be coded in python?:

• FLX

Should do our best to document what we are doing and what we have accomplished. Whether this is in the form suggested in #7 or some other format, it should be done.

Convert sedtran.m to Python

Pythonize totalsedimenterosionMUDsine.m