sediment.in: Difference between revisions
From WikiROMS
Jump to navigationJump to search
No edit summary (change visibility) |
No edit summary (change visibility) |
||
| (7 intermediate revisions by the same user 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 <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 [[ | 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. | ||
| Line 13: | Line 13: | ||
* 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> | * 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> | ||
== | ==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== | ==Non-cohesive Sediment Parameters== | ||
| Line 71: | Line 54: | ||
* Surface erosion rate (kg/m<sup>2</sup>/s).<div class="box"> [[Variables#Erate|SAND_ERATE]] == 5.0d-4</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. | * 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> | * 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> | ||
| Line 77: | Line 60: | ||
* 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> | * 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> | ||
* 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 ! | * 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> | |||
* Nudging/relaxation time scales, inverse scales will be computed internally.<div class="box"> [[Variables#Tnudg|SAND_TNUDG]] == 0.0d0 ! days</div> | * Nudging/relaxation time scales, inverse scales will be computed internally.<div class="box"> [[Variables#Tnudg|SAND_TNUDG]] == 0.0d0 ! days</div> | ||
| Line 83: | Line 68: | ||
* 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> | * 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> | ||
* Logical switch to specify which variables to consider on tracers point Sources/Sinks (like river runoff).<div class="box"> | * 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> | |||
*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> | |||
* 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> | |||
* Logical switches to activate writing of non-cohesive sediment into output history file. | * Logical switches to activate writing of non-cohesive sediment into output history file. | ||
:<div class="box"> [[Variables#Hout|Hout(idsand) | :<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> | ||
* 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> | |||
* Logical switches to activate writing of time-averaged, non-cohesive sediment fields into averages output file. | * Logical switches to activate writing of time-averaged, non-cohesive sediment fields into averages output file. | ||
:<div class="box"> [[Variables#Aout|Aout(idsand)]] == T | :<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> | ||
* Logical switches to activate writing of time-averaged, non-cohesive sediment diagnostic terms into output diagnostics file. | * Logical switches to activate writing of time-averaged, non-cohesive sediment diagnostic terms into output diagnostics file. | ||
:<div class="box"> [[Variables#Dout|Dout(STrate) | :<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> | ||
==Suspended Cohesive Sediment Parameters== | |||
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. | |||
* Median sediment grain diameter (mm).<div class="box"> [[Variables#sd50|MUD_SD50]] == 0.01d0 0.005d0</div> | |||
* Sediment concentration (kg/m<sup>3</sup>).<div class="box"> [[Variables#Csed|MUD_CSED]] == 0.0d0 0.0d0</div> | |||
* Sediment grain density (kg/m<sup>3</sup>).<div class="box"> [[Variables#Srho|MUD_SRHO]] == 2650.0d0 2400.0d0</div> | |||
* Particle settling velocity (mm/s).<div class="box"> [[Variables#Wsed|MUD_WSED]] == 0.1d0 0.01d0</div> | |||
* Surface erosion rate (kg/m<sup>2</sup>/s).<div class="box"> [[Variables#Erate|MUD_ERATE]] == 5.0d-4 5.0d-4</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> | |||
* 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> | |||
* 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> | |||
== | * 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> | ||
* 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> | |||
* Nudging/relaxation time scales, inverse scales will be computed internally.<div class="box"> [[Variables#Tnudg|MUD_TNUDG]] == 0.0d0 0.0d0 ! days</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> | |||
* 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> | |||
*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> | |||
* 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> | |||
* 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 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> | |||
* 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 to activate writing of | * 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.
Notice: Detailed information about ROMS input script file syntax can be found here.
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 AlgorithmThe 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.
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)
- Adjoint-based algorithms can have different horizontal and vertical schemes for biological tracers.[1:NBT,Ngrids] values are expected.
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.
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.! 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.
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 load transport rate coefficient, [1:Ngrids] values are expected.BEDLOAD_COEFF == 0.05d0
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
- Critical shear for erosion and deposition (N/m2).
- 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.
- 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.
- 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).
- 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.
- 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.
- 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
