/[MITgcm]/MITgcm_contrib/torge/itd/code/SEAICE_OPTIONS.h

Diff of /MITgcm_contrib/torge/itd/code/SEAICE_OPTIONS.h

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

-revision 1.3 by torge,
Wed Sep 26 16:59:29 2012 UTC
+revision 1.4 by torge,
Tue Oct 30 21:51:44 2012 UTC
 Line 1
  C $Header$
+ C $Name$
  C     *==========================================================*
  C     | SEAICE_OPTIONS.h
-Line 33 
 C     in a separate external package, fo
+Line 34 
 C     in a separate external package, fo
  C     modified for sea-ice effects by pkg/seaice.
  #define SEAICE_EXTERNAL_FLUXES
- C--   By default, the sea-ice package uses 2-category thermodynamics.
+ C--   The actual number of ice categories used to solve for seaice flux is
- C     When this flag is set, an 8-category calculation of ice
+ C     now a run-time parameter (SEAICE_multDim).
- C     thermodynamics is carried out in groatb.F
+ C     This CPP-flag will be completely removed soon (no longer in main code);
- C     Note the pickup_seaice.* generated by this option differ
+ C     it is just used to set default number of categories, i.e., =1 if undef,
- C     from those generated with the default 2-category model.
+ C     or =MULTDIM if defined (MULTDIM=7 in default SEAICE_SIZE.h).
- C     Therefore it is not possible to switch between the two
+ C     Note: be aware of pickup_seaice.* compatibility issues when restarting
- C     in the middle of an integration.
+ C     a simulation with a different number of categories.
- #define SEAICE_MULTICATEGORY
+ c#define SEAICE_MULTICATEGORY
- CToM<<<
- C--   run with sea ice thickness distribution;
- C     set number of categories (nITD) in SIZE.h
- #define SEAICE_ITD
- C
- C--   SEAICE_ITD replaces SEAICE_MULTICATEGORY
- #ifdef SEAICE_ITD
- #undef SEAICE_MULTICATEGORY
- #endif
- C>>>ToM
  C--   Use the Old version of seaice_growth (close to cvs version 1.70)
  C     otherwise, use the merged version (with some of Ian Fenty s code)
-Line 61 
 C--   Since the missing sublimation term
+Line 51 
 C--   Since the missing sublimation term
  C     this flag is needed for backward compatibility
  #undef SEAICE_DISABLE_SUBLIM
- C--   Suspected missing term in coupled ocn-ice heat budget (to be further tested)
+ C--   Suspected missing term in coupled ocn-ice heat budget (to be confirmed)
  #undef SEAICE_DISABLE_HEATCONSFIX
- C--   options only available in the merged version
+ C--   Default is constant seaice salinity (SEAICE_salt0); Define the following
- ccc#ifndef SEAICE_GROWTH_LEGACY
+ C     flag to consider (space & time) variable salinity: advected and forming
- C-    to ensure heat conservation in the coupled ocean-seaice system
+ C     seaice with a fraction (=SEAICE_saltFrac) of freezing seawater salinity.
- ccc# undef SEAICE_HEAT_CONSERV_FIX
+ C- Note: SItracer also offers an alternative way to handle variable salinity.
- C-    to preclude infinitesimal ice concentrations:
+ #undef SEAICE_VARIABLE_SALINITY
- ccc# undef ALLOW_PRECLUDE_INFINITESIMAL_AREA
- ccc#endif
- CToM<<<
- C--   The freezing point of water is a simple linear dependence to salinity.
- C     The flag below is now obsolete, but kept to ease the transition
- C     from the old constant freezing point default (see seaice_check.F).
- #define SEAICE_VARIABLE_FREEZING_POINT
- C>>>ToM
- C--   By default sea ice is fresh.  Set following flag for salty ice.
- #define SEAICE_VARIABLE_SALINITY
  C--   Tracers of ice and/or ice cover.
- #define ALLOW_SITRACER
+ #undef ALLOW_SITRACER
  #ifdef ALLOW_SITRACER
  C--   To try avoid 'spontaneous generation' of tracer maxima by advdiff.
  # define ALLOW_SITRACER_ADVCAP
-Line 95 
 C     (not thoroughly) test version on a
+Line 74 
 C     (not thoroughly) test version on a
  C--   Only for the C-grid version it is possible to
  #ifdef SEAICE_CGRID
+ C     enable JFNK code by defining the following flag
+ # define SEAICE_ALLOW_JFNK
  C     enable EVP code by defining the following flag
  # define SEAICE_ALLOW_EVP
  # ifdef SEAICE_ALLOW_EVP
- C--   When set use SEAICE_zetaMin and SEAICE_evpDampC to limit
+ C--   When set use SEAICE_zetaMin and SEAICE_evpDampC to limit viscosities
- C--   viscosities from below and above in seaice_evp
+ C     from below and above in seaice_evp: not necessary, and not recommended
- C--   not necessary, and not recommended
  #  undef SEAICE_ALLOW_CLIPZETA
  # endif /* SEAICE_ALLOW_EVP */
+ C     regularize zeta to zmax with a smooth tanh-function instead
+ C     of a min(zeta,zmax). This improves convergence of iterative
+ C     solvers (Lemieux and Tremblay 2009, JGR). No effect on EVP
+ # undef SEAICE_ZETA_SMOOTHREG
  C     allow the truncated ellipse rheology (runtime flag SEAICEuseTEM)
  # undef SEAICE_ALLOW_TEM
  #else /* not SEAICE_CGRID, but old B-grid */
-Line 115 
 C--   By default for B-grid dynamics sol
+Line 99 
 C--   By default for B-grid dynamics sol
  C     indirectly via geostrophic velocities. Define following CPP
  C     in order to use ETAN instead.
  # define EXPLICIT_SSH_SLOPE
- C--   Defining this flag turns on a FV-discretization of the B-grid LSOR
+ C--   Defining this flag turns on FV-discretization of the B-grid LSOR solver.
- C     solver. It is smoother and includes all metric terms, similar to the
+ C     It is smoother and includes all metric terms, similar to C-grid solvers.
- C     C-grid solver. It is here for completeness, but its usefulness is
+ C     It is here for completeness, but its usefulness is unclear.
- C     unclear.
  # undef SEAICE_LSRBNEW
  #endif /* SEAICE_CGRID */
- C--   When set use MAX_HEFF to cap sea ice thickness in seaice_growth
+ C--   When set use MAX_HEFF to cap seaice thickness in seaice_growth;
+ C     currently only relevant for SEAICE_GROWTH_LEGACY
  #undef SEAICE_CAP_HEFF
  C--   When set limit the Ice-Loading to mass of 1/5 of Surface ocean grid-box
  #undef SEAICE_CAP_ICELOAD
  C--   When set use SEAICE_clipVelocties = .true., to clip U/VICE at 40cm/s,
- C--   not recommended
+ C     not recommended
  #undef SEAICE_ALLOW_CLIPVELS
- C--   When set cap the sublimation latent heat flux in solve4temp
+ C--   When set cap the sublimation latent heat flux in solve4temp according
- C     according to the available amount of ice+snow. Otherwise this
+ C     to the available amount of ice+snow. Otherwise this term is treated
- C     term is treated like all of the others -- residuals heat and
+ C     like all of the others -- residuals heat and fw stocks are passed to
- C     fw stocks are passed to the ocean at the end of seaice_growth
+ C     the ocean at the end of seaice_growth in a conservative manner.
- C     in a conservative manner. SEAICE_CAP_SUBLIM is not needed
+ C     SEAICE_CAP_SUBLIM is not needed as of now, but kept just in case.
- C     as of now, but kept just in case.
  #undef SEAICE_CAP_SUBLIM
- C     enable free drift code
+ C--   Enable free drift code
  #undef SEAICE_ALLOW_FREEDRIFT
+ CToM<<<
+ C--   run with sea ice thickness distribution;
+ C     set number of categories (nITD) in SIZE.h
+ #undef SEAICE_ITD
+ C
+ C--   SEAICE_ITD replaces SEAICE_MULTICATEGORY
+ #ifdef SEAICE_ITD
+ #undef SEAICE_MULTICATEGORY
+ #endif
+ C>>>ToM
  #endif /* ALLOW_SEAICE */
  #endif /* SEAICE_OPTIONS_H */

 Legend:



Removed from v.1.3
 


changed lines


 
Added in v.1.4
 Legend:



Removed from v.1.3
 


changed lines


 
Added in v.1.4
-Removed from v.1.3
+Added in v.1.4

	ViewVC Help
Powered by ViewVC 1.1.22