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

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

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

revision 1.2.2.1 by dimitri, Sat Feb 15 12:40:30 2003 UTC revision 1.34 by mlosch, Fri Apr 22 09:13:50 2016 UTC
# Line 1  Line 1 
1  C  C $Header$
2    C $Name$
 C     /==========================================================\  
 C     | SEAICE_OPTIONS.h                                         |  
 C     | o CPP options file for sea ice package.                  |  
 C     |==========================================================|  
 C     | Use this file for selecting options within the sea ice   |  
 C     | package.  Sea ice model is enabled with ALLOW_SEAICE in  |  
 C     | CPP_OPTIONS.h                                            |  
 C     \==========================================================/  
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"  #include "CPP_OPTIONS.h"
16    
17  #ifdef ALLOW_SEAICE  #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
# 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  
   
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
33  C     in a separate external package, for example, pkg/exf, and then  C     in a separate external package, for example, pkg/exf, and then
34  C     modified for sea-ice effects by pkg/seaice.  C     modified for sea-ice effects by pkg/seaice.
35  #undef SEAICE_EXTERNAL_FLUXES  #define SEAICE_EXTERNAL_FLUXES
36    
37    C--   This CPP flag has been retired.  The number of ice categories
38    C     used to solve for seaice flux is now specified by run-time
39    C     parameter SEAICE_multDim.
40    C     Note: be aware of pickup_seaice.* compatibility issues when
41    C     restarting a simulation with a different number of categories.
42    c#define SEAICE_MULTICATEGORY
43    
44    C--   run with sea Ice Thickness Distribution (ITD);
45    C     set number of categories (nITD) in SEAICE_SIZE.h
46    #undef SEAICE_ITD
47    
48    C--   Since the missing sublimation term is now included
49    C     this flag is needed for backward compatibility
50    #undef SEAICE_DISABLE_SUBLIM
51    
52    C--   Suspected missing term in coupled ocn-ice heat budget (to be confirmed)
53    #undef SEAICE_DISABLE_HEATCONSFIX
54    
55    C--   Default is constant seaice salinity (SEAICE_salt0); Define the following
56    C     flag to consider (space & time) variable salinity: advected and forming
57    C     seaice with a fraction (=SEAICE_saltFrac) of freezing seawater salinity.
58    C- Note: SItracer also offers an alternative way to handle variable salinity.
59    #undef SEAICE_VARIABLE_SALINITY
60    
61    C--   Tracers of ice and/or ice cover.
62    #define ALLOW_SITRACER
63    #ifdef ALLOW_SITRACER
64    C--   To try avoid 'spontaneous generation' of tracer maxima by advdiff.
65    # define ALLOW_SITRACER_ADVCAP
66    #endif
67    
68    C--   By default the seaice model is discretized on a B-Grid (for
69    C     historical reasons). Define the following flag to use a new
70    C     (not thoroughly) test version on a C-grid
71    #define SEAICE_CGRID
72    
73    C--   Only for the C-grid version it is possible to
74    #ifdef SEAICE_CGRID
75    C     enable JFNK code by defining the following flag
76    # undef SEAICE_ALLOW_JFNK
77    C     enable LSR to use global (multi-tile) tri-diagonal solver
78    # undef SEAICE_GLOBAL_3DIAG_SOLVER
79    C     enable EVP code by defining the following flag
80    # define SEAICE_ALLOW_EVP
81    # ifdef SEAICE_ALLOW_EVP
82    C--   When set use SEAICE_zetaMin and SEAICE_evpDampC to limit viscosities
83    C     from below and above in seaice_evp: not necessary, and not recommended
84    #  undef SEAICE_ALLOW_CLIPZETA
85    # endif /* SEAICE_ALLOW_EVP */
86    C     regularize zeta to zmax with a smooth tanh-function instead
87    C     of a min(zeta,zmax). This improves convergence of iterative
88    C     solvers (Lemieux and Tremblay 2009, JGR). No effect on EVP
89    # undef SEAICE_ZETA_SMOOTHREG
90    C     allow the truncated ellipse rheology (runtime flag SEAICEuseTEM)
91    # undef SEAICE_ALLOW_TEM
92    C     Use LSR vector code; not useful on non-vector machines, because it
93    C     slows down convergence considerably, but the extra iterations are
94    C     more than made up by the much faster code on vector machines. For
95    C     the only regularly test vector machine these flags a specified
96    C     in the build options file SUPER-UX_SX-8_sxf90_awi, so that we comment
97    C     them out here.
98    # undef SEAICE_VECTORIZE_LSR
99    C     Use zebra-method (alternate lines) for line-successive-relaxation
100    C     This modification improves the convergence of the vector code
101    C     dramatically, so that is may actually be useful in general, but
102    C     that needs to be tested. Can be used without vectorization options.
103    # define SEAICE_LSR_ZEBRA
104    C     Use parameterisation of grounding ice for a better representation
105    C     of fastice in shallow seas
106    # define SEAICE_ALLOW_BOTTOMDRAG
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 FV-discretization of the B-grid LSOR solver.
118    C     It is smoother and includes all metric terms, similar to C-grid solvers.
119    C     It is here for completeness, but its usefulness is unclear.
120    # undef SEAICE_LSRBNEW
121    #endif /* SEAICE_CGRID */
122    
123    C--   When set limit the Ice-Loading to mass of 1/5 of Surface ocean grid-box
124    #undef SEAICE_CAP_ICELOAD
125    C--   When set use SEAICE_clipVelocties = .true., to clip U/VICE at 40cm/s,
126    C     not recommended
127    #undef SEAICE_ALLOW_CLIPVELS
128    C--   When set cap the sublimation latent heat flux in solve4temp according
129    C     to the available amount of ice+snow. Otherwise this term is treated
130    C     like all of the others -- residuals heat and fw stocks are passed to
131    C     the ocean at the end of seaice_growth in a conservative manner.
132    C     SEAICE_CAP_SUBLIM is not needed as of now, but kept just in case.
133    #undef SEAICE_CAP_SUBLIM
134    
135    C--   Enable free drift code
136    #define SEAICE_ALLOW_FREEDRIFT
137    
138  #endif /* ALLOW_SEAICE */  #endif /* ALLOW_SEAICE */
139    #endif /* SEAICE_OPTIONS_H */
140    
141    CEH3 ;;; Local Variables: ***
142    CEH3 ;;; mode:fortran ***
143    CEH3 ;;; End: ***

Legend:
Removed from v.1.2.2.1  
changed lines
  Added in v.1.34

  ViewVC Help
Powered by ViewVC 1.1.22