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

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

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

Revision 1.3 - (show annotations) (download)
Wed May 25 19:36:02 2011 UTC (12 years, 11 months ago) by gforget
Branch: MAIN
CVS Tags: checkpoint62y
Changes since 1.2: +3 -3 lines
File MIME type: text/plain 
-seaice_growth.F :
1) in EVOLUTION branch, always include a_QSWbyATM_cover in qnet (as it should be).
2) remove the Reproduce_v1p109 bloc that is essentially the same as the reference code.
3) remove duplicated blocs of code. After the a_QSWbyATM_cover fix (thanks to I. Fenty)
  FENTY_DELTA_HEFF_OPEN_WATER_FLUXES is redundant with SEAICE_DO_OPEN_WATER_GROWTH and
  FENTY_OPEN_WATER_FLUXES_MELT_ICE is redundant with SEAICE_DO_OPEN_WATER_MELT.
- seaice_check.F : FENTY_AREA_EXPANSION_CONTRACTION requires SEAICE_DO_OPEN_WATER_GROWTH.
- verif/1D_ocean_ice_column : update SEAICE_OPTIONS.h according to item 3).
- verif/global_ocean.cs32x15 : update results according to item 1)

note: it turns out that FENTY_OPEN_WATER_FLUXES_MELT_ICE is not actually
tested in 1D_ocean_ice_column over 10 time steps, so I did tests over
a year to check that SEAICE_DO_OPEN_WATER_MELT was indeed equivalent to
FENTY_OPEN_WATER_FLUXES_MELT_ICE. It is.
1 C $Header: /u/gcmpack/MITgcm/verification/1D_ocean_ice_column/code/SEAICE_OPTIONS.h,v 1.2 2011/04/28 02:09:00 ifenty 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 C-- Write "text-plots" of certain fields in STDOUT for debugging.
18 #undef SEAICE_DEBUG
19
20 C-- Allow sea-ice dynamic code.
21 C This option is provided to allow use of TAMC
22 C on the thermodynamics component of the code only.
23 C Sea-ice dynamics can also be turned off at runtime
24 C using variable SEAICEuseDYNAMICS.
25 #define SEAICE_ALLOW_DYNAMICS
26
27 C-- By default, the sea-ice package uses its own integrated bulk
28 C formulae to compute fluxes (fu, fv, EmPmR, Qnet, and Qsw) over
29 C open-ocean. When this flag is set, these variables are computed
30 C in a separate external package, for example, pkg/exf, and then
31 C modified for sea-ice effects by pkg/seaice.
32 #define SEAICE_EXTERNAL_FLUXES
33
34 C-- By default, the sea-ice package uses 2-category thermodynamics.
35 C When this flag is set, an 8-category calculation of ice
36 C thermodynamics is carried out in groatb.F
37 C Note the pickup_seaice.* generated by this option differ
38 C from those generated with the default 2-category model.
39 C Therefore it is not possible to switch between the two
40 C in the middle of an integration.
41 #undef SEAICE_MULTICATEGORY
42 #undef SEAICE_MODIFY_GROWTH_ADJ
43 #undef SEAICE_SIMPLIFY_GROWTH_ADJ
44
45 c Use the McPhee Ocean-Ice heat flux
46 #define MCPHEE_OCEAN_ICE_HEAT_FLUX
47 #define GRADIENT_MIXED_LAYER_TURBULENCE_FACTOR
48
49 #define SEAICE_DO_OPEN_WATER_GROWTH
50 #define SEAICE_DO_OPEN_WATER_MELT
51
52 #define FENTY_AREA_EXPANSION_CONTRACTION
53
54 #undef SEAICE_ADD_SUBLIMATION_TO_FWBUDGET
55
56
57
58 C-- By default cdm wind stress under sea-ice is set to the
59 C same value as it would be if there was no sea-ice.
60 C This is a placeholder until a more physically-sensible
61 C stable solution is developed. Apparently ocean stress
62 C that is based on sea-ice velocity can cause trouble in
63 C Arctic Ocean configurations. This temporary solution
64 C also corresponds to that adopted by CCSM-2. For the
65 C adventurous user who requires ice-modified stress,
66 C following options are available.
67 #undef SEAICE_ORIGINAL_BAD_ICE_STRESS
68 #undef SEAICE_TEST_ICE_STRESS_1
69
70 C-- By default the freezing point of water is set to the value of
71 C the parameter SEAICE_freeze (=-1.96 by default). To use a
72 C simple linear dependence of the freezing point on salinity,
73 C set the following flag (pressure is assumed to have no effect,
74 C which is a good assumption for the top 20 meters). With this
75 C option defined the parameter SEAICE_freeze has no effect.
76 #undef SEAICE_VARIABLE_FREEZING_POINT
77
78 C-- By default the seaice model is discretized on a B-Grid (for
79 C historical reasons). Define the following flag to use a new
80 C (not thoroughly) test version on a C-grid
81 #define SEAICE_CGRID
82
83 C-- Only for the C-grid version it is possible to enable EVP code by
84 C defining the following flag
85 #ifdef SEAICE_CGRID
86 #define SEAICE_ALLOW_EVP
87 #endif /* SEAICE_CGRID */
88
89 C-- Seaice flooding
90 #define ALLOW_SEAICE_FLOODING
91
92 C-- By default sea ice is fresh. Set following flag for salty ice.
93 #define SEAICE_VARIABLE_SALINITY
94
95 c-- Enable modified sea-ice thermodynamic by Ian Fenty
96 c-- Should improve adjoint, but doesnt work so far.
97 #undef SEAICE_ALLOW_TD_IF
98
99 #endif /* SEAICE_OPTIONS_H */
100
101 CEH3 ;;; Local Variables: ***
102 CEH3 ;;; mode:fortran ***
103 CEH3 ;;; End: ***
  ViewVC Help
Powered by ViewVC 1.1.22