/[MITgcm]/MITgcm/pkg/thsice/thsice_readparms.F
ViewVC logotype

Diff of /MITgcm/pkg/thsice/thsice_readparms.F

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

revision 1.11 by jmc, Wed Apr 4 02:11:35 2007 UTC revision 1.13 by jscott, Mon Apr 9 17:44:13 2007 UTC
# Line 62  C--   Th-Sea-ICE parameter Line 62  C--   Th-Sea-ICE parameter
62    
63        NAMELIST /THSICE_PARM01/        NAMELIST /THSICE_PARM01/
64       &     startIceModel, stepFwd_oceMxL,       &     startIceModel, stepFwd_oceMxL,
65       &     thSIce_deltaT, ocean_deltaT, tauRelax_MxL,       &     thSIce_deltaT, thSIce_dtTemp,
66         &     ocean_deltaT, tauRelax_MxL,
67       &     hMxL_default, sMxL_default, vMxL_default,       &     hMxL_default, sMxL_default, vMxL_default,
68       &     thSIce_diffK, thSIceAdvScheme, stressReduction,       &     thSIce_diffK, thSIceAdvScheme, stressReduction,
69       &     thSIce_taveFreq, thSIce_diagFreq, thSIce_monFreq,       &     thSIce_taveFreq, thSIce_diagFreq, thSIce_monFreq,
# Line 126  C--   Default values (parameters) Line 127  C--   Default values (parameters)
127        stepFwd_oceMxL  = .FALSE.        stepFwd_oceMxL  = .FALSE.
128        startIceModel   = 0        startIceModel   = 0
129        thSIce_deltaT   = dTtracerLev(1)        thSIce_deltaT   = dTtracerLev(1)
130          thSIce_dtTemp  = UNSET_RL
131        ocean_deltaT    = dTtracerLev(1)        ocean_deltaT    = dTtracerLev(1)
132        tauRelax_MxL    = 0. _d 0        tauRelax_MxL    = 0. _d 0
133        hMxL_default    = 50. _d 0        hMxL_default    = 50. _d 0
# Line 173  C--   Read parameters from open data fil Line 175  C--   Read parameters from open data fil
175  C--   Close the open data file  C--   Close the open data file
176        CLOSE(iUnit)        CLOSE(iUnit)
177    
178    C     If using the same time step for both icetop temp solver
179    C     and ice thickness/growth, use thSIce_deltaT value
180          IF ( thSIce_dtTemp .EQ. UNSET_RL ) thSIce_dtTemp=thSIce_deltaT
181    
182  C-    Define other constants (from previous ones):  C-    Define other constants (from previous ones):
183        Tmlt1=-mu_Tf*S_winton        Tmlt1=-mu_Tf*S_winton
184        rhoiw = rhosw - rhoi        rhoiw = rhosw - rhoi
# Line 197  C     Set I/O parameters Line 203  C     Set I/O parameters
203        ENDIF        ENDIF
204  #endif  #endif
205    
206    C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
207  C--   Check parameter consistency:  C--   Check parameter consistency:
208          IF ( thSIceAdvScheme.EQ.0 .AND. thSIce_diffK.NE.0. ) THEN
209            WRITE(msgBuf,'(2A)')
210         &   'THSICE_READPARMS: to use thSIce_diffK, needs to select',
211         &   ' one advection scheme (thSIceAdvScheme<>0)'
212            CALL PRINT_ERROR( msgBuf , myThid )
213            STOP 'ABNORMAL END: THSICE_READPARMS'
214          ENDIF
215    #ifndef ALLOW_GENERIC_ADVDIFF
216          IF ( thSIceAdvScheme.NE.0 ) THEN
217            WRITE(msgBuf,'(2A)')
218         &   'THSICE_READPARMS: Need to compile "generic_advdiff" pkg',
219         &   ' in order to use thSIceAdvScheme'
220            CALL PRINT_ERROR( msgBuf , myThid )
221            STOP 'ABNORMAL END: THSICE_READPARMS'
222          ENDIF
223    #endif /* ndef ALLOW_GENERIC_ADVDIFF */
224    
225        IF ( useSEAICE .AND. stressReduction.NE.0. _d 0 ) THEN        IF ( useSEAICE .AND. stressReduction.NE.0. _d 0 ) THEN
226  C--     If useSEAICE=.true., the stress is computed in seaice_model,  C--     If useSEAICE=.true., the stress is computed in seaice_model,
227  C--     so that it does not need any further reduction  C--     so that it does not need any further reduction
# Line 282  C--   Everyone else must wait for the pa Line 306  C--   Everyone else must wait for the pa
306    
307        RETURN        RETURN
308        END        END
   
Legend:
Removed from v.1.11  
changed lines
  Added in v.1.13
  ViewVC Help
Powered by ViewVC 1.1.22