:============================================================================= : NAMELIST FOR OBSTRUCTIONS VARIABLE PARAMETERS :============================================================================= : obst_var_main : main variable parameters : obst_var_option : variable option for it behaviour : obst_var_init : parameters relative with variable initialization : obst_var_flexibility : parameters relative to the flexible obctructions : obst_var_roughdrag : parameters relative to the roughness and drag behaviour : obst_var_fracxy : parameters relative to small-scale patchiness correction : obst_var_bstress : parameters relative to the bottom shear stress : :------------------------------------------------------------------------------ : obst_var_main : main variable parameters :------------------------------------------------------------------------------ : r_obst_varname : Name (identifier) of the variable : r_obst_type : choice between "UP", "DO", "3D" : "UP" : if variable start from the bed : "DO" : if variable hang from sea-surface (down) : "3D" : if variable is full 3D (based on vertical density variation) : r_l_obst_cylinder : If variable representation is a cylinder (if False : parallelepiped) :------------------------------------------------------------------------------ &obst_var_main r_obst_varname = 'Obstruct' r_obst_type = "UP" r_l_obst_cylinder = .FALSE. / : :------------------------------------------------------------------------------ : obst_var_option : variable option for it behaviour :------------------------------------------------------------------------------ : r_l_obst_flexible : If variable is flexible (if False : rigid) : r_l_obst_noturb : If variable should be represented as macro-roughness (instead of turbulent approach) : r_l_obst_filetimeserie : To use a time-series (not spatial) of obstructions characteristics : !WARNING ! NOT AVAILABLE with r_l_obst_init_spatial = .TRUE. : r_obst_fn_timeserie : File containing temporal obstructions characteristics (if r_l_obst_filetimeserie = .TRUE.) : r_l_obst_filedistri : To use a file describing the vertical distribution of obstruction density : r_obst_fn_distrib : File containing the vertical distribution of obstruction density (if r_l_obst_filedistri = .TRUE.) :------------------------------------------------------------------------------ &obst_var_option r_l_obst_flexible = .TRUE. r_l_obst_noturb = .FALSE. r_l_obst_filetimeserie = .FALSE. r_obst_fn_timeserie = "" r_l_obst_filedistri = .TRUE. r_obst_fn_distrib = 'TEST_CASES/obstruction_seagrass_distri.txt'/ : :------------------------------------------------------------------------------ : obst_var_init : parameters relative with variable initialization :------------------------------------------------------------------------------ : r_l_obst_init_spatial : To use spatially variable file (not temporal) of obstructions charcateristics : !WARNING ! NOT AVAILABLE with r_l_obst_filetimeserie = .TRUE. : r_obst_fn_initspatial : File containing spatial obstructions characteristics (if r_l_obst_init_spatial = .TRUE.) : (variables : height,width,thick and dens) : r_obst_i_height : Initial height (unbent, eg. leaf-length for segrasses) of obstructions, : used if not spatial initialisation (r_l_obst_init_spatial = .FALSE.) : r_obst_i_width : Initial width (or diameter for cylindric obstructions) of obstructions (perpendicular to flow), : used if not spatial initialisation (r_l_obst_init_spatial = .FALSE.) : r_obst_i_thick : Initial thick of obstructions (along the flow), : used if not spatial initialisation (r_l_obst_init_spatial = .FALSE.) : r_obst_i_dens : Initial density of obstructions (maximum density if using a vertical distribution file), : used if not spatial initialisation (r_l_obst_init_spatial = .FALSE.) :------------------------------------------------------------------------------ &obst_var_init r_l_obst_init_spatial = .FALSE. r_obst_fn_initspatial = "" r_obst_i_height = 0.073 r_obst_i_width = 0.00082 r_obst_i_thick = 0.0002 r_obst_i_dens = 80940.0 / : :------------------------------------------------------------------------------ : obst_var_flexibility : parameters relative to the flexible obstructions (r_l_obst_flexible = .TRUE.) :------------------------------------------------------------------------------ : r_l_obst_abdelposture : To use Abdhelhrman's (2007) procedure to compute bending : r_l_obst_param_height : To use empirical (exponential desrease) formulation to compute bending : r_obst_c_abdel_nmax : Number of segments for Abdhelhrman's (2007) procedure : r_obst_c_rho : Volumic mass of obstructions for Abdhelhrman's (2007) procedure : r_obst_c_lift : Lift coefficient for Abdhelhrman's (2007) procedure : r_obst_c_shelter : Sheltering coefficient Afor bdhelhrman's (2007) procedure : r_obst_c_height_x0 : First parameter for empirical formulation : r_obst_c_height_x1 : Second parameter for empirical formulation :------------------------------------------------------------------------------ &obst_var_flexibility r_l_obst_abdelposture = .TRUE. r_l_obst_param_height = .FALSE. r_obst_c_abdel_nmax = 5 r_obst_c_rho = 650.0 r_obst_c_lift = 0.0 r_obst_c_shelter = 100.0 r_obst_c_height_x0 = 0.99 r_obst_c_height_x1 = -3.92 / : :------------------------------------------------------------------------------ : obst_var_roughdrag : parameters relative to the roughness and drag behaviour :------------------------------------------------------------------------------ : r_l_obst_drag_cste : To use a constant drag coefficient for obstructions in hydrodynamics : (if false, drag varies depending on the bending angle) : Only works for FLEXIBLE obstructions : r_l_obst_abdelrough_cste : To use a constant drag during reconstruction of velocity profiles : using Abdelhrman's (2003) procedure (variable drag if false) : r_obst_c_crough_x0 : First coefficient for drag coefficient during Abdelhrman's (2003) procedure : r_obst_c_crough_x1 : Second coefficient for drag coefficient during Abdelhrman's (2003) procedure : r_obst_c_drag : Drag coefficient (maximum value if not constant) for obstructions elements : r_obst_c_lz : Coefficient for turbulent dissipation time-scale between obstructions elements :------------------------------------------------------------------------------ &obst_var_roughdrag r_l_obst_drag_cste = .TRUE. r_l_obst_abdelrough_cste = .FALSE. r_obst_c_crough_x0 = 2.6 r_obst_c_crough_x1 = 0.0 r_obst_c_drag = 1.5 r_obst_c_lz = 0.10 / : :------------------------------------------------------------------------------ : obst_var_fracxy : parameters relative to small-scale patchiness correction :------------------------------------------------------------------------------ : r_l_obst_fracxy : To take account for patchiness correction (if false, no correction) : r_obst_fracxy_type : Kind of correction method : : 0 : patchiness correction is equal to the fraction of cell occupied by obstructions : 1 : patchiness correction is equal to an exponential of the fraction of cell occupied : by obstructions with one coefficient (r_obst_c_fracxy_k0) : 2 : patchiness correction is equal to an exponential of the fraction of cell occupied : by obstructions with several coefficients (r_obst_c_fracxy_k0, r_obst_c_fracxy_k1 and : r_obst_c_fracxy_l) : 3 : patchiness correction is equal to the product of the fraction of cell occupied by obstructions : and r_obst_c_fracxy_k0 : r_obst_c_fracxy_k0 : Coefficient for the corrections (type 1 and 2 and 3) : r_obst_c_fracxy_k1 : First parameter for correction of the exponential coefficient (type 2) : r_obst_c_fracxy_l : Second parameter for correction of the exponential coefficient (type 2) :------------------------------------------------------------------------------ &obst_var_fracxy r_l_obst_fracxy = .FALSE. r_obst_fracxy_type = 2 r_obst_c_fracxy_k0 = 6.0 r_obst_c_fracxy_k1 = 12.0 r_obst_c_fracxy_l = 3.5 / : :------------------------------------------------------------------------------ : obst_var_bstress : parameters relative to the bottom shear stress (used only if #MUSTANG) :------------------------------------------------------------------------------ : r_l_obst_z0bstress : To activate the impact of obstruction on roughness length : used to compute the bottom shear stress (only for UP type) : r_obst_z0bstress_option : Option to compute the obstruction induced roughness length: : 0 : constant z0 (r_obst_c_z0bstress) : 1 : parameterization : r_obst_c_z0bstress : Constant (uncorrected value of roughness length) : r_obst_c_z0bstress_x0 : First parameter for rouhgness length computation (in 3D) : r_obst_c_z0bstress_x1 : Second parameter for rouhgness length computation (in 3D) :------------------------------------------------------------------------------ &obst_var_bstress r_l_obst_z0bstress = .TRUE. r_obst_z0bstress_option = 0 r_obst_c_z0bstress = 0.00012 r_obst_c_z0bstress_x0 = 0.00055 r_obst_c_z0bstress_x1 = -0.397 /