:=============================================================================
:         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 /
