pkg/thsice/thsice_readparms.F

C $Header: /u/gcmpack/MITgcm/pkg/thsice/thsice_readparms.F,v 1.8 2006/05/25 18:03:24 jmc Exp $
C $Name:  $

#include "THSICE_OPTIONS.h"

CBOP
C     !ROUTINE: THSICE_READPARMS
C     !INTERFACE:
      SUBROUTINE THSICE_READPARMS( myThid )

C     !DESCRIPTION: \bv
C     *==========================================================*
C     | S/R THSICE_READPARMS
C     | o Routine to initialize THSICE parameters and constants
C     *==========================================================*
C     | Initialize Th-Sea-ICE parameters, read in data.ice
C     *==========================================================*
C     \ev

C     !USES:
      IMPLICIT NONE

C     === Global variables ===
#include "SIZE.h"
#include "EEPARAMS.h"
#include "PARAMS.h"
#include "GRID.h"
#include "THSICE_PARAMS.h"
#ifdef ALLOW_MNC
#include "MNC_PARAMS.h"
#endif

C     !INPUT/OUTPUT PARAMETERS:
C     === Routine arguments ===
      INTEGER myThid
CEOP

#ifdef ALLOW_THSICE

C     === Local variables ===
C     msgBuf      - Informational/error message buffer
C     iUnit       - Work variable for IO unit number
      CHARACTER*(MAX_LEN_MBUF) msgBuf
      INTEGER iUnit

C--   Th-Sea-ICE parameter
      NAMELIST /THSICE_CONST/
     &  rhos, rhoi, rhosw, rhofw,
     &  cpice, cpwater,
     &  kice, ksnow,
     &  transcoef, Lfresh, qsnow,
     &  albColdSnow, albWarmSnow, albOldSnow,
     &  hNewSnowAge, snowAgTime,
     &  albIceMax, albIceMin, hAlbIce, hAlbSnow,
     &  i0, ksolar,
     &  saltice, S_winton, mu_Tf,
     &  Tf0kel,
     &  himin, Terrmax, nitMaxTsf, hiMax, hsMax,
     &  iceMaskmax, iceMaskmin, himin0,
     &  frac_energy, hihig

      NAMELIST /THSICE_PARM01/
     &     startIceModel, stepFwd_oceMxL,
     &     thSIce_deltaT, ocean_deltaT, tauRelax_MxL,
     &     hMxL_default, sMxL_default, vMxL_default,
     &     stressReduction,
     &     thSIce_taveFreq, thSIce_diagFreq, thSIce_monFreq,
     &     thSIce_tave_mnc, thSIce_snapshot_mnc, thSIce_mon_mnc,
     &     thSIce_pickup_read_mnc, thSIce_pickup_write_mnc,
     &     thSIceFract_InitFile, thSIceThick_InitFile,
     &     thSIceSnowH_InitFile, thSIceSnowA_InitFile,
     &     thSIceEnthp_InitFile, thSIceTsurf_InitFile

      _BEGIN_MASTER(myThid)

      WRITE(msgBuf,'(A)') ' THSICE_READPARMS: opening data.ice'
      CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
     &                    SQUEEZE_RIGHT , 1)

      CALL OPEN_COPY_DATA_FILE(
     I                          'data.ice', 'THSICE_READPARMS',
     O                          iUnit,
     I                          myThid )

C--   Default values (constants)
      rhos     = 330. _d 0
      rhoi     = 900. _d 0
      rhosw    = rhoConst
      rhofw    = rhoConstFresh
      cpice    = 2106. _d 0
      cpwater  = HeatCapacity_Cp
      kice     = 2.03 _d 0
      ksnow    = 0.30 _d 0
      transcoef=0.006 _d 0
      Lfresh   = 3.34 _d 5
      qsnow    = Lfresh
      albColdSnow= 0.85 _d 0
      albWarmSnow= 0.70 _d 0
      albOldSnow = 0.55 _d 0
      albIceMax  = 0.65 _d 0
      albIceMin  = 0.20 _d 0
      hAlbIce    = 0.50 _d 0
      hAlbSnow   = 0.30 _d 0
      hNewSnowAge= 2. _d -3
      snowAgTime = 50. _d 0 * 86400. _d 0
      i0       = 0.3 _d 0
      ksolar   = 1.5 _d 0
      saltice  = 4. _d 0
      S_winton = 1. _d 0
      mu_Tf    = 0.054 _d 0
      Tf0kel   = celsius2K
      himin    = 0.01 _d 0
      Terrmax  = 5.0 _d -1
      nitMaxTsf= 20
      hiMax      = 10. _d 0
      hsMax      = 10. _d 0
      iceMaskmax = 1.  _d 0
      iceMaskmin =  .1 _d 0
      himin0     = 0.2 _d 0
      frac_energy=  .4 _d 0
      hihig      = 2.5 _d 0

C--   Default values (parameters)
      stepFwd_oceMxL  = .FALSE.
      startIceModel   = 0
      thSIce_deltaT   = dTtracerLev(1)
      ocean_deltaT    = dTtracerLev(1)
      tauRelax_MxL    = 0. _d 0
      hMxL_default    = 50. _d 0
      sMxL_default    = 35. _d 0
      vMxL_default    = 5. _d -2
      stressReduction = 1. _d 0
      IF ( useSEAICE ) stressReduction = 0. _d 0
      thSIce_taveFreq = taveFreq
      thSIce_diagFreq = dumpFreq
      thSIce_monFreq  = monitorFreq
#ifdef ALLOW_MNC
      thSIce_tave_mnc     = timeave_mnc
      thSIce_snapshot_mnc = snapshot_mnc
      thSIce_mon_mnc      = monitor_mnc
      thSIce_pickup_read_mnc  = pickup_read_mnc
      thSIce_pickup_write_mnc = pickup_write_mnc
#else
      thSIce_tave_mnc     = .FALSE.
      thSIce_snapshot_mnc = .FALSE.
      thSIce_mon_mnc      = .FALSE.
      thSIce_pickup_read_mnc  = .FALSE.
      thSIce_pickup_write_mnc = .FALSE.
#endif
      thSIceFract_InitFile = ' '
      thSIceThick_InitFile = ' '
      thSIceSnowH_InitFile = ' '
      thSIceSnowA_InitFile = ' '
      thSIceEnthp_InitFile = ' '
      thSIceTsurf_InitFile = ' '


C--   Read parameters from open data file
      READ(UNIT=iUnit,NML=THSICE_CONST)
      WRITE(msgBuf,'(A)') ' THSICE_READPARMS: read THSICE_CONST'
      CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
     &                    SQUEEZE_RIGHT , 1)

      READ(UNIT=iUnit,NML=THSICE_PARM01)
      WRITE(msgBuf,'(A)') ' THSICE_READPARMS: read THSICE_PARM01'
      CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
     &                    SQUEEZE_RIGHT , 1)

C--   Close the open data file
      CLOSE(iUnit)

C-    Define other constants (from previous ones):
      Tmlt1=-mu_Tf*S_winton
      rhoiw = rhosw - rhoi

C     Set I/O parameters
      thSIce_tave_mdsio     = .TRUE.
      thSIce_snapshot_mdsio = .TRUE.
      thSIce_mon_stdio      = .TRUE.
      thSIce_pickup_write_mdsio = .TRUE.
#ifdef ALLOW_MNC
      IF (useMNC) THEN
        IF ( .NOT.outputTypesInclusive
     &       .AND. thSIce_tave_mnc ) thSIce_tave_mdsio = .FALSE.
        IF ( .NOT.outputTypesInclusive
     &       .AND. thSIce_snapshot_mnc )
     &       thSIce_snapshot_mdsio = .FALSE.
        IF ( .NOT.outputTypesInclusive
     &       .AND. thSIce_mon_mnc  ) thSIce_mon_stdio  = .FALSE.
        IF ( .NOT.outputTypesInclusive
     &       .AND. thSIce_pickup_write_mnc  )
     &       thSIce_pickup_write_mdsio = .FALSE.
      ENDIF
#endif

C--   Check parameter consistency:
      IF ( useSEAICE .AND. stressReduction.NE.0. _d 0 ) THEN
C--     If useSEAICE=.true., the stress is computed in seaice_model,
C--     so that it does not need any further reduction
        WRITE(msgBuf,'(2A)')
     &   'THSICE_READPARMS: if useSEAICE, stress will be computed',
     &   ' by SEAICE pkg => no reduction'
        CALL PRINT_MESSAGE( msgBuf, errorMessageUnit,
     &                    SQUEEZE_RIGHT , myThid)
        WRITE(msgBuf,'(A)')
     &   'THSICE_READPARMS: WARNING: reset stressReduction to zero'
        CALL PRINT_MESSAGE( msgBuf, errorMessageUnit,
     &                    SQUEEZE_RIGHT , myThid)
        stressReduction = 0. _d 0
      ENDIF

C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
      iUnit = standardMessageUnit
c     CALL MDSFINDUNIT( iUnit, mythid )
c     OPEN(iUnit,file='thsice_check_params',status='unknown')
      WRITE(iUnit,*) 'ThSI: rhos    =',rhos
      WRITE(iUnit,*) 'ThSI: rhoi    =',rhoi
      WRITE(iUnit,*) 'ThSI: rhosw   =',rhosw
      WRITE(iUnit,*) 'ThSI: rhofw   =',rhofw
      WRITE(iUnit,*) 'ThSI: rhoiw   =',rhoiw
      WRITE(iUnit,*) 'ThSI: cpice   =',cpice
      WRITE(iUnit,*) 'ThSI: cpwater =',cpwater
      WRITE(iUnit,*) 'ThSI: kice    =',kice
      WRITE(iUnit,*) 'ThSI: ksnow   =',ksnow
      WRITE(iUnit,*) 'ThSI: transcoef=',transcoef
      WRITE(iUnit,*) 'ThSI: Lfresh  =',Lfresh
      WRITE(iUnit,*) 'ThSI: qsnow   =',qsnow
      WRITE(iUnit,*) 'ThSI: albColdSnow=',albColdSnow
      WRITE(iUnit,*) 'ThSI: albWarmSnow=',albWarmSnow
      WRITE(iUnit,*) 'ThSI: albOldSnow =',albOldSnow
      WRITE(iUnit,*) 'ThSI: albIceMax =',albIceMax
      WRITE(iUnit,*) 'ThSI: albIceMin =',albIceMin
      WRITE(iUnit,*) 'ThSI: hAlbIce   =',hAlbIce
      WRITE(iUnit,*) 'ThSI: hAlbSnow  =',hAlbSnow
      WRITE(iUnit,*) 'ThSI: hNewSnowAge=',hNewSnowAge
      WRITE(iUnit,*) 'ThSI: snowAgTime =',snowAgTime
      WRITE(iUnit,*) 'ThSI: i0      =',i0
      WRITE(iUnit,*) 'ThSI: ksolar  =',ksolar
      WRITE(iUnit,*) 'ThSI: saltice =',saltice
      WRITE(iUnit,*) 'ThSI: S_winton=',S_winton
      WRITE(iUnit,*) 'ThSI: mu_Tf   =',mu_Tf
      WRITE(iUnit,*) 'ThSI: Tf0kel  =',Tf0kel
      WRITE(iUnit,*) 'ThSI: Tmlt1   =',Tmlt1
      WRITE(iUnit,*) 'ThSI: himin   =',himin
      WRITE(iUnit,*) 'ThSI: Terrmax =',Terrmax
      WRITE(iUnit,*) 'ThSI: nitMaxTsf=',nitMaxTsf
      WRITE(iUnit,*) 'ThSI: hiMax   =',hiMax
      WRITE(iUnit,*) 'ThSI: hsMax   =',hsMax
      WRITE(iUnit,*) 'ThSI: iceMaskmax=',iceMaskmax
      WRITE(iUnit,*) 'ThSI: iceMaskmin=',iceMaskmin
      WRITE(iUnit,*) 'ThSI: himin0  =',himin0
      WRITE(iUnit,*) 'ThSI: frac_energy',frac_energy
      WRITE(iUnit,*) 'ThSI: hihig   =',hihig
      WRITE(iUnit,*) 'ThSI: stressReduction =',stressReduction
      WRITE(iUnit,*) 'ThSI: thSIce_deltaT =',thSIce_deltaT
      WRITE(iUnit,*) 'ThSI: ocean_deltaT  =',ocean_deltaT
      WRITE(iUnit,*) 'ThSI: stepFwd_oceMxL=',stepFwd_oceMxL
      WRITE(iUnit,*) 'ThSI: tauRelax_MxL  =',tauRelax_MxL
      WRITE(iUnit,*) 'ThSI: hMxL_default  =',hMxL_default
      WRITE(iUnit,*) 'ThSI: sMxL_default  =',sMxL_default
      WRITE(iUnit,*) 'ThSI: vMxL_default  =',vMxL_default
      WRITE(iUnit,*) 'ThSI: thSIce_taveFreq=',thSIce_taveFreq
      WRITE(iUnit,*) 'ThSI: thSIce_diagFreq=',thSIce_diagFreq
      WRITE(iUnit,*) 'ThSI: thSIce_monFreq =',thSIce_monFreq
      WRITE(iUnit,*) 'ThSI: startIceModel =',startIceModel
      IF (iUnit.NE.standardMessageUnit) CLOSE(iUnit)
C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|

      _END_MASTER(myThid)

C--   Everyone else must wait for the parameters to be loaded
      _BARRIER

#endif /* ALLOW_THSICE */

      RETURN
      END
1	jmc	1.9	C $Header: /u/gcmpack/MITgcm/pkg/thsice/thsice_readparms.F,v 1.8 2006/05/25 18:03:24 jmc Exp $
2	jmc	1.1	C $Name: $
3
4			#include "THSICE_OPTIONS.h"
5
6	jmc	1.3	CBOP
7			C !ROUTINE: THSICE_READPARMS
8			C !INTERFACE:
9	jmc	1.1	SUBROUTINE THSICE_READPARMS( myThid )
10	jmc	1.3
11			C !DESCRIPTION: \bv
12			C ==========================================================
13			C \| S/R THSICE_READPARMS
14	jmc	1.1	C \| o Routine to initialize THSICE parameters and constants
15	jmc	1.3	C ==========================================================
16	jmc	1.1	C \| Initialize Th-Sea-ICE parameters, read in data.ice
17	jmc	1.3	C ==========================================================
18			C \ev
19
20			C !USES:
21	jmc	1.1	IMPLICIT NONE
22
23			C === Global variables ===
24			#include "SIZE.h"
25			#include "EEPARAMS.h"
26			#include "PARAMS.h"
27			#include "GRID.h"
28			#include "THSICE_PARAMS.h"
29	jmc	1.5	#ifdef ALLOW_MNC
30			#include "MNC_PARAMS.h"
31			#endif
32	jmc	1.1
33	jmc	1.3	C !INPUT/OUTPUT PARAMETERS:
34	jmc	1.1	C === Routine arguments ===
35			INTEGER myThid
36	jmc	1.3	CEOP
37	jmc	1.1
38			#ifdef ALLOW_THSICE
39
40			C === Local variables ===
41	jmc	1.3	C msgBuf - Informational/error message buffer
42	jmc	1.1	C iUnit - Work variable for IO unit number
43			CHARACTER*(MAX_LEN_MBUF) msgBuf
44			INTEGER iUnit
45
46			C-- Th-Sea-ICE parameter
47			NAMELIST /THSICE_CONST/
48			& rhos, rhoi, rhosw, rhofw,
49			& cpice, cpwater,
50			& kice, ksnow,
51			& transcoef, Lfresh, qsnow,
52	jmc	1.8	& albColdSnow, albWarmSnow, albOldSnow,
53			& hNewSnowAge, snowAgTime,
54	jmc	1.2	& albIceMax, albIceMin, hAlbIce, hAlbSnow,
55	jmc	1.1	& i0, ksolar,
56			& saltice, S_winton, mu_Tf,
57			& Tf0kel,
58			& himin, Terrmax, nitMaxTsf, hiMax, hsMax,
59			& iceMaskmax, iceMaskmin, himin0,
60			& frac_energy, hihig
61	jmc	1.8
62	jmc	1.1	NAMELIST /THSICE_PARM01/
63	edhill	1.7	& startIceModel, stepFwd_oceMxL,
64			& thSIce_deltaT, ocean_deltaT, tauRelax_MxL,
65			& hMxL_default, sMxL_default, vMxL_default,
66			& stressReduction,
67			& thSIce_taveFreq, thSIce_diagFreq, thSIce_monFreq,
68			& thSIce_tave_mnc, thSIce_snapshot_mnc, thSIce_mon_mnc,
69			& thSIce_pickup_read_mnc, thSIce_pickup_write_mnc,
70			& thSIceFract_InitFile, thSIceThick_InitFile,
71			& thSIceSnowH_InitFile, thSIceSnowA_InitFile,
72			& thSIceEnthp_InitFile, thSIceTsurf_InitFile
73	jmc	1.1
74			_BEGIN_MASTER(myThid)
75
76			WRITE(msgBuf,'(A)') ' THSICE_READPARMS: opening data.ice'
77			CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
78			& SQUEEZE_RIGHT , 1)
79	jmc	1.8
80	jmc	1.1	CALL OPEN_COPY_DATA_FILE(
81			I 'data.ice', 'THSICE_READPARMS',
82			O iUnit,
83			I myThid )
84
85			C-- Default values (constants)
86			rhos = 330. _d 0
87			rhoi = 900. _d 0
88			rhosw = rhoConst
89			rhofw = rhoConstFresh
90			cpice = 2106. _d 0
91			cpwater = HeatCapacity_Cp
92			kice = 2.03 _d 0
93			ksnow = 0.30 _d 0
94			transcoef=0.006 _d 0
95			Lfresh = 3.34 _d 5
96			qsnow = Lfresh
97	jmc	1.2	albColdSnow= 0.85 _d 0
98			albWarmSnow= 0.70 _d 0
99			albOldSnow = 0.55 _d 0
100			albIceMax = 0.65 _d 0
101			albIceMin = 0.20 _d 0
102			hAlbIce = 0.50 _d 0
103			hAlbSnow = 0.30 _d 0
104			hNewSnowAge= 2. _d -3
105	jmc	1.8	snowAgTime = 50. _d 0 * 86400. _d 0
106	jmc	1.1	i0 = 0.3 _d 0
107			ksolar = 1.5 _d 0
108			saltice = 4. _d 0
109			S_winton = 1. _d 0
110			mu_Tf = 0.054 _d 0
111			Tf0kel = celsius2K
112			himin = 0.01 _d 0
113			Terrmax = 5.0 _d -1
114			nitMaxTsf= 20
115			hiMax = 10. _d 0
116			hsMax = 10. _d 0
117			iceMaskmax = 1. _d 0
118			iceMaskmin = .1 _d 0
119			himin0 = 0.2 _d 0
120			frac_energy= .4 _d 0
121			hihig = 2.5 _d 0
122
123			C-- Default values (parameters)
124	jmc	1.3	stepFwd_oceMxL = .FALSE.
125			startIceModel = 0
126	jmc	1.4	thSIce_deltaT = dTtracerLev(1)
127			ocean_deltaT = dTtracerLev(1)
128	jmc	1.3	tauRelax_MxL = 0. _d 0
129			hMxL_default = 50. _d 0
130			sMxL_default = 35. _d 0
131			vMxL_default = 5. _d -2
132	jmc	1.1	stressReduction = 1. _d 0
133	jmc	1.9	IF ( useSEAICE ) stressReduction = 0. _d 0
134	jmc	1.8	thSIce_taveFreq = taveFreq
135	jmc	1.1	thSIce_diagFreq = dumpFreq
136			thSIce_monFreq = monitorFreq
137	jmc	1.5	#ifdef ALLOW_MNC
138	edhill	1.7	thSIce_tave_mnc = timeave_mnc
139			thSIce_snapshot_mnc = snapshot_mnc
140			thSIce_mon_mnc = monitor_mnc
141			thSIce_pickup_read_mnc = pickup_read_mnc
142			thSIce_pickup_write_mnc = pickup_write_mnc
143	jmc	1.5	#else
144	edhill	1.7	thSIce_tave_mnc = .FALSE.
145			thSIce_snapshot_mnc = .FALSE.
146			thSIce_mon_mnc = .FALSE.
147			thSIce_pickup_read_mnc = .FALSE.
148			thSIce_pickup_write_mnc = .FALSE.
149	jmc	1.5	#endif
150			thSIceFract_InitFile = ' '
151			thSIceThick_InitFile = ' '
152			thSIceSnowH_InitFile = ' '
153			thSIceSnowA_InitFile = ' '
154			thSIceEnthp_InitFile = ' '
155			thSIceTsurf_InitFile = ' '
156
157	jmc	1.1
158			C-- Read parameters from open data file
159			READ(UNIT=iUnit,NML=THSICE_CONST)
160			WRITE(msgBuf,'(A)') ' THSICE_READPARMS: read THSICE_CONST'
161			CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
162			& SQUEEZE_RIGHT , 1)
163
164			READ(UNIT=iUnit,NML=THSICE_PARM01)
165			WRITE(msgBuf,'(A)') ' THSICE_READPARMS: read THSICE_PARM01'
166			CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
167			& SQUEEZE_RIGHT , 1)
168
169			C-- Close the open data file
170			CLOSE(iUnit)
171
172			C- Define other constants (from previous ones):
173			Tmlt1=-mu_Tf*S_winton
174			rhoiw = rhosw - rhoi
175
176	edhill	1.7	C Set I/O parameters
177			thSIce_tave_mdsio = .TRUE.
178			thSIce_snapshot_mdsio = .TRUE.
179			thSIce_mon_stdio = .TRUE.
180			thSIce_pickup_write_mdsio = .TRUE.
181	jmc	1.5	#ifdef ALLOW_MNC
182			IF (useMNC) THEN
183			IF ( .NOT.outputTypesInclusive
184			& .AND. thSIce_tave_mnc ) thSIce_tave_mdsio = .FALSE.
185	jmc	1.8	IF ( .NOT.outputTypesInclusive
186			& .AND. thSIce_snapshot_mnc )
187	edhill	1.7	& thSIce_snapshot_mdsio = .FALSE.
188	jmc	1.5	IF ( .NOT.outputTypesInclusive
189			& .AND. thSIce_mon_mnc ) thSIce_mon_stdio = .FALSE.
190	edhill	1.7	IF ( .NOT.outputTypesInclusive
191	jmc	1.8	& .AND. thSIce_pickup_write_mnc )
192	edhill	1.7	& thSIce_pickup_write_mdsio = .FALSE.
193	jmc	1.5	ENDIF
194			#endif
195
196	jmc	1.9	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	jmc	1.1	C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-\|--+----\|
213			iUnit = standardMessageUnit
214	jmc	1.5	c CALL MDSFINDUNIT( iUnit, mythid )
215	jmc	1.1	c OPEN(iUnit,file='thsice_check_params',status='unknown')
216			WRITE(iUnit,*) 'ThSI: rhos =',rhos
217			WRITE(iUnit,*) 'ThSI: rhoi =',rhoi
218			WRITE(iUnit,*) 'ThSI: rhosw =',rhosw
219			WRITE(iUnit,*) 'ThSI: rhofw =',rhofw
220			WRITE(iUnit,*) 'ThSI: rhoiw =',rhoiw
221			WRITE(iUnit,*) 'ThSI: cpice =',cpice
222			WRITE(iUnit,*) 'ThSI: cpwater =',cpwater
223			WRITE(iUnit,*) 'ThSI: kice =',kice
224			WRITE(iUnit,*) 'ThSI: ksnow =',ksnow
225			WRITE(iUnit,*) 'ThSI: transcoef=',transcoef
226			WRITE(iUnit,*) 'ThSI: Lfresh =',Lfresh
227			WRITE(iUnit,*) 'ThSI: qsnow =',qsnow
228	jmc	1.2	WRITE(iUnit,*) 'ThSI: albColdSnow=',albColdSnow
229			WRITE(iUnit,*) 'ThSI: albWarmSnow=',albWarmSnow
230			WRITE(iUnit,*) 'ThSI: albOldSnow =',albOldSnow
231			WRITE(iUnit,*) 'ThSI: albIceMax =',albIceMax
232			WRITE(iUnit,*) 'ThSI: albIceMin =',albIceMin
233			WRITE(iUnit,*) 'ThSI: hAlbIce =',hAlbIce
234			WRITE(iUnit,*) 'ThSI: hAlbSnow =',hAlbSnow
235	jmc	1.8	WRITE(iUnit,*) 'ThSI: hNewSnowAge=',hNewSnowAge
236			WRITE(iUnit,*) 'ThSI: snowAgTime =',snowAgTime
237	jmc	1.1	WRITE(iUnit,*) 'ThSI: i0 =',i0
238			WRITE(iUnit,*) 'ThSI: ksolar =',ksolar
239			WRITE(iUnit,*) 'ThSI: saltice =',saltice
240			WRITE(iUnit,*) 'ThSI: S_winton=',S_winton
241			WRITE(iUnit,*) 'ThSI: mu_Tf =',mu_Tf
242			WRITE(iUnit,*) 'ThSI: Tf0kel =',Tf0kel
243			WRITE(iUnit,*) 'ThSI: Tmlt1 =',Tmlt1
244			WRITE(iUnit,*) 'ThSI: himin =',himin
245			WRITE(iUnit,*) 'ThSI: Terrmax =',Terrmax
246			WRITE(iUnit,*) 'ThSI: nitMaxTsf=',nitMaxTsf
247			WRITE(iUnit,*) 'ThSI: hiMax =',hiMax
248			WRITE(iUnit,*) 'ThSI: hsMax =',hsMax
249			WRITE(iUnit,*) 'ThSI: iceMaskmax=',iceMaskmax
250			WRITE(iUnit,*) 'ThSI: iceMaskmin=',iceMaskmin
251			WRITE(iUnit,*) 'ThSI: himin0 =',himin0
252			WRITE(iUnit,*) 'ThSI: frac_energy',frac_energy
253			WRITE(iUnit,*) 'ThSI: hihig =',hihig
254			WRITE(iUnit,*) 'ThSI: stressReduction =',stressReduction
255	jmc	1.3	WRITE(iUnit,*) 'ThSI: thSIce_deltaT =',thSIce_deltaT
256			WRITE(iUnit,*) 'ThSI: ocean_deltaT =',ocean_deltaT
257			WRITE(iUnit,*) 'ThSI: stepFwd_oceMxL=',stepFwd_oceMxL
258			WRITE(iUnit,*) 'ThSI: tauRelax_MxL =',tauRelax_MxL
259			WRITE(iUnit,*) 'ThSI: hMxL_default =',hMxL_default
260			WRITE(iUnit,*) 'ThSI: sMxL_default =',sMxL_default
261			WRITE(iUnit,*) 'ThSI: vMxL_default =',vMxL_default
262	jmc	1.1	WRITE(iUnit,*) 'ThSI: thSIce_taveFreq=',thSIce_taveFreq
263			WRITE(iUnit,*) 'ThSI: thSIce_diagFreq=',thSIce_diagFreq
264			WRITE(iUnit,*) 'ThSI: thSIce_monFreq =',thSIce_monFreq
265			WRITE(iUnit,*) 'ThSI: startIceModel =',startIceModel
266	jmc	1.5	IF (iUnit.NE.standardMessageUnit) CLOSE(iUnit)
267	jmc	1.1	C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-\|--+----\|
268
269			_END_MASTER(myThid)
270
271			C-- Everyone else must wait for the parameters to be loaded
272			_BARRIER
273
274			#endif /* ALLOW_THSICE */
275
276			RETURN
277			END