I am considering this blog as a location for my notes on major model advancements and as an additional means to track what has been incorporated into the releases.
Subversion version control.
The idea is that we start with a distributed version of ROMS, roms_v2.2.
This was our starting point and we entered roms_v2.2 into SVN version control. As we make modifications for the sediment advancements, we will make new (internal) releases and use the second dot (which is the third number), such as
roms_sed_v2.2.1, v2.2.2, v2.2.3, etc etc. As more features become available, and as we gain confidence in the third number, we will merge back to the main roms trunk and make roms_v2.3.
*******************************************************************
Starting with roms_v2.2 we have added capabilites for:
– ssw_bbl : bottom boundary layer model to account for increased bottom stress due to wave+currents
– suspended sediment transport : minor additions to account for multiple sediment classes and new input file structures.
– bedload transport : included 2 methods 1) Meyer-Peter Mueller and 2) Soulsby +Damgaard methodologies. Both use bottom stress for waves + currents but the SD method is more applicable to wave environemnts. The MPM method was developed just for currents.
– bed model : 3D bed strucutre that tracks stratigraphy, multiple sed class distributions, resuspension + deposition interaction with suspended sed.
– morphology : as bed evolves, the changing sea floor elevation is used as boundary condition to omega. Entire model knows about changing bed elevation (if activated).
– wetting/drying : this is based on cell face blocking methodology. The user specifies a ‘critical depth.’ During computation, if the total water depth (h+zeta) at a rho point is less than the critical depth, then transport is prevented out of that cell. That is it. Water is always allowed to flow into a cell. Cells with an initial rho_mask = 0 will always be dry.
–surface tke fluxes : there are 2 methods to account for flux of turbulent kinetic energy from breaking surface waves 1) Craig/Banner with Charnok coefficient that is based on a surface wind stress and 2) wave energy dissipation that is based on wave breaking. Both methods are incorporated into the GLS turbulence closure model.
– wave/current interaction : this is currently based on Mellor 2003; 2005 JPO and adds ‘radiation stress terms’ to the momentum equations. These terms add forcing due to momentum flux from waves. The formulation requires information of wave height, length, and direction.
– model coupling : we use the Model Coupling Toolkit to couple ROMS to SWAN (v40.41AB). Users need to install and compile MCT to make the librarires for linking. The model is driven with the ROMS/External/coupling*.in file.
These features now comprise our version roms_sed_v2.2.2.
*******************************************************************
roms_sed_v2.2.3
– minor updates to several test cases
*******************************************************************
roms_sed_v2.2.4
– update to new directory structure of
— src
—-Compilers
—-Lib
—-Master
—-ROMS
——Adjoint
——Bin
——Drivers
——External
——Include
——Modules
——Nonlinear
——Programs
——Representer
——SeaIce
——Tangent
——Utility
—-SWAN