Stony Brook Parallel Ocean Model

What's new

  • 10-Sep-2011: New sbPOM version including support for large files (greater than 2 Gb) and 4th order pressure gradient scheme (see parameter npg in pom.nml). In addition, results are now identical (to the last bit) to the sequential POM. You can download the new version for the three test cases (seamount, shallowtide, and ecoast) from the codes page.
  • 28-May-2010: Elimination of the optimization option (-O3) in the makefile. We have experienced problems with optimization in some systems. You can use the optimization option, but be careful about using an aggressive optimization.
  • 25-May-2010: Update of the US east coast test case. The call to subroutine ztosig (in initialize.f) was wrong, although it did not cause any problem because nz=kb in this tets case.
  • 25-May-2010: Update of the shallow tide test case. Boundary conditions are specified in bounds_forcing.f. Initial temperature and salinity files are in z coordinate. A Matlab© function (ellipse.m) is included to plot ellipses for tidal analysis.
  • 24-May-2010: Update of the seamount test case. Boundary conditions are specified in bounds_forcing.f. Some variables in subroutine profq are updated at the border of local domains. Initial temperature and salinity files are in z coordinate.
  • 24-May-2010: In the Matlab© files (in prep and proc folders), command "ncmex" have been changed to "mexnc" because of the new mexcdf package. This modification has been also applied to the pomviz package. Download the new pomviz package HERE. Thanks to Yanping Feng.
  • 7-May-2010: New sbPOM test case for the East coast of U.S.
  • 05-Apr-2010: Several bugs corrected in the shallow tides test case. Beacuse density is a function of pressure, the corresponding line in subroutine profq is un-commented. Boundary flags are initialized at subroutine initial_conditions (thanks to Josh). The variable edge(3) is set to k instead of kb in subroutines read_clim_ts_pnetcdf and read_initial_ts_pnetcdf.
  • 21-Mar-2010: New sbPOM test case for tides in a shallow basin. In addition, several parts of the code are updated. Some variables in subroutine profq are updated at the border of local domains. Momentum fluxes at the bottom are also updated at the border of local domains (subroutines advave, profu and profv). The sponge layer in lateral viscosity (subroutine lateral_viscosity) is supressed. Vertical kinematic viscosity and diffusivity is included in the output (subroutine write_output_pnetcdf).
  • 05-Mar-2010: Update of subroutine ztosig and bcond. Several errors associated with the definition of variables in subroutine ztosig caused that temperature and salinity were zero. With correct values of temperature and salinity, the model blows up due to the boundary conditions. So, boundary conditions are also modified (subrountine bcond) to make the code stable.
  • 19-Oct-2009: Update of subroutine dens. Density was only a function of temperature. Density is now a function of temperature and salinity according to the UNESCO equation of state.
  • 06-Oct-2009: Update of subroutine ztosig. The interpolation is now performed in the (global) boundaries.
  • 01-Oct-2009: Update of matlab function pom_vertical_coords (in prep directory). When a uniform spaced vertical grid was generated, the length of z was kb+1 instead of kb.
  • 22-Jul-2009: Update of subroutine advt2. The values for net advective fluxes are now updated at the border of local domains. Prior to this update, numerical 'noise' in the model results was present when the Smolarkiewicz iterative upstream scheme (advection scheme = 2) was used. Thanks to Raffaele Bernardello for the update.
  • 17-Jun-2009: Update of subroutine ztosig. The interpolation was wrong at the border of local domains.
  • 10-Mar-2009: First realese of sbPOM test case for a seamount.
  •  

    Antoni JordiIMEDEA (UIB-CSIC)