/[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.8 by jmc, Thu May 25 18:03:24 2006 UTC revision 1.9 by jmc, Sun Jun 25 22:31:02 2006 UTC
# Line 130  C--   Default values (parameters) Line 130  C--   Default values (parameters)
130        sMxL_default    = 35. _d 0        sMxL_default    = 35. _d 0
131        vMxL_default    = 5. _d -2        vMxL_default    = 5. _d -2
132        stressReduction = 1. _d 0        stressReduction = 1. _d 0
133          IF ( useSEAICE ) stressReduction = 0. _d 0
134        thSIce_taveFreq = taveFreq        thSIce_taveFreq = taveFreq
135        thSIce_diagFreq = dumpFreq        thSIce_diagFreq = dumpFreq
136        thSIce_monFreq  = monitorFreq        thSIce_monFreq  = monitorFreq
# Line 192  C     Set I/O parameters Line 193  C     Set I/O parameters
193        ENDIF        ENDIF
194  #endif  #endif
195    
196    C--   Check parameter consistency:
197          IF ( useSEAICE .AND. stressReduction.NE.0. _d 0 ) THEN
198    C--     If useSEAICE=.true., the stress is computed in seaice_model,
199    C--     so that it does not need any further reduction
200            WRITE(msgBuf,'(2A)')
201         &   'THSICE_READPARMS: if useSEAICE, stress will be computed',
202         &   ' by SEAICE pkg => no reduction'
203            CALL PRINT_MESSAGE( msgBuf, errorMessageUnit,
204         &                    SQUEEZE_RIGHT , myThid)
205            WRITE(msgBuf,'(A)')
206         &   'THSICE_READPARMS: WARNING: reset stressReduction to zero'
207            CALL PRINT_MESSAGE( msgBuf, errorMessageUnit,
208         &                    SQUEEZE_RIGHT , myThid)
209            stressReduction = 0. _d 0
210          ENDIF
211    
212  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
213        iUnit = standardMessageUnit        iUnit = standardMessageUnit
214  c     CALL MDSFINDUNIT( iUnit, mythid )  c     CALL MDSFINDUNIT( iUnit, mythid )
Legend:
Removed from v.1.8  
changed lines
  Added in v.1.9
  ViewVC Help
Powered by ViewVC 1.1.22