seaice_itd/code/SEAICE_OPTIONS.h

C $Header: /u/gcmpack/MITgcm/verification/offline_exf_seaice/code/SEAICE_OPTIONS.h,v 1.24 2013/02/13 08:56:14 mlosch Exp $
C $Name:  $

C     *==========================================================*
C     | SEAICE_OPTIONS.h
C     | o CPP options file for sea ice package.
C     *==========================================================*
C     | Use this file for selecting options within the sea ice
C     | package.
C     *==========================================================*

#ifndef SEAICE_OPTIONS_H
#define SEAICE_OPTIONS_H
#include "PACKAGES_CONFIG.h"
#include "CPP_OPTIONS.h"

#ifdef ALLOW_SEAICE
C     Package-specific Options & Macros go here

C--   Write "text-plots" of certain fields in STDOUT for debugging.
#undef SEAICE_DEBUG

C--   Allow sea-ice dynamic code.
C     This option is provided to allow use of TAMC
C     on the thermodynamics component of the code only.
C     Sea-ice dynamics can also be turned off at runtime
C     using variable SEAICEuseDYNAMICS.
#define SEAICE_ALLOW_DYNAMICS

C--   By default, the sea-ice package uses its own integrated bulk
C     formulae to compute fluxes (fu, fv, EmPmR, Qnet, and Qsw) over
C     open-ocean.  When this flag is set, these variables are computed
C     in a separate external package, for example, pkg/exf, and then
C     modified for sea-ice effects by pkg/seaice.
#define SEAICE_EXTERNAL_FLUXES

C--   This CPP flag has been retired.  The number of ice categories
C     used to solve for seaice flux is now specified by run-time
C     parameter SEAICE_multDim.
C     Note: be aware of pickup_seaice.* compatibility issues when
C     restarting a simulation with a different number of categories.
c#define SEAICE_MULTICATEGORY

C--   run with sea Ice Thickness Distribution (ITD);
C     set number of categories (nITD) in SEAICE_SIZE.h
#define SEAICE_ITD

C--   Since the missing sublimation term is now included
C     this flag is needed for backward compatibility
#undef SEAICE_DISABLE_SUBLIM

C--   Suspected missing term in coupled ocn-ice heat budget (to be confirmed)
#undef SEAICE_DISABLE_HEATCONSFIX

C--   Default is constant seaice salinity (SEAICE_salt0); Define the following
C     flag to consider (space & time) variable salinity: advected and forming
C     seaice with a fraction (=SEAICE_saltFrac) of freezing seawater salinity.
C- Note: SItracer also offers an alternative way to handle variable salinity.
#undef SEAICE_VARIABLE_SALINITY

C--   Tracers of ice and/or ice cover.
#undef ALLOW_SITRACER
#ifdef ALLOW_SITRACER
C--   To try avoid 'spontaneous generation' of tracer maxima by advdiff.
# define ALLOW_SITRACER_ADVCAP
#endif

C--   By default the seaice model is discretized on a B-Grid (for
C     historical reasons). Define the following flag to use a new
C     (not thoroughly) test version on a C-grid
#define SEAICE_CGRID

C--   Only for the C-grid version it is possible to
#ifdef SEAICE_CGRID
C     enable JFNK code by defining the following flag
# define SEAICE_ALLOW_JFNK
C     enable this flag to reproduce old verification results for JFNK
# define SEAICE_PRECOND_EXTRA_EXCHANGE
C     enable LSR to use global (multi-tile) tri-diagonal solver
# define SEAICE_GLOBAL_3DIAG_SOLVER
C     enable EVP code by defining the following flag
# define SEAICE_ALLOW_EVP
# ifdef SEAICE_ALLOW_EVP
C--   When set use SEAICE_zetaMin and SEAICE_evpDampC to limit viscosities
C     from below and above in seaice_evp: not necessary, and not recommended
#  undef SEAICE_ALLOW_CLIPZETA
# endif /* SEAICE_ALLOW_EVP */
C     regularize zeta to zmax with a smooth tanh-function instead
C     of a min(zeta,zmax). This improves convergence of iterative
C     solvers (Lemieux and Tremblay 2009, JGR). No effect on EVP
# define SEAICE_ZETA_SMOOTHREG
C     allow the truncated ellipse rheology (runtime flag SEAICEuseTEM)
# undef SEAICE_ALLOW_TEM
#else /* not SEAICE_CGRID, but old B-grid */
C--   By default for B-grid dynamics solver wind stress under sea-ice is
C     set to the same value as it would be if there was no sea-ice.
C     Define following CPP flag for B-grid ice-ocean stress coupling.
# define SEAICE_BICE_STRESS

C--   By default for B-grid dynamics solver surface tilt is obtained
C     indirectly via geostrophic velocities. Define following CPP
C     in order to use ETAN instead.
# define EXPLICIT_SSH_SLOPE
C--   Defining this flag turns on FV-discretization of the B-grid LSOR solver.
C     It is smoother and includes all metric terms, similar to C-grid solvers.
C     It is here for completeness, but its usefulness is unclear.
# undef SEAICE_LSRBNEW
#endif /* SEAICE_CGRID */

C--   When set limit the Ice-Loading to mass of 1/5 of Surface ocean grid-box
#undef SEAICE_CAP_ICELOAD
C--   When set use SEAICE_clipVelocties = .true., to clip U/VICE at 40cm/s,
C     not recommended
#undef SEAICE_ALLOW_CLIPVELS
C--   When set cap the sublimation latent heat flux in solve4temp according
C     to the available amount of ice+snow. Otherwise this term is treated
C     like all of the others -- residuals heat and fw stocks are passed to
C     the ocean at the end of seaice_growth in a conservative manner.
C     SEAICE_CAP_SUBLIM is not needed as of now, but kept just in case.
#undef SEAICE_CAP_SUBLIM

C--   Enable free drift code
#define SEAICE_ALLOW_FREEDRIFT

#endif /* ALLOW_SEAICE */
#endif /* SEAICE_OPTIONS_H */

CEH3 ;;; Local Variables: ***
CEH3 ;;; mode:fortran ***
CEH3 ;;; End: ***
1	C $Header: /u/gcmpack/MITgcm/verification/offline_exf_seaice/code/SEAICE_OPTIONS.h,v 1.24 2013/02/13 08:56:14 mlosch Exp $
2	C $Name: $
3
4	C ==========================================================
5	C \| SEAICE_OPTIONS.h
6	C \| o CPP options file for sea ice package.
7	C ==========================================================
8	C \| Use this file for selecting options within the sea ice
9	C \| package.
10	C ==========================================================
11
12	#ifndef SEAICE_OPTIONS_H
13	#define SEAICE_OPTIONS_H
14	#include "PACKAGES_CONFIG.h"
15	#include "CPP_OPTIONS.h"
16
17	#ifdef ALLOW_SEAICE
18	C Package-specific Options & Macros go here
19
20	C-- Write "text-plots" of certain fields in STDOUT for debugging.
21	#undef SEAICE_DEBUG
22
23	C-- Allow sea-ice dynamic code.
24	C This option is provided to allow use of TAMC
25	C on the thermodynamics component of the code only.
26	C Sea-ice dynamics can also be turned off at runtime
27	C using variable SEAICEuseDYNAMICS.
28	#define SEAICE_ALLOW_DYNAMICS
29
30	C-- By default, the sea-ice package uses its own integrated bulk
31	C formulae to compute fluxes (fu, fv, EmPmR, Qnet, and Qsw) over
32	C open-ocean. When this flag is set, these variables are computed
33	C in a separate external package, for example, pkg/exf, and then
34	C modified for sea-ice effects by pkg/seaice.
35	#define SEAICE_EXTERNAL_FLUXES
36
37	C-- This CPP flag has been retired. The number of ice categories
38	C used to solve for seaice flux is now specified by run-time
39	C parameter SEAICE_multDim.
40	C Note: be aware of pickup_seaice.* compatibility issues when
41	C restarting a simulation with a different number of categories.
42	c#define SEAICE_MULTICATEGORY
43
44	C-- run with sea Ice Thickness Distribution (ITD);
45	C set number of categories (nITD) in SEAICE_SIZE.h
46	#define SEAICE_ITD
47
48	C-- Since the missing sublimation term is now included
49	C this flag is needed for backward compatibility
50	#undef SEAICE_DISABLE_SUBLIM
51
52	C-- Suspected missing term in coupled ocn-ice heat budget (to be confirmed)
53	#undef SEAICE_DISABLE_HEATCONSFIX
54
55	C-- Default is constant seaice salinity (SEAICE_salt0); Define the following
56	C flag to consider (space & time) variable salinity: advected and forming
57	C seaice with a fraction (=SEAICE_saltFrac) of freezing seawater salinity.
58	C- Note: SItracer also offers an alternative way to handle variable salinity.
59	#undef SEAICE_VARIABLE_SALINITY
60
61	C-- Tracers of ice and/or ice cover.
62	#undef ALLOW_SITRACER
63	#ifdef ALLOW_SITRACER
64	C-- To try avoid 'spontaneous generation' of tracer maxima by advdiff.
65	# define ALLOW_SITRACER_ADVCAP
66	#endif
67
68	C-- By default the seaice model is discretized on a B-Grid (for
69	C historical reasons). Define the following flag to use a new
70	C (not thoroughly) test version on a C-grid
71	#define SEAICE_CGRID
72
73	C-- Only for the C-grid version it is possible to
74	#ifdef SEAICE_CGRID
75	C enable JFNK code by defining the following flag
76	# define SEAICE_ALLOW_JFNK
77	C enable this flag to reproduce old verification results for JFNK
78	# define SEAICE_PRECOND_EXTRA_EXCHANGE
79	C enable LSR to use global (multi-tile) tri-diagonal solver
80	# define SEAICE_GLOBAL_3DIAG_SOLVER
81	C enable EVP code by defining the following flag
82	# define SEAICE_ALLOW_EVP
83	# ifdef SEAICE_ALLOW_EVP
84	C-- When set use SEAICE_zetaMin and SEAICE_evpDampC to limit viscosities
85	C from below and above in seaice_evp: not necessary, and not recommended
86	# undef SEAICE_ALLOW_CLIPZETA
87	# endif /* SEAICE_ALLOW_EVP */
88	C regularize zeta to zmax with a smooth tanh-function instead
89	C of a min(zeta,zmax). This improves convergence of iterative
90	C solvers (Lemieux and Tremblay 2009, JGR). No effect on EVP
91	# define SEAICE_ZETA_SMOOTHREG
92	C allow the truncated ellipse rheology (runtime flag SEAICEuseTEM)
93	# undef SEAICE_ALLOW_TEM
94	#else /* not SEAICE_CGRID, but old B-grid */
95	C-- By default for B-grid dynamics solver wind stress under sea-ice is
96	C set to the same value as it would be if there was no sea-ice.
97	C Define following CPP flag for B-grid ice-ocean stress coupling.
98	# define SEAICE_BICE_STRESS
99
100	C-- By default for B-grid dynamics solver surface tilt is obtained
101	C indirectly via geostrophic velocities. Define following CPP
102	C in order to use ETAN instead.
103	# define EXPLICIT_SSH_SLOPE
104	C-- Defining this flag turns on FV-discretization of the B-grid LSOR solver.
105	C It is smoother and includes all metric terms, similar to C-grid solvers.
106	C It is here for completeness, but its usefulness is unclear.
107	# undef SEAICE_LSRBNEW
108	#endif /* SEAICE_CGRID */
109
110	C-- When set limit the Ice-Loading to mass of 1/5 of Surface ocean grid-box
111	#undef SEAICE_CAP_ICELOAD
112	C-- When set use SEAICE_clipVelocties = .true., to clip U/VICE at 40cm/s,
113	C not recommended
114	#undef SEAICE_ALLOW_CLIPVELS
115	C-- When set cap the sublimation latent heat flux in solve4temp according
116	C to the available amount of ice+snow. Otherwise this term is treated
117	C like all of the others -- residuals heat and fw stocks are passed to
118	C the ocean at the end of seaice_growth in a conservative manner.
119	C SEAICE_CAP_SUBLIM is not needed as of now, but kept just in case.
120	#undef SEAICE_CAP_SUBLIM
121
122	C-- Enable free drift code
123	#define SEAICE_ALLOW_FREEDRIFT
124
125	#endif /* ALLOW_SEAICE */
126	#endif /* SEAICE_OPTIONS_H */
127
128	CEH3 ;;; Local Variables: ***
129	CEH3 ;;; mode:fortran ***
130	CEH3 ;;; End: ***