/[MITgcm]/MITgcm/pkg/seaice/SEAICE_OPTIONS.h
ViewVC logotype

Contents of /MITgcm/pkg/seaice/SEAICE_OPTIONS.h

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

Revision 1.58 - (show annotations) (download)
Tue Feb 7 18:41:09 2012 UTC (12 years, 3 months ago) by gforget
Branch: MAIN
Changes since 1.57: +1 -4 lines
File MIME type: text/plain 
pkg/seaice : turbulent flux term specification via run time params.
- replace CPP switches with run time switches to facilitate testing
- retired CPPs : MCPHEE_OCEAN_ICE_HEAT_FLUX, GRADIENT_MIXED_LAYER_TURBULENCE_FACTOR
- added runtime parameter : SEAICEturbFluxFormula
- move Mcphee constants to SEAICE_PARAMS.h
- seaice_readparms.F : read turbulent flux term params (SEAICEturbFluxFormula;
  SEAICE_availHeatFrac/Frz, or SEAICE_gamma_t/_frz if apply); remove
  seaice_summary call (moved to seaice_init_fixed)
- seaice_init_fixed.F : set SEAICE_availHeatFrac/Frz depending on
  SEAICEturbFluxFormula, and seaice_summary call (moved from seaice_readparms)
- seaice_summary.F : add missing report of read turbulent flux term set up (and SITRACER).
- seaice_check.F : check range; inform and stop if retired CPPs are used (incl. SEAICE_AGE).
- seaice_growth.F : unified use of SEAICE_availHeatFrac/Frz, without
	CPP switch, same computation as before for the different options.
1 C $Header: /u/gcmpack/MITgcm/pkg/seaice/SEAICE_OPTIONS.h,v 1.57 2012/02/07 02:19:14 jmc 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-- By default, the sea-ice package uses 2-category thermodynamics.
38 C When this flag is set, an 8-category calculation of ice
39 C thermodynamics is carried out in groatb.F
40 C Note the pickup_seaice.* generated by this option differ
41 C from those generated with the default 2-category model.
42 C Therefore it is not possible to switch between the two
43 C in the middle of an integration.
44 #undef SEAICE_MULTICATEGORY
45
46 C-- Use the Old version of seaice_growth (close to cvs version 1.70)
47 C otherwise, use the merged version (with some of Ian Fenty s code)
48 #define SEAICE_GROWTH_LEGACY
49
50 C-- options only available in the merged version (from Ian Fenty s code)
51 #ifndef SEAICE_GROWTH_LEGACY
52 C- to ensure heat conservation in the coupled ocean-seaice system
53 # undef SEAICE_HEAT_CONSERV_FIX
54 C- to switch on/off open-water freezing contribution to thickness tendency:
55 # define SEAICE_DO_OPEN_WATER_GROWTH
56 C- ifdef SEAICE_DO_OPEN_WATER_GROWTH then define SEAICE_DO_OPEN_WATER_MELT
57 C to also allow open-water air-sea heat fluxes melt ice
58 # undef SEAICE_DO_OPEN_WATER_MELT
59 C- to switch on/off ocean heat contribution to seaice cover reduction:
60 # define SEAICE_OCN_MELT_ACT_ON_AREA
61 C- to preclude infinitesimal ice concentrations:
62 # undef ALLOW_PRECLUDE_INFINITESIMAL_AREA
63 #endif
64
65 C-- By default the freezing point of water is set to the value of
66 C the parameter SEAICE_freeze (=-1.96 by default). To use a
67 C simple linear dependence of the freezing point on salinity,
68 C set the following flag (pressure is assumed to have no effect,
69 C which is a good assumption for the top 20 meters). With this
70 C option defined the parameter SEAICE_freeze has no effect.
71 #undef SEAICE_VARIABLE_FREEZING_POINT
72
73 C-- Allow SEAICEuseFlooding, which converts snow to ice if submerged.
74 #define ALLOW_SEAICE_FLOODING
75
76 C-- By default sea ice is fresh. Set following flag for salty ice.
77 #undef SEAICE_VARIABLE_SALINITY
78
79 C-- Tracers of ice and/or ice cover.
80 #undef ALLOW_SITRACER
81 #ifdef ALLOW_SITRACER
82 C-- To try avoid 'spontaneous generation' of tracer maxima by advdiff.
83 # define ALLOW_SITRACER_ADVCAP
84 #endif
85
86 C-- By default the seaice model is discretized on a B-Grid (for
87 C historical reasons). Define the following flag to use a new
88 C (not thoroughly) test version on a C-grid
89 #define SEAICE_CGRID
90
91 C-- Only for the C-grid version it is possible to
92 #ifdef SEAICE_CGRID
93 C enable EVP code by defining the following flag
94 # define SEAICE_ALLOW_EVP
95 # ifdef SEAICE_ALLOW_EVP
96 C-- When set use SEAICE_zetaMin and SEAICE_evpDampC to limit
97 C-- viscosities from below and above in seaice_evp
98 C-- not necessary, and not recommended
99 # undef SEAICE_ALLOW_CLIPZETA
100 # endif /* SEAICE_ALLOW_EVP */
101 C allow the truncated ellipse rheology (runtime flag SEAICEuseTEM)
102 # undef SEAICE_ALLOW_TEM
103 #else /* not SEAICE_CGRID, but old B-grid */
104 C-- By default for B-grid dynamics solver wind stress under sea-ice is
105 C set to the same value as it would be if there was no sea-ice.
106 C Define following CPP flag for B-grid ice-ocean stress coupling.
107 # define SEAICE_BICE_STRESS
108
109 C-- By default for B-grid dynamics solver surface tilt is obtained
110 C indirectly via geostrophic velocities. Define following CPP
111 C in order to use ETAN instead.
112 # define EXPLICIT_SSH_SLOPE
113 C-- Defining this flag turns on a FV-discretization of the B-grid LSOR
114 C solver. It is smoother and includes all metric terms, similar to the
115 C C-grid solver. It is here for completeness, but its usefulness is
116 C unclear.
117 # undef SEAICE_LSRBNEW
118 #endif /* SEAICE_CGRID */
119
120 C-- When set use MAX_HEFF to cap sea ice thickness in seaice_growth
121 #undef SEAICE_CAP_HEFF
122 C-- When set limit the Ice-Loading to mass of 1/5 of Surface ocean grid-box
123 #undef SEAICE_CAP_ICELOAD
124 C-- When set use SEAICE_clipVelocties = .true., to clip U/VICE at 40cm/s,
125 C-- not recommended
126 #undef SEAICE_ALLOW_CLIPVELS
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