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

Annotation 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 - (hide 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 gforget 1.3 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 ifenty 1.1 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 gforget 1.3 #define SEAICE_DO_OPEN_WATER_GROWTH
50     #define SEAICE_DO_OPEN_WATER_MELT
51 ifenty 1.1
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 ifenty 1.2 #define SEAICE_VARIABLE_SALINITY
94 ifenty 1.1
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