bset.F

subroutine bset

Description

Initial temperature and salinity fields are set and the initial buoyancy is calculated (through a call to subroutine bcalc).

Firstly, uniform temperature and salinity fields are assigned using the values of tmp_init and sal_init in parameters.dat. If an analytical initial temperature (or buoyancy) condition is required (analytic_init=.true.) then it is calculated from the formula in tmpfield.F. Initial temperature and/or salinity fields may be read in from a file (read_initial_ts=.true.). The format of the file, ts_form, is set in parameters.dat and is either Fortran format, eg '(10f8.2)', or 'unf' for a Fortran unformatted file. The input arrays are of size and shape (l,m,n). A temperature array is expected unless no_tmp is true and a salinity array is expected unless no_sal is true.

The initial temperature and salinity fields are multiplied by the mask of scalar points ipexb so that values at redundant points are set equal to zero.

Subroutine Arguments

none

Local variables

i, j, k local indices.

zd depth below surface (used in tmpfield.F).

temp3d temporary array with dimensions (l,m,n) used in reading in initial temperature and/or salinity fields.

work3d temporary array with dimensions (n,l,m) - transpose of temp3d.

Global variables changed

tmp, sal

Logical units

16 Formatted or unformatted file containing initial temperature and/or salinity arrays.

Order of Things

1. Set initial values of temperature and salinity.
2. Calculate initial buoyancy (through a call to bcalc).

Calls

dist3d

gsynchwrk

bcalc

Called By

b3dinit

Options - Logical

read_initial_ts - initial temperature and/or salinity fields are to be read in from a file (unit number 16)

analytic_init - an analytical initial condition is calculated.

no_tmp - temperature is constant in time and space - initial temperature data is not read from file.

no_sal - salinity is constant in time and space - initial salinity data is not read from file.

Options - Compiler

none

Known Issues

none