/[MITgcm]/MITgcm/verification/lab_sea/code_ad/SEAICE_OPTIONS.h
ViewVC logotype

Diff of /MITgcm/verification/lab_sea/code_ad/SEAICE_OPTIONS.h

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

revision 1.3 by heimbach, Wed Sep 26 04:38:20 2007 UTC revision 1.16 by gforget, Tue Feb 7 18:51:04 2012 UTC
# Line 1  Line 1 
1  C $Header$  C $Header$
2    C $Name$
3    
4  C     /==========================================================\  C     *==========================================================*
5  C     | SEAICE_OPTIONS.h                                         |  C     | SEAICE_OPTIONS.h
6  C     | o CPP options file for sea ice package.                  |  C     | o CPP options file for sea ice package.
7  C     |==========================================================|  C     *==========================================================*
8  C     | Use this file for selecting options within the sea ice   |  C     | Use this file for selecting options within the sea ice
9  C     | package.                                                 |  C     | package.
10  C     \==========================================================/  C     *==========================================================*
11    
12  #ifndef SEAICE_OPTIONS_H  #ifndef SEAICE_OPTIONS_H
13  #define SEAICE_OPTIONS_H  #define SEAICE_OPTIONS_H
14  #include "PACKAGES_CONFIG.h"  #include "PACKAGES_CONFIG.h"
15  #include "CPP_OPTIONS.h"  #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.  C--   Write "text-plots" of certain fields in STDOUT for debugging.
21  #undef SEAICE_DEBUG  #undef SEAICE_DEBUG
22    
# Line 23  C     Sea-ice dynamics can also be turne Line 27  C     Sea-ice dynamics can also be turne
27  C     using variable SEAICEuseDYNAMICS.  C     using variable SEAICEuseDYNAMICS.
28  #define SEAICE_ALLOW_DYNAMICS  #define SEAICE_ALLOW_DYNAMICS
29    
 C--   By default, pkg/seaice defines forcing variable internally, in  
 C     SEAICE_FFIELDS.h, and reads them in from files.  When CPP option  
 C     SEAICE_EXTERNAL_FORCING is defined, the above forcing variables  
 C     are defined and provided by an external package.  At present this  
 C     option is hardwired for pkg/exf and the variables are passed using  
 C     include file exf_fields.h.  
 #define SEAICE_EXTERNAL_FORCING  
 #ifdef SEAICE_EXTERNAL_FORCING  
 cph#include "EXF_OPTIONS.h"  
 #endif /* SEAICE_EXTERNAL_FORCING */  
   
30  C--   By default, the sea-ice package uses its own integrated bulk  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  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  C     open-ocean.  When this flag is set, these variables are computed
# Line 48  C     Note the pickup_seaice.* generated Line 41  C     Note the pickup_seaice.* generated
41  C     from those generated with the default 2-category model.  C     from those generated with the default 2-category model.
42  C     Therefore it is not possible to switch between the two  C     Therefore it is not possible to switch between the two
43  C     in the middle of an integration.  C     in the middle of an integration.
44  #undef SEAICE_MULTILEVEL  #define SEAICE_MULTICATEGORY
45    
46  C--   By default cdm wind stress under sea-ice is set to the  C--   Use the Old version of seaice_growth (close to cvs version 1.70)
47  C     same value as it would be if there was no sea-ice.  C     otherwise, use the merged version (with some of Ian Fenty s code)
48  C     This is a placeholder until a more physically-sensible  #define SEAICE_GROWTH_LEGACY
49  C     stable solution is developed.  Apparently ocean stress  
50  C     that is based on sea-ice velocity can cause trouble in  C--   options only available in the merged version (from Ian Fenty s code)
51  C     Arctic Ocean configurations.  This temporary solution  #ifndef SEAICE_GROWTH_LEGACY
52  C     also corresponds to that adopted by CCSM-2.  For the  C-    to ensure heat conservation in the coupled ocean-seaice system
53  C     adventurous user who requires ice-modified stress,  # undef SEAICE_HEAT_CONSERV_FIX
54  C     following options are available.  C-    to switch on/off open-water freezing contribution to thickness tendency:
55  #undef SEAICE_ORIGINAL_BAD_ICE_STRESS  # define SEAICE_DO_OPEN_WATER_GROWTH
56  #undef SEAICE_TEST_ICE_STRESS_1  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--   Use the Old version of seaice_solve4temp (formerly seaice_budget_ice)
66    C     otherwise, use Ian Fenty s version
67    #define SEAICE_SOLVE4TEMP_LEGACY
68    
69  C--   By default the freezing point of water is set to the value of  C--   By default the freezing point of water is set to the value of
70  C     the parameter SEAICE_freeze (=-1.96 by default). To use a  C     the parameter SEAICE_freeze (=-1.96 by default). To use a
71  C     simple linear dependence of the freezing point on salinity,  C     simple linear dependence of the freezing point on salinity,
72  C     set the following flag (pressure is assumed to have no effect,  C     set the following flag (pressure is assumed to have no effect,
73  C     which is a good assumption for the top 20 meters). With this  C     which is a good assumption for the top 20 meters). With this
74  C     option defined the parameter SEAICE_freeze has no effect.  C     option defined the parameter SEAICE_freeze has no effect.
75  #undef SEAICE_VARIABLE_FREEZING_POINT  #undef SEAICE_VARIABLE_FREEZING_POINT
76    
77  C--   By default the seaice model is discretized on a B-Grid (for  C--   Allow SEAICEuseFlooding, which converts snow to ice if submerged.
78    #define ALLOW_SEAICE_FLOODING
79    
80    C--   By default sea ice is fresh.  Set following flag for salty ice.
81    #define SEAICE_VARIABLE_SALINITY
82    
83    C--   Tracers of ice and/or ice cover.
84    #undef ALLOW_SITRACER
85    #ifdef ALLOW_SITRACER
86    C--   To try avoid 'spontaneous generation' of tracer maxima by advdiff.
87    # define ALLOW_SITRACER_ADVCAP
88    #endif
89    
90    C--   By default the seaice model is discretized on a B-Grid (for
91  C     historical reasons). Define the following flag to use a new  C     historical reasons). Define the following flag to use a new
92  C     (not thoroughly) test version on a C-grid  C     (not thoroughly) test version on a C-grid
93  #define SEAICE_CGRID  #define SEAICE_CGRID
94    
95  C--   Only for the C-grid version it is possible to enable EVP code by  C--   Only for the C-grid version it is possible to
 C     defining the following flag  
96  #ifdef SEAICE_CGRID  #ifdef SEAICE_CGRID
97  #undef SEAICE_ALLOW_EVP  C     enable EVP code by defining the following flag
98    # define SEAICE_ALLOW_EVP
99    # ifdef SEAICE_ALLOW_EVP
100    C--   When set use SEAICE_zetaMin and SEAICE_evpDampC to limit
101    C--   viscosities from below and above in seaice_evp
102    C--   not necessary, and not recommended
103    #  undef SEAICE_ALLOW_CLIPZETA
104    # endif /* SEAICE_ALLOW_EVP */
105    C     allow the truncated ellipse rheology (runtime flag SEAICEuseTEM)
106    # undef SEAICE_ALLOW_TEM
107    #else /* not SEAICE_CGRID, but old B-grid */
108    C--   By default for B-grid dynamics solver wind stress under sea-ice is
109    C     set to the same value as it would be if there was no sea-ice.
110    C     Define following CPP flag for B-grid ice-ocean stress coupling.
111    # define SEAICE_BICE_STRESS
112    
113    C--   By default for B-grid dynamics solver surface tilt is obtained
114    C     indirectly via geostrophic velocities. Define following CPP
115    C     in order to use ETAN instead.
116    # define EXPLICIT_SSH_SLOPE
117    C--   Defining this flag turns on a FV-discretization of the B-grid LSOR
118    C     solver. It is smoother and includes all metric terms, similar to the
119    C     C-grid solver. It is here for completeness, but its usefulness is
120    C     unclear.
121    # undef SEAICE_LSRBNEW
122  #endif /* SEAICE_CGRID */  #endif /* SEAICE_CGRID */
123    
124  C--   Seaice flooding  C--   When set use MAX_HEFF to cap sea ice thickness in seaice_growth
125  #define ALLOW_SEAICE_FLOODING  #undef SEAICE_CAP_HEFF
126    C--   When set limit the Ice-Loading to mass of 1/5 of Surface ocean grid-box
127  C--   By default sea ice is fresh.  Set following flag for salty ice.  #undef SEAICE_CAP_ICELOAD
128  #define SEAICE_SALINITY  C--   When set use SEAICE_clipVelocties = .true., to clip U/VICE at 40cm/s,
129    C--   not recommended
130    #undef SEAICE_ALLOW_CLIPVELS
131    
132    C     enable free drift code
133    #define SEAICE_ALLOW_FREEDRIFT
134    
135    #endif /* ALLOW_SEAICE */
136  #endif /* SEAICE_OPTIONS_H */  #endif /* SEAICE_OPTIONS_H */
137    
138  CEH3 ;;; Local Variables: ***  CEH3 ;;; Local Variables: ***

Legend:
Removed from v.1.3  
changed lines
  Added in v.1.16

  ViewVC Help
Powered by ViewVC 1.1.22