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

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

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

-revision 1.2 by dimitri,
Fri Apr 27 22:25:23 2012 UTC
+revision 1.7 by torge,
Fri May  3 18:59:39 2013 UTC
 Line 7 
 C     | o Basic parameter header for sea
  C     *==========================================================*
  C--   COMMON /SEAICE_PARM_L/ Logical parameters of sea ice model.
- C
+ C - dynamics:
- C     SEAICEwriteState  :: If true, write sea ice state to file;
- C                          default is false.
  C     SEAICEuseDYNAMICS :: If false, do not use dynamics;
  C                          default is to use dynamics.
- C     SEAICEuseTEM      :: to use truncated ellipse method (see Geiger et al.
- C                          1998) set this parameter to true, default is false
- C     SEAICEuseEVP      :: If false, use Zhangs LSR solver for VP equations
- C                          if true use elastic viscous plastic solver
  C     SEAICEuseFREEDRIFT :: If True use free drift velocity instead of EVP
  C                           or LSR
- C     SEAICEheatConsFix  :: If true then fix ocn<->seaice advective heat flux.
+ C     SEAICEuseEVP      :: If true use elastic viscous plastic solver
+ C     SEAICEuseEVPstar  :: If true use modified elastic viscous plastic
+ C                          solver (EVP*) by Lemieux et al (2012)
  C     SEAICEuseEVPpickup :: Set to false in order to start EVP solver with
  C                          non-EVP pickup files.  Default is true.
  C                          Applied only if SEAICEuseEVP=.TRUE.
- C     SEAICEuseFluxForm :: use flux form for advection and diffusion
+ C     SEAICEuseMultiTileSolver :: in LSR, use full domain tri-diagonal solver
- C                          of seaice
+ C     SEAICEuseJFNK      :: If true, use Jacobi-free Newton-Krylov solver
+ C                           instead of LSR (default: false)
+ C     SEAICEuseTEM      :: to use truncated ellipse method (see Geiger et al.
+ C                          1998) set this parameter to true, default is false
+ C     SEAICEuseTilt      :: If true then include surface tilt term in dynamics
  C     SEAICEuseMetricTerms :: use metric terms for dynamics solver
  C                          (default = .true. )
- C     SEAICEuseFlooding :: turn on scheme to convert submerged snow into ice
+ C     SEAICE_no_slip    :: apply no slip boundary conditions to seaice velocity
+ C     SEAICE_maskRHS    :: mask the RHS of the solver where there is no ice
+ C     SEAICE_clipVelocities :: clip velocities to +/- 40cm/s
+ C     useHB87stressCoupling :: use an intergral over ice and ocean surface
+ C                          layer to define surface stresses on ocean
+ C                          following Hibler and Bryan (1987, JPO)
+ C - advection:
+ C     SEAICEuseFluxForm :: use flux form for advection and diffusion
+ C                          of seaice
  C     SEAICEadvHeff     :: turn on advection of effective thickness
  C                          (default = .true.)
  C     SEAICEadvArea     :: turn on advection of fraction area
-Line 35 
 C     SEAICEadvSnow     :: turn on advec
+Line 43 
 C     SEAICEadvSnow     :: turn on advec
  C                          non-default Leap-frog scheme for advection)
  C     SEAICEadvSalt     :: turn on advection of salt (does not work with
  C                          non-default Leap-frog scheme for advection)
- C     useHB87stressCoupling :: use an intergral over ice and ocean surface
+ C - thermodynamics:
- C                          layer to define surface stresses on ocean
- C                          following Hibler and Bryan (1987, JPO)
  C     usePW79thermodynamics :: use "0-layer" thermodynamics as described in
  C                           Parkinson and Washington (1979) and Hibler (1979)
+ C     SEAICE_useMultDimSnow :: use same fixed pdf for snow as for
+ C                              MULITCATEGORY ice
+ C     SEAICEuseFlooding  :: turn on scheme to convert submerged snow into ice
+ C     SEAICEheatConsFix  :: If true then fix ocn<->seaice advective heat flux.
  C     useMaykutSatVapPoly :: use Maykut Polynomial for saturation vapor pressure
  C                         instead of extended temp-range exponential law; def=F.
- C     SEAICE_mcPheeStepFunc    :: use step function (not linear tapering) in
+ C     SEAICE_mcPheeStepFunc :: use step function (not linear tapering) in
  C                           ocean-ice turbulent flux
  C     SEAICE_doOpenWaterGrowth :: use open water heat flux directly to grow ice
  C                           (when false cool ocean, and grow later if needed)
-Line 50 
 C     SEAICE_doOpenWaterMelt   :: use op
+Line 60 
 C     SEAICE_doOpenWaterMelt   :: use op
  C                           (when false warm ocean, and melt later if needed)
  C     SEAICE_salinityTracer    :: use SItracer to exchange and trace ocean
  C                           salt in ice
- C     SEAICE_age Tracer        :: use SItracer to trace the age of ice
+ C     SEAICE_ageTracer         :: use SItracer to trace the age of ice
  C     SEAICErestoreUnderIce :: restore surface T/S also underneath ice
  C                          ( default is false )
- C     SEAICE_no_slip    :: apply no slip boundary conditions to seaice velocity
+ C - other (I/O, ...):
- C     SEAICE_clipVelocities :: clip velocities to +/- 40cm/s
+ C     SEAICEwriteState  :: If true, write sea ice state to file;
- C     SEAICE_maskRHS    :: mask the RHS of the solver where there is no ice
+ C                          default is false.
  C     SEAICE_tave_mdsio :: write TimeAverage output using MDSIO
  C     SEAICE_dump_mdsio :: write snap-shot output   using MDSIO
  C     SEAICE_mon_stdio  :: write monitor to std-outp
-Line 63 
 C     SEAICE_tave_mnc   :: write TimeAve
+Line 73 
 C     SEAICE_tave_mnc   :: write TimeAve
  C     SEAICE_dump_mnc   :: write snap-shot output   using MNC
  C     SEAICE_mon_mnc    :: write monitor to netcdf file
        LOGICAL
-      &     SEAICEwriteState, SEAICEuseDYNAMICS, SEAICEuseEVP,
+      &     SEAICEuseDYNAMICS, SEAICEuseFREEDRIFT,
-      &     SEAICEuseFREEDRIFT, SEAICEuseTEM,
+      &     SEAICEuseEVP, SEAICEuseEVPstar, SEAICEuseEVPpickup,
-      &     SEAICEheatConsFix,
+      &     SEAICEuseMultiTileSolver,
-      &     SEAICEuseMetricTerms,
+      &     SEAICEuseJFNK,
-      &     SEAICEuseEVPpickup, SEAICEuseFlooding,
+      &     SEAICEuseTEM, SEAICEuseTilt, SEAICEuseMetricTerms,
-      &     SEAICEadvHeff, SEAICEadvArea,
+      &     SEAICE_no_slip, SEAICE_maskRHS,
+      &     SEAICE_clipVelocities, useHB87stressCoupling,
+      &     SEAICEuseFluxForm, SEAICEadvHeff, SEAICEadvArea,
       &     SEAICEadvSnow, SEAICEadvSalt,
-      &     SEAICEuseFluxForm, useHB87stressCoupling,
+      &     usePW79thermodynamics,
-      &     usePW79thermodynamics, useMaykutSatVapPoly,
+      &     SEAICE_useMultDimSnow, SEAICEuseFlooding, SEAICEheatConsFix,
-      &     SEAICE_mcPheeStepFunc,
+      &     useMaykutSatVapPoly, SEAICE_mcPheeStepFunc,
       &     SEAICE_doOpenWaterGrowth, SEAICE_doOpenWaterMelt,
       &     SEAICE_salinityTracer, SEAICE_ageTracer,
       &     SEAICErestoreUnderIce,
-      &     SEAICE_no_slip, SEAICE_clipVelocities, SEAICE_maskRHS,
+      &     SEAICEwriteState,
       &     SEAICE_tave_mdsio, SEAICE_dump_mdsio, SEAICE_mon_stdio,
       &     SEAICE_tave_mnc,   SEAICE_dump_mnc,   SEAICE_mon_mnc
        COMMON /SEAICE_PARM_L/
-      &     SEAICEwriteState, SEAICEuseDYNAMICS, SEAICEuseEVP,
+      &     SEAICEuseDYNAMICS, SEAICEuseFREEDRIFT,
-      &     SEAICEuseFREEDRIFT, SEAICEuseTEM,
+      &     SEAICEuseEVP, SEAICEuseEVPstar, SEAICEuseEVPpickup,
-      &     SEAICEheatConsFix,
+      &     SEAICEuseMultiTileSolver,
-      &     SEAICEuseMetricTerms,
+      &     SEAICEuseJFNK,
-      &     SEAICEuseEVPpickup, SEAICEuseFlooding,
+      &     SEAICEuseTEM, SEAICEuseTilt, SEAICEuseMetricTerms,
-      &     SEAICEadvHeff, SEAICEadvArea,
+      &     SEAICE_no_slip, SEAICE_maskRHS,
+      &     SEAICE_clipVelocities, useHB87stressCoupling,
+      &     SEAICEuseFluxForm, SEAICEadvHeff, SEAICEadvArea,
       &     SEAICEadvSnow, SEAICEadvSalt,
-      &     SEAICEuseFluxForm, useHB87stressCoupling,
+      &     usePW79thermodynamics,
-      &     usePW79thermodynamics, useMaykutSatVapPoly,
+      &     SEAICE_useMultDimSnow, SEAICEuseFlooding, SEAICEheatConsFix,
-      &     SEAICE_mcPheeStepFunc,
+      &     useMaykutSatVapPoly, SEAICE_mcPheeStepFunc,
       &     SEAICE_doOpenWaterGrowth, SEAICE_doOpenWaterMelt,
       &     SEAICE_salinityTracer, SEAICE_ageTracer,
       &     SEAICErestoreUnderIce,
-      &     SEAICE_no_slip, SEAICE_clipVelocities, SEAICE_maskRHS,
+      &     SEAICEwriteState,
       &     SEAICE_tave_mdsio, SEAICE_dump_mdsio, SEAICE_mon_stdio,
       &     SEAICE_tave_mnc,   SEAICE_dump_mnc,   SEAICE_mon_mnc
  C--   COMMON /SEAICE_PARM_I/ Integer valued parameters of sea ice model.
- C     IMAX_TICE         :: number of iterations for ice surface temp (default=10)
+ C     IMAX_TICE         :: number of iterations for ice surface temp
- C     postSolvTempIter :: select flux calculation after surf. temp solver iteration
+ C                          (default=10)
- C                         0 = none, i.e., from last iter ; 2 = full non-lin form
+ C     postSolvTempIter :: select flux calculation after surf. temp solver
- C                         1 = use linearized approx (consistent with tsurf finding)
+ C                         iteration
+ C                         0 = none, i.e., from last iter
+ C                         1 = use linearized approx (consistent with tsurf
+ C                             finding)
+ C                         2 = full non-lin form
  C     SOLV_MAX_ITERS    :: maximum number of allowed LSR-solver iterations
  C     SOLV_NCHECK       :: iteration interval for solver convergence test
  C     NPSEUDOTIMESTEPS  :: number of extra pseudo time steps (>= 2)
- C     LSR_mixIniGuess   :: control mixing of free-drift sol. into LSR initial guess
+ C     SEAICEnEVPstarSteps :: number of evp*-steps
- C                       :: =0 : no mix ; =2,4 : mix with (1/err)^2,4 factor
+ C     SEAICEnewtonIterMax :: maximum number of allowed Newton iterations
+ C                          in JFNK-solver
+ C     SEAICEkrylovIterMax :: maximum number of allowed Krylov iterations
+ C                          in JFNK-solver
+ C     SEAICE_JFNK_lsIter  :: number of Newton iterations after which the
+ C                            line search is started
+ C     SEAICE_JFNK_tolIter :: number of Newton iterations after which the
+ C                            the tolerance is relaxed again (default = 100)
+ C     SEAICE_OLx/y      :: overlaps for LSR-preconditioner in JFNK solver;
+ C                          for 0 < SEAICE_OLx/y 0 <= OLx/y-2
+ C                          the preconditioner is a restricted additive
+ C                           Schwarz method (default = OLx/y-2).
+ C     LSR_mixIniGuess   :: control mixing of free-drift sol. into LSR initial
+ C                          guess
+ C                       :: =0 : nothing; =1 : no mix, but print free-drift
+ C                          resid.;
+ C                       :: =2,4 : mix with (1/local-residual)^2,4 factor
+ C     SEAICEpresPow0    :: HEFF exponent for ice strength below SEAICEpresH0
+ C     SEAICEpresPow1    :: HEFF exponent for ice strength above SEAICEpresH0
  C     SEAICEadvScheme   :: sets the advection scheme for thickness and area
  C     SEAICEadvSchArea  :: sets the advection scheme for area
  C     SEAICEadvSchHeff  :: sets the advection scheme for effective thickness
-Line 118 
 C     SEAICE_areaLossFormula :: selects
+Line 154 
 C     SEAICE_areaLossFormula :: selects
  C                        :: 1=from all but only melt conributions by ATM and OCN
  C                        :: 2=from net melt-growth>0 by ATM and OCN
  C                        :: 3=from predicted melt by ATM
- C     SEAICE_areaGainFormula :: selects formula for ice cover gain from open water growth
+ C     SEAICE_areaGainFormula :: selects formula for ice cover gain from open
+ C                               water growth
  C                        :: 1=from growth by ATM
  C                        :: 2=from predicted growth by ATM
+ C     SEAICEetaZmethod   :: determines how shear-viscosity eta is computed at
+ C                           Z-points
+ C                           0=simple averaging from C-points (default and old)
+ C                           1=linear averaging of strain rates to Z-points
+ C                           2=averaging of squares of strain rates
  C     SEAICE_multDim     :: number of ice categories
  C     SEAICE_debugPointI :: I,J index for seaice-specific debuggin
  C     SEAICE_debugPointJ
-Line 129 
 C
+Line 171 
 C
        INTEGER SOLV_MAX_ITERS, SOLV_NCHECK
        INTEGER NPSEUDOTIMESTEPS
        INTEGER LSR_mixIniGuess
+       INTEGER SEAICEnEVPstarSteps
+       INTEGER SEAICEnewtonIterMax, SEAICEkrylovIterMax
+       INTEGER SEAICE_JFNK_lsIter, SEAICE_JFNK_tolIter
+       INTEGER SEAICE_OLx, SEAICE_OLy
        INTEGER SEAICEadvScheme
        INTEGER SEAICEadvSchArea
        INTEGER SEAICEadvSchHeff
-Line 137 
 C
+Line 183 
 C
        INTEGER SEAICEadjMODE
        INTEGER SEAICE_areaLossFormula
        INTEGER SEAICE_areaGainFormula
+       INTEGER SEAICEetaZmethod
        INTEGER SEAICE_multDim
        INTEGER SEAICE_debugPointI
        INTEGER SEAICE_debugPointJ
+       INTEGER SEAICEpresPow0, SEAICEpresPow1
        COMMON /SEAICE_PARM_I/
       &     IMAX_TICE, postSolvTempIter,
       &     SOLV_MAX_ITERS, SOLV_NCHECK,
       &     NPSEUDOTIMESTEPS,
       &     LSR_mixIniGuess,
+      &     SEAICEnEVPstarSteps,
+      &     SEAICEnewtonIterMax, SEAICEkrylovIterMax,
+      &     SEAICE_JFNK_lsIter, SEAICE_OLx, SEAICE_OLy,
+      &     SEAICE_JFNK_tolIter,
+      &     SEAICEpresPow0, SEAICEpresPow1,
       &     SEAICEadvScheme,
       &     SEAICEadvSchArea,
       &     SEAICEadvSchHeff,
-Line 154 
 C
+Line 207 
 C
       &     SEAICE_areaLossFormula,
       &     SEAICE_areaGainFormula,
       &     SEAICE_multDim,
+      &     SEAICEetaZmethod,
       &     SEAICE_debugPointI,
       &     SEAICE_debugPointJ
-Line 182 
 C
+Line 236 
 C
  C--   COMMON /SEAICE_PARM_RL/ Real valued parameters of sea ice model.
  C     SEAICE_deltaTtherm :: Seaice timestep for thermodynamic equations (s)
  C     SEAICE_deltaTdyn   :: Seaice timestep for dynamic solver          (s)
+ C     SEAICE_LSRrelaxU/V :: relaxation parameter for LSR-solver: U/V-component
  C     SEAICE_deltaTevp   :: Seaice timestep for EVP solver              (s)
  C     SEAICE_elasticParm :: parameter that sets relaxation timescale
  C                           tau = SEAICE_elasticParm * SEAICE_deltaTdyn
  C     SEAICE_evpTauRelax :: relaxation timescale tau                    (s)
  C     SEAICE_evpDampC    :: evp damping constant (Hunke,JCP,2001)       (kg/m^2)
+ C     SEAICE_evpAlpha    :: dimensionless parameter 2*evpTauRelax/deltaTevp
+ C     SEAICE_evpBeta     :: dimensionless parameter deltaTdyn/deltaTevp
+ C     JFNKgamma_nonlin   :: non-linear tolerance parameter for JFNK solver
+ C     JFNKgamma_lin_min/max :: tolerance parameters for linear JFNK solver
+ C     JFNKres_t          :: tolerance parameter for FGMRES residual
+ C     JFNKres_tFac       :: if set, JFNKres_t=JFNKres_tFac*(initial residual)
+ C     SEAICE_JFNKepsilon :: step size for the FD-gradient in s/r seaice_jacvec
  C     SEAICE_zetaMaxFac  :: factor determining the maximum viscosity    (s)
  C                          (default = 5.e+12/2.e4 = 2.5e8)
  C     SEAICE_zetaMin     :: lower bound for viscosity (default = 0)    (N s/m^2)
+ C     SEAICEpresH0       :: HEFF threshold for ice strength            (m)
  C     SEAICE_monFreq     :: SEAICE monitor frequency.                   (s)
  C     SEAICE_dumpFreq    :: SEAICE dump frequency.                      (s)
  C     SEAICE_taveFreq    :: SEAICE time-averaging frequency.            (s)
-Line 276 
 C     SEAICEdiffKhSnow   :: sets the dif
+Line 339 
 C     SEAICEdiffKhSnow   :: sets the dif
  C     SEAICEdiffKhSalt   :: sets the diffusivity for sea ice salinity (m^2/s)
  C     SEAICE_airTurnAngle   :: turning angles of air-ice interfacial stress
  C     SEAICE_waterTurnAngle :: and ice-water interfacial stress (in degrees)
  C     SEAICE_tauAreaObsRelax :: Timescale of relaxation to observed
  C                               sea ice concentration (s), default=unset
  C
        _RL SEAICE_deltaTtherm, SEAICE_deltaTdyn, SEAICE_deltaTevp
+       _RL SEAICE_LSRrelaxU, SEAICE_LSRrelaxV
        _RL SEAICE_monFreq, SEAICE_dumpFreq, SEAICE_taveFreq
        _RL SEAICE_initialHEFF
        _RL SEAICE_rhoAir, SEAICE_rhoIce, SEAICE_rhoSnow, ICE2WATR
-Line 303 
 C
+Line 367 
 C
        _RL facOpenGrow, facOpenMelt
        _RL SEAICE_tempFrz0, SEAICE_dTempFrz_dS
        _RL OCEAN_drag, LSR_ERROR, DIFF1
+       _RL JFNKgamma_nonlin, JFNKres_t, JFNKres_tFac
+       _RL JFNKgamma_lin_min, JFNKgamma_lin_max, SEAICE_JFNKepsilon
        _RL SEAICE_area_reg, SEAICE_hice_reg
        _RL SEAICE_area_floor, SEAICE_area_max
        _RL SEAICE_airTurnAngle, SEAICE_waterTurnAngle
        _RL SEAICE_elasticParm, SEAICE_evpTauRelax
+       _RL SEAICE_evpAlpha, SEAICE_evpBeta
        _RL SEAICE_evpDampC, SEAICE_zetaMin, SEAICE_zetaMaxFac
+       _RL SEAICEpresH0
        _RL SEAICEdiffKhArea, SEAICEdiffKhHeff, SEAICEdiffKhSnow
        _RL SEAICEdiffKhSalt
        _RL SEAICE_tauAreaObsRelax
        COMMON /SEAICE_PARM_RL/
       &    SEAICE_deltaTtherm, SEAICE_deltaTdyn,
+      &    SEAICE_LSRrelaxU, SEAICE_LSRrelaxV,
       &    SEAICE_deltaTevp, SEAICE_elasticParm, SEAICE_evpTauRelax,
+      &    SEAICE_evpAlpha, SEAICE_evpBeta,
       &    SEAICE_evpDampC, SEAICE_zetaMin, SEAICE_zetaMaxFac,
+      &    SEAICEpresH0,
       &    SEAICE_monFreq, SEAICE_dumpFreq, SEAICE_taveFreq,
       &    SEAICE_initialHEFF,
       &    SEAICE_rhoAir, SEAICE_rhoIce, SEAICE_rhoSnow, ICE2WATR,
-Line 338 
 C
+Line 409 
 C
       &    facOpenGrow, facOpenMelt,
       &    SEAICE_tempFrz0, SEAICE_dTempFrz_dS,
       &    OCEAN_drag, LSR_ERROR, DIFF1,
+      &    JFNKgamma_nonlin, JFNKres_t, JFNKres_tFac,
+      &    JFNKgamma_lin_min, JFNKgamma_lin_max, SEAICE_JFNKepsilon,
       &    SEAICE_area_reg, SEAICE_hice_reg,
       &    SEAICE_area_floor, SEAICE_area_max,
       &    SEAICEdiffKhArea, SEAICEdiffKhHeff, SEAICEdiffKhSnow,
-Line 345 
 C
+Line 418 
 C
       &    SEAICE_airTurnAngle, SEAICE_waterTurnAngle
  C--   COMMON /SEAICE_BOUND_RL/ Various bounding values
- C     MAX_HEFF          :: maximum ice thickness     (m)
  C     MIN_ATEMP         :: minimum air temperature   (deg C)
  C     MIN_LWDOWN        :: minimum downward longwave (W/m^2)
  C     MIN_TICE          :: minimum ice temperature   (deg C)
  C     SEAICE_EPS        :: small number used to reduce derivative singularities
  C     SEAICE_EPS_SQ     :: small number square
- CToM<<<
- C     Hlimit            :: ice thickness category limits (m), array of size nITD+1
- C     Hlimit_c1,_c2,_c3 :: coefficients set in seaice_readparams.F to calculate Hlimit in seaice_init_fixed.F
- C>>>ToM
  C
-       _RL MAX_HEFF, MIN_ATEMP, MIN_LWDOWN, MIN_TICE
+       _RL MIN_ATEMP, MIN_LWDOWN, MIN_TICE
        _RL SEAICE_EPS, SEAICE_EPS_SQ
- CToM<<<
+       COMMON /SEAICE_BOUND_RL/
+      &     MIN_ATEMP, MIN_LWDOWN, MIN_TICE,
+      &     SEAICE_EPS, SEAICE_EPS_SQ
+ #ifdef SEAICE_ITD
+ C     Hlimit            :: ice thickness category limits (m), array of
+ C                          size nITD+1
+ C     Hlimit_c1,_c2,_c3 :: coefficients set in seaice_readparams.F to
+ C                          calculate Hlimit in seaice_init_fixed.F
        _RL Hlimit(0:nITD)
        _RL Hlimit_c1, Hlimit_c2, Hlimit_c3
- C>>>ToM
+       COMMON /SEAICE_BOUND_ITD_RL/
-       COMMON /SEAICE_BOUND_RL/
-      &     MAX_HEFF, MIN_ATEMP, MIN_LWDOWN, MIN_TICE,
- CToM<<<
- C     &     SEAICE_EPS, SEAICE_EPS_SQ
-      &     SEAICE_EPS, SEAICE_EPS_SQ,
       &     Hlimit,
       &     Hlimit_c1,Hlimit_c2,Hlimit_c3
- C>>>ToM
+ #endif /* SEAICE_ITD */
  C--   Constants used by sea-ice model
        _RL         ZERO           , ONE           , TWO

 Legend:



Removed from v.1.2
 


changed lines


 
Added in v.1.7
 Legend:



Removed from v.1.2
 


changed lines


 
Added in v.1.7
-Removed from v.1.2
+Added in v.1.7

	ViewVC Help
Powered by ViewVC 1.1.22