The TELEMAC-MASCARET system  trunk
Functions/Subroutines
suspension_main.f File Reference

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)
 

Function/Subroutine Documentation

◆ suspension_main()

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 
)
Parameters
[in,out]slvtra[in, out] ENTET Logical, if yes information is given on mass conservation
[in]ENTETSLogical, if yes information is given on mass conservation for suspension
[in,out]ESThickness of each layer (m)
[in,out]ES_VASEThickness of the mud layer (m)
[in,out]ES_SABLEThickness of the sand layer (m)
[in]FDMGrain size per sediment class
[in,out]FLBORTRAFluxes at boundaries tracer
[in,out]FLBOR_SISFluxes at boundaries sisyphe
[in]FLBOR_TELFluxes at boundaries telemac
[in,out]FLUDPDeposition flux (m/s)
[in,out]FLUDPTDeposition flux (implicit)
[in,out]FLUERErosion flux (m/s)
[in,out]FLUER_VASEErosion flux (m/s)for mixed sediments
[in]GRAVAcceleration of gravity
[in]HMINMinimum value of water depth (m)
[in]HNWater depth (m)
[in]HN_TELWater depth sent by telemac or calling code
[in,out]HPROPPropagation depth (done in cvdftr)
[in]ICQFlag for reference concentration formula
[in]IELMTNumber of elements
[in]IFAMASA modified ifabor when elements are masked
[in]IMP_INFLOW_CLogical, imposed equilibrium concentration at the inflow or not
[in,out]IT1Integer work array in a bief_obj structure
[in,out]IT2Integer work array in a bief_obj structure
[in,out]IT3Integer work array in a bief_obj structure
[in,out]IT4Integer work array in a bief_obj structure
[in]KARMANVon karman constant
[in]KDDLConvention for degree of freedom
[in]KDIRConvention for dirichlet point
[in]KENTConvention for liquid input with prescribed value
[in]KLOGConvention for solid boundary
[in]KNEUConvention for neumann condition
[in]KSTotal bed roughness
[in]KSORTConvention for free output
[in]KSPSkin bed roughness
[in]KSRRipple bed roughness
[in,out]KXCoefficients of the dispersion tensor (dim. npoin)
[in,out]KYCoefficients of the dispersion tensor (dim. npoin)
[in,out]KZCoefficients of the dispersion tensor (dim. npoin)
[in]LICBORBoundary conditions for sediment
[in,out]LIMDIFBoundary conditions for diffusion
[in]LTIteration
[out]MASDEPTotal deposited mass
[out]MASDEPTDeposited mass during the time step
[in,out]MASED0Suspended mass balance
[out]MASFINMass at the end
[in,out]MASINIInitial mass
[in]MASKELMasking of elements
[in]MASKPTMasking per point
[in,out]MASKTRMasking for tracers, per point
[out]MASSOUMass of tracer added by source term SEE DIFSOU
[out]MASTCP??? ne sert a rien, a supprimer
[in,out]MASTENMass entered through liquid boundary
[in,out]MASTOUMass created by source term
[in]MAXADVMaximum number of iterations for advection schemes
[in,out]MBORMatrix object
[in,out]MESHMesh structure
[in]MIXTELogical, mixte sediment or not
[in]MSKIf yes, there is masked elements
[in,out]MS_SABLEMass of sand per layer (kg/m2)
[in,out]MS_VASEMass of mud per layerv (kg/m2)
[in]MUCorrection factor for bed roughness
[in]NOMBLAYNumber of layers for consolidation
[in]NFRLIQNumber of liquid boundaries
[in]NITTotal number of iterations
[in]NPOINNumber of points
[in]NPTFRNumber of boundary points
[in]NSICLANumber of size classes for bed materials
[in]NUMLIQLiquid boundary number of boundary points
[in]OPDTRAOption for the diffusion of tracers
[in]OPTBANOption for the treatment of tidal flats
[in]OPTDIFOption for the dispersion
[in]OPTADVScheme option for advection
[in]PARTHENIADESConstant of the krone and partheniades erosion law (m/s)
[in,out]PASSIn fact pass_susp in sisyphe.f, arrives as .true. AT FIRST CALL AND IS CHANGED INTO .FALSE. BELOW
[in,out]QSCLXSSuspended load transport rate for each class x-direction
[in,out]QSCLYSSuspended load transport rate for each class y-direction
[in,out]QSCL_SSuspended load transport rate
[in,out]QSXSSolid discharge x (suspension)
[in,out]QSYSSolid discharge y (suspension)
[in]QS_CBedload transport rate
[in,out]QS_SSuspended load transport rate
[in]RESOLChoice of advection scheme
[in,out]SVoid structure
[in]SEDCOLogical, sediment cohesive or not
[in,out]SLVTRASlvcfg structure
[in]SOLSYSSlvcfg structure
[in,out]T1Work bief_obj structure
[in,out]T10Work bief_obj structure
[in,out]T11Work bief_obj structure
[in,out]T12Work bief_obj structure
[in,out]T2Work bief_obj structure
[in,out]T3Work bief_obj structure
[in,out]T4Work bief_obj structure
[in,out]T5Work bief_obj structure
[in,out]T6Work bief_obj structure
[in,out]T7Work bief_obj structure
[in,out]T8Work bief_obj structure
[in,out]T9Work bief_obj structure
[in]TBBlock of working arrays
[in,out]TE1Working array for elements
[in,out]TE2Working array for elements
[in,out]TE3Working array for elements
[in,out]TETA_SUSPImplicitation factor for the deposition flux and diffusion
[in]TOBBed shear stress (total friction)
[in,out]TOCE_SABLECritical shear stress for sand (n/m2)
[in,out]TOCE_MIXTECritical shear stress for mixed sediments (n/m2)
[in,out]TOCE_VASECritical erosion shear stress of the mud per layer (n/m2)
[in]U2DMean flow velocity x-direction
[in,out]UCONVX-component advection field (sisyphe)
[in]UCONV_TELX-component advection field (telemac)
[in]UNSV2DInverse of integrals of test functions
[in]V2DMean flow velocity y-direction
[in]V2DPARIntegral of test functions, assembled in parallel
[in]VCEFlow viscosity
[in,out]VCONVY-component advection field (sisyphe)
[in]VCONV_TELY-component advection field (telemac)
[in]VISC_TELVelocity diffusivity (telemac)
[in]VITCDCritical shear velocity for mud deposition
[in]VOLU2DIntegral of bases
[in]XKXCoefficient used for computing the dispersion DEPENDS OF OPTIONS
[in]XKYCoefficient used for computing the dispersion DEPENDS OF OPTIONS
[in]XMVEFluid density
[in]XMVSSediment density
[in]XWCSettling velocities per class of sediment
[in]ZCONVThe piece-wise constant part of advection field IS DM1*GRAD(ZCONV)
[in]ZEROZero
[in]ZFElevation of bottom
[in,out]ZFCL_SBed evolution per class, due to suspended sediment
[in,out]ZF_SAccumulated bed evolution due to suspended sediment
[in]ZREFReference elevation

Definition at line 25 of file suspension_main.f.

+ Here is the call graph for this function:
+ Here is the caller graph for this function: