/*
 * ** svn $Id: icetest.h 1307 2008-01-10 00:22:36Z bgalton $
 * *******************************************************************************
 * ** Copyright (c) 2002-2008 The ROMS/TOMS Group                               **
 * **   Licensed under a MIT/X style license                                    **
 * **   See License_ROMS.txt                                                    **
 * *******************************************************************************
 * **
 * ** Options for Weddell Idealized simulations
 * ** Solid walls in north, open boundries east and west
 * ** Input scripts determine resolution, timestep, 
 * ** grid, forcing and initialization
 * **
 * ** See ROMS/Include/cppdefs.h for documentation of flags (don't trust ROMSwiki)
* */


/* Define basic model configuration */
#define WEDDKYST
#define SOLVE3D
#define NONLIN_EOS
#define MASKING
#define CURVGRID

/* Iceshelf */
#define ICESHELF
#define ICESHELF_3EQN_VBC
#undef LIMIT_BSTRESS /* kaitlyn uses this *//* limit bottom stress */
#undef LIMIT_ICESTRESS /* kaitlyn uses this *//* limit ice shelf basal stress */

/* Momentum equations */
#define UV_ADV
#define UV_COR
#define UV_QDRAG

#undef  UV_LOGDRAG
#undef  UV_LDRAG
#undef  UV_PSOURCE
#undef  Q_PSOURCE
#undef  UV_VIS4 /* kaitlyn uses this */

#undef RI_SPLINES /* kaitlyn uses this...but I don't know what it is! */

/* TIDES */

#define WS_TIDES

#ifdef WS_TIDES
# define SSH_TIDES
# define UV_TIDES
# define ADD_FSOBC
# define ADD_M2OBC
# define RAMP_TIDES
#endif

/* Horizontal mixing or momentum */
#undef  VISC_GRID /* kaitlyn uses this */
#define MIX_S_UV
#undef MIX_GEO_UV/* kaitlyn uses this */

/* Tracer equations */
#define SALINITY
#define TS_U3HADVECTION
#define TS_C4VADVECTION /* This is the default, 4th-order Akima vertical advection */
#undef TS_DIF2 /* kaitlyn uses this */
#undef SSFLUX_EXTRA /* kaitlyn uses this.*/
#undef COOL_SKIN/* kaitlyn uses this */
#undef LONGWAVE /* kaitlyn uses this */
#undef LIMIT_STFLX_COOLING
#define UPWIND_LIMITER
/* ****  Look into SSFLUX_EXTRA ******* */
 
/* 3D-advection of tracers:
use default: 3rd order upstream horizontal,
centered 4th order vertical,
add no explicit diffusion*/

#define MIX_S_TS
#undef  MIX_GEO_TS
#undef  MIX_ISO_TS /* kaitlyn uses this */

#undef  TS_U3ADV_SPLIT
#undef  TS_A4HADVECTION/* kaitlyn uses this */
#undef  TS_C2HADVECTION
#undef  TS_MPDATA
#undef  TS_A4VADVECTION /* kaitlyn uses this */
#undef  TS_C2VADVECTION
#undef  TS_C4VADVECTION
#undef  TS_SVADVECTION
#undef  TS_DIF4
#undef  TS_SMAGORINSKY
#undef  TS_FIXED
#undef  T_PASSIVE

#undef  DIFF_GRID /* kaitlyn uses this */

/* pressure gradient */
#define DJ_GRADPS

/* Vertical mixing parameterization:
LMD-mixing may be modified by activation of ICETEST flag */

#define LMD_MIXING /* kaitlyn uses this */

#undef  BVF_MIXING
#undef  GLS_MIXING
#undef  MY25_MIXING
#undef  ANA_VMIX

#ifdef MY25_MIXING
#define N2S2_HORAVG
#define K_C4ADVECTION
#endif

#ifdef  LMD_MIXING
#define LMD_CONVEC /* kaitlyn uses this */
#define LMD_DDMIX /* kaitlyn uses this */
#define LMD_RIMIX/* kaitlyn uses this */
#define LMD_SKPP/* kaitlyn uses this */
#define LMD_BKPP/* kaitlyn uses this */
#undef LMD_NONLOCAL /* kaitlyn uses this */
#endif


/* Define output */
#define PERFECT_RESTART
#define  AVERAGES
#undef  AVERAGES_QUADRATIC
#undef  STATIONS

#undef  DIAGNOSTICS_UV
#undef  DIAGNOSTICS_TS

/*sponge, climatology and nudging is now set in ocean.in file */

/* surface heat and salt from file*/
#define QCORRECTION
#define SRELAXATION
#define TRELAXATION_WK
#define WEDDKYST1DAY
#undef SCORRECTION_MASK
/* boundary conditions are defined in ocean.in*/

#define ANA_STFLUX
#define ANA_SSFLUX
#define ANA_BSFLUX/* kaitlyn uses this */
#define ANA_BTFLUX/* kaitlyn uses this */
#define ANA_SRFLUX/* kaitlyn uses this */

#undef  ANA_GRID
#undef  ANA_SEAICE
#undef  ANA_INITIAL /* to read initial data from clm file when perfect restart is applied
#undef  ANA_SMFLUX/* need to turn that on to skip forcing files */
#undef ANA_FSOBC
#undef ANA_TOBC
#undef ANA_M3OBC
#undef ANA_M2OBC

#undef ANA_SPONGE /* kaitlyn uses this: sponge layer at open boundary */

/* Internal flags that are set automatically */
/*obc inputs*/
#undef SPONGE