COAWST sediment example results inquiry

[Honolulu_Wang]

Hi All,

I set to run the sediment transport example in COAWST/Projects/Sed_floc_toy/
It contains 20 layers in the sea-bed as in the *.in file.
But when I finished the run then looked the results e.g,
The output variable: bed_porosity seemed unchanged, it is the same as the input file which is 0.5.
So is there sth I needed to turn on? I used the default CPP flags: sed_floc_toy.h in the folder.
Additionally, when I turned on the SED_BIODIFF, the results of bed_biodiff is also zero, I checked the code, it seemed it is still under developing?

Thanks.

Code: Select all

** svn $Id: sed_toy.h 429 2009-12-20 17:30:26Z arango $
*******************************************************************************
** Copyright (c) 2002-2010 The ROMS/TOMS Group                               **
**   Licensed under a MIT/X style license                                    **
**   See License_ROMS.txt                                                    **
*******************************************************************************
**
** Options for One-Dimensional (vertical) Sediment Toy.
**
** Application flag:   SED_TOY
** Input scripts:      ocean_sed_toy.in
**                     sediment_floc_toy.in
*/

#define ROMS_MODEL
#define BODYFORCE
#undef  LOG_PROFILE
#define DJ_GRADPS
#undef  TS_U3HADVECTION
#undef  TS_C2VADVECTION
#define TS_MPDATA
#undef  SALINITY
#undef  SPLINES
#define OUT_DOUBLE
#define ANA_GRID
#undef  ANA_INITIAL
#define ANA_SMFLUX
#define ATM_PRESS
#undef  ANA_PAIR
#define SOLVE3D
#ifdef SOLVE3D
# undef  ANA_SEDIMENT
# define ANA_BPFLUX
# define ANA_BSFLUX
# define ANA_BTFLUX
# define ANA_SPFLUX
# define ANA_SRFLUX
# define ANA_SSFLUX
# define ANA_STFLUX
#endif
#undef   ANA_VMIX
#undef   ANA_WWAVE

#undef  UV_LOGDRAG
#undef  UV_LDRAG
#undef  UV_QDRAG


#undef  SG_BBL
#ifdef SG_BBL
# undef  SG_CALC_ZNOT
# undef  SG_LOGINT
#endif

#undef  MB_BBL
#ifdef MB_BBL
# undef  MB_CALC_ZNOT
# undef  MB_Z0BIO
# undef  MB_Z0BL
# undef  MB_Z0RIP
#endif

#define SSW_BBL
#ifdef SSW_BBL
# undef  SSW_CALC_UB
# define SSW_CALC_ZNOT
# define SSW_LOGINT
#endif

#define GLS_MIXING
#ifdef GLS_MIXING
# define KANTHA_CLAYSON
# define N2S2_HORAVG
# undef  CRAIG_BANNER
# undef  CHARNOK
# undef  ZOS_HSIG
# undef  TKE_WAVEDISS
#endif

#define SEDIMENT
#ifdef SEDIMENT
# define SUSPLOAD
# undef  BEDLOAD_SOULSBY
# undef  BEDLOAD_MPM
# define SED_DENS
# undef  COHESIVE_BED
# undef  BF_TCR
# undef  LINEAR_TCR
# undef  POWERLAW_TCR
# define MIXED_BED
# undef  SED_MORPH
# define SED_FLOCS
# undef  FLOC_TURB_DISS
# define FLOC_BBL_DISS
# define SED_DEFLOC
# define SED_TAU_CD_CONST
# undef  SED_TAU_CD_LIN
# undef  SED_BIODIFF
#endif

[ckharris]

Our group is running a modified version of sed_floc_toy also.

The main thing you should see with this test case is sediment flocculation and disaggregation in the water column. So the sediment size distributions in the water column should be changing with flocculation processes.

This test case is not very focused on the sediment bed. The sediment bed porosity would not depend on the flocculation processes. I am not sure if the biodiffusion is a part of this test case, but the version of COAWST that I downloaded recently does have the routine sed_biodiff.F in Nonlinear/Sediment .

Hope that helps - Courtney Harris