/[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.7 by torge, Wed Mar 27 18:59:53 2013 UTC revision 1.8 by torge, Wed Mar 27 19:13:51 2013 UTC
# Line 142  C--   SEAICE parameters Line 142  C--   SEAICE parameters
142       & SEAICE_monFreq, SEAICE_dumpFreq, SEAICE_taveFreq,       & SEAICE_monFreq, SEAICE_dumpFreq, SEAICE_taveFreq,
143       & SEAICE_tave_mnc, SEAICE_dump_mnc, SEAICE_mon_mnc,       & SEAICE_tave_mnc, SEAICE_dump_mnc, SEAICE_mon_mnc,
144  #ifdef SEAICE_ITD  #ifdef SEAICE_ITD
145       & Hlimit_c1, Hlimit_c2, Hlimit_c3       & Hlimit_c1, Hlimit_c2, Hlimit_c3,
146  #endif  #endif
147       & SEAICE_debugPointI, SEAICE_debugPointJ       & SEAICE_debugPointI, SEAICE_debugPointJ
148    
# Line 249  C     Hunke, JCP, 2001 use 615 kg/m^2 fo Line 249  C     Hunke, JCP, 2001 use 615 kg/m^2 fo
249        SEAICEnEVPstarSteps = UNSET_I        SEAICEnEVPstarSteps = UNSET_I
250        SEAICE_initialHEFF = ZERO        SEAICE_initialHEFF = ZERO
251  #ifdef SEAICE_ITD  #ifdef SEAICE_ITD
 CToM<<<  
252  C     Coefficients used to calculate sea ice thickness category limits  C     Coefficients used to calculate sea ice thickness category limits
253  C     after Lipscomb et al. (2001, JGR), Equ. 22  C     after Lipscomb et al. (2001, JGR), Equ. 22
254  C     choose between  C     choose between
# Line 260  C        c1=1.5/N, c2=42*c1, c3=3.3 Line 259  C        c1=1.5/N, c2=42*c1, c3=3.3
259        Hlimit_c1          = 3.0        Hlimit_c1          = 3.0
260        Hlimit_c2          = 15.        Hlimit_c2          = 15.
261        Hlimit_c3          = 3.0        Hlimit_c3          = 3.0
 C>>>ToM  
262  #endif  #endif
263        SEAICE_rhoIce      = 0.91   _d +03        SEAICE_rhoIce      = 0.91   _d +03
264        SEAICE_rhoSnow     = 330.   _d 0        SEAICE_rhoSnow     = 330.   _d 0
# Line 322  C     old default value of 0.97001763668 Line 320  C     old default value of 0.97001763668
320        SEAICE_shortwave   = 0.30   _d 0        SEAICE_shortwave   = 0.30   _d 0
321        SEAICE_salt0       = 0.0    _d 0        SEAICE_salt0       = 0.0    _d 0
322        SEAICE_saltFrac    = 0.0    _d 0        SEAICE_saltFrac    = 0.0    _d 0
 CToM<<<  
323  #ifdef SEAICE_ITD  #ifdef SEAICE_ITD
324  C in case defined(SEAICE_ITD) MULTDIM = nITD (see SEAICE_SIZE.h)  C in case defined(SEAICE_ITD) MULTDIM = nITD (see SEAICE_SIZE.h)
325  C>>>ToM  c      SEAICE_multDim     = MULTDIM
326        SEAICE_multDim     = MULTDIM  C the switch MULTICATEGORY (and with it parameter MULTDIM) has been retired
327    C  and SEAICE_multDim is now a runtime parameter;
328    C  in case SEAICE_multDim is given in data.seaice it needs to be overwritten
329    C  after PARM01 was read (see below)
330          SEAICE_multDim     = nITD
331  #else  #else
332        SEAICE_multDim     = 1        SEAICE_multDim     = 1
333  #endif  #endif
# Line 514  C--   Read settings from model parameter Line 515  C--   Read settings from model parameter
515    
516        CLOSE(iUnit)        CLOSE(iUnit)
517    
518    #ifdef SEAICE_ITD
519    C SEAICE_multDim has become a runtime parameter but if SEAICE_ITD is defined
520    C  it needs to equal nITD because of shared code (mostly in seaice_growth.F).
521    C nITD is set in SEAICE_SIZE.h
522          SEAICE_multDim     = nITD
523    #endif
524    
525        WRITE(msgBuf,'(A)')        WRITE(msgBuf,'(A)')
526       &     ' SEAICE_READPARMS: finished reading data.seaice'       &     ' SEAICE_READPARMS: finished reading data.seaice'
527        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
Legend:
Removed from v.1.7  
changed lines
  Added in v.1.8
  ViewVC Help
Powered by ViewVC 1.1.22