/[MITgcm]/MITgcm/verification/1D_ocean_ice_column/code_ad/SEAICE_OPTIONS.h
ViewVC logotype

Contents of /MITgcm/verification/1D_ocean_ice_column/code_ad/SEAICE_OPTIONS.h

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph


Revision 1.13 - (show annotations) (download)
Sat Dec 22 00:52:05 2012 UTC (11 years, 2 months ago) by dimitri
Branch: MAIN
CVS Tags: checkpoint64b
Changes since 1.12: +77 -27 lines
File MIME type: text/plain
Replace SEAICE_MULTICATEGORY check with "warning" instead of "error+stop".
Update SEAICE_OPTIONS.h for main code and "almost-up-to-date" verification experiments.

1 C $Header: /u/gcmpack/MITgcm/verification/1D_ocean_ice_column/code_ad/SEAICE_OPTIONS.h,v 1.12 2012/03/08 01:17:43 gforget 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 #undef 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 #undef SEAICE_ITD
47
48 C-- Use the Old version of seaice_growth (close to cvs version 1.70)
49 C otherwise, use the merged version (with some of Ian Fenty s code)
50 #undef SEAICE_GROWTH_LEGACY
51 #undef SEAICE_MODIFY_GROWTH_ADJ
52
53 C-- Since the missing sublimation term is now included
54 C this flag is needed for backward compatibility
55 #define SEAICE_DISABLE_SUBLIM
56
57 C-- Suspected missing term in coupled ocn-ice heat budget (to be confirmed)
58 #define SEAICE_DISABLE_HEATCONSFIX
59
60 C-- Default is constant seaice salinity (SEAICE_salt0); Define the following
61 C flag to consider (space & time) variable salinity: advected and forming
62 C seaice with a fraction (=SEAICE_saltFrac) of freezing seawater salinity.
63 C- Note: SItracer also offers an alternative way to handle variable salinity.
64 #define SEAICE_VARIABLE_SALINITY
65
66 C-- Tracers of ice and/or ice cover.
67 #undef ALLOW_SITRACER
68 #ifdef ALLOW_SITRACER
69 C-- To try avoid 'spontaneous generation' of tracer maxima by advdiff.
70 # define ALLOW_SITRACER_ADVCAP
71 #endif
72
73 C-- By default the seaice model is discretized on a B-Grid (for
74 C historical reasons). Define the following flag to use a new
75 C (not thoroughly) test version on a C-grid
76 #define SEAICE_CGRID
77
78 C-- Only for the C-grid version it is possible to
79 #ifdef SEAICE_CGRID
80 C enable JFNK code by defining the following flag
81 # undef SEAICE_ALLOW_JFNK
82 C enable LSR to use global (multi-tile) tri-diagonal solver
83 # undef SEAICE_GLOBAL_3DIAG_SOLVER
84 C enable EVP code by defining the following flag
85 # define SEAICE_ALLOW_EVP
86 # ifdef SEAICE_ALLOW_EVP
87 C-- When set use SEAICE_zetaMin and SEAICE_evpDampC to limit viscosities
88 C from below and above in seaice_evp: not necessary, and not recommended
89 # undef SEAICE_ALLOW_CLIPZETA
90 # endif /* SEAICE_ALLOW_EVP */
91 C regularize zeta to zmax with a smooth tanh-function instead
92 C of a min(zeta,zmax). This improves convergence of iterative
93 C solvers (Lemieux and Tremblay 2009, JGR). No effect on EVP
94 # undef SEAICE_ZETA_SMOOTHREG
95 C allow the truncated ellipse rheology (runtime flag SEAICEuseTEM)
96 # undef SEAICE_ALLOW_TEM
97 #else /* not SEAICE_CGRID, but old B-grid */
98 C-- By default for B-grid dynamics solver wind stress under sea-ice is
99 C set to the same value as it would be if there was no sea-ice.
100 C Define following CPP flag for B-grid ice-ocean stress coupling.
101 # define SEAICE_BICE_STRESS
102
103 C-- By default for B-grid dynamics solver surface tilt is obtained
104 C indirectly via geostrophic velocities. Define following CPP
105 C in order to use ETAN instead.
106 # define EXPLICIT_SSH_SLOPE
107 C-- Defining this flag turns on FV-discretization of the B-grid LSOR solver.
108 C It is smoother and includes all metric terms, similar to C-grid solvers.
109 C It is here for completeness, but its usefulness is unclear.
110 # undef SEAICE_LSRBNEW
111 #endif /* SEAICE_CGRID */
112
113 C-- When set use MAX_HEFF to cap seaice thickness in seaice_growth;
114 C currently only relevant for SEAICE_GROWTH_LEGACY
115 #undef SEAICE_CAP_HEFF
116 C-- When set limit the Ice-Loading to mass of 1/5 of Surface ocean grid-box
117 #undef SEAICE_CAP_ICELOAD
118 C-- When set use SEAICE_clipVelocties = .true., to clip U/VICE at 40cm/s,
119 C not recommended
120 #undef SEAICE_ALLOW_CLIPVELS
121 C-- When set cap the sublimation latent heat flux in solve4temp according
122 C to the available amount of ice+snow. Otherwise this term is treated
123 C like all of the others -- residuals heat and fw stocks are passed to
124 C the ocean at the end of seaice_growth in a conservative manner.
125 C SEAICE_CAP_SUBLIM is not needed as of now, but kept just in case.
126 #undef SEAICE_CAP_SUBLIM
127
128 C-- Enable free drift code
129 #undef SEAICE_ALLOW_FREEDRIFT
130
131 #endif /* ALLOW_SEAICE */
132 #endif /* SEAICE_OPTIONS_H */
133
134 CEH3 ;;; Local Variables: ***
135 CEH3 ;;; mode:fortran ***
136 CEH3 ;;; End: ***

  ViewVC Help
Powered by ViewVC 1.1.22