pkg/thsice/THSICE_PARAMS.h

C $Header: /u/gcmpack/MITgcm/pkg/thsice/THSICE_PARAMS.h,v 1.10 2007/04/16 22:38:24 heimbach Exp $
C $Name:  $

#ifdef ALLOW_THSICE
C     *==========================================================*
C     | THSICE_PARAMS.h
C     | o Header file for Therm_SeaIce package parameters:
C     |   - basic parameter ( I/O frequency, etc ...)
C     |   - physical constants (used in therm_SeaIce pkg)
C     *==========================================================*

C----------------------------------------------------------------------------
C.. Common blocks for almost everything that the sea ice model passes around.
C----------------------------------------------------------------------------

C--   COMMON / THSICE_PHYSPAR_R / physical (real) parameter
C.. densities
C     rhos      ::   density of snow [kg/m^3]
C     rhoi      ::   density of ice [kg/m^3]
C     rhosw     ::   density of seawater [kg/m^3]
C     rhofw     ::   density of fresh water [kg/m^3]
C     rhoiw     ::   ocean-ice density difference [kg/m^3]
C.. specific heats
C     cpIce     ::   specific heat of fresh ice [J/kg/K]
C     cpWater   ::   specific heat of water [J/kg/K]
C .. thermal conductivity. QQ check units
C     kIce      ::   thermal conductivity of pure ice [W/m/K]
C     kSnow     ::   thermal conductivity of snow [W/m/K]
C .. heat transfer coefficient
C     bMeltCoef ::   base-melting heat transfer coefficient
C                    (between ice & water) [no unit]
C .. latent heat
C     Lfresh    ::   latent heat of melting of pure ice [J/kg]
C .. Enthalpy
C     qsnow     ::   snow enthalpy [J/kg]
C .. Albedo
C     albColdSnow :: albedo of cold (=dry) new snow (Tsfc < tempSnowAlb)
C     albWarmSnow :: albedo of warm (=wet) new snow (Tsfc = 0)
C     tempSnowAlb :: temperature transition from ColdSnow to WarmSnow Alb. [oC]
C     albOldSnow  :: albedo of old snow (snowAge > 35.d)
C     albIceMax   :: max albedo of bare ice (thick ice)
C     albIceMin   :: minimum ice albedo (very thin ice)
C     hAlbIce     :: ice thickness for albedo transition: thin/thick ice albedo
C     hAlbSnow    :: snow thickness for albedo transition: snow/ice albedo
C     hNewSnowAge :: new snow thickness that refresh the snow-age (by 1/e)
C     snowAgTime  :: snow aging time scale [s]
C .. Solar parameters
C     i0swFrac  ::   fraction of penetrating solar rad
C     ksolar    ::   bulk solar abs coeff of sea ice [m^-1]
C     dhSnowLin ::   half slope of linear distribution of snow thickness within
C                    the grid-cell (from hSnow-dhSnow to hSnow+dhSnow, if full
C                    ice & snow cover) [m] ; (only used for SW radiation).
C .. Salinity
C     saltIce   ::   salinity of ice [g/kg]
C     S_winton  ::   Winton salinity of ice [g/kg]
C .. melting
C     Tf0kel    ::   Freezing temp of fresh water in Kelvin = 273.15
C     mu_Tf     ::   linear dependance of melting temperature on Salinity [oC/psu]
C                     Tf(sea-water) = -mu_Tf * S
C     Tmlt1     ::   Winton melting temperature: Tmlt1 = -mu_Tf * S_winton
C     Terrmax   ::   Temperature convergence criteria [oC]
C .. Min/Max
C     hIceMin   ::   Minimum ice  thickness [m]
C     hiMax     ::   Maximum ice  thickness [m]
C     hsMax     ::   Maximum snow thickness [m]
C .. for fractional ice
C     iceMaskMax  :: maximum Ice fraction (=1 for no fractional ice)
C     iceMaskMin  :: mimimum Ice fraction (=1 for no fractional ice)
C     fracEnFreez :: fraction of energy going to lateral freezing (vs height increase)
C     fracEnMelt  :: fraction of energy going to lateral melting (vs height decrease)
C                                         (=0 for no fract. ice)
C     hThinIce    :: ice height above which fracEnMelt/Freez are applied [m]
C                                         (=hIceMin for no fractional ice)
C     hThickIce   :: ice height below which fracEnMelt/Freez are applied [m]
C                                         (=large for no fractional ice)
C     hNewIceMax  :: new ice maximum thickness [m]
C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|

      COMMON / THSICE_PHYSPAR_R /
     &  rhos, rhoi, rhosw, rhofw, rhoiw,
     &  cpIce, cpWater,
     &  kIce, kSnow,
     &  bMeltCoef, Lfresh, qsnow,
     &  albColdSnow, albWarmSnow, tempSnowAlb,
     &  albOldSnow, hNewSnowAge, snowAgTime,
     &  albIceMax, albIceMin, hAlbIce, hAlbSnow,
     &  i0swFrac, ksolar, dhSnowLin,
     &  saltIce, S_winton, mu_Tf,
     &  Tf0kel, Tmlt1, Terrmax,
     &  hIceMin, hiMax, hsMax,
     &  iceMaskMax, iceMaskMin,
     &  fracEnMelt, fracEnFreez,
     &  hThinIce, hThickIce, hNewIceMax

      _RL  rhos
      _RL  rhoi
      _RL  rhosw
      _RL  rhofw
      _RL  rhoiw
      _RL  cpIce
      _RL  cpWater
      _RL  kIce
      _RL  kSnow
      _RL  bMeltCoef
      _RL  Lfresh
      _RL  qsnow
      _RL  albColdSnow
      _RL  albWarmSnow
      _RL  tempSnowAlb
      _RL  albOldSnow
      _RL  hNewSnowAge
      _RL  snowAgTime
      _RL  albIceMax
      _RL  albIceMin
      _RL  hAlbIce
      _RL  hAlbSnow
      _RL  i0swFrac
      _RL  ksolar
      _RL  dhSnowLin
      _RL  saltIce
      _RL  S_winton
      _RL  mu_Tf
      _RL  Tf0kel
      _RL  Tmlt1
      _RL  Terrmax
      _RL  hIceMin
      _RL  hiMax
      _RL  hsMax
      _RL iceMaskMax
      _RL iceMaskMin
      _RL fracEnMelt
      _RL fracEnFreez
      _RL hThinIce
      _RL hThickIce
      _RL hNewIceMax

C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|

C--   COMMON / THSICE_PAR_L / ice model (logical) parameters
C     stepFwd_oceMxL        :: step forward mixed-layer T & S (slab-ocean)
C     thSIce_tave_mdsio     :: write TimeAverage output using MDSIO
C     thSIce_snapshot_mdsio :: write snap-shot output   using MDSIO
C     thSIce_mon_stdio      :: write monitor to std-outp
C     thSIce_tave_mnc       :: write TimeAverage output using MNC
C     thSIce_snapshot_mnc   :: write snap-shot output   using MNC
C     thSIce_mon_mnc        :: write monitor to netcdf file
C     thSIce_pickup_read_mnc    :: pickup read w/ MNC
C     thSIce_pickup_write_mnc   :: pickup write w/ MNC
C     thSIce_pickup_write_mdsio :: pickup write w/ MDSIO
      COMMON / THSICE_PAR_L /
     &     stepFwd_oceMxL,
     &     thSIce_tave_mdsio, thSIce_snapshot_mdsio, thSIce_mon_stdio,
     &     thSIce_tave_mnc,   thSIce_snapshot_mnc,   thSIce_mon_mnc,
     &     thSIce_pickup_read_mnc,
     &     thSIce_pickup_write_mdsio,
     &     thSIce_pickup_write_mnc

      LOGICAL stepFwd_oceMxL
      LOGICAL thSIce_tave_mdsio, thSIce_snapshot_mdsio, thSIce_mon_stdio
      LOGICAL thSIce_tave_mnc,   thSIce_snapshot_mnc,   thSIce_mon_mnc
      LOGICAL thSIce_pickup_read_mnc
      LOGICAL thSIce_pickup_write_mdsio
      LOGICAL thSIce_pickup_write_mnc

C--   COMMON / THSICE_PAR_I / ice model (integer) parameters
C     startIceModel   :: =1 : start ice model at nIter0 ; =0 : use pickup files
C                     :: -1 : start from a small pickup (without Mix.Layer)
C     nitMaxTsf       :: maximum Nb of iter to find Surface Temp (Trsf)
C     thSIceAdvScheme :: thSIce Advection scheme selector
      COMMON / THSICE_PAR_I /
     &  startIceModel, nitMaxTsf, thSIceAdvScheme

      INTEGER startIceModel
      INTEGER nitMaxTsf
      INTEGER thSIceAdvScheme

C--   COMMON / THSICE_PAR_R / ice model (real) parameters
C     thSIce_deltaT   :: ice model time-step, seaice thicken/extend [s]
C     thSIce_dtTemp   :: ice model time-step, solve4temp [s]
C     ocean_deltaT    :: ocean mixed-layer time-step [s]
C     tauRelax_MxL    :: Relaxation time scale for MixLayer T & S [s]
C     hMxL_default    :: default value for ocean MixLayer thickness [m]
C     sMxL_default    :: default value for salinity in MixLayer [psu]
C     vMxL_default    :: default value for ocean current velocity in MxL [m/s]
C     thSIce_diffK    :: thickness (horizontal) diffusivity [m^2/s]
C     stressReduction :: reduction factor for wind-stress under sea-ice [0-1]
C     thSIce_taveFreq :: Frequency^-1 for time-Aver. output [s]
C     thSIce_diagFreq :: Frequency^-1 for diagnostic output [s]
C     thSIce_monFreq  :: Frequency^-1 for monitor    output [s]
      COMMON / THSICE_PAR_R /
     &  thSIce_deltaT,  thSIce_dtTemp, ocean_deltaT, tauRelax_MxL,
     &  hMxL_default,  sMxL_default, vMxL_default,
     &  thSIce_diffK,  stressReduction,
     &  thSIce_taveFreq, thSIce_diagFreq, thSIce_monFreq

      _RL  thSIce_deltaT, thSIce_dtTemp, ocean_deltaT, tauRelax_MxL
      _RL  hMxL_default, sMxL_default, vMxL_default
      _RL  thSIce_diffK,  stressReduction
      _RL  thSIce_taveFreq, thSIce_diagFreq, thSIce_monFreq

C--   COMMON / THSICE_PAR_C / ice model (character) parameters
C     thSIceFract_InitFile :: File name for initial ice fraction
C     thSIceThick_InitFile :: File name for initial ice thickness
C     thSIceSnowH_InitFile :: File name for initial snow thickness
C     thSIceSnowA_InitFile :: File name for initial snow Age
C     thSIceEnthp_InitFile :: File name for initial ice enthalpy
C     thSIceTsurf_InitFile :: File name for initial surf. temp
      COMMON / THSICE_PAR_C /
     &  thSIceFract_InitFile,
     &  thSIceThick_InitFile,
     &  thSIceSnowH_InitFile,
     &  thSIceSnowA_InitFile,
     &  thSIceEnthp_InitFile,
     &  thSIceTsurf_InitFile
      CHARACTER*(MAX_LEN_FNAM) thSIceFract_InitFile
      CHARACTER*(MAX_LEN_FNAM) thSIceThick_InitFile
      CHARACTER*(MAX_LEN_FNAM) thSIceSnowH_InitFile
      CHARACTER*(MAX_LEN_FNAM) thSIceSnowA_InitFile
      CHARACTER*(MAX_LEN_FNAM) thSIceEnthp_InitFile
      CHARACTER*(MAX_LEN_FNAM) thSIceTsurf_InitFile

#ifdef ALLOW_AUTODIFF_TAMC
      INTEGER iicekey
      INTEGER ikey_1, ikey_2, ikey_3, ikey_4
      INTEGER MaxTsf
      PARAMETER ( MaxTsf = 20 )
#endif

C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|

#endif /* ALLOW_THSICE */

CEH3 ;;; Local Variables: ***
CEH3 ;;; mode:fortran ***
CEH3 ;;; End: ***
1	jmc	1.11	C $Header: /u/gcmpack/MITgcm/pkg/thsice/THSICE_PARAMS.h,v 1.10 2007/04/16 22:38:24 heimbach Exp $
2	jmc	1.1	C $Name: $
3
4			#ifdef ALLOW_THSICE
5			C ==========================================================
6			C \| THSICE_PARAMS.h
7			C \| o Header file for Therm_SeaIce package parameters:
8			C \| - basic parameter ( I/O frequency, etc ...)
9			C \| - physical constants (used in therm_SeaIce pkg)
10			C ==========================================================
11
12			C----------------------------------------------------------------------------
13			C.. Common blocks for almost everything that the sea ice model passes around.
14			C----------------------------------------------------------------------------
15
16			C-- COMMON / THSICE_PHYSPAR_R / physical (real) parameter
17			C.. densities
18	jmc	1.11	C rhos :: density of snow [kg/m^3]
19			C rhoi :: density of ice [kg/m^3]
20			C rhosw :: density of seawater [kg/m^3]
21			C rhofw :: density of fresh water [kg/m^3]
22			C rhoiw :: ocean-ice density difference [kg/m^3]
23	jmc	1.1	C.. specific heats
24	jmc	1.11	C cpIce :: specific heat of fresh ice [J/kg/K]
25			C cpWater :: specific heat of water [J/kg/K]
26	jmc	1.1	C .. thermal conductivity. QQ check units
27	jmc	1.11	C kIce :: thermal conductivity of pure ice [W/m/K]
28			C kSnow :: thermal conductivity of snow [W/m/K]
29	jmc	1.1	C .. heat transfer coefficient
30	jmc	1.11	C bMeltCoef :: base-melting heat transfer coefficient
31			C (between ice & water) [no unit]
32	jmc	1.1	C .. latent heat
33	jmc	1.11	C Lfresh :: latent heat of melting of pure ice [J/kg]
34	jmc	1.2	C .. Enthalpy
35	jmc	1.11	C qsnow :: snow enthalpy [J/kg]
36	jmc	1.1	C .. Albedo
37	jmc	1.7	C albColdSnow :: albedo of cold (=dry) new snow (Tsfc < tempSnowAlb)
38	jmc	1.2	C albWarmSnow :: albedo of warm (=wet) new snow (Tsfc = 0)
39	jmc	1.11	C tempSnowAlb :: temperature transition from ColdSnow to WarmSnow Alb. [oC]
40	jmc	1.2	C albOldSnow :: albedo of old snow (snowAge > 35.d)
41			C albIceMax :: max albedo of bare ice (thick ice)
42			C albIceMin :: minimum ice albedo (very thin ice)
43			C hAlbIce :: ice thickness for albedo transition: thin/thick ice albedo
44			C hAlbSnow :: snow thickness for albedo transition: snow/ice albedo
45			C hNewSnowAge :: new snow thickness that refresh the snow-age (by 1/e)
46	jmc	1.11	C snowAgTime :: snow aging time scale [s]
47	jmc	1.2	C .. Solar parameters
48	jmc	1.11	C i0swFrac :: fraction of penetrating solar rad
49			C ksolar :: bulk solar abs coeff of sea ice [m^-1]
50			C dhSnowLin :: half slope of linear distribution of snow thickness within
51			C the grid-cell (from hSnow-dhSnow to hSnow+dhSnow, if full
52			C ice & snow cover) [m] ; (only used for SW radiation).
53	jmc	1.1	C .. Salinity
54	jmc	1.11	C saltIce :: salinity of ice [g/kg]
55			C S_winton :: Winton salinity of ice [g/kg]
56	jmc	1.1	C .. melting
57	jmc	1.11	C Tf0kel :: Freezing temp of fresh water in Kelvin = 273.15
58			C mu_Tf :: linear dependance of melting temperature on Salinity [oC/psu]
59			C Tf(sea-water) = -mu_Tf * S
60			C Tmlt1 :: Winton melting temperature: Tmlt1 = -mu_Tf * S_winton
61			C Terrmax :: Temperature convergence criteria [oC]
62	jmc	1.1	C .. Min/Max
63	jmc	1.11	C hIceMin :: Minimum ice thickness [m]
64			C hiMax :: Maximum ice thickness [m]
65			C hsMax :: Maximum snow thickness [m]
66	jmc	1.1	C .. for fractional ice
67	jmc	1.11	C iceMaskMax :: maximum Ice fraction (=1 for no fractional ice)
68			C iceMaskMin :: mimimum Ice fraction (=1 for no fractional ice)
69			C fracEnFreez :: fraction of energy going to lateral freezing (vs height increase)
70			C fracEnMelt :: fraction of energy going to lateral melting (vs height decrease)
71			C (=0 for no fract. ice)
72			C hThinIce :: ice height above which fracEnMelt/Freez are applied [m]
73			C (=hIceMin for no fractional ice)
74			C hThickIce :: ice height below which fracEnMelt/Freez are applied [m]
75			C (=large for no fractional ice)
76			C hNewIceMax :: new ice maximum thickness [m]
77	jmc	1.1	C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-\|--+----\|
78
79	jmc	1.6	COMMON / THSICE_PHYSPAR_R /
80	jmc	1.1	& rhos, rhoi, rhosw, rhofw, rhoiw,
81	jmc	1.11	& cpIce, cpWater,
82			& kIce, kSnow,
83			& bMeltCoef, Lfresh, qsnow,
84	jmc	1.7	& albColdSnow, albWarmSnow, tempSnowAlb,
85			& albOldSnow, hNewSnowAge, snowAgTime,
86	jmc	1.2	& albIceMax, albIceMin, hAlbIce, hAlbSnow,
87	jmc	1.11	& i0swFrac, ksolar, dhSnowLin,
88			& saltIce, S_winton, mu_Tf,
89			& Tf0kel, Tmlt1, Terrmax,
90			& hIceMin, hiMax, hsMax,
91			& iceMaskMax, iceMaskMin,
92			& fracEnMelt, fracEnFreez,
93			& hThinIce, hThickIce, hNewIceMax
94	jmc	1.1
95			_RL rhos
96			_RL rhoi
97			_RL rhosw
98			_RL rhofw
99			_RL rhoiw
100	jmc	1.11	_RL cpIce
101			_RL cpWater
102			_RL kIce
103			_RL kSnow
104			_RL bMeltCoef
105	jmc	1.1	_RL Lfresh
106			_RL qsnow
107	jmc	1.2	_RL albColdSnow
108			_RL albWarmSnow
109	jmc	1.7	_RL tempSnowAlb
110	jmc	1.2	_RL albOldSnow
111	jmc	1.7	_RL hNewSnowAge
112			_RL snowAgTime
113	jmc	1.2	_RL albIceMax
114			_RL albIceMin
115			_RL hAlbIce
116			_RL hAlbSnow
117	jmc	1.11	_RL i0swFrac
118	jmc	1.1	_RL ksolar
119	jmc	1.11	_RL dhSnowLin
120			_RL saltIce
121	jmc	1.1	_RL S_winton
122			_RL mu_Tf
123			_RL Tf0kel
124			_RL Tmlt1
125			_RL Terrmax
126	jmc	1.11	_RL hIceMin
127	jmc	1.1	_RL hiMax
128			_RL hsMax
129	jmc	1.11	_RL iceMaskMax
130			_RL iceMaskMin
131			_RL fracEnMelt
132			_RL fracEnFreez
133			_RL hThinIce
134			_RL hThickIce
135			_RL hNewIceMax
136	jmc	1.1
137			C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-\|--+----\|
138
139	jmc	1.3	C-- COMMON / THSICE_PAR_L / ice model (logical) parameters
140	edhill	1.5	C stepFwd_oceMxL :: step forward mixed-layer T & S (slab-ocean)
141			C thSIce_tave_mdsio :: write TimeAverage output using MDSIO
142			C thSIce_snapshot_mdsio :: write snap-shot output using MDSIO
143			C thSIce_mon_stdio :: write monitor to std-outp
144			C thSIce_tave_mnc :: write TimeAverage output using MNC
145			C thSIce_snapshot_mnc :: write snap-shot output using MNC
146			C thSIce_mon_mnc :: write monitor to netcdf file
147			C thSIce_pickup_read_mnc :: pickup read w/ MNC
148			C thSIce_pickup_write_mnc :: pickup write w/ MNC
149			C thSIce_pickup_write_mdsio :: pickup write w/ MDSIO
150	jmc	1.6	COMMON / THSICE_PAR_L /
151	edhill	1.5	& stepFwd_oceMxL,
152			& thSIce_tave_mdsio, thSIce_snapshot_mdsio, thSIce_mon_stdio,
153			& thSIce_tave_mnc, thSIce_snapshot_mnc, thSIce_mon_mnc,
154	jmc	1.6	& thSIce_pickup_read_mnc,
155	edhill	1.5	& thSIce_pickup_write_mdsio,
156			& thSIce_pickup_write_mnc
157	jmc	1.3
158			LOGICAL stepFwd_oceMxL
159	edhill	1.5	LOGICAL thSIce_tave_mdsio, thSIce_snapshot_mdsio, thSIce_mon_stdio
160			LOGICAL thSIce_tave_mnc, thSIce_snapshot_mnc, thSIce_mon_mnc
161			LOGICAL thSIce_pickup_read_mnc
162			LOGICAL thSIce_pickup_write_mdsio
163			LOGICAL thSIce_pickup_write_mnc
164	jmc	1.3
165	jmc	1.1	C-- COMMON / THSICE_PAR_I / ice model (integer) parameters
166	jmc	1.8	C startIceModel :: =1 : start ice model at nIter0 ; =0 : use pickup files
167			C :: -1 : start from a small pickup (without Mix.Layer)
168			C nitMaxTsf :: maximum Nb of iter to find Surface Temp (Trsf)
169			C thSIceAdvScheme :: thSIce Advection scheme selector
170	jmc	1.6	COMMON / THSICE_PAR_I /
171	jmc	1.8	& startIceModel, nitMaxTsf, thSIceAdvScheme
172	jmc	1.1
173			INTEGER startIceModel
174			INTEGER nitMaxTsf
175	jmc	1.8	INTEGER thSIceAdvScheme
176	jmc	1.1
177			C-- COMMON / THSICE_PAR_R / ice model (real) parameters
178	jscott	1.9	C thSIce_deltaT :: ice model time-step, seaice thicken/extend [s]
179			C thSIce_dtTemp :: ice model time-step, solve4temp [s]
180	jmc	1.8	C ocean_deltaT :: ocean mixed-layer time-step [s]
181	jmc	1.3	C tauRelax_MxL :: Relaxation time scale for MixLayer T & S [s]
182			C hMxL_default :: default value for ocean MixLayer thickness [m]
183			C sMxL_default :: default value for salinity in MixLayer [psu]
184			C vMxL_default :: default value for ocean current velocity in MxL [m/s]
185	jmc	1.8	C thSIce_diffK :: thickness (horizontal) diffusivity [m^2/s]
186	jmc	1.3	C stressReduction :: reduction factor for wind-stress under sea-ice [0-1]
187			C thSIce_taveFreq :: Frequency^-1 for time-Aver. output [s]
188			C thSIce_diagFreq :: Frequency^-1 for diagnostic output [s]
189			C thSIce_monFreq :: Frequency^-1 for monitor output [s]
190	jmc	1.1	COMMON / THSICE_PAR_R /
191	jscott	1.9	& thSIce_deltaT, thSIce_dtTemp, ocean_deltaT, tauRelax_MxL,
192	jmc	1.3	& hMxL_default, sMxL_default, vMxL_default,
193	jmc	1.8	& thSIce_diffK, stressReduction,
194	jmc	1.1	& thSIce_taveFreq, thSIce_diagFreq, thSIce_monFreq
195
196	jscott	1.9	_RL thSIce_deltaT, thSIce_dtTemp, ocean_deltaT, tauRelax_MxL
197	jmc	1.6	_RL hMxL_default, sMxL_default, vMxL_default
198	jmc	1.8	_RL thSIce_diffK, stressReduction
199	jmc	1.1	_RL thSIce_taveFreq, thSIce_diagFreq, thSIce_monFreq
200
201	jmc	1.4	C-- COMMON / THSICE_PAR_C / ice model (character) parameters
202			C thSIceFract_InitFile :: File name for initial ice fraction
203			C thSIceThick_InitFile :: File name for initial ice thickness
204			C thSIceSnowH_InitFile :: File name for initial snow thickness
205			C thSIceSnowA_InitFile :: File name for initial snow Age
206			C thSIceEnthp_InitFile :: File name for initial ice enthalpy
207			C thSIceTsurf_InitFile :: File name for initial surf. temp
208			COMMON / THSICE_PAR_C /
209			& thSIceFract_InitFile,
210			& thSIceThick_InitFile,
211			& thSIceSnowH_InitFile,
212			& thSIceSnowA_InitFile,
213			& thSIceEnthp_InitFile,
214	jmc	1.6	& thSIceTsurf_InitFile
215	jmc	1.4	CHARACTER*(MAX_LEN_FNAM) thSIceFract_InitFile
216			CHARACTER*(MAX_LEN_FNAM) thSIceThick_InitFile
217			CHARACTER*(MAX_LEN_FNAM) thSIceSnowH_InitFile
218			CHARACTER*(MAX_LEN_FNAM) thSIceSnowA_InitFile
219			CHARACTER*(MAX_LEN_FNAM) thSIceEnthp_InitFile
220			CHARACTER*(MAX_LEN_FNAM) thSIceTsurf_InitFile
221
222	heimbach	1.10	#ifdef ALLOW_AUTODIFF_TAMC
223			INTEGER iicekey
224			INTEGER ikey_1, ikey_2, ikey_3, ikey_4
225			INTEGER MaxTsf
226			PARAMETER ( MaxTsf = 20 )
227			#endif
228
229	jmc	1.1	C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-\|--+----\|
230
231			#endif /* ALLOW_THSICE */
232	edhill	1.5
233			CEH3 ;;; Local Variables: ***
234			CEH3 ;;; mode:fortran ***
235			CEH3 ;;; End: ***