Difference between revisions of "sediment.in"

From WikiROMS
Jump to navigationJump to search
 
(17 intermediate revisions by 3 users not shown)
Line 1: Line 1:
<div class="title">Sediment Model Input Script - <span class="red">sediment.in</span></div>
<div class="title">Sediment Model Input Script - <span class="red">sediment.in</span></div>


The sediment model input ASCII files <span class="red">sediment.in</span> is organized in several sections as shown below, with links to more detailed explanation where required. The format is keyword == value(s); the order in which keywords appear in the file is not important.
The <span class="red">sediment.in</span> file sets the parameters for sediment output. The name of this file is set by the <span class="blue">SPARNAM</span> keyword in the [[roms.in]] file. This standard input ASCII file is organized in several sections as shown below, with links to more detailed explanation where required.


{{note}} '''Notice:''' A default <span class="red">sediment.in</span> input script is provided in the '''User/External''' subdirectory. Also there are several sediment input scripts in the '''ROMS/External''' subdirectory which are used in the distributed [[Test Cases|test cases]]. They are usually named <span class="red">sediment_</span><span class="blue">app</span><span class="red">.in</span> where <span class="blue">app</span> is the lowercase of the test case [[Options|cpp option]].


==Sediment model control switches==
{{note}} '''Notice:''' Detailed information about ROMS input script file syntax can be found [[Input Parameter Files|here]].


* Switch is used to control sediment model computation within nested and/or multiple connected grids, [1:[[Ngrids]]] values are expected.<div class="box">     [[Lsediment]] == T</div>
{{note}} '''Notice:''' A default <span class="red">sediment.in</span> input script is provided in the '''User/External''' subdirectory. Also there are several sediment input scripts in the '''ROMS/External''' subdirectory which are used in the distributed [[Test Cases|test cases]]. They are usually named <span class="red">sediment_</span><span class="blue">app</span><span class="red">.in</span> where <span class="blue">app</span> is the lowercase of the test case [[C Preprocessor|CPP option]].


==General sediment bed model controls==


* Depositional bed layer thickness (m) criteria to create a new layer. If deposition exceeds this value, then a new layer is created, [1:[[Ngrids]]] values are expected.<div class="box">[[newlayer_thick|NEWLAYER_THICK]] == 0.01d0</div>
==Sediment Model Control Switches==


* Bed load transport rate coefficient, [1:[[Ngrids]]] values are expected.<div class="box"> [[bedload_coeff|BEDLOAD_COEFF]] == 0.05d0</div>
* Switch is used to control sediment model computation within nested and/or multiple connected grids, [1:[[Variables#Ngrids|Ngrids]]] values are expected.<div class="box">     [[Variables#Lsediment|Lsediment]] == T</div>


==Suspended Cohesive Sediment Parameters==
==Tracer Advection Schemes==
 
*Set horizontal and vertical advection schemes for biological tracers. A different advection scheme is allowed for each tracer. For example, a positive-definite (monotonic) algorithm can be activated for salinity and biological tracers, while a different one is set for temperature.<div class="box">Keyword    Advection Algorithm<br /><br />A4        4th-order Akima (horizontal/vertical)<br />C2        2nd-order centered differences (horizontal/vertical)<br />C4        4th-order centered differences (horizontal/vertical)<br />HSIMT      3th-order HSIMT with TVD limiter (horizontal/vertical)<br />MPDATA    recursive flux corrected MPDATA (horizontal/vertical)<br />SPLINES    parabolic splines reconstruction (only vertical)<br />SU3        split third-order upstream (horizontal/vertical)<br />U3        3rd-order upstresm-bias (only horizontal)</div>The user has the option of specifying the full Keyword or the first two letters, regardless if using uppercase or lowercase. If nested grids, specify values for each grid. For more details please read the [[roms.in#Tracer_Advection_Schemes|roms.in Tracer Advection Schemes]] section. [1:[[Variables#NBT|NBT]],[[Variables#Ngrids|Ngrids]]] values are expected.<div class="box">    [[Variables#Hadvection|Hadvection]] == HSIMT                          ! [[Variables#idsed|idsed]](:), compact<br /><br />    [[Variables#Vadvection|Vadvection]] == HSIMT                          ! [[Variables#idsed|idsed]](:), compact</div>
 
*Adjoint-based algorithms can have different horizontal and vertical schemes for biological tracers.[1:[[Variables#NBT|NBT]],[[Variables#Ngrids|Ngrids]]] values are expected.<div class="box"> [[Variables#ad_Hadvection|ad_Hadvection]] == U3                            ! [[Variables#idsed|idsed]](:), compact<br /><br /> [[Variables#ad_Vadvection|ad_Vadvection]] == C4                            ! [[Variables#idsed|idsed]](:), compact</div>
 
==Lateral Open Boundary Conditions Parameters==
 
*The lateral boundary conditions are entered with a keyword. A value is expected for each boundary segment per nested grid for each state variable. The sediment tracer variables (cohesive and non-cohesive) require [1:4,1:[[Variables#NCS|NCS]]+[[Variables#NNS|NNS]],[[Variables#Ngrids|Ngrids]]] values. If specifying every tracer, enter first 1:[[Variables#NCS|NCS]] cohesive sediment entries followed by 1:[[Variables#NNS|NNS]] non-cohesive sediment entries. The boundary order is: '''1'''=west, '''2'''=south, '''3'''=east, and '''4'''=north. That is, anticlockwise starting at the western boundary.<br /><br />The keyword is case insensitive and usually has three characters. However, it is possible to have compound keywords, if applicable. For example, the keyword '''RadNud''' implies radiation boundary condition with nudging. This combination is usually used in active/passive radiation conditions.<br /><br />{{note}}'''Notice:''' It is possible to specify the lateral boundary conditions for all biological tracers in a compact form with a single entry. If so, all the biological tracers are assumed to have the same boundary condition as the single entry.<div class="box"><span class="twilightBlue">!  Keyword    Lateral Boundary Condition Type<br />!<br />!  Cla        Clamped                                _____N_____    j=Mm<br />!  Clo        Closed                                |    4    |<br />!  Gra        Gradient                              |          |<br />!  Nes        Nested                              1 W          E 3<br />!  Nud        Nudging                              |          |<br />!  Per        Periodic                              |_____S_____|<br />!  Rad        Radiation                                  2          j=1<br />!                                                  i=1        i=Lm<br />!                  W      S      E      N<br />!                  e      o      a      o<br />!                  s      u      s      r<br />!                  t      t      t      t<br />!                          h              h<br />!<br />!                  1      2      3      4</span><br /><br />  [[Variables#LBC|LBC]](isTvar) ==  Per    Clo    Per    Clo      ! [[Variables#idsed|idsed]](:), compact</div>
 
*Adjoint-based algorithms can have different lateral boundary conditions keywords.<div class="box">[[Variables#ad_LBC|ad_LBC]](isTvar) ==  Per    Clo    Per    Clo      ! [[Variables#idsed|idsed]](:), compact</div>
 
==General Sediment Bed Model Controls==
 
* Depositional bed layer thickness (m) criteria to create a new layer. If deposition exceeds this value, then a new layer is created, [1:[[Variables#Ngrids|Ngrids]]] values are expected.<div class="box">[[Variables#newlayer_thick|NEWLAYER_THICK]] == 0.01d0</div>
 
* Bed load transport rate coefficient, [1:[[Variables#Ngrids|Ngrids]]] values are expected.<div class="box"> [[Variables#bedload_coeff|BEDLOAD_COEFF]] == 0.05d0</div>
 
==Bed Layer and Bottom Sediment Parameters==
 
* Logical switches to activate writing of bed layer and bottom sediment fields into history output NetCDF file(s), [1:[[Variables#Ngrids|Ngrids]]] values are expected.
:<div class="box">  [[Variables#Hout|Hout]](ithck) == T      ! bed_thickness          sediment layer thickness<br />  [[Variables#Hout|Hout]](iaged) == T      ! bed_age                sediment layer age<br />  [[Variables#Hout|Hout]](iporo) == T      ! bed_porosity            sediment layer porosity<br />  [[Variables#Hout|Hout]](idiff) == F      ! bed_biodiff            biodiffusivity<br/><br />  [[Variables#Hout|Hout]](ids50) == T      ! grain_diameter          mean grain diameter<br />  [[Variables#Hout|Hout]](idens) == T      ! grain_density          mean grain density<br />  [[Variables#Hout|Hout]](iwsed) == T      ! settling_vel            mean settling velocity<br />  [[Variables#Hout|Hout]](itauc) == T      ! erosion_stress          critical erosion stress<br />  [[Variables#Hout|Hout]](irlen) == T      ! ripple_length          ripple length<br />  [[Variables#Hout|Hout]](irhgt) == T      ! ripple_height          ripple height<br />  [[Variables#Hout|Hout]](ibwav) == T      ! bed_wave_amp            wave excursion amplitude<br />  [[Variables#Hout|Hout]](izdef) == T      ! Zo_def                  default bottom roughness<br />  [[Variables#Hout|Hout]](izapp) == T      ! Zo_app                  apparent bottom roughness<br />  [[Variables#Hout|Hout]](izNik) == T      ! Zo_Nik                  Nikuradse bottom roughness<br />  [[Variables#Hout|Hout]](izbio) == T      ! Zo_bio                  biological bottom roughness<br />  [[Variables#Hout|Hout]](izbfm) == T      ! Zo_bedform              bed form bottom roughness<br />  [[Variables#Hout|Hout]](izbld) == T      ! Zo_bedload              bed load bottom roughness<br />  [[Variables#Hout|Hout]](izwbl) == T      ! Zo_wbl                  wave bottom roughness<br />  [[Variables#Hout|Hout]](iactv) == T      ! active_layer_thickness  active layer thickness<br />  [[Variables#Hout|Hout]](ishgt) == T      ! saltation              saltation height<br /></div>
 
* Logical switches to activate writing of bed layer and bottom sediment fields into quicksave output NetCDF file(s), [1:[[Variables#Ngrids|Ngrids]]] values are expected.
:<div class="box">  [[Variables#Qout|Qout]](ithck) == F      ! bed_thickness          sediment layer thickness<br />  [[Variables#Qout|Qout]](iaged) == F      ! bed_age                sediment layer age<br />  [[Variables#Qout|Qout]](iporo) == F      ! bed_porosity            sediment layer porosity<br />  [[Variables#Qout|Qout]](idiff) == F      ! bed_biodiff            biodiffusivity<br/><br />  [[Variables#Qout|Qout]](ids50) == F      ! grain_diameter          mean grain diameter<br />  [[Variables#Qout|Qout]](idens) == F      ! grain_density          mean grain density<br />  [[Variables#Qout|Qout]](iwsed) == F      ! settling_vel            mean settling velocity<br />  [[Variables#Qout|Qout]](itauc) == F      ! erosion_stress          critical erosion stress<br />  [[Variables#Qout|Qout]](irlen) == F      ! ripple_length          ripple length<br />  [[Variables#Qout|Qout]](irhgt) == F      ! ripple_height          ripple height<br />  [[Variables#Qout|Qout]](ibwav) == F      ! bed_wave_amp            wave excursion amplitude<br />  [[Variables#Qout|Qout]](izdef) == F      ! Zo_def                  default bottom roughness<br />  [[Variables#Qout|Qout]](izapp) == F      ! Zo_app                  apparent bottom roughness<br />  [[Variables#Qout|Qout]](izNik) == F      ! Zo_Nik                  Nikuradse bottom roughness<br />  [[Variables#Qout|Qout]](izbio) == F      ! Zo_bio                  biological bottom roughness<br />  [[Variables#Qout|Qout]](izbfm) == F      ! Zo_bedform              bed form bottom roughness<br />  [[Variables#Qout|Qout]](izbld) == F      ! Zo_bedload              bed load bottom roughness<br />  [[Variables#Qout|Qout]](izwbl) == F      ! Zo_wbl                  wave bottom roughness<br />  [[Variables#Qout|Qout]](iactv) == F      ! active_layer_thickness  active layer thickness<br />  [[Variables#Qout|Qout]](ishgt) == F      ! saltation              saltation height<br /></div>
 
 
==Non-cohesive Sediment Parameters==
 
There are [[Variables#NNS|NNS]] non-cohesive (sand) sediment tracers per nested grid in ROMS. Therefore, [1:[[Variables#NNS|NNS]],1:[[Variables#Ngrids|Ngrids]]] values are expected for these parameters.
 
* Median sediment grain diameter (mm).<div class="box">    [[Variables#sd50|SAND_SD50]] == 1.0d0</div>
 
* Sediment concentration (kg/m<sup>3</sup>).<div class="box">    [[Variables#Csed|SAND_CSED]] == 0.0d0</div>
 
* Sediment grain density (kg/m<sup>3</sup>).<div class="box">    [[Variables#Srho|SAND_SRHO]] == 2650.0d0</div>
 
* Particle settling velocity (mm/s).<div class="box">    [[Variables#Wsed|SAND_WSED]] == 1.0d0</div>
 
* Surface erosion rate (kg/m<sup>2</sup>/s).<div class="box">    [[Variables#Erate|SAND_ERATE]] == 5.0d-4</div>
 
* Critical shear for erosion and deposition (N/m<sup>2</sup>).<div class="box">  [[Variables#tau_ce|SAND_TAU_CE]] == 0.1d0<br />  [[Variables#tau_cd|SAND_TAU_CD]] == 0.1d0</div>
 
* Porosity (nondimensional): <span class="red">Vwater/(Vwater+Vsed)</span>. Its values range between '''0''' and '''1'''.<div class="box">    [[Variables#poros|SAND_POROS]] == 0.5d0</div>
 
* Lateral, constant, harmonic/biharmonic horizontal diffusion of tracer for the nonlinear model and adjoint-based algorithms.<div class="box">    [[Variables#tnu2|SAND_TNU2]] == 0.0d0                          ! m<sup>2</sup>/s<br />    [[Variables#tnu4|SAND_TNU4]] == 0.0d0                          ! m<sup>4</sup>/s<br /><br />  [[Variables#ad_tnu2|ad_SAND_TNU2]] == 0.0d0                          ! m<sup>2</sup>/s<br />  [[Variables#ad_tnu4|ad_SAND_TNU4]] == 0.0d0                          ! m<sup>4</sup>/s</div>
 
* Logical switches to increase horizontal diffusivity of noncohesive sediment trace in specific areas of the application domain (like sponge areas) for the desired grid.<div class="box">    [[Variables#sand_sponge|SAND_Sponge]] == F</div>
 
* Vertical mixing coefficients for tracers in nonlinear model and basic state scale factor in adjoint-based algorithms.<div class="box">  [[Variables#Akt_bak|SAND_AKT_BAK]] == 5.0d-6                          ! m<sup>2</sup>/s<br />  [[Variables#ad_Akt_fac|SAND_AKT_fac]] == 1.0d0                          ! nondimensional</div>


There are [[NCS]] cohesive (mud) sediment tracers per nested grid in ROMS. Therefore, [1:[[NCS]],1:[[Ngrids]]] values are expected for these parameters.
* Nudging/relaxation time scales, inverse scales will be computed internally.<div class="box">    [[Variables#Tnudg|SAND_TNUDG]] == 0.0d0                          ! days</div>


* Median sediment grain diameter (mm).<div class="box">     [[sd50|MUD_SD50]] == 0.01d0  0.06d0  0.10d0</div>
* Morphological time scale factor (greater than or equal to 1.0). A value of 1.0 has no scale effect.<div class="box">[[Variables#morph_fac|SAND_MORPH_FAC]] == 1.0d0                          ! nondimensional</div>


* Sediment concentration (kg/m3).<div class="box">     [[Csed|MUD_CSED]] == 0.0d0  0.0d0  0.0d0</div>
* Logical switch to specify which variables to consider on tracers point Sources/Sinks (like river runoff). [1:[[Variables#NCS|NCS]],[[Variables#Ngrids|Ngrids]]] values expected.
:<div class="box">   [[Variables#LtracerSrc|SAND_Ltsrc]] == F</div>


* Sediment grain density (kg/m3).<div class="box">     [[Srho|MUD_SRHO]] == 2650.0d0  2650.0d0  2650.0d0</div>
*Logical switches to specify which variables to process for tracers climatology: [1:[[Variables#NCS|NCS]],[[Variables#Ngrids|Ngrids]]] values expected.
:<div class="box">   [[Variables#SAND_Ltclm|SAND_Ltclm]] == F</div>


* Particle settling velocity (mm/s).<div class="box">     [[Wsed|MUD_WSED]] == 0.1d0  1.0d0  2.0d0</div>
* Logical switches to nudge the desired noncohesive sediment tracer climatology field. If not analytical climatology fields, users need to turn on the logical switches above to process the fields from the climatology NetCDF file that are needed for nudging. <div class="box">   [[Variables#SAND_Tnudge|SAND_Tnudge]] == F</div>


* Surface erosion rate (kg/m2/s).<div class="box">     [[Erate|MUD_ERATE]] == 5.0d-5 5.0d-5 5.0d-5</div>
* Logical switches to activate writing of non-cohesive sediment into output history file.
:<div class="box"> [[Variables#Hout|Hout]](idsand) == T      ! sand_01, ...            suspended concentration<br />  [[Variables#Hout|Hout]](iSfrac) == T      ! sandfrac_01, ...       bed layer fraction<br />  [[Variables#Hout|Hout]](iSmass) == T      ! sandmass_01, ...        bed layer mass<br /> [[Variables#Hout|Hout]](iSUbld) == T      ! bedload_Usand_01, ...  bed load at U-points<br /> [[Variables#Hout|Hout]](iSVbld) == T      ! bedload_Vsand_01, ...   bed load at V-points</div>


* Critical shear for erosion (N/m2).<div class="box">   [[tau_ce|MUD_TAU_CE]] == 0.05d0 0.05d0 0.05d0</div>
* Logical switches to activate writing of non-cohesive sediment into output quicksave file.
:<div class="box"> [[Variables#Qout|Qout]](idsand) == F      ! sand_01, ...            suspended concentration<br />  [[Variables#Qout|Qout]](iSsand) == F      ! sand_01_sur, ...        surface concentration<br />  [[Variables#Qout|Qout]](iSfrac) == F      ! sandfrac_01, ...        bed layer fraction<br />  [[Variables#Qout|Qout]](iSmass) == F      ! sandmass_01, ...        bed layer mass<br /> [[Variables#Qout|Qout]](iSUbld) == F      ! bedload_Usand_01, ...  bed load at U-points<br /> [[Variables#Qout|Qout]](iSVbld) == F      ! bedload_Vsand_01, ...   bed load at V-points</div>


* Critical stress for deposition (N/m2), not used.<div class="box">   [[tau_cd|MUD_TAU_CD]] == 0.01d0 0.01d0 0.01d0</div>
* Logical switches to activate writing of time-averaged, non-cohesive sediment fields into averages output file.
:<div class="box"> [[Variables#Aout|Aout]](idsand) == T      ! sand_01, ...            suspended concentration<br /><br />  [[Variables#Aout|Aout]](iSTTav) == T      ! sand_01_2, ...          quadratic <t*t> tracer terms<br />  [[Variables#Aout|Aout]](iSUTav) == T      ! u_sand_01, ...          quadratic <u*t> tracer terms<br />  [[Variables#Aout|Aout]](iSVTav) == T      ! v_sand_01, ...          quadratic <v*t> tracer terms<br />  [[Variables#Aout|Aout]](SHUTav) == T      ! Huon_sand_01, ...      tracer volume flux, <Huon*t><br /> [[Variables#Aout|Aout]](SHVTav) == T      ! Hvom_sand_01, ...      tracer volume flux, <Hvom*t><br /><br /> [[Variables#Aout|Aout]](iSUbld) == T      ! bedload_Usand_01, ...  bed load at U-points<br />  [[Variables#Aout|Aout]](iSVbld) == T      ! bedload_Vsand_01, ...   bed load at V-points</div>


* Porosity (nondimensional): <span class="red">Vwater/(Vwater+Vsed)</span>. Its values range between '''0''' and '''1'''.<div class="box">     [[poros|MUD_POROS]] == 0.0d0 0.0d0 0.0d0</div>
* Logical switches to activate writing of time-averaged, non-cohesive sediment diagnostic terms into output diagnostics file.
:<div class="box"> [[Variables#Dout|Dout]](STrate) == T      ! sand_01_rate, ...      time rate of change<br />  [[Variables#Dout|Dout]](SThadv) == T      ! sand_01_hadv, ...      horizontal total advection<br />  [[Variables#Dout|Dout]](STxadv) == T      ! sand_01_xadv, ...      horizontal XI-advection<br />  [[Variables#Dout|Dout]](STyadv) == T      ! sand_01_yadv, ...      horizontal ETA-advection<br />  [[Variables#Dout|Dout]](STvadv) == T      ! sand_01_vadv, ...      vertical advection<br /> [[Variables#Dout|Dout]](SThdif) == T      ! sand_01_hdiff, ...     horizontal total diffusion<br />  [[Variables#Dout|Dout]](STxdif) == T      ! sand_01_xdiff, ...      horizontal XI-diffusion<br /> [[Variables#Dout|Dout]](STydif) == T      ! sand_01_ydiff, ...     horizontal ETA-diffusion<br /> [[Variables#Dout|Dout]](STsdif) == T      ! sand_01_sdiff, ...      horizontal S-diffusion<br /> [[Variables#Dout|Dout]](STvdif) == T      ! sand_01_vdiff, ...     vertical diffusion</div>


* Lateral, harmonic and biharmonic, constant, diffusion coefficient.<div class="box">      [[tnu2|MUD_TNU2]] == 0.0d0  0.0d0  0.0d0            ! m2/s<br />      [[tnu4|MUD_TNU4]] == 0.0d0  0.0d0  0.0d0            ! m4/s</div>
==Suspended Cohesive Sediment Parameters==


* Vertical diffusion coefficients.<div class="box">  [[Akt_bak|MUD_AKT_BAK]] == 5.0d-6  5.0d-6  5.0d-6          ! m2/s</div>
There are [[Variables#NCS|NCS]] cohesive (mud) sediment tracers per nested grid in ROMS. Therefore, [1:[[Variables#NCS|NCS]],1:[[Variables#Ngrids|Ngrids]]] values are expected for these parameters.


* Nudging/relaxation time scales, inverse scales will be computed internally.<div class="box">     [[Tnudg|MUD_TNUDG]] == 0.0d0 0.0d0 0.0d0              ! days</div>
* Median sediment grain diameter (mm).<div class="box">     [[Variables#sd50|MUD_SD50]] == 0.01d0  0.005d0</div>


* Morphological time scale factor (greater than or equal to 1.0). A value of 1.0 has no scale effect.<div class="box"> [[morph_fac|MUD_MORPH_FAC]] == 1.0d0 1.0d0 1.0d0               ! nondimensional</div>
* Sediment concentration (kg/m<sup>3</sup>).<div class="box">     [[Variables#Csed|MUD_CSED]] == 0.0d0 0.0d0</div>


* Logical switch ('''TRUE'''/'''FALSE''') to activate writing of cohesive sediment into output NetCDF file(s).<div class="box"> [[Hout|Hout(idmud)]] == T                              ! suspended concentration<br /> [[Hout|Hout(iMfrac)]] == T                              ! bed layer fraction<br />  [[Hout|Hout(iMmass)]] == T                              ! bed layer mass<br />  [[Hout|Hout(iMUbld)]] == T                              ! bed load at U-points<br />  [[Hout|Hout(iMVbld)]] == T                              ! bed load at V-points</div>
* Sediment grain density (kg/m<sup>3</sup>).<div class="box">     [[Variables#Srho|MUD_SRHO]] == 2650.0d0 2400.0d0</div>


==Non-cohesive Sediment Parameters==
* Particle settling velocity (mm/s).<div class="box">      [[Variables#Wsed|MUD_WSED]] == 0.1d0  0.01d0</div>


There are [[NNS]] non-cohesive (sand) sediment tracers per nested grid in ROMS. Therefore, [1:[[NNS]],1:[[Ngrids]]] values are expected for these parameters.
* Surface erosion rate (kg/m<sup>2</sup>/s).<div class="box">    [[Variables#Erate|MUD_ERATE]] == 5.0d-4  5.0d-4</div>


* Median sediment grain diameter (mm).<div class="box">     [[sd50|SAND_SD50]] == 1.0d0</div>
* Critical shear for erosion and deposition (N/m<sup>2</sup>).<div class="box">   [[Variables#tau_ce|MUD_TAU_CE]] == 0.1d0  0.1d0<br />    [[Variables#tau_cd|MUD_TAU_CD]] == 0.1d0  0.1d0</div>


* Sediment concentration (kg/m3).<div class="box">    [[Csed|SAND_CSED]] == 0.0d0</div>
* Porosity (nondimensional): <span class="red">Vwater/(Vwater+Vsed)</span>. Its values range between '''0''' and '''1'''.<div class="box">    [[Variables#poros|MUD_POROS]] == 0.9d0  0.9d0</div>


* Sediment grain density (kg/m3).<div class="box">     [[Srho|SAND_SRHO]] == 2650.0d0</div>
* Lateral, constant, harmonic/biharmonic horizontal diffusion of tracer for the nonlinear model and adjoint-based algorithms.<div class="box">     [[Variables#tnu2|MUD_TNU2]] == 0.0d0 0.0d0                    ! m<sup>2</sup>/s<br />      [[Variables#tnu4|MUD_TNU4]] == 0.0d0  0.0d0                    ! m<sup>4</sup>/s<br /><br />  [[Variables#ad_tnu2|ad_MUD_TNU2]] == 0.0d0  0.0d0                    ! m<sup>2</sup>/s<br />  [[Variables#ad_tnu4|ad_MUD_TNU4]] == 0.0d0  0.0d0                    ! m<sup>4</sup>/s</div>


* Particle settling velocity (mm/s).<div class="box">     [[Wsed|SAND_WSED]] == 1.0d0</div>
* Logical switches to increase horizontal diffusivity of cohesive sediment trace in specific areas of the application domain (like sponge areas) for the desired grid.<div class="box">   [[Variables#MUD_sponge|MUD_Sponge]] == F F</div>


* Surface erosion rate (kg/m2/s).<div class="box">   [[Erate|SAND_ERATE]] == 1.0d-5</div>
* Vertical mixing coefficients for tracers in nonlinear model and basic state scale factor in adjoint-based algorithms.<div class="box">   [[Variables#Akt_bak|MUD_AKT_BAK]] == 5.0d-6  5.0d-6                  ! m<sup>2</sup>/s<br />  [[Variables#ad_Akt_fac|MUD_AKT_fac]] == 2*1.0d0                        ! nondimensional</div>


* Critical shear for erosion and deposition (N/m2).<div class="box">   [[tau_ce|SAND_TAU_CE]] == 0.07d0<br />  [[tau_cd|SAND_TAU_CD]] == 0.02d0</div>
* Nudging/relaxation time scales, inverse scales will be computed internally.<div class="box">     [[Variables#Tnudg|MUD_TNUDG]] == 0.0d0 0.0d0                    ! days</div>


* Porosity (nondimensional): <span class="red">Vwater/(Vwater+Vsed)</span>. Its values range between '''0''' and '''1'''.<div class="box">   [[poros|SAND_POROS]] == 0.5d0</div>
* Morphological time scale factor (greater than or equal to 1.0). A value of 1.0 has no scale effect.<div class="box"> [[Variables#morph_fac|MUD_MORPH_FAC]] == 1.0d0 1.0d0                    ! nondimensional</div>


* Lateral, harmonic and biharmonic, constant, diffusion coefficient.<div class="box">    [[tnu2|SAND_TNU2]] == 0.0d0                          ! m2/s<br />    [[tnu4|SAND_TNU4]] == 0.0d0                          ! m4/s</div>
* Logical switch to specify which variables to consider on tracers point Sources/Sinks (like river runoff). [1:[[Variables#NCS|NCS]],[[Variables#Ngrids|Ngrids]]] values expected.
:<div class="box">    [[Variables#LtracerSrc|MUD_Ltsrc]] == F F</div>


* Vertical diffusion coefficients.<div class="box"> [[Akt_bak|SAND_AKT_BAK]] == 5.0d-6                          ! m2/s</div>
*Logical switches to specify which variables to process for tracers climatology: [1:[[Variables#NCS|NCS]],[[Variables#Ngrids|Ngrids]]] values expected.
:<div class="box">     [[Variables#MUD_Ltclm|MUD_Ltclm]] == F F</div>


* Nudging/relaxation time scales, inverse scales will be computed internally.<div class="box">    [[Tnudg|SAND_TNUDG]] == 0.0d0                          ! days</div>
* Logical switches to nudge the desired cohesive sediment tracer climatology field. If not analytical climatology fields, users need to turn on the logical switches above to process the fields from the climatology NetCDF file that are needed for nudging. <div class="box">    [[Variables#MUD_Tnudge|MUD_Tnudge]] == F F</div>


* Morphological time scale factor (greater than or equal to 1.0). A value of 1.0 has no scale effect.<div class="box">[[morph_fac|SAND_MORPH_FAC]] == 1.0d0                          ! nondimensional</div>
* Logical switches to activate writing of cohesive sediment into output history file.
:<div class="box"> [[Variables#Hout|Hout]](idmud)  == T      ! mud_01, ...             suspended concentration<br />  [[Variables#Hout|Hout]](iMfrac) == T      ! mudfrac_01, ...        bed layer fraction<br />  [[Variables#Hout|Hout]](iMmass) == T      ! mudmass_01, ...        bed layer mass<br />  [[Variables#Hout|Hout]](iMUbld) == T      ! bedload_Umud_01, ...    bed load at U-points<br />  [[Variables#Hout|Hout]](iMVbld) == T      ! bedload_Vmud_01, ...    bed load at V-points</div>


* Logical switches ('''TRUE'''/'''FALSE''') to activate writing of non-cohesive sediment fields into output NetCDF file(s).<div class="box">  [[Hout|Hout(isand)]] == T                              ! suspended concentration<br />  [[Hout|Hout(iSfrac)]] == T                              ! bed layer fraction<br />  [[Hout|Hout(iSmass)]] == T                              ! bed layer mass<br />  [[Hout|Hout(iSUbld)]] == T                              ! bed load at U-points<br />  [[Hout|Hout(iSVbld)]] == T                              ! bed load at V-points</div>
* Logical switches to activate writing of cohesive sediment into output quicksave file.
:<div class="box">  [[Variables#Qout|Qout]](idmud)  == F      ! mud_01, ...            suspended concentration<br />  [[Variables#Qout|Qout]](iSmud) == F      ! mud_01_sur, ...        surface concentration<br />  [[Variables#Qout|Qout]](iMfrac) == F      ! mudfrac_01, ...        bed layer fraction<br />  [[Variables#Qout|Qout]](iMmass) == F      ! mudmass_01, ...        bed layer mass<br />  [[Variables#Qout|Qout]](iMUbld) == F      ! bedload_Umud_01, ...    bed load at U-points<br />  [[Variables#Qout|Qout]](iMVbld) == F      ! bedload_Vmud_01, ...    bed load at V-points</div>


==Bed layer and bottom sediment parameters==
* Logical switches to activate writing of time-averaged, cohesive sediment fields into averages output file.
:<div class="box">  [[Variables#Aout|Aout]](idmud)  == T      ! mud_01, ...            suspended concentration<br /><br />  [[Variables#Aout|Aout]](iMTTav) == T      ! mud_01_2, ...          quadratic <t*t> tracer terms<br />  [[Variables#Aout|Aout]](iMUTav) == T      ! u_mud_01, ...          quadratic <u*t> tracer terms<br />  [[Variables#Aout|Aout]](iMVTav) == T      ! v_mud_01, ...          quadratic <v*t> tracer terms<br />  [[Variables#Aout|Aout]](MHUTav) == T      ! Huon_mud_01, ...        tracer volume flux, <Huon*t><br />  [[Variables#Aout|Aout]](MHVTav) == T      ! Hvom_mud_01, ...        tracer volume flux, <Hvom*t><br /><br />  [[Variables#Aout|Aout]](iMUbld) == T      ! bedload_Umud_01, ...    bed load at U-points<br />  [[Variables#Aout|Aout]](iMVbld) == T      ! bedload_Vmud_01, ...    bed load at V-points</div>


* Logical switches ('''TRUE'''/'''FALSE''') to activate writing of bed layer and bottom sediment fields into output NetCDF file(s), [1:[[Ngrids]]] values are expected.<div class="box">   [[Hout|Hout(ithck)]] == T                               ! sediment layer thickness<br />   [[Hout|Hout(iaged)]] == T                               ! sediment layer age<br />   [[Hout|Hout(iporo)]] == T                               ! sediment layer porosity<br />   [[Hout|Hout(idiff)]] == F                              ! biodiffusivity</div>
* Logical switches to activate writing of time-averaged, cohesive sediment diagnostic terms into output diagnostics file.
:<div class="box">  [[Variables#Dout|Dout]](MTrate) == T      ! mud_01_rate, ...        time rate of change<br />  [[Variables#Dout|Dout]](MThadv) == T      ! mud_01_hadv, ...       horizontal total advection<br /> [[Variables#Dout|Dout]](MTxadv) == T      ! mud_01_xadv, ...        horizontal XI-advection<br />  [[Variables#Dout|Dout]](MTyadv) == T       ! mud_01_yadv, ...        horizontal ETA-advection<br /> [[Variables#Dout|Dout]](MTvadv) == T      ! mud_01_vadv, ...        vertical advection<br />  [[Variables#Dout|Dout]](MThdif) == T       ! mud_01_hdiff, ...      horizontal total diffusion<br /> [[Variables#Dout|Dout]](MTxdif) == T      ! mud_01_xdiff, ...      horizontal XI-diffusion<br />  [[Variables#Dout|Dout]](MTydif) == T       ! mud_01_ydiff, ...      horizontal ETA-diffusion<br /> [[Variables#Dout|Dout]](MTsdif) == T      ! mud_01_sdiff, ...      horizontal S-diffusion<br />  [[Variables#Dout|Dout]](MTvdif) == T      ! mud_01_vdiff, ...      vertical diffusion</div>

Latest revision as of 17:51, 16 January 2020

Sediment Model Input Script - sediment.in

The sediment.in file sets the parameters for sediment output. The name of this file is set by the SPARNAM keyword in the roms.in file. This standard input ASCII file is organized in several sections as shown below, with links to more detailed explanation where required.


Note Notice: Detailed information about ROMS input script file syntax can be found here.

Note Notice: A default sediment.in input script is provided in the User/External subdirectory. Also there are several sediment input scripts in the ROMS/External subdirectory which are used in the distributed test cases. They are usually named sediment_app.in where app is the lowercase of the test case CPP option.


Sediment Model Control Switches

  • Switch is used to control sediment model computation within nested and/or multiple connected grids, [1:Ngrids] values are expected.
    Lsediment == T

Tracer Advection Schemes

  • Set horizontal and vertical advection schemes for biological tracers. A different advection scheme is allowed for each tracer. For example, a positive-definite (monotonic) algorithm can be activated for salinity and biological tracers, while a different one is set for temperature.
    Keyword Advection Algorithm

    A4 4th-order Akima (horizontal/vertical)
    C2 2nd-order centered differences (horizontal/vertical)
    C4 4th-order centered differences (horizontal/vertical)
    HSIMT 3th-order HSIMT with TVD limiter (horizontal/vertical)
    MPDATA recursive flux corrected MPDATA (horizontal/vertical)
    SPLINES parabolic splines reconstruction (only vertical)
    SU3 split third-order upstream (horizontal/vertical)
    U3 3rd-order upstresm-bias (only horizontal)
    The user has the option of specifying the full Keyword or the first two letters, regardless if using uppercase or lowercase. If nested grids, specify values for each grid. For more details please read the roms.in Tracer Advection Schemes section. [1:NBT,Ngrids] values are expected.
    Hadvection == HSIMT  ! idsed(:), compact

    Vadvection == HSIMT  ! idsed(:), compact

Lateral Open Boundary Conditions Parameters

  • The lateral boundary conditions are entered with a keyword. A value is expected for each boundary segment per nested grid for each state variable. The sediment tracer variables (cohesive and non-cohesive) require [1:4,1:NCS+NNS,Ngrids] values. If specifying every tracer, enter first 1:NCS cohesive sediment entries followed by 1:NNS non-cohesive sediment entries. The boundary order is: 1=west, 2=south, 3=east, and 4=north. That is, anticlockwise starting at the western boundary.

    The keyword is case insensitive and usually has three characters. However, it is possible to have compound keywords, if applicable. For example, the keyword RadNud implies radiation boundary condition with nudging. This combination is usually used in active/passive radiation conditions.

    NoteNotice: It is possible to specify the lateral boundary conditions for all biological tracers in a compact form with a single entry. If so, all the biological tracers are assumed to have the same boundary condition as the single entry.
    ! Keyword Lateral Boundary Condition Type
    !
    ! Cla Clamped _____N_____ j=Mm
    ! Clo Closed | 4 |
    ! Gra Gradient | |
    ! Nes Nested 1 W E 3
    ! Nud Nudging | |
    ! Per Periodic |_____S_____|
    ! Rad Radiation 2 j=1
    ! i=1 i=Lm
    ! W S E N
    ! e o a o
    ! s u s r
    ! t t t t
    ! h h
    !
    ! 1 2 3 4

    LBC(isTvar) == Per Clo Per Clo  ! idsed(:), compact
  • Adjoint-based algorithms can have different lateral boundary conditions keywords.
    ad_LBC(isTvar) == Per Clo Per Clo  ! idsed(:), compact

General Sediment Bed Model Controls

  • Depositional bed layer thickness (m) criteria to create a new layer. If deposition exceeds this value, then a new layer is created, [1:Ngrids] values are expected.
    NEWLAYER_THICK == 0.01d0

Bed Layer and Bottom Sediment Parameters

  • Logical switches to activate writing of bed layer and bottom sediment fields into history output NetCDF file(s), [1:Ngrids] values are expected.
Hout(ithck) == T  ! bed_thickness sediment layer thickness
Hout(iaged) == T  ! bed_age sediment layer age
Hout(iporo) == T  ! bed_porosity sediment layer porosity
Hout(idiff) == F  ! bed_biodiff biodiffusivity

Hout(ids50) == T  ! grain_diameter mean grain diameter
Hout(idens) == T  ! grain_density mean grain density
Hout(iwsed) == T  ! settling_vel mean settling velocity
Hout(itauc) == T  ! erosion_stress critical erosion stress
Hout(irlen) == T  ! ripple_length ripple length
Hout(irhgt) == T  ! ripple_height ripple height
Hout(ibwav) == T  ! bed_wave_amp wave excursion amplitude
Hout(izdef) == T  ! Zo_def default bottom roughness
Hout(izapp) == T  ! Zo_app apparent bottom roughness
Hout(izNik) == T  ! Zo_Nik Nikuradse bottom roughness
Hout(izbio) == T  ! Zo_bio biological bottom roughness
Hout(izbfm) == T  ! Zo_bedform bed form bottom roughness
Hout(izbld) == T  ! Zo_bedload bed load bottom roughness
Hout(izwbl) == T  ! Zo_wbl wave bottom roughness
Hout(iactv) == T  ! active_layer_thickness active layer thickness
Hout(ishgt) == T  ! saltation saltation height
  • Logical switches to activate writing of bed layer and bottom sediment fields into quicksave output NetCDF file(s), [1:Ngrids] values are expected.
Qout(ithck) == F  ! bed_thickness sediment layer thickness
Qout(iaged) == F  ! bed_age sediment layer age
Qout(iporo) == F  ! bed_porosity sediment layer porosity
Qout(idiff) == F  ! bed_biodiff biodiffusivity

Qout(ids50) == F  ! grain_diameter mean grain diameter
Qout(idens) == F  ! grain_density mean grain density
Qout(iwsed) == F  ! settling_vel mean settling velocity
Qout(itauc) == F  ! erosion_stress critical erosion stress
Qout(irlen) == F  ! ripple_length ripple length
Qout(irhgt) == F  ! ripple_height ripple height
Qout(ibwav) == F  ! bed_wave_amp wave excursion amplitude
Qout(izdef) == F  ! Zo_def default bottom roughness
Qout(izapp) == F  ! Zo_app apparent bottom roughness
Qout(izNik) == F  ! Zo_Nik Nikuradse bottom roughness
Qout(izbio) == F  ! Zo_bio biological bottom roughness
Qout(izbfm) == F  ! Zo_bedform bed form bottom roughness
Qout(izbld) == F  ! Zo_bedload bed load bottom roughness
Qout(izwbl) == F  ! Zo_wbl wave bottom roughness
Qout(iactv) == F  ! active_layer_thickness active layer thickness
Qout(ishgt) == F  ! saltation saltation height


Non-cohesive Sediment Parameters

There are NNS non-cohesive (sand) sediment tracers per nested grid in ROMS. Therefore, [1:NNS,1:Ngrids] values are expected for these parameters.

  • Median sediment grain diameter (mm).
    SAND_SD50 == 1.0d0
  • Sediment concentration (kg/m3).
    SAND_CSED == 0.0d0
  • Sediment grain density (kg/m3).
    SAND_SRHO == 2650.0d0
  • Particle settling velocity (mm/s).
    SAND_WSED == 1.0d0
  • Surface erosion rate (kg/m2/s).
    SAND_ERATE == 5.0d-4
  • Porosity (nondimensional): Vwater/(Vwater+Vsed). Its values range between 0 and 1.
    SAND_POROS == 0.5d0
  • Lateral, constant, harmonic/biharmonic horizontal diffusion of tracer for the nonlinear model and adjoint-based algorithms.
    SAND_TNU2 == 0.0d0  ! m2/s
    SAND_TNU4 == 0.0d0  ! m4/s

    ad_SAND_TNU2 == 0.0d0  ! m2/s
    ad_SAND_TNU4 == 0.0d0  ! m4/s
  • Logical switches to increase horizontal diffusivity of noncohesive sediment trace in specific areas of the application domain (like sponge areas) for the desired grid.
    SAND_Sponge == F
  • Vertical mixing coefficients for tracers in nonlinear model and basic state scale factor in adjoint-based algorithms.
    SAND_AKT_BAK == 5.0d-6  ! m2/s
    SAND_AKT_fac == 1.0d0  ! nondimensional
  • Nudging/relaxation time scales, inverse scales will be computed internally.
    SAND_TNUDG == 0.0d0  ! days
  • Morphological time scale factor (greater than or equal to 1.0). A value of 1.0 has no scale effect.
    SAND_MORPH_FAC == 1.0d0  ! nondimensional
  • Logical switch to specify which variables to consider on tracers point Sources/Sinks (like river runoff). [1:NCS,Ngrids] values expected.
SAND_Ltsrc == F
  • Logical switches to specify which variables to process for tracers climatology: [1:NCS,Ngrids] values expected.
SAND_Ltclm == F
  • Logical switches to nudge the desired noncohesive sediment tracer climatology field. If not analytical climatology fields, users need to turn on the logical switches above to process the fields from the climatology NetCDF file that are needed for nudging.
    SAND_Tnudge == F
  • Logical switches to activate writing of non-cohesive sediment into output history file.
Hout(idsand) == T  ! sand_01, ... suspended concentration
Hout(iSfrac) == T  ! sandfrac_01, ... bed layer fraction
Hout(iSmass) == T  ! sandmass_01, ... bed layer mass
Hout(iSUbld) == T  ! bedload_Usand_01, ... bed load at U-points
Hout(iSVbld) == T  ! bedload_Vsand_01, ... bed load at V-points
  • Logical switches to activate writing of non-cohesive sediment into output quicksave file.
Qout(idsand) == F  ! sand_01, ... suspended concentration
Qout(iSsand) == F  ! sand_01_sur, ... surface concentration
Qout(iSfrac) == F  ! sandfrac_01, ... bed layer fraction
Qout(iSmass) == F  ! sandmass_01, ... bed layer mass
Qout(iSUbld) == F  ! bedload_Usand_01, ... bed load at U-points
Qout(iSVbld) == F  ! bedload_Vsand_01, ... bed load at V-points
  • Logical switches to activate writing of time-averaged, non-cohesive sediment fields into averages output file.
Aout(idsand) == T  ! sand_01, ... suspended concentration

Aout(iSTTav) == T  ! sand_01_2, ... quadratic <t*t> tracer terms
Aout(iSUTav) == T  ! u_sand_01, ... quadratic <u*t> tracer terms
Aout(iSVTav) == T  ! v_sand_01, ... quadratic <v*t> tracer terms
Aout(SHUTav) == T  ! Huon_sand_01, ... tracer volume flux, <Huon*t>
Aout(SHVTav) == T  ! Hvom_sand_01, ... tracer volume flux, <Hvom*t>

Aout(iSUbld) == T  ! bedload_Usand_01, ... bed load at U-points
Aout(iSVbld) == T  ! bedload_Vsand_01, ... bed load at V-points
  • Logical switches to activate writing of time-averaged, non-cohesive sediment diagnostic terms into output diagnostics file.
Dout(STrate) == T  ! sand_01_rate, ... time rate of change
Dout(SThadv) == T  ! sand_01_hadv, ... horizontal total advection
Dout(STxadv) == T  ! sand_01_xadv, ... horizontal XI-advection
Dout(STyadv) == T  ! sand_01_yadv, ... horizontal ETA-advection
Dout(STvadv) == T  ! sand_01_vadv, ... vertical advection
Dout(SThdif) == T  ! sand_01_hdiff, ... horizontal total diffusion
Dout(STxdif) == T  ! sand_01_xdiff, ... horizontal XI-diffusion
Dout(STydif) == T  ! sand_01_ydiff, ... horizontal ETA-diffusion
Dout(STsdif) == T  ! sand_01_sdiff, ... horizontal S-diffusion
Dout(STvdif) == T  ! sand_01_vdiff, ... vertical diffusion

Suspended Cohesive Sediment Parameters

There are NCS cohesive (mud) sediment tracers per nested grid in ROMS. Therefore, [1:NCS,1:Ngrids] values are expected for these parameters.

  • Median sediment grain diameter (mm).
    MUD_SD50 == 0.01d0 0.005d0
  • Sediment concentration (kg/m3).
    MUD_CSED == 0.0d0 0.0d0
  • Sediment grain density (kg/m3).
    MUD_SRHO == 2650.0d0 2400.0d0
  • Particle settling velocity (mm/s).
    MUD_WSED == 0.1d0 0.01d0
  • Surface erosion rate (kg/m2/s).
    MUD_ERATE == 5.0d-4 5.0d-4
  • Critical shear for erosion and deposition (N/m2).
    MUD_TAU_CE == 0.1d0 0.1d0
    MUD_TAU_CD == 0.1d0 0.1d0
  • Porosity (nondimensional): Vwater/(Vwater+Vsed). Its values range between 0 and 1.
    MUD_POROS == 0.9d0 0.9d0
  • Lateral, constant, harmonic/biharmonic horizontal diffusion of tracer for the nonlinear model and adjoint-based algorithms.
    MUD_TNU2 == 0.0d0 0.0d0  ! m2/s
    MUD_TNU4 == 0.0d0 0.0d0  ! m4/s

    ad_MUD_TNU2 == 0.0d0 0.0d0  ! m2/s
    ad_MUD_TNU4 == 0.0d0 0.0d0  ! m4/s
  • Logical switches to increase horizontal diffusivity of cohesive sediment trace in specific areas of the application domain (like sponge areas) for the desired grid.
    MUD_Sponge == F F
  • Vertical mixing coefficients for tracers in nonlinear model and basic state scale factor in adjoint-based algorithms.
    MUD_AKT_BAK == 5.0d-6 5.0d-6  ! m2/s
    MUD_AKT_fac == 2*1.0d0  ! nondimensional
  • Nudging/relaxation time scales, inverse scales will be computed internally.
    MUD_TNUDG == 0.0d0 0.0d0  ! days
  • Morphological time scale factor (greater than or equal to 1.0). A value of 1.0 has no scale effect.
    MUD_MORPH_FAC == 1.0d0 1.0d0  ! nondimensional
  • Logical switch to specify which variables to consider on tracers point Sources/Sinks (like river runoff). [1:NCS,Ngrids] values expected.
MUD_Ltsrc == F F
  • Logical switches to specify which variables to process for tracers climatology: [1:NCS,Ngrids] values expected.
MUD_Ltclm == F F
  • Logical switches to nudge the desired cohesive sediment tracer climatology field. If not analytical climatology fields, users need to turn on the logical switches above to process the fields from the climatology NetCDF file that are needed for nudging.
    MUD_Tnudge == F F
  • Logical switches to activate writing of cohesive sediment into output history file.
Hout(idmud) == T  ! mud_01, ... suspended concentration
Hout(iMfrac) == T  ! mudfrac_01, ... bed layer fraction
Hout(iMmass) == T  ! mudmass_01, ... bed layer mass
Hout(iMUbld) == T  ! bedload_Umud_01, ... bed load at U-points
Hout(iMVbld) == T  ! bedload_Vmud_01, ... bed load at V-points
  • Logical switches to activate writing of cohesive sediment into output quicksave file.
Qout(idmud) == F  ! mud_01, ... suspended concentration
Qout(iSmud) == F  ! mud_01_sur, ... surface concentration
Qout(iMfrac) == F  ! mudfrac_01, ... bed layer fraction
Qout(iMmass) == F  ! mudmass_01, ... bed layer mass
Qout(iMUbld) == F  ! bedload_Umud_01, ... bed load at U-points
Qout(iMVbld) == F  ! bedload_Vmud_01, ... bed load at V-points
  • Logical switches to activate writing of time-averaged, cohesive sediment fields into averages output file.
Aout(idmud) == T  ! mud_01, ... suspended concentration

Aout(iMTTav) == T  ! mud_01_2, ... quadratic <t*t> tracer terms
Aout(iMUTav) == T  ! u_mud_01, ... quadratic <u*t> tracer terms
Aout(iMVTav) == T  ! v_mud_01, ... quadratic <v*t> tracer terms
Aout(MHUTav) == T  ! Huon_mud_01, ... tracer volume flux, <Huon*t>
Aout(MHVTav) == T  ! Hvom_mud_01, ... tracer volume flux, <Hvom*t>

Aout(iMUbld) == T  ! bedload_Umud_01, ... bed load at U-points
Aout(iMVbld) == T  ! bedload_Vmud_01, ... bed load at V-points
  • Logical switches to activate writing of time-averaged, cohesive sediment diagnostic terms into output diagnostics file.
Dout(MTrate) == T  ! mud_01_rate, ... time rate of change
Dout(MThadv) == T  ! mud_01_hadv, ... horizontal total advection
Dout(MTxadv) == T  ! mud_01_xadv, ... horizontal XI-advection
Dout(MTyadv) == T  ! mud_01_yadv, ... horizontal ETA-advection
Dout(MTvadv) == T  ! mud_01_vadv, ... vertical advection
Dout(MThdif) == T  ! mud_01_hdiff, ... horizontal total diffusion
Dout(MTxdif) == T  ! mud_01_xdiff, ... horizontal XI-diffusion
Dout(MTydif) == T  ! mud_01_ydiff, ... horizontal ETA-diffusion
Dout(MTsdif) == T  ! mud_01_sdiff, ... horizontal S-diffusion
Dout(MTvdif) == T  ! mud_01_vdiff, ... vertical diffusion
Retrieved from "https://www.myroms.org/wiki/index.php?title=sediment.in&oldid=5761"