Go to the source code of this file.
Functions/Subroutines | |
| subroutine | suspension_computation (SLVTRA, HN, HN_TEL, UCONV, VCONV, MU, TOB, FDM, FD90, KSP, KSR, KS, ELAY, AVA, AFBOR, BFBOR, LIMDIF, CLT, MASKEL, MASKTR, MASKPT, IFAMAS, NPOIN, IELM, NPTFR, ITRA, LT, NIT, RESOL, OPTBAN, KENT, KDDL, KDIR, KSORT, KLOG, KNEU, OPTADV, OPDTRA, DEBUG, CSF_SABLE, TETA_SUSP, DT, MASED0, ZERO, XWC, KARMAN, XMVE, XMVS, VCE, GRAV, HMIN, VITCD, PARTHENIADES, ENTETS, BILMA, MSK, CHARR, IMP_INFLOW_C, MESH, ZF, CS, CST, CTILD, CBOR, DISP, IT1, IT2, IT3, IT4, TB, T1, T2, T3, T4, T8, T9, T10, T11, T12, T14, TE1, TE2, TE3, S, AM1_S, AM2_S, MBOR, MASTEN, MASTOU, MASINI, AC, ZFCL_S, FLUDPT, FLUDP, FLUER, HPROP, DISP_C, CSTAEQ, CSRATIO, MASFIN, MASDEPT, MASDEP, MASSOU, QS_C, ICQ, ZREF, CORR_CONV, U2D, V2D, SEDCO, DIFT, DM1, ZCONV, UCONV_TEL, VCONV_TEL, SOLSYS, FLBOR_TEL, FLBOR_SIS, FLBORTRA, CODE, VOLU2D, V2DPAR, UNSV2D, NUMLIQ, NFRLIQ, LICBOR, MIXTE, AVAIL, NSICLA, ES, ES_SABLE, ES_VASE, NOMBLAY, CONC, TOCE_VASE, TOCE_SABLE, FLUER_VASE, TOCE_MIXTE, MS_SABLE, MS_VASE, DIRFLU, QSCLXS, QSCLYS, MAXADV) |
Function/Subroutine Documentation
◆ suspension_computation()
| subroutine suspension_computation | ( | type (slvcfg), intent(inout) | SLVTRA, |
| type (bief_obj), intent(in), target | HN, | ||
| type (bief_obj), intent(in), target | HN_TEL, | ||
| type (bief_obj), intent(inout) | UCONV, | ||
| type (bief_obj), intent(inout) | VCONV, | ||
| type (bief_obj), intent(in) | MU, | ||
| type (bief_obj), intent(in) | TOB, | ||
| double precision, intent(in) | FDM, | ||
| double precision, 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(inout) | ELAY, | ||
| double precision, dimension(npoin), intent(in) | AVA, | ||
| type (bief_obj), intent(in) | AFBOR, | ||
| type (bief_obj), intent(in) | BFBOR, | ||
| type (bief_obj), intent(inout) | LIMDIF, | ||
| type (bief_obj), intent(inout) | CLT, | ||
| type (bief_obj), intent(in) | MASKEL, | ||
| type (bief_obj), intent(inout) | MASKTR, | ||
| type (bief_obj), intent(in) | MASKPT, | ||
| type (bief_obj), intent(in) | IFAMAS, | ||
| integer, intent(in) | NPOIN, | ||
| integer, intent(in) | IELM, | ||
| integer, intent(in) | NPTFR, | ||
| integer, intent(in) | ITRA, | ||
| integer, intent(in) | LT, | ||
| integer, intent(in) | NIT, | ||
| integer, intent(in) | RESOL, | ||
| integer, intent(in) | OPTBAN, | ||
| integer, intent(in) | KENT, | ||
| integer, intent(in) | KDDL, | ||
| integer, intent(in) | KDIR, | ||
| integer, intent(in) | KSORT, | ||
| integer, intent(in) | KLOG, | ||
| integer, intent(in) | KNEU, | ||
| integer, intent(in) | OPTADV, | ||
| integer, intent(in) | OPDTRA, | ||
| integer, intent(in) | DEBUG, | ||
| double precision, intent(in) | CSF_SABLE, | ||
| double precision, intent(in) | TETA_SUSP, | ||
| double precision, intent(in) | DT, | ||
| double precision, intent(in) | MASED0, | ||
| double precision, intent(in) | ZERO, | ||
| double precision, intent(in) | XWC, | ||
| double precision, intent(in) | KARMAN, | ||
| double precision, intent(in) | XMVE, | ||
| double precision, intent(in) | XMVS, | ||
| double precision, intent(in) | VCE, | ||
| double precision, intent(in) | GRAV, | ||
| double precision, intent(in) | HMIN, | ||
| double precision, intent(in) | VITCD, | ||
| double precision, intent(in) | PARTHENIADES, | ||
| logical, intent(in) | ENTETS, | ||
| 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(in) | ZF, | ||
| 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), target | 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) | 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, | ||
| double precision, intent(inout) | MASTEN, | ||
| double precision, intent(inout) | MASTOU, | ||
| double precision, intent(inout) | MASINI, | ||
| double precision, intent(inout) | AC, | ||
| type (bief_obj), intent(inout) | ZFCL_S, | ||
| type (bief_obj), intent(inout) | FLUDPT, | ||
| type (bief_obj), intent(inout) | FLUDP, | ||
| type (bief_obj), intent(inout) | FLUER, | ||
| type (bief_obj), intent(inout) | HPROP, | ||
| type (bief_obj), intent(inout) | DISP_C, | ||
| type (bief_obj), intent(inout) | CSTAEQ, | ||
| type (bief_obj), intent(inout) | CSRATIO, | ||
| double precision, intent(inout) | MASFIN, | ||
| double precision, intent(inout) | MASDEPT, | ||
| double precision, intent(inout) | MASDEP, | ||
| double precision, intent(inout) | MASSOU, | ||
| type (bief_obj), intent(in) | QS_C, | ||
| integer, intent(in) | ICQ, | ||
| type (bief_obj), intent(inout) | ZREF, | ||
| logical, intent(in) | CORR_CONV, | ||
| type (bief_obj), intent(in) | U2D, | ||
| type (bief_obj), intent(in) | V2D, | ||
| logical, intent(in) | SEDCO, | ||
| logical, intent(in) | DIFT, | ||
| type (bief_obj), intent(in) | DM1, | ||
| type (bief_obj), intent(in) | ZCONV, | ||
| type (bief_obj), intent(in) | UCONV_TEL, | ||
| type (bief_obj), intent(in) | VCONV_TEL, | ||
| integer, intent(in) | SOLSYS, | ||
| type (bief_obj), intent(in) | FLBOR_TEL, | ||
| type (bief_obj), intent(inout) | FLBOR_SIS, | ||
| type (bief_obj), intent(inout) | FLBORTRA, | ||
| character(len=24), intent(in) | CODE, | ||
| type (bief_obj), intent(in) | VOLU2D, | ||
| type (bief_obj), intent(in) | V2DPAR, | ||
| type (bief_obj), intent(in) | UNSV2D, | ||
| integer, dimension(*), intent(in) | NUMLIQ, | ||
| integer, intent(in) | NFRLIQ, | ||
| type (bief_obj), intent(in) | LICBOR, | ||
| logical, intent(in) | MIXTE, | ||
| double precision, dimension(npoin,nomblay,nsicla), intent(inout) | AVAIL, | ||
| integer, intent(in) | NSICLA, | ||
| double precision, dimension(npoin,nomblay), intent(inout) | ES, | ||
| double precision, dimension(*), intent(inout) | ES_SABLE, | ||
| double precision, dimension(*), intent(inout) | ES_VASE, | ||
| 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, | ||
| double precision, dimension(*), intent(inout) | MS_SABLE, | ||
| double precision, dimension(*), intent(inout) | MS_VASE, | ||
| integer, intent(in) | DIRFLU, | ||
| type (bief_obj), intent(inout) | QSCLXS, | ||
| type (bief_obj), intent(inout) | QSCLYS, | ||
| integer, intent(in) | MAXADV | ||
| ) |
- Parameters
-
[in,out] slvtra [in, out] AC Critical shields parameter [in] ACLADM Mean diameter of sediment [in] AFBOR Boundary condition on f: nu*df/dn=afbor*f+bfbor [in,out] AM1_S Matrix object [in,out] AM2_S Matrix object [in,out] AVAIL Volume percent of each class [in] BFBOR Boundary condition on f: nu*df/dn=afbor*f+bfbor [in] BILMA Mass balance [in,out] CBOR Imposed suspended sand concentration at the boundary [in] CF Quadratic friction coefficient [in] CHARR Logical, bedload or not [in,out] CLT Boundary conditions for tracer (modified litbor) [in] CODE Hydrodynamic code in case of coupling [in,out] CONC_VASE Mud concentration for each layer [in] CORR_CONV Correction on convection velocity [in,out] CS Concentration at time n [in] CSF_SABLE Volume concentration of the sand bed [in,out] CSRATIO Equilibrium concentration for soulsby-van rijn eq. [in,out] CST Concentration at time t(n+1) [in,out] CSTAEQ Equilibrium concentration [in,out] CTILD Concentration after advection [in] DEBUG Flag for debugging [in] DIFT Diffusion of suspended sediment concentration [in] DISP Viscosity coefficients along x,y and z . IF P0 : PER ELEMENT IF P1 : PERR POINT [in,out] DISP_S Work array for saving dispc [in] DM1 The piece-wise constant part of advection field IS DM1*GRAD(ZCONV) [in] DTS Time step for suspension [in,out] ELAY Thickness of each layer [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 Layer thicknesses as double precision [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 [in,out] FLUDPT Deposition flux (implicit) [in,out] FLUER Erosion flux [in,out] FLUER_VASE For mixed sediments [in] GRAV Acceleration of gravity [in] HMIN Minimum value of water depth [in] HN Water depth [in] HN_TEL Water depth as sent by telemac or calling code [in,out] HPROP Propagation depth (done in cvdftr) [in] ICQ Reference concentration formula [in] IELM Type of element [in] IFAMAS A modified ifabor when elements are masked [in] IMP_INFLOW_C Imposed concentration in inflow [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] KINC Convention for incident wave boundary condition [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,out] MAXADV Maximum number of iterations of advection schemes [in,out] MBOR Matrix object [in,out] MESH Mesh structure [in] MIXTE Mixture of cohesive and non cohesive sediment [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 layer (kg/m2) [in] MU Correction factor for bed roughness [in] NOMBLAY Number of layers [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 Transport rate for each class x-direction [in,out] QSCLYS 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 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] TASS Consolidation taken into account [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 bed shear stress of sand [in,out] TOCE_MIXTE Critical bed shear stress of the mixed seduiment per layer [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,out] W1 Working array [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 [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,out] ZREF Reference elevation
Definition at line 25 of file suspension_computation.f.
Here is the call graph for this function: Here is the caller graph for this function: