!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
!
! ICE namelist :  ice model option and parameter input
! -------------
!
!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
!
!'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
!       namicerun   Share parameters for dynamics/advection/thermo
!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
!
!  ln_limdyn        : switch for ice dynamics (true) or not (false)
!  acrit(1/2)  : minimum fraction for leads in the Northern (Southern) Hemisphere
!  hnsdif      : computation of temperature in snow (=0.0) or not (=9999.0)
!  hicdif      : computation of temperature in ice  (=0.0) or not (=9999.0)
!  cai         : atmospheric drag over sea ice
!  cao         : atmospheric drag over ocean
!  ln_nicep  : Ice points output for debug (yes or no)
!
&namicerun
   ln_limdyn    = .true. 
   acrit   = 1.0e-02 , 1.0e-02
   hsndif  = 0.0 
   hicdif  = 0.0
   cai     = 1.40e-3
   cao     = 1.00e-3
   ln_nicep  = .true.
/
!'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
!       namiceini   parameters for ice initialisation
!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
!
!  ttest   : threshold water temperature for initial sea ice
!  hninn   : initial snow thickness in the north
!  hginn_u : initial undeformed ice thickness in the north
!  aginn_u : initial undeformed ice concentration in the north
!  hginn_d : initial  deformed  ice thickness in the north
!  aginn_d : initial  deformed  ice concentration in the north
!  hnins   : initial snow thickness in the south
!  hgins_u : initial undeformed ice thickness in the north
!  agins_u : initial undeformed ice concentration in the north
!  hgins_d : initial  deformed  ice thickness in the north
!  agins_d : initial  deformed  ice concentration in the north
!  sinn    : initial salinity in the north
!  sins    : initial salinity in the south
!
&namiceini
   ttest   = 2.0
   hninn   = 0.3
   hginn_u = 3.50
   aginn_u = 0.95
   hginn_d = 0.0
   aginn_d = 0.00
   hnins = 0.1
   hgins_u = 1.0
   agins_u = 0.9
   hgins_d = 0.0
   agins_d = 0.00
   sinn  = 6.301
   sins  = 6.301
/
!'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
!       namicedia   parameters for ice diagnostics
!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
!
!  fmtinf : format of the output values 
!  nfrinf : number of variables written in one line
!  ntmoy  : instantaneous values of ice evolution or averaging
!  ninfo  : frequency of ouputs on file ice_evolu in case of averaging
!
&namicedia
   fmtinf = '1PE13.5 '
   nfrinf = 4
   ntmoy  = 1
   ninfo  = 1
/
!'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
!       namicedyn   parameters for ice dynamic
!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
!
!  epsd   : tolerance parameter
!  alpha  : coefficient for semi-implicit coriolis
!  bound  : boundary conditions (=0.0 no-slip, =1.0 free-slip)
!  dm     : diffusion constant for dynamics.
!  nbiter : number of sub-time steps for relaxation
!  nbitdr : maximum number of iterations for relaxation
!  om     : relaxation constant 
!  resl   : maximum value for the residual of relaxation 
!  cw     : drag coefficient for oceanic stress
!  angvg  : turning angle for oceanic stress
!  pstar  : first bulk-rheology parameter
!  c_rhg  : second bulk-rhelogy parameter
!  etamn  : minimun value for viscosity
!  creepl : creep limit
!  ecc    : eccentricity of the elliptical yield curve
!  ahi0   : horizontal eddy diffusivity coefficient for sea-ice (m2/s)
!  nevp   : number of iterations for subcycling in EVP, SB, 400
!  telest : timescale for elastic waves, SB, 720.0
!  alphaevp : coefficient for the solution of internal ice stresses
!
&namicedyn 
   epsd   = 1.0e-20
   alpha  = 0.5
   dm     = 0.6e+03
   nbiter = 1
   nbitdr = 500
   om     = 0.5
   resl   = 5.0e-05
   cw     = 5.0e-03
   angvg  = 0.0
   pstar  = 4.00e+04 
   c_rhg  = 20.0
   etamn  = 0.0e+07
   creepl = 1.0e-12
   ecc    = 2.0
   ahi0   = 350.e0
   nevp   = 300
   telast = 9600.0
   alphaevp = 1.0
/
!
!'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
!       namicetrp   parameters for ice advection
!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
!  bound  : boundary conditions (=0.0 no-slip, =1.0 free-slip)
&namicetrp 
   bound  = 0.
/
!'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
!       namicethd   parameters for thermodynamic computation
!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
!
!  hmelt       : maximum melting at the bottom
!  hiccrit(1/2): ice thickness for lateral accretion in the Northern (Southern) Hemisphere
!                caution 1.0, 1.0 best value to be used!!! (gilles G.)
!  fraz_swi    : use of frazil ice collection thickness in function of wind (1.0) or not (0.0)
!  maxfrazb    : maximum portion of frazil ice collecting at the ice bottom
!  vfrazb      : thresold drift speed for frazil ice collecting at the ice bottom
!  Cfrazb      : squeezing coefficient for frazil ice collecting at the ice bottom
!  hicmin      : ice thickness corr. to max. energy stored in brine pocket
!  hiclim      : minimum ice thickness
!  amax        : maximum lead fraction
!  sbeta       : numerical caracteritic of the scheme for diffusion in ice
!               Cranck-Nicholson (=0.5), implicit (=1), explicit (=0)
!  parlat      : percentage of energy used for lateral ablation
!  hakspl      : slope of distr. for Hakkinen-Mellor's lateral melting
!  hibspl      : slope of distribution for Hibler's lateral melting
!  exld        : exponent for leads-closure rate
!  hakdif      : coefficient for diffusions of ice and snow
!  thth         : threshold thickness for comp. of eq. thermal conductivity
!  hnzst       : thickness of the surf. layer in temp. computation
!  parsub      : switch for snow sublimation or not
!  alphs       : coefficient for snow density when snow ice formation
!  betas       : exponent in lead-ice fractionation of snow precipitation 0.66
!                betas = 1 -> equipartition, betas < 1 -> more on leads
!  kappa_i     : extinction radiation parameter in sea ice (1.0)
!  nconv_i_thd : maximal number of iterations for heat diffusion computation
!  maxer_i_thd : maximal error in temperature for heat diffusion computation
!  thcon_i_swi : switch for computation of thermal conductivity in the ice
!                (0) Untersteiner (1964), (1) Pringle et al. (2007)
! 
&namicethd
   hmelt   = -0.15
   hiccrit = 0.1 , 0.1
   fraz_swi= 0.0
   maxfrazb= 0.0
   vfrazb  = 0.4166667
   Cfrazb  = 5.0
   hicmin  = 0.2
   hiclim  = 0.10
   amax    = 0.999 
   sbeta   = 1.
   parlat  = 0.0
   hakspl  = 0.25
   hibspl  = 0.5
   exld    = 2.0
   hakdif  = 1.0
   thth    = 0.2
   hnzst   = 0.1
   parsub  = 1.0
   alphs   = 1.0
   betas   = 0.6
   kappa_i = 1.0 
   nconv_i_thd = 50
   maxer_i_thd = 0.0001
   thcon_i_swi = 1
/
!'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
!       namicesal   parameters for salinity calculation
!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
!
!   num_sal = number of salinity option
!      1 -> S = bulk_sal
!      2 -> S = S(z,t) with a simple parameterization
!      3 -> S = S(z) profile of Scwharzacher [1959]
!      4 -> S = S(h) Cox and Weeks [1974]
!   bulk_sal = if 1 is used, it represents the ice salinity
!   sal_G  : restoring salinity for GD 
!   time_G : restoring time for GD
!   sal_F  : restoring salinity for flushing
!   time_F : restoring time for flushing 
!   s_i_max : Maximum salinity     
!   s_i_min : Minimum tolerated ice salinity   
!   s_i_0   : 1st salinity for salinity profile 
!   s_i_1   : 2nd salinity for salinity profile 

!
&namicesal
   num_sal  = 2
   bulk_sal = 4.0
   sal_G    = 5.00
   time_G   = 1.728e+06
   sal_F    = 2.00
   time_F   = 8.640e+05
   s_i_max  = 20.0
   s_i_min  =  0.1
   s_i_0    =  3.5
   s_i_1    =  4.5
/
!'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
!       namiceitdme   parameters for mechanical redistribution of ice 
!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
!
!   ridge_scheme_swi ! which ridging scheme using (1=Rothrock,else=Hibler79)
!   Cs               ! shearing energy contribution to ridging 
!   Cf               ! ratio of ridging work to PE change in ridging        
!   fsnowrdg         ! snow fraction that survives in ridging
!   fsnowrdg         ! snow fraction that survives in rafting
!   Gstar            ! fractional area of thin ice being ridged
!   astar            ! equivalent of gstar (0.05 for TH75 and 0.03 for weaker ice)
!   Hstar            ! parameter determining the maximum thickness of ridged ice
!   raftswi          ! rafting or not
!   hparmeter        ! threshold thickness for rafting or not
!   Craft            ! coefficient used in the rafting function
!   ridge_por        ! intial porosity of the ridged ice (typically 0.30)
!   sal_max_ridge    ! maximum ridged ice salinity
!   partfun_swi      ! participation function linear, TH75 (0) or exponential Letal07 (1)
!   transfun_swi     ! transfer function uniform of H80 (0) or exponential Letal07 (1)
!   brinstren_swi    ! (1) use brine volume to diminish ice strength
&namiceitdme
   ridge_scheme_swi = 0
   Cs = 0.50
   Cf = 17.0
   fsnowrdg  = 0.5
   fsnowrft  = 0.5
   Gstar     = 0.15
   astar     = 0.05
   Hstar     = 100.0
   raftswi   = 1
   hparmeter = 0.75
   Craft     = 5.0
   ridge_por = 0.3
   sal_max_ridge = 15.0
   partfun_swi   = 1
   transfun_swi  = 0
   brinstren_swi = 0
/
!'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
!       namiceout   parameters for outputs
!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
!
!   noumef      : number of fields
!   field
!        title  : title of the field
!        name   : name of the field
!        unit   : unit of the field
!        nc     : switch for saving field ( = 1 ) or not ( = 0 )
!        cmulti : multiplicative constant
!        cadd   : additive constant
!  add_diag_swi : 1 -> diagnose distribution in thickness space
!                 0 -> only simple diagnostics
!
&namiceout
   noumef   = 37
   field_1  = 'Ice concentration                  ', 'iiceconc', '%       ', 1, 1.0, 0.0
   field_2  = 'Ice thickness                      ', 'iicethic', 'm       ', 1, 1.0, 0.0
   field_3  = 'Snow thickness                     ', 'isnowthi', 'm       ', 1, 1.0, 0.0
   field_4  = 'Daily bottom thermo ice production ', 'iicebopr', 'cm/day  ', 1, 100., 0.0
   field_5  = 'Daily dynamic ice production       ', 'iicedypr', 'cm/day  ', 1, 100., 0.0
   field_6  = 'Oceanic flux at the ice base       ', 'ioceflxb', 'w/m2    ', 1, 1.0, 0.0
   field_7  = 'Ice velocity u                     ', 'iicevelu', 'm/s     ', 1, 1.0, 0.0
   field_8  = 'Ice velocity v                     ', 'iicevelv', 'm/s     ', 1, 1.0, 0.0
   field_9  = 'Sea surface temperature            ', 'isstempe', 'C       ', 1, 1.0, -273.15
   field_10 = 'Sea surface salinity               ', 'isssalin', 'PSU     ', 1, 1.0, 0.0
   field_11 = 'Total flux at ocean surface        ', 'iocetflx', 'w/m2    ', 1, 1.0, 0.0
   field_12 = 'Solar flux at ocean surface        ', 'iocesflx', 'w/m2    ', 1, 1.0, 0.0
   field_13 = 'Non-solar flux at ocean surface    ', 'iocwnsfl', 'w/m2    ', 1, 1.0, 0.0
   field_14 = 'Heat flux due to brine release     ', 'iocehebr', 'w/m2    ', 1, 1.0, 0.0
   field_15 = 'Wind stress u                      ', 'iocestru', 'Pa      ', 1, 1.0, 0.0
   field_16 = 'Wind stress v                      ', 'iocestrv', 'Pa      ', 1, 1.0, 0.0 
   field_17 = 'Solar flux at ice/ocean surface    ', 'iicesflx', 'w/m2    ', 1, 1.0, 0.0
   field_18 = 'Non-solar flux at ice/ocean surface', 'iicenflx', 'w/m2    ', 1, 1.0, 0.0
   field_19 = 'Snow precipitation                 ', 'isnowpre', 'kg/day  ', 1, 1.0, 0.0
   field_20 = 'Mean ice salinity                  ', 'iicesali', 'psu     ', 1, 1.0, 0.0
   field_21 = 'Mean ice age                       ', 'iiceages', 'years   ', 1, 0.002739, 0.0
   field_22 = 'Daily lateral thermo ice prod.     ', 'iicelapr', 'cm/day  ', 1, 100., 0.0
   field_23 = 'Daily snowice ice production       ', 'iicesipr', 'cm/day  ', 1, 100., 0.0
   field_24 = 'Mean ice temperature               ', 'iicetemp', 'C       ', 1, 1.0, -273.15
   field_25 = 'Ice total heat content             ', 'iiceheco', '10^9 J  ', 1, 1.0, 0.0
   field_26 = 'Ice surface temperature            ', 'iicesurt', 'C       ', 1, 1.0, -273.15
   field_27 = 'Snow temperature                   ', 'isnotem2', 'C       ', 1, 1.0, -273.15
   field_28 = 'Fsbri - brine salt flux            ', 'iicfsbri', 'kg/m2/s ', 1, 1.0, 0.0
   field_29 = 'Fseqv - equivalent FW salt flux    ', 'iicfseqv', 'kg/m2/s ', 1, 1.0, 0.0
   field_30 = 'Brine volume                       ', 'ibrinvol', '%       ', 1, 100.0, 0.0
   field_31 = 'Frazil ice collection thickness    ', 'iicecolf', 'm       ', 1, 1.0, 0.0
   field_32 = 'Ice strength                       ', 'iicestre', 'N/m     ', 1, 0.001, 0.0
   field_33 = 'Ice velocity                       ', 'iicevelo', 'm/s     ', 1, 1.0, 0.0
   field_34 = 'Surface melt                       ', 'iicesume', 'cm/day  ', 1, 100., 0.0
   field_35 = 'Bottom melt                        ', 'iicebome', 'cm/day  ', 1, 100., 0.0
   field_36 = 'Divergence                         ', 'iicedive', '10-8s-1 ', 1, 1.0e8, 0.0
   field_37 = 'Shear                              ', 'iiceshea', '10-8s-1 ', 1, 1.0e8, 0.0
   add_diag_swi = 1
/