/[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.16 by jmc, Wed Jul 8 23:35:05 2009 UTC revision 1.17 by dfer, Wed Sep 23 20:24:47 2009 UTC
# Line 62  C--   Th-Sea-ICE parameter Line 62  C--   Th-Sea-ICE parameter
62        NAMELIST /THSICE_PARM01/        NAMELIST /THSICE_PARM01/
63       &     startIceModel, stepFwd_oceMxL, thSIce_calc_albNIR,       &     startIceModel, stepFwd_oceMxL, thSIce_calc_albNIR,
64       &     thSIce_deltaT, thSIce_dtTemp,       &     thSIce_deltaT, thSIce_dtTemp,
65       &     ocean_deltaT, tauRelax_MxL,       &     ocean_deltaT, tauRelax_MxL, tauRelax_MxL_salt,
66       &     hMxL_default, sMxL_default, vMxL_default,       &     hMxL_default, sMxL_default, vMxL_default,
67       &     thSIce_diffK, thSIceAdvScheme, stressReduction,       &     thSIce_diffK, thSIceAdvScheme, stressReduction,
68       &     thSIce_taveFreq, thSIce_diagFreq, thSIce_monFreq,       &     thSIce_taveFreq, thSIce_diagFreq, thSIce_monFreq,
# Line 133  C--   Default values (parameters) Line 133  C--   Default values (parameters)
133        thSIce_dtTemp  = UNSET_RL        thSIce_dtTemp  = UNSET_RL
134        ocean_deltaT    = dTtracerLev(1)        ocean_deltaT    = dTtracerLev(1)
135        tauRelax_MxL    = 0. _d 0        tauRelax_MxL    = 0. _d 0
136          tauRelax_MxL_salt = UNSET_RL
137        hMxL_default    = 50. _d 0        hMxL_default    = 50. _d 0
138        sMxL_default    = 35. _d 0        sMxL_default    = 35. _d 0
139        vMxL_default    = 5. _d -2        vMxL_default    = 5. _d -2
# Line 185  C     If using the same time step for bo Line 186  C     If using the same time step for bo
186  C     and ice thickness/growth, use thSIce_deltaT value  C     and ice thickness/growth, use thSIce_deltaT value
187        IF ( thSIce_dtTemp .EQ. UNSET_RL ) thSIce_dtTemp=thSIce_deltaT        IF ( thSIce_dtTemp .EQ. UNSET_RL ) thSIce_dtTemp=thSIce_deltaT
188    
189    C-    If undef, set salt relax to temperature relax
190          IF ( tauRelax_MxL_salt .EQ. UNSET_RL ) THEN
191               tauRelax_MxL_salt = tauRelax_MxL
192          ENDIF
193    
194  C-    Define other constants (from previous ones):  C-    Define other constants (from previous ones):
195        Tmlt1=-mu_Tf*S_winton        Tmlt1=-mu_Tf*S_winton
196        floodFac = (rhosw - rhoi)/rhos        floodFac = (rhosw - rhoi)/rhos
# Line 296  c     OPEN(iUnit,file='thsice_check_para Line 302  c     OPEN(iUnit,file='thsice_check_para
302        WRITE(iUnit,*) 'ThSI: ocean_deltaT  =',ocean_deltaT        WRITE(iUnit,*) 'ThSI: ocean_deltaT  =',ocean_deltaT
303        WRITE(iUnit,*) 'ThSI: stepFwd_oceMxL=',stepFwd_oceMxL        WRITE(iUnit,*) 'ThSI: stepFwd_oceMxL=',stepFwd_oceMxL
304        WRITE(iUnit,*) 'ThSI: tauRelax_MxL  =',tauRelax_MxL        WRITE(iUnit,*) 'ThSI: tauRelax_MxL  =',tauRelax_MxL
305          WRITE(iUnit,*) 'ThSI: tauRelax_MxL_salt  =',tauRelax_MxL_salt
306        WRITE(iUnit,*) 'ThSI: hMxL_default  =',hMxL_default        WRITE(iUnit,*) 'ThSI: hMxL_default  =',hMxL_default
307        WRITE(iUnit,*) 'ThSI: sMxL_default  =',sMxL_default        WRITE(iUnit,*) 'ThSI: sMxL_default  =',sMxL_default
308        WRITE(iUnit,*) 'ThSI: vMxL_default  =',vMxL_default        WRITE(iUnit,*) 'ThSI: vMxL_default  =',vMxL_default
Legend:
Removed from v.1.16  
changed lines
  Added in v.1.17
  ViewVC Help
Powered by ViewVC 1.1.22