/[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.4 by jmc, Sat Dec 4 00:19:32 2004 UTC revision 1.5 by jmc, Fri Dec 17 04:29:07 2004 UTC
# Line 26  C     === Global variables === Line 26  C     === Global variables ===
26  #include "PARAMS.h"  #include "PARAMS.h"
27  #include "GRID.h"  #include "GRID.h"
28  #include "THSICE_PARAMS.h"  #include "THSICE_PARAMS.h"
29    #ifdef ALLOW_MNC
30    #include "MNC_PARAMS.h"
31    #endif
32    
33  C     !INPUT/OUTPUT PARAMETERS:  C     !INPUT/OUTPUT PARAMETERS:
34  C     === Routine arguments ===  C     === Routine arguments ===
# Line 60  C--   Th-Sea-ICE parameter Line 63  C--   Th-Sea-ICE parameter
63       &         thSIce_deltaT, ocean_deltaT, tauRelax_MxL,       &         thSIce_deltaT, ocean_deltaT, tauRelax_MxL,
64       &         hMxL_default, sMxL_default, vMxL_default,       &         hMxL_default, sMxL_default, vMxL_default,
65       &         stressReduction,       &         stressReduction,
66       &         thSIce_taveFreq, thSIce_diagFreq, thSIce_monFreq       &         thSIce_taveFreq, thSIce_diagFreq, thSIce_monFreq,
67    c    &         thSIce_tave_mnc, thSIce_diag_mnc, thSIce_mon_mnc,
68         &         thSIceFract_InitFile, thSIceThick_InitFile,
69         &         thSIceSnowH_InitFile, thSIceSnowA_InitFile,
70         &         thSIceEnthp_InitFile, thSIceTsurf_InitFile
71    
72        _BEGIN_MASTER(myThid)        _BEGIN_MASTER(myThid)
73    
# Line 123  C--   Default values (parameters) Line 130  C--   Default values (parameters)
130        thSIce_taveFreq = taveFreq        thSIce_taveFreq = taveFreq
131        thSIce_diagFreq = dumpFreq        thSIce_diagFreq = dumpFreq
132        thSIce_monFreq  = monitorFreq        thSIce_monFreq  = monitorFreq
133    #ifdef ALLOW_MNC
134          thSIce_tave_mnc = timeave_mnc
135          thSIce_diag_mnc = snapshot_mnc
136          thSIce_mon_mnc  = monitor_mnc
137    #else
138          thSIce_tave_mnc = .FALSE.
139          thSIce_diag_mnc = .FALSE.
140          thSIce_mon_mnc  = .FALSE.
141    #endif
142          thSIceFract_InitFile = ' '
143          thSIceThick_InitFile = ' '
144          thSIceSnowH_InitFile = ' '
145          thSIceSnowA_InitFile = ' '
146          thSIceEnthp_InitFile = ' '
147          thSIceTsurf_InitFile = ' '
148    
149    
150  C--   Read parameters from open data file  C--   Read parameters from open data file
151        READ(UNIT=iUnit,NML=THSICE_CONST)        READ(UNIT=iUnit,NML=THSICE_CONST)
# Line 142  C-    Define other constants (from previ Line 165  C-    Define other constants (from previ
165        Tmlt1=-mu_Tf*S_winton        Tmlt1=-mu_Tf*S_winton
166        rhoiw = rhosw - rhoi        rhoiw = rhosw - rhoi
167    
168    C-    Set other parameters (from previous ones):
169          thSIce_tave_mdsio = .TRUE.
170          thSIce_diag_mdsio = .TRUE.
171          thSIce_mon_stdio  = .TRUE.
172    #ifdef ALLOW_MNC
173          IF (useMNC) THEN
174    C-------
175    C-    MNC not yet available in this package => turn off MNC flags for now
176            thSIce_tave_mnc = .FALSE.
177            thSIce_diag_mnc = .FALSE.
178            thSIce_mon_mnc  = .FALSE.
179    C-------
180            IF ( .NOT.outputTypesInclusive
181         &       .AND. thSIce_tave_mnc ) thSIce_tave_mdsio = .FALSE.
182            IF ( .NOT.outputTypesInclusive)
183         &       .AND. thSIce_diag_mnc ) thSIce_diag_mdsio = .FALSE.
184            IF ( .NOT.outputTypesInclusive
185         &       .AND. thSIce_mon_mnc  ) thSIce_mon_stdio  = .FALSE.
186          ENDIF
187    #endif
188    
189  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
190        iUnit = standardMessageUnit        iUnit = standardMessageUnit
191  c     iUnit = 88  c     CALL MDSFINDUNIT( iUnit, mythid )
192  c     OPEN(iUnit,file='thsice_check_params',status='unknown')  c     OPEN(iUnit,file='thsice_check_params',status='unknown')
193        WRITE(iUnit,*) 'ThSI: rhos    =',rhos        WRITE(iUnit,*) 'ThSI: rhos    =',rhos
194        WRITE(iUnit,*) 'ThSI: rhoi    =',rhoi        WRITE(iUnit,*) 'ThSI: rhoi    =',rhoi
# Line 195  c     OPEN(iUnit,file='thsice_check_para Line 239  c     OPEN(iUnit,file='thsice_check_para
239        WRITE(iUnit,*) 'ThSI: thSIce_diagFreq=',thSIce_diagFreq        WRITE(iUnit,*) 'ThSI: thSIce_diagFreq=',thSIce_diagFreq
240        WRITE(iUnit,*) 'ThSI: thSIce_monFreq =',thSIce_monFreq        WRITE(iUnit,*) 'ThSI: thSIce_monFreq =',thSIce_monFreq
241        WRITE(iUnit,*) 'ThSI: startIceModel =',startIceModel        WRITE(iUnit,*) 'ThSI: startIceModel =',startIceModel
242        IF (iUnit.EQ.88) CLOSE(iUnit)        IF (iUnit.NE.standardMessageUnit) CLOSE(iUnit)
243  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
244    
245        _END_MASTER(myThid)        _END_MASTER(myThid)
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22