c		****************************************
c		Parameters for thermodynamic computation
c		****************************************
c
c		hmelt, Maximum melting at the bottom
      -0.15
c		acrit(1), Minimum fraction for leads in the NH
      1.0e-06
c		acrit(2), Minimum fraction for leads in the SH
      1.0e-06
c		hgcrit(1), Ice thickness for lateral accretion in the NH
      0.3
c               hgcrit(2), Ice thickness for lateral accretion in the SH
      0.3
c               hgmin, Ice thickness corr. to max. energy stored in brine pocket
      0.2
c		hndif, Computation of temp. in snow (=0.0) or not (=9999.0)
      0.0      9999.0   
c		hgdif, Computation of temp. in ice (=0.0) or not (=9999.0)
      0.0     9999.0
c		hglim, Minimum ice thickness
      0.05
c		amax, Maximum lead fraction
      0.999   0.85
c		swiqst, Energy stored in brine pocket (=1) or not (=0)	
      1
c		beta, Numerical caracteritic of the scheme for diffusion in ice
c		      Cranck-Nicholson (=0.5), implicit (=1), explicit (=0)
      1
c		ddtb, Time step in thermodynamic calculation
      5760.00
c		nbits, Number of time steps in Newton -Raphson procedure
      1
c               parlat, Percentage of energy used for lateral ablation
      0.0
c               hakspl, Slope of distr. for Hakkinen-Mellor's lateral melting
      0.5
c               hibspl, Slope of distribution for Hibler's lateral melting
      0.5
c               exld, Exponent for leads-closure rate
      2.0
c               hakdif, Coefficient for diffusions of ice and snow
      1.0
c               hth, Threshold thickness for comp. of eq. thermal conductivity
      0.2
c               hnzst, Thickness of the surf. layer in temp. computation
      0.1
c               parsub, Switch for snow sublimation or not
      1.0
c               alphs, Coefficient for snow density when snow ice formation
      1.0       1.1236363636