pkg/thsice/THSICE_PARAMS.h

C $Header: /u/gcmpack/MITgcm/pkg/thsice/THSICE_PARAMS.h,v 1.8 2007/04/04 02:11:34 jmc 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     transcoef ::   transfer coef between ice and water (unitless)
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     i0        ::   fraction of penetrating solar rad
C     ksolar    ::   bulk solar abs coeff of sea ice (m-1)
C .. Salinity
C     saltice   ::   salinity of ice (o/oo)
C     S_winton  ::   winton salinity of ice (o/oo)
C     mu_Tf     ::   Tf:brine salinity ratio (C/ppt)
C .. melting
C     Tf0kel    ::   Freezing temp of fresh ice in Kelvin = 273.15
C     Tmlt1     ::    melting temp; depends on S (C)
C .. Min/Max
C     himin     ::   minimum thickness for ice (m)
C     Terrmax   ::   temperature convergence precision (C)
C     hiMax     ::   Maximum thickness for ice  (m)
C     hsMax     ::   Maximum thickness for snow (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     himin0     :: minimum ice height   (=himin for no fractional ice)
C     frac_energy:: fraction of energy going to melting (=0 for no fract. ice)
C     hihig      :: ice height above which freezing only occurs over open ocean
C                                        (=large for no fractional ice)
C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|

      COMMON / THSICE_PHYSPAR_R /
     &  rhos, rhoi, rhosw, rhofw, rhoiw,
     &  cpice, cpwater,
     &  kice, ksnow,
     &  transcoef, Lfresh, qsnow,
     &  albColdSnow, albWarmSnow, tempSnowAlb,
     &  albOldSnow, hNewSnowAge, snowAgTime,
     &  albIceMax, albIceMin, hAlbIce, hAlbSnow,
     &  i0, ksolar,
     &  saltice, S_winton, mu_Tf,
     &  Tf0kel, Tmlt1,
     &  himin, Terrmax, hiMax, hsMax,
     &  iceMaskmax, iceMaskmin, himin0,
     &  frac_energy, hihig

      _RL  rhos
      _RL  rhoi
      _RL  rhosw
      _RL  rhofw
      _RL  rhoiw
      _RL  cpice
      _RL  cpwater
      _RL  kice
      _RL  ksnow
      _RL  transcoef
      _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  i0
      _RL  ksolar
      _RL  saltice
      _RL  S_winton
      _RL  mu_Tf
      _RL  Tf0kel
      _RL  Tmlt1
      _RL  himin
      _RL  Terrmax
      _RL  hiMax
      _RL  hsMax
      _RL iceMaskmax
      _RL iceMaskmin
      _RL himin0
      _RL frac_energy
      _RL hihig

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

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

#endif /* ALLOW_THSICE */

CEH3 ;;; Local Variables: ***
CEH3 ;;; mode:fortran ***
CEH3 ;;; End: ***
1	jscott	1.9	C $Header: /u/gcmpack/MITgcm/pkg/thsice/THSICE_PARAMS.h,v 1.8 2007/04/04 02:11:34 jmc 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			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			C.. specific heats
24			C cpice :: specific heat of fresh ice (J/kg/K)
25			C cpwater :: specific heat of water (J/kg/K)
26			C .. thermal conductivity. QQ check units
27			C kice :: thermal conductivity of pure ice (W/m/K)
28			C ksnow :: thermal conductivity of snow (W/m/K)
29			C .. heat transfer coefficient
30			C transcoef :: transfer coef between ice and water (unitless)
31			C .. latent heat
32			C Lfresh :: latent heat of melting of pure ice (J/kg)
33	jmc	1.2	C .. Enthalpy
34	jmc	1.1	C qsnow :: snow enthalpy (J/kg)
35			C .. Albedo
36	jmc	1.7	C albColdSnow :: albedo of cold (=dry) new snow (Tsfc < tempSnowAlb)
37	jmc	1.2	C albWarmSnow :: albedo of warm (=wet) new snow (Tsfc = 0)
38	jmc	1.7	C tempSnowAlb :: temperature transition from ColdSnow to WarmSnow Alb. (oC)
39	jmc	1.2	C albOldSnow :: albedo of old snow (snowAge > 35.d)
40			C albIceMax :: max albedo of bare ice (thick ice)
41			C albIceMin :: minimum ice albedo (very thin ice)
42			C hAlbIce :: ice thickness for albedo transition: thin/thick ice albedo
43			C hAlbSnow :: snow thickness for albedo transition: snow/ice albedo
44			C hNewSnowAge :: new snow thickness that refresh the snow-age (by 1/e)
45	jmc	1.6	C snowAgTime :: snow aging time scale (s)
46	jmc	1.2	C .. Solar parameters
47	jmc	1.1	C i0 :: fraction of penetrating solar rad
48			C ksolar :: bulk solar abs coeff of sea ice (m-1)
49			C .. Salinity
50			C saltice :: salinity of ice (o/oo)
51			C S_winton :: winton salinity of ice (o/oo)
52			C mu_Tf :: Tf:brine salinity ratio (C/ppt)
53			C .. melting
54			C Tf0kel :: Freezing temp of fresh ice in Kelvin = 273.15
55			C Tmlt1 :: melting temp; depends on S (C)
56			C .. Min/Max
57			C himin :: minimum thickness for ice (m)
58			C Terrmax :: temperature convergence precision (C)
59			C hiMax :: Maximum thickness for ice (m)
60			C hsMax :: Maximum thickness for snow (m)
61			C .. for fractional ice
62			C iceMaskmax :: maximum Ice fraction (=1 for no fractional ice)
63			C iceMaskmin :: mimimum Ice fraction (=1 for no fractional ice)
64			C himin0 :: minimum ice height (=himin for no fractional ice)
65			C frac_energy:: fraction of energy going to melting (=0 for no fract. ice)
66			C hihig :: ice height above which freezing only occurs over open ocean
67			C (=large for no fractional ice)
68			C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-\|--+----\|
69
70	jmc	1.6	COMMON / THSICE_PHYSPAR_R /
71	jmc	1.1	& rhos, rhoi, rhosw, rhofw, rhoiw,
72			& cpice, cpwater,
73			& kice, ksnow,
74			& transcoef, Lfresh, qsnow,
75	jmc	1.7	& albColdSnow, albWarmSnow, tempSnowAlb,
76			& albOldSnow, hNewSnowAge, snowAgTime,
77	jmc	1.2	& albIceMax, albIceMin, hAlbIce, hAlbSnow,
78	jmc	1.1	& i0, ksolar,
79			& saltice, S_winton, mu_Tf,
80			& Tf0kel, Tmlt1,
81			& himin, Terrmax, hiMax, hsMax,
82			& iceMaskmax, iceMaskmin, himin0,
83			& frac_energy, hihig
84
85			_RL rhos
86			_RL rhoi
87			_RL rhosw
88			_RL rhofw
89			_RL rhoiw
90			_RL cpice
91			_RL cpwater
92			_RL kice
93			_RL ksnow
94			_RL transcoef
95			_RL Lfresh
96			_RL qsnow
97	jmc	1.2	_RL albColdSnow
98			_RL albWarmSnow
99	jmc	1.7	_RL tempSnowAlb
100	jmc	1.2	_RL albOldSnow
101	jmc	1.7	_RL hNewSnowAge
102			_RL snowAgTime
103	jmc	1.2	_RL albIceMax
104			_RL albIceMin
105			_RL hAlbIce
106			_RL hAlbSnow
107	jmc	1.1	_RL i0
108			_RL ksolar
109			_RL saltice
110			_RL S_winton
111			_RL mu_Tf
112			_RL Tf0kel
113			_RL Tmlt1
114			_RL himin
115			_RL Terrmax
116			_RL hiMax
117			_RL hsMax
118			_RL iceMaskmax
119			_RL iceMaskmin
120			_RL himin0
121			_RL frac_energy
122			_RL hihig
123
124			C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-\|--+----\|
125
126	jmc	1.3	C-- COMMON / THSICE_PAR_L / ice model (logical) parameters
127	edhill	1.5	C stepFwd_oceMxL :: step forward mixed-layer T & S (slab-ocean)
128			C thSIce_tave_mdsio :: write TimeAverage output using MDSIO
129			C thSIce_snapshot_mdsio :: write snap-shot output using MDSIO
130			C thSIce_mon_stdio :: write monitor to std-outp
131			C thSIce_tave_mnc :: write TimeAverage output using MNC
132			C thSIce_snapshot_mnc :: write snap-shot output using MNC
133			C thSIce_mon_mnc :: write monitor to netcdf file
134			C thSIce_pickup_read_mnc :: pickup read w/ MNC
135			C thSIce_pickup_write_mnc :: pickup write w/ MNC
136			C thSIce_pickup_write_mdsio :: pickup write w/ MDSIO
137	jmc	1.6	COMMON / THSICE_PAR_L /
138	edhill	1.5	& stepFwd_oceMxL,
139			& thSIce_tave_mdsio, thSIce_snapshot_mdsio, thSIce_mon_stdio,
140			& thSIce_tave_mnc, thSIce_snapshot_mnc, thSIce_mon_mnc,
141	jmc	1.6	& thSIce_pickup_read_mnc,
142	edhill	1.5	& thSIce_pickup_write_mdsio,
143			& thSIce_pickup_write_mnc
144	jmc	1.3
145			LOGICAL stepFwd_oceMxL
146	edhill	1.5	LOGICAL thSIce_tave_mdsio, thSIce_snapshot_mdsio, thSIce_mon_stdio
147			LOGICAL thSIce_tave_mnc, thSIce_snapshot_mnc, thSIce_mon_mnc
148			LOGICAL thSIce_pickup_read_mnc
149			LOGICAL thSIce_pickup_write_mdsio
150			LOGICAL thSIce_pickup_write_mnc
151	jmc	1.3
152	jmc	1.1	C-- COMMON / THSICE_PAR_I / ice model (integer) parameters
153	jmc	1.8	C startIceModel :: =1 : start ice model at nIter0 ; =0 : use pickup files
154			C :: -1 : start from a small pickup (without Mix.Layer)
155			C nitMaxTsf :: maximum Nb of iter to find Surface Temp (Trsf)
156			C thSIceAdvScheme :: thSIce Advection scheme selector
157	jmc	1.6	COMMON / THSICE_PAR_I /
158	jmc	1.8	& startIceModel, nitMaxTsf, thSIceAdvScheme
159	jmc	1.1
160			INTEGER startIceModel
161			INTEGER nitMaxTsf
162	jmc	1.8	INTEGER thSIceAdvScheme
163	jmc	1.1
164			C-- COMMON / THSICE_PAR_R / ice model (real) parameters
165	jscott	1.9	C thSIce_deltaT :: ice model time-step, seaice thicken/extend [s]
166			C thSIce_dtTemp :: ice model time-step, solve4temp [s]
167	jmc	1.8	C ocean_deltaT :: ocean mixed-layer time-step [s]
168	jmc	1.3	C tauRelax_MxL :: Relaxation time scale for MixLayer T & S [s]
169			C hMxL_default :: default value for ocean MixLayer thickness [m]
170			C sMxL_default :: default value for salinity in MixLayer [psu]
171			C vMxL_default :: default value for ocean current velocity in MxL [m/s]
172	jmc	1.8	C thSIce_diffK :: thickness (horizontal) diffusivity [m^2/s]
173	jmc	1.3	C stressReduction :: reduction factor for wind-stress under sea-ice [0-1]
174			C thSIce_taveFreq :: Frequency^-1 for time-Aver. output [s]
175			C thSIce_diagFreq :: Frequency^-1 for diagnostic output [s]
176			C thSIce_monFreq :: Frequency^-1 for monitor output [s]
177	jmc	1.1	COMMON / THSICE_PAR_R /
178	jscott	1.9	& thSIce_deltaT, thSIce_dtTemp, ocean_deltaT, tauRelax_MxL,
179	jmc	1.3	& hMxL_default, sMxL_default, vMxL_default,
180	jmc	1.8	& thSIce_diffK, stressReduction,
181	jmc	1.1	& thSIce_taveFreq, thSIce_diagFreq, thSIce_monFreq
182
183	jscott	1.9	_RL thSIce_deltaT, thSIce_dtTemp, ocean_deltaT, tauRelax_MxL
184	jmc	1.6	_RL hMxL_default, sMxL_default, vMxL_default
185	jmc	1.8	_RL thSIce_diffK, stressReduction
186	jmc	1.1	_RL thSIce_taveFreq, thSIce_diagFreq, thSIce_monFreq
187
188	jmc	1.4	C-- COMMON / THSICE_PAR_C / ice model (character) parameters
189			C thSIceFract_InitFile :: File name for initial ice fraction
190			C thSIceThick_InitFile :: File name for initial ice thickness
191			C thSIceSnowH_InitFile :: File name for initial snow thickness
192			C thSIceSnowA_InitFile :: File name for initial snow Age
193			C thSIceEnthp_InitFile :: File name for initial ice enthalpy
194			C thSIceTsurf_InitFile :: File name for initial surf. temp
195			COMMON / THSICE_PAR_C /
196			& thSIceFract_InitFile,
197			& thSIceThick_InitFile,
198			& thSIceSnowH_InitFile,
199			& thSIceSnowA_InitFile,
200			& thSIceEnthp_InitFile,
201	jmc	1.6	& thSIceTsurf_InitFile
202	jmc	1.4	CHARACTER*(MAX_LEN_FNAM) thSIceFract_InitFile
203			CHARACTER*(MAX_LEN_FNAM) thSIceThick_InitFile
204			CHARACTER*(MAX_LEN_FNAM) thSIceSnowH_InitFile
205			CHARACTER*(MAX_LEN_FNAM) thSIceSnowA_InitFile
206			CHARACTER*(MAX_LEN_FNAM) thSIceEnthp_InitFile
207			CHARACTER*(MAX_LEN_FNAM) thSIceTsurf_InitFile
208
209	jmc	1.1	C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-\|--+----\|
210
211			#endif /* ALLOW_THSICE */
212	edhill	1.5
213			CEH3 ;;; Local Variables: ***
214			CEH3 ;;; mode:fortran ***
215			CEH3 ;;; End: ***