Download Sedtrans05

Transcript
21
Direct access to the computation routines from other programs
It is possible to access the single-point model Sedtrans05 from other program, either to
develop another user interface, or to integrate Sedtrans05 in larger 2D/3D
hydrodynamic/sedimentary models.
All the computations in Sedtrans05 is accessed through the subroutine SEDTRANS05 (in
the file "sedtrans-core.f"), with following two exceptions:
- the text of the bedform predictions is generated in the subroutine BEDFORM (in the
file "bedform.f"), and
- density and viscosity of water is computed from salinity and temperature in the
subroutine DENSVISC (in the file "cohesive.f").
Before the first call to SEDTRANS05, it is necessary to make one call to the subroutine
INICONST.
The input parameters of Sedtrans05 are transmitted as arguments to SEDTRANS05. The
main results are returned as arguments from SEDTRANS05, additional results are
returned in the common blocks OSCL, THRE, FRICF, TIMIN, PROF, PROFIL
(definition of common blocks in file "sed_param.h").
For explanations of the input/output argument of SEDTRANS05, see the beginning of the
SEDTRANS05 subroutine. For explanations of the additional results see the comments in
the file "sed_param.h".
NO input argument check is done in the subroutine SEDTRANS05
Specific for the non-cohesive methods (call to SEDTRANS05 with IOPT<>7):
+ If IOPT1.NE.7, then none of the input arguments for the Cohesive Sediment
Algorithm will be accessed; the arguments BEDCHA and CONC do not have to be
arrays in this case.
Specific for the Cohesive Sediment Algorithm (call to SEDTRANS05 with IOPT=7):
+ Various constants and parameters of the Cohesive Sediment Algorithm are specified in
the common block CCONST, which is initialized by the subroutine INICONST.
Expert users can change these constants/parameters by declaring the common block
CCONST (using "INCLUDE sed_param.h") in the calling program and modifying
some values AFTER the call to INICONST.
+ The subroutine BEDINI (in file "cohesive.f") is NOT tested so far.
+ The subroutine COMPACT (in file "cohesive.f") is NOT tested so far.
+ The array variable BEDCHA (describing the variations of the cohesive bed with
depth) must have at least 3 row and may have up to 50 rows. The function
CHECKBEDCHA can be used to check that information in BEDCHA is valid.