Latest

Main page
Online manual
Applications
Latest
Publications
User group
Useful links

Dep.Mat.UoB
CMU
IMR

Downloads:

BOM 5.0 ALPHA4:
download the code
BOM 5.0 ALPHA3:
download the code
BOM 5.0 ALPHA2:
download the code
BOM 4.1:
download the code
userguide: PDF -->

2011-07-15: BOM 5 alpha4

• Initial, but functional, support for the new TEOS-10 equation of state.
• Add ability to specify river concentrations, and support for rivers when advecting density
• Use a global shapiro filter instead of TANH transitions to smooth edges of sponge regions
• Add ability to override river inlet geometries on the grid level for full control.
• Add support for dived river inlets
• Add support for reading binary data in "binary stream" format, more compatible with the C world.
• Some minor cleanups, bug fixes and improvements.

• The automatic conversion between absolute salinity and PSU is not 100% complete.
• 1way nesting is in the pipeline
• Documentation should be transferred to LaTeX and PDF

2010-12-16: BOM 5 alpha3

• More support for reading lateral and surface BC's from file
• New keyword for permanently increasing viscosity in specific regions ("add_sponge")
• Improved the info available via the command "./bomhelp.py"
• Support for settling velocity in the particle model
• Improved robustness of timeseries dumping technique, did not work well on > 256 processors.
• A lot of minor cleanup, bug fixes and improvements.

• More tuning of boundary conditions
• Almost done; support for easily doing 1way nesting using the new BC code above.
• If a river with code 3 or 4 is located at a processor boundary the ghost zones are not correctly updated. the trivial fix is not good for performance, so this will be fixed if it is really needed only.
• Sync documentation with code. Mostly ok, except coordinate change and some new keywords.
• Documentation should be transferred to LaTeX and PDF

2010-07-16: BOM 5 alpha2

• A few major and minor bug fixes w.r. binreg averaging BC, plotting.
• Overhauled the advection of tracers.
• Added capability to plot array slices "as is" in pgplot, optionally each time "dump_state" is called. this means changes in most arrays can be visualized during a timestep, which is nice for debugging.
• The switch from specification of horisontal ranges using grid indices into x,y coordinates in meters has begun and is close to complete.

• More tuning of boundary conditions
• More methods for forcing will come eventually, on a needed by project basis.
• We have some known performance issues on more than 100processors with the binreg boundary conditions that can be fixed by using communicators instead of global collective calls. will be done if needed.
• If a river with code 3 or 4 is located at a processor boundary the ghost zones are not correctly updated. the trivial fix is not good for performance, so this will be fixed if it is really needed only.
• Sync documentation with code. Mostly ok, except coordinate change and some new keywords.
• Documentation should be transferred to LaTeX and PDF

2010-06-04: BOM 5 alpha

• Non hydrostatic pressure solver
• 2D vertical slice mode via setting JM=1 is integrated
• MPI parallelisation
• Checkpoint/restart facility
• New configuration system, moving almost all case specific config into setup file. Idealized winds, rivers, point sources, plotting, initialization, some simple types of boundary conditions and more are working.
• User configurable dumping of data to ascii files
• Run time plotting via PGPLOT

It is a goal to let all features be turned on or off via the setup file.

There will be more features added on a "needed by project" basis. There is a framework in place for offering user selectable boundary conditions. As of today FRS, Dirichlet, Neuman=0, and periodic is implemented. More will come for realistic forcing, nesting etc.

Scalability for certain choices of open boundary conditions need optimization for use on more than 30-40 processors depending on case size.

A particle tracking module is under development and testing.

One upcoming change is that we will switch from specifying coordinates in the plotting commands with indices and "grid coordinates" into using meters as the unit. this will make config files more generic, and enable very quick changes of experiment resolution.

The code is tested with Portland (pgf90) and the GNU Fortran compiler (gfortran) with MPICH2. Intel ifort does not yet work well with optimization on for all files.

2005-05-06: Internal pressure and rotated grids

This new algorithm can easily be implemented in the BOM 4.1 version by following the below steps:

In the file state.f90 add the following line

     REAL, DIMENSION(IM,JM) :: DCORNER
     REAL, DIMENSION(IM,JM,KB) :: DELTAUROT, DELTAVROT

In the BOM source file Main.f90, glue in the below lines around each call to the INTERNAL subroutine:

      ...
      CALL INTERNALROT
      DELTAUROT = DELTAU
      DELTAVROT = DELTAV

      CALL INTERNAL

      DELTAU = (2.5*DELTAU+DELTAUROT)/3.5
      DELTAV = (2.5*DELTAV+DELTAVROT)/3.5
      ...

Download the file at this link into your BOM source directory, call it e.g. internalrot.f90. Then, in the Makefile add internalrot.o to the OBJECTS variable and recompile.

2004-09-01: BOM 4.1 released

It is very hard to see differences in the results after rerunning experiements and comparing with version 4.0, but this may be case dependant. The 2D solution seems to be non-affected, while e.g. the stratification seems somewhat improved. Generally the bug has lead to a damping of the solution, so it will generally be necessary to increase the viscosity up to levels corresponding to 3.x, likewise reduce the time step, in order to get stable results.

2003-11-18: BOM 4.0 released

The alternative internal pressures scheme (due to McCalpin) with fourth order discretization along sigma layers has been made a runtime option. By applying the integration by parts rule to the internal pressure term it is possible to do the vertical integral in a more robust way. This scheme has been proven effective to reduce internal pressure error in the infamous Seamount case. It is not, however, clear how beneficial the effects are in more realistic cases.

The numerical kernel of BOM seems to stabilize, so it is not expected that there will be major changes to the core in the future. A non-hydrostatic version of the modesplit routine is under testing, and may be added soon, and depending on funding there may come a separate MPI version of 4.0 next year.

New animations of the contintental shelf on the examplespage

2003-03-03: BOM 3.2 released

Also updated Makefile with flags for Intel Fortran and Pentium4.

2002-09-06: Small bug in 3.0 fixed

2002-08-20: BOM 3.0 released

New features:

• Predictor/corrector time-stepping scheme, with the theta-method on the corrector step. Hopefully this approach will be useful for studies with low viscosity, and prevent/reduce unphysical grovth of certain modes, e.g. from errors in internal pressure calculations.
• Additional recent advection scheme due to Shchepetkin. Gives sharp fronts, but may have some overshoot.