| 1 | dimitri | 1.1 | C $Header: /u/gcmpack/MITgcm/verification/1D_ocean_ice_column/code/SEAICE_OPTIONS.h,v 1.4 2011/06/19 00:06:28 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 |  |  | #define 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 |  |  | cph#define MCPHEE_OCEAN_ICE_HEAT_FLUX | 
| 47 |  |  | cph#define GRADIENT_MIXED_LAYER_TURBULENCE_FACTOR | 
| 48 |  |  |  | 
| 49 |  |  | cph#define SEAICE_DO_OPEN_WATER_GROWTH | 
| 50 |  |  | cph#define SEAICE_DO_OPEN_WATER_MELT | 
| 51 |  |  |  | 
| 52 |  |  | cph#define FENTY_AREA_EXPANSION_CONTRACTION | 
| 53 |  |  |  | 
| 54 |  |  | #define 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 |  |  | #define 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 |  |  | #undef 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 |  |  | #endif /* SEAICE_OPTIONS_H */ | 
| 96 |  |  |  | 
| 97 |  |  | CEH3 ;;; Local Variables: *** | 
| 98 |  |  | CEH3 ;;; mode:fortran *** | 
| 99 |  |  | CEH3 ;;; End: *** |