The TELEMAC-MASCARET system
trunk
|
Go to the source code of this file.
Functions/Subroutines | |
subroutine | suspension_main (SLVTRA, HN, HN_TEL, MU, TOB, FDM, FD90, KSP, KSR, KS, VOLU2D, V2DPAR, UNSV2D, AFBOR, BFBOR, ZF, LICBOR, IFAMAS, MASKEL, MASKPT, U2D, V2D, NSICLA, NPOIN, NPTFR, IELMT, OPTDIF, RESOL, LT, NIT, OPTBAN, OPTADV, OPDTRA, KENT, KSORT, KLOG, KNEU, KDIR, KDDL, DEBUG, DTS, CSF_SABLE, ZERO, GRAV, XKX, XKY, KARMAN, XMVE, XMVS, VCE, HMIN, XWC, VITCD, PARTHENIADES, BILMA, MSK, CHARR, IMP_INFLOW_C, MESH, ZF_S, CS, CST, CTILD, CBOR, DISP, IT1, IT2, IT3, IT4, TB, T1, T2, T3, T4, T8, T9, T10, T11, T12, T14, TE1, CLT, TE2, TE3, S, AM1_S, AM2_S, MBOR, ELAY, LIMDIF, MASKTR, TETA_SUSP, AC, MASED0, MASINI, MASTEN, MASTOU, ES, ES_SABLE, ES_VASE, AVAIL, ENTETS, PASS, ZFCL_S, HPROP, FLUDPT, FLUDP, FLUER, DISP_C, KX, KY, KZ, UCONV, VCONV, QSXS, QSYS, QSCLXS, QSCLYS, QSCL_S, QS_S, QS_C, CSTAEQ, CSRATIO, ICQ, MASTCP, MASFIN, MASDEPT, MASDEP, MASSOU, CORR_CONV, ZREF, SEDCO, VISC_TEL, CODE, DIFT, DM1, UCONV_TEL, VCONV_TEL, ZCONV, SOLSYS, FLBOR_TEL, FLBOR_SIS, FLBORTRA, NUMLIQ, NFRLIQ, MIXTE, NOMBLAY, CONC, TOCE_VASE, TOCE_SABLE, FLUER_VASE, TOCE_MIXTE, MS_SABLE, MS_VASE, DIRFLU, MAXADV) |
subroutine suspension_main | ( | type (slvcfg), intent(inout) | SLVTRA, |
type (bief_obj), intent(in) | HN, | ||
type (bief_obj), intent(in) | HN_TEL, | ||
type (bief_obj), intent(in) | MU, | ||
type (bief_obj), intent(in) | TOB, | ||
double precision, dimension(nsicla), intent(in) | FDM, | ||
double precision, dimension(nsicla), intent(in) | FD90, | ||
type (bief_obj), intent(in) | KSP, | ||
type (bief_obj), intent(in) | KSR, | ||
type (bief_obj), intent(in) | KS, | ||
type (bief_obj), intent(in) | VOLU2D, | ||
type (bief_obj), intent(in) | V2DPAR, | ||
type (bief_obj), intent(in) | UNSV2D, | ||
type (bief_obj), intent(in) | AFBOR, | ||
type (bief_obj), intent(in) | BFBOR, | ||
type (bief_obj), intent(in) | ZF, | ||
type (bief_obj), intent(in) | LICBOR, | ||
type (bief_obj), intent(in) | IFAMAS, | ||
type (bief_obj), intent(in) | MASKEL, | ||
type (bief_obj), intent(in) | MASKPT, | ||
type (bief_obj), intent(in) | U2D, | ||
type (bief_obj), intent(in) | V2D, | ||
integer, intent(in) | NSICLA, | ||
integer, intent(in) | NPOIN, | ||
integer, intent(in) | NPTFR, | ||
integer, intent(in) | IELMT, | ||
integer, intent(in) | OPTDIF, | ||
integer, intent(in) | RESOL, | ||
integer, intent(in) | LT, | ||
integer, intent(in) | NIT, | ||
integer, intent(in) | OPTBAN, | ||
integer, intent(in) | OPTADV, | ||
integer, intent(in) | OPDTRA, | ||
integer, intent(in) | KENT, | ||
integer, intent(in) | KSORT, | ||
integer, intent(in) | KLOG, | ||
integer, intent(in) | KNEU, | ||
integer, intent(in) | KDIR, | ||
integer, intent(in) | KDDL, | ||
integer, intent(in) | DEBUG, | ||
double precision, intent(in) | DTS, | ||
double precision, intent(in) | CSF_SABLE, | ||
double precision, intent(in) | ZERO, | ||
double precision, intent(in) | GRAV, | ||
double precision, intent(in) | XKX, | ||
double precision, intent(in) | XKY, | ||
double precision, intent(in) | KARMAN, | ||
double precision, intent(in) | XMVE, | ||
double precision, intent(in) | XMVS, | ||
double precision, intent(in) | VCE, | ||
double precision, intent(in) | HMIN, | ||
double precision, dimension(nsicla), intent(in) | XWC, | ||
double precision, intent(in) | VITCD, | ||
double precision, intent(in) | PARTHENIADES, | ||
logical, intent(in) | BILMA, | ||
logical, intent(in) | MSK, | ||
logical, intent(in) | CHARR, | ||
logical, intent(in) | IMP_INFLOW_C, | ||
type(bief_mesh), intent(inout) | MESH, | ||
type (bief_obj), intent(inout) | ZF_S, | ||
type (bief_obj), intent(inout) | CS, | ||
type (bief_obj), intent(inout) | CST, | ||
type (bief_obj), intent(inout) | CTILD, | ||
type (bief_obj), intent(inout) | CBOR, | ||
type (bief_obj), intent(inout) | DISP, | ||
type (bief_obj), intent(inout) | IT1, | ||
type (bief_obj), intent(inout) | IT2, | ||
type (bief_obj), intent(inout) | IT3, | ||
type (bief_obj), intent(inout) | IT4, | ||
type (bief_obj), intent(inout) | TB, | ||
type (bief_obj), intent(inout) | T1, | ||
type (bief_obj), intent(inout) | T2, | ||
type (bief_obj), intent(inout) | T3, | ||
type (bief_obj), intent(inout) | T4, | ||
type (bief_obj), intent(inout) | T8, | ||
type (bief_obj), intent(inout) | T9, | ||
type (bief_obj), intent(inout) | T10, | ||
type (bief_obj), intent(inout) | T11, | ||
type (bief_obj), intent(inout) | T12, | ||
type (bief_obj), intent(inout) | T14, | ||
type (bief_obj), intent(inout) | TE1, | ||
type (bief_obj), intent(inout) | CLT, | ||
type (bief_obj), intent(inout) | TE2, | ||
type (bief_obj), intent(inout) | TE3, | ||
type (bief_obj), intent(inout) | S, | ||
type (bief_obj), intent(inout) | AM1_S, | ||
type (bief_obj), intent(inout) | AM2_S, | ||
type (bief_obj), intent(inout) | MBOR, | ||
type (bief_obj), intent(inout) | ELAY, | ||
type (bief_obj), intent(inout) | LIMDIF, | ||
type (bief_obj), intent(inout) | MASKTR, | ||
double precision, intent(inout) | TETA_SUSP, | ||
double precision, dimension(nsicla), intent(inout) | AC, | ||
double precision, dimension(nsicla), intent(inout) | MASED0, | ||
double precision, dimension(nsicla), intent(inout) | MASINI, | ||
double precision, dimension(nsicla), intent(inout) | MASTEN, | ||
double precision, dimension(nsicla), intent(inout) | MASTOU, | ||
double precision, dimension(npoin,nomblay), intent(inout) | ES, | ||
double precision, dimension(npoin,nomblay), intent(inout) | ES_SABLE, | ||
double precision, dimension(npoin,nomblay), intent(inout) | ES_VASE, | ||
double precision, dimension(npoin,nomblay,nsicla), intent(inout) | AVAIL, | ||
logical, intent(inout) | ENTETS, | ||
logical, intent(inout) | PASS, | ||
type (bief_obj), intent(inout) | ZFCL_S, | ||
type (bief_obj), intent(inout) | HPROP, | ||
type (bief_obj), intent(inout) | FLUDPT, | ||
type (bief_obj), intent(inout) | FLUDP, | ||
type (bief_obj), intent(inout) | FLUER, | ||
type (bief_obj), intent(inout) | DISP_C, | ||
type (bief_obj), intent(inout) | KX, | ||
type (bief_obj), intent(inout) | KY, | ||
type (bief_obj), intent(inout) | KZ, | ||
type (bief_obj), intent(inout) | UCONV, | ||
type (bief_obj), intent(inout) | VCONV, | ||
type (bief_obj), intent(inout) | QSXS, | ||
type (bief_obj), intent(inout) | QSYS, | ||
type (bief_obj), intent(inout) | QSCLXS, | ||
type (bief_obj), intent(inout) | QSCLYS, | ||
type (bief_obj), intent(inout) | QSCL_S, | ||
type (bief_obj), intent(inout) | QS_S, | ||
type (bief_obj), intent(in) | QS_C, | ||
type (bief_obj), intent(inout) | CSTAEQ, | ||
type (bief_obj), intent(inout) | CSRATIO, | ||
integer, intent(in) | ICQ, | ||
double precision, dimension(nsicla), intent(out) | MASTCP, | ||
double precision, dimension(nsicla), intent(out) | MASFIN, | ||
double precision, dimension(nsicla), intent(out) | MASDEPT, | ||
double precision, dimension(nsicla), intent(out) | MASDEP, | ||
double precision, intent(out) | MASSOU, | ||
logical, intent(in) | CORR_CONV, | ||
type (bief_obj), intent(inout) | ZREF, | ||
logical, dimension(nsicla), intent(in) | SEDCO, | ||
type (bief_obj), intent(in) | VISC_TEL, | ||
character(len=24), intent(in) | CODE, | ||
logical, intent(in) | DIFT, | ||
type (bief_obj), intent(in) | DM1, | ||
type (bief_obj), intent(in) | UCONV_TEL, | ||
type (bief_obj), intent(in) | VCONV_TEL, | ||
type (bief_obj), intent(in) | ZCONV, | ||
integer, intent(in) | SOLSYS, | ||
type (bief_obj), intent(in) | FLBOR_TEL, | ||
type (bief_obj), intent(inout) | FLBOR_SIS, | ||
type (bief_obj), intent(inout) | FLBORTRA, | ||
integer, dimension(nptfr), intent(in) | NUMLIQ, | ||
integer, intent(in) | NFRLIQ, | ||
logical, intent(in) | MIXTE, | ||
integer, intent(in) | NOMBLAY, | ||
double precision, dimension(npoin,nomblay), intent(inout) | CONC, | ||
double precision, dimension(nomblay), intent(inout) | TOCE_VASE, | ||
double precision, intent(inout) | TOCE_SABLE, | ||
type (bief_obj), intent(inout) | FLUER_VASE, | ||
type (bief_obj), intent(inout) | TOCE_MIXTE, | ||
type (bief_obj), intent(inout) | MS_SABLE, | ||
type (bief_obj), intent(inout) | MS_VASE, | ||
integer, intent(in) | DIRFLU, | ||
integer, intent(in) | MAXADV | ||
) |
[in,out] | slvtra | [in, out] ENTET Logical, if yes information is given on mass conservation |
[in] | ENTETS | Logical, if yes information is given on mass conservation for suspension |
[in,out] | ES | Thickness of each layer (m) |
[in,out] | ES_VASE | Thickness of the mud layer (m) |
[in,out] | ES_SABLE | Thickness of the sand layer (m) |
[in] | FDM | Grain size per sediment class |
[in,out] | FLBORTRA | Fluxes at boundaries tracer |
[in,out] | FLBOR_SIS | Fluxes at boundaries sisyphe |
[in] | FLBOR_TEL | Fluxes at boundaries telemac |
[in,out] | FLUDP | Deposition flux (m/s) |
[in,out] | FLUDPT | Deposition flux (implicit) |
[in,out] | FLUER | Erosion flux (m/s) |
[in,out] | FLUER_VASE | Erosion flux (m/s)for mixed sediments |
[in] | GRAV | Acceleration of gravity |
[in] | HMIN | Minimum value of water depth (m) |
[in] | HN | Water depth (m) |
[in] | HN_TEL | Water depth sent by telemac or calling code |
[in,out] | HPROP | Propagation depth (done in cvdftr) |
[in] | ICQ | Flag for reference concentration formula |
[in] | IELMT | Number of elements |
[in] | IFAMAS | A modified ifabor when elements are masked |
[in] | IMP_INFLOW_C | Logical, imposed equilibrium concentration at the inflow or not |
[in,out] | IT1 | Integer work array in a bief_obj structure |
[in,out] | IT2 | Integer work array in a bief_obj structure |
[in,out] | IT3 | Integer work array in a bief_obj structure |
[in,out] | IT4 | Integer work array in a bief_obj structure |
[in] | KARMAN | Von karman constant |
[in] | KDDL | Convention for degree of freedom |
[in] | KDIR | Convention for dirichlet point |
[in] | KENT | Convention for liquid input with prescribed value |
[in] | KLOG | Convention for solid boundary |
[in] | KNEU | Convention for neumann condition |
[in] | KS | Total bed roughness |
[in] | KSORT | Convention for free output |
[in] | KSP | Skin bed roughness |
[in] | KSR | Ripple bed roughness |
[in,out] | KX | Coefficients of the dispersion tensor (dim. npoin) |
[in,out] | KY | Coefficients of the dispersion tensor (dim. npoin) |
[in,out] | KZ | Coefficients of the dispersion tensor (dim. npoin) |
[in] | LICBOR | Boundary conditions for sediment |
[in,out] | LIMDIF | Boundary conditions for diffusion |
[in] | LT | Iteration |
[out] | MASDEP | Total deposited mass |
[out] | MASDEPT | Deposited mass during the time step |
[in,out] | MASED0 | Suspended mass balance |
[out] | MASFIN | Mass at the end |
[in,out] | MASINI | Initial mass |
[in] | MASKEL | Masking of elements |
[in] | MASKPT | Masking per point |
[in,out] | MASKTR | Masking for tracers, per point |
[out] | MASSOU | Mass of tracer added by source term SEE DIFSOU |
[out] | MASTCP | ??? ne sert a rien, a supprimer |
[in,out] | MASTEN | Mass entered through liquid boundary |
[in,out] | MASTOU | Mass created by source term |
[in] | MAXADV | Maximum number of iterations for advection schemes |
[in,out] | MBOR | Matrix object |
[in,out] | MESH | Mesh structure |
[in] | MIXTE | Logical, mixte sediment or not |
[in] | MSK | If yes, there is masked elements |
[in,out] | MS_SABLE | Mass of sand per layer (kg/m2) |
[in,out] | MS_VASE | Mass of mud per layerv (kg/m2) |
[in] | MU | Correction factor for bed roughness |
[in] | NOMBLAY | Number of layers for consolidation |
[in] | NFRLIQ | Number of liquid boundaries |
[in] | NIT | Total number of iterations |
[in] | NPOIN | Number of points |
[in] | NPTFR | Number of boundary points |
[in] | NSICLA | Number of size classes for bed materials |
[in] | NUMLIQ | Liquid boundary number of boundary points |
[in] | OPDTRA | Option for the diffusion of tracers |
[in] | OPTBAN | Option for the treatment of tidal flats |
[in] | OPTDIF | Option for the dispersion |
[in] | OPTADV | Scheme option for advection |
[in] | PARTHENIADES | Constant of the krone and partheniades erosion law (m/s) |
[in,out] | PASS | In fact pass_susp in sisyphe.f, arrives as .true. AT FIRST CALL AND IS CHANGED INTO .FALSE. BELOW |
[in,out] | QSCLXS | Suspended load transport rate for each class x-direction |
[in,out] | QSCLYS | Suspended load transport rate for each class y-direction |
[in,out] | QSCL_S | Suspended load transport rate |
[in,out] | QSXS | Solid discharge x (suspension) |
[in,out] | QSYS | Solid discharge y (suspension) |
[in] | QS_C | Bedload transport rate |
[in,out] | QS_S | Suspended load transport rate |
[in] | RESOL | Choice of advection scheme |
[in,out] | S | Void structure |
[in] | SEDCO | Logical, sediment cohesive or not |
[in,out] | SLVTRA | Slvcfg structure |
[in] | SOLSYS | Slvcfg structure |
[in,out] | T1 | Work bief_obj structure |
[in,out] | T10 | Work bief_obj structure |
[in,out] | T11 | Work bief_obj structure |
[in,out] | T12 | Work bief_obj structure |
[in,out] | T2 | Work bief_obj structure |
[in,out] | T3 | Work bief_obj structure |
[in,out] | T4 | Work bief_obj structure |
[in,out] | T5 | Work bief_obj structure |
[in,out] | T6 | Work bief_obj structure |
[in,out] | T7 | Work bief_obj structure |
[in,out] | T8 | Work bief_obj structure |
[in,out] | T9 | Work bief_obj structure |
[in] | TB | Block of working arrays |
[in,out] | TE1 | Working array for elements |
[in,out] | TE2 | Working array for elements |
[in,out] | TE3 | Working array for elements |
[in,out] | TETA_SUSP | Implicitation factor for the deposition flux and diffusion |
[in] | TOB | Bed shear stress (total friction) |
[in,out] | TOCE_SABLE | Critical shear stress for sand (n/m2) |
[in,out] | TOCE_MIXTE | Critical shear stress for mixed sediments (n/m2) |
[in,out] | TOCE_VASE | Critical erosion shear stress of the mud per layer (n/m2) |
[in] | U2D | Mean flow velocity x-direction |
[in,out] | UCONV | X-component advection field (sisyphe) |
[in] | UCONV_TEL | X-component advection field (telemac) |
[in] | UNSV2D | Inverse of integrals of test functions |
[in] | V2D | Mean flow velocity y-direction |
[in] | V2DPAR | Integral of test functions, assembled in parallel |
[in] | VCE | Flow viscosity |
[in,out] | VCONV | Y-component advection field (sisyphe) |
[in] | VCONV_TEL | Y-component advection field (telemac) |
[in] | VISC_TEL | Velocity diffusivity (telemac) |
[in] | VITCD | Critical shear velocity for mud deposition |
[in] | VOLU2D | Integral of bases |
[in] | XKX | Coefficient used for computing the dispersion DEPENDS OF OPTIONS |
[in] | XKY | Coefficient used for computing the dispersion DEPENDS OF OPTIONS |
[in] | XMVE | Fluid density |
[in] | XMVS | Sediment density |
[in] | XWC | Settling velocities per class of sediment |
[in] | ZCONV | The piece-wise constant part of advection field IS DM1*GRAD(ZCONV) |
[in] | ZERO | Zero |
[in] | ZF | Elevation of bottom |
[in,out] | ZFCL_S | Bed evolution per class, due to suspended sediment |
[in,out] | ZF_S | Accumulated bed evolution due to suspended sediment |
[in] | ZREF | Reference elevation |
Definition at line 25 of file suspension_main.f.