/[MITgcm]/MITgcm_contrib/torge/itd/code/seaice_readparms.F
ViewVC logotype

Diff of /MITgcm_contrib/torge/itd/code/seaice_readparms.F

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

revision 1.1 by dimitri, Fri Apr 27 22:22:17 2012 UTC revision 1.3 by torge, Wed Sep 26 17:03:29 2012 UTC
# Line 131  C--   SEAICE parameters Line 131  C--   SEAICE parameters
131       & MAX_HEFF, MIN_ATEMP, MIN_LWDOWN, MAX_TICE, MIN_TICE,       & MAX_HEFF, MIN_ATEMP, MIN_LWDOWN, MAX_TICE, MIN_TICE,
132       & SEAICE_EPS, SEAICE_EPS_SQ,       & SEAICE_EPS, SEAICE_EPS_SQ,
133       & SEAICE_tave_mnc, SEAICE_dump_mnc, SEAICE_mon_mnc,       & SEAICE_tave_mnc, SEAICE_dump_mnc, SEAICE_mon_mnc,
134       & SEAICE_debugPointI, SEAICE_debugPointJ  CToM<<<
135    C     & SEAICE_debugPointI, SEAICE_debugPointJ
136         & SEAICE_debugPointI, SEAICE_debugPointJ,
137         & Hlimit_c1, Hlimit_c2, Hlimit_c3
138    C>>>ToM
139    
140  #ifdef ALLOW_COST  #ifdef ALLOW_COST
141        NAMELIST /SEAICE_PARM02/        NAMELIST /SEAICE_PARM02/
# Line 224  C     Hunke, JCP, 2001 use 615 kg/m^2 fo Line 228  C     Hunke, JCP, 2001 use 615 kg/m^2 fo
228        SEAICE_mon_mnc  = .FALSE.        SEAICE_mon_mnc  = .FALSE.
229  #endif  #endif
230        SEAICE_initialHEFF = ZERO        SEAICE_initialHEFF = ZERO
231    CToM<<<
232    C     Coefficients used to calculate sea ice thickness category limits
233    C     after Lipscomb et al. (2001, JGR), Equ. 22
234    C     choose between
235    C      - original parameters of Lipscomb et al. (2001):
236    C        c1=3.0/N, c2=15*c1, c3=3.0
237    C      - and a higher resolution of thin end of ITD:
238    C        c1=1.5/N, c2=42*c1, c3=3.3
239          Hlimit_c1          = 3.0
240          Hlimit_c2          = 15.
241          Hlimit_c3          = 3.0
242    C>>>ToM
243        SEAICE_rhoIce      = 0.91   _d +03        SEAICE_rhoIce      = 0.91   _d +03
244        SEAICE_rhoSnow     = 330.   _d 0        SEAICE_rhoSnow     = 330.   _d 0
245        ICE2WATR           = UNSET_RL        ICE2WATR           = UNSET_RL
# Line 284  C     old default value of 0.97001763668 Line 300  C     old default value of 0.97001763668
300        SEAICE_shortwave   = 0.30   _d 0        SEAICE_shortwave   = 0.30   _d 0
301        SEAICE_salt0       = 0.0    _d 0        SEAICE_salt0       = 0.0    _d 0
302        SEAICE_saltFrac    = 0.0    _d 0        SEAICE_saltFrac    = 0.0    _d 0
303  #ifdef SEAICE_MULTICATEGORY  CToM<<<
304    C#ifdef SEAICE_MULTICATEGORY
305    #if defined(SEAICE_MULTICATEGORY) || defined(SEAICE_ITD)
306    C in case defined(SEAICE_ITD) MULTDIM = nITD (see SEAICE_SIZE.h)
307    C because SEAICE_ITD uses code related to SEAICE_MULTICATEGORY
308    C>>>ToM
309        SEAICE_multDim     = MULTDIM        SEAICE_multDim     = MULTDIM
310  #else  #else
311        SEAICE_multDim     = 1        SEAICE_multDim     = 1
# Line 348  c     DIFF1      = .004 _d 0 Line 369  c     DIFF1      = .004 _d 0
369        SEAICE_area_reg   = siEPS        SEAICE_area_reg   = siEPS
370  #endif  #endif
371        SEAICE_hice_reg   = 0.05 _d 0        SEAICE_hice_reg   = 0.05 _d 0
372        SEAICE_area_max    = 1.00 _d 0        SEAICE_area_max   = 1.00 _d 0
373    
374        SEAICE_airTurnAngle   = 0.0 _d 0        SEAICE_airTurnAngle   = 0.0 _d 0
375        SEAICE_waterTurnAngle = 0.0 _d 0        SEAICE_waterTurnAngle = 0.0 _d 0
# Line 620  C     to set SEAICE_mcPheePiston once dr Line 641  C     to set SEAICE_mcPheePiston once dr
641          STOP 'ABNORMAL END: S/R SEAICE_READPARMS'          STOP 'ABNORMAL END: S/R SEAICE_READPARMS'
642        ENDIF        ENDIF
643    
644    CToM<<<
645    #ifdef SEAICE_ITD
646    C     The ice thickness distribution module can only be used with
647    C      package seaice calling
648    C      SEAICE_ADVDIFF and
649    C      SEAICE_GROWTH, i.e. needs usePW79thermodynamics = .TRUE.
650          useTHSice = .FALSE.
651    #endif
652    C>>>ToM
653        IF ( useThSice ) THEN        IF ( useThSice ) THEN
654  C     If the thsice package with the Winton thermodynamics is used  C     If the thsice package with the Winton thermodynamics is used
655  C     is does not make sense to have the following parameters defined,  C     is does not make sense to have the following parameters defined,
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22