c  Physical constants first
      parameter (hfusion=2.5e6,rhoocean=1026.)
      parameter (cpice=2090., rhoice=910., rhowater=1026.,
     +           hfusionice=3.34e5)

c  Iteration parameter for icethermo   
      parameter (dqmax=0.1, tconvin=0.01, tconvgr=1.30, tconvmax=.3)

c  Default values for common block parameter YOU HAVE TO SET THEM !!??

      parameter (albedooceandef=0.07,albedoicedef=0.7,albedofdef=1.)
      parameter (tfreezedef=-1.8+273.15,cicemaxdef=0.99,
     +           hsnowdef=0.06,sicedef=6.)
      parameter (itermaxdef=50,ssticemaxdef=-1+273.15)
      parameter (hicemindef=0.05)
      parameter (tksnowdef=0.33 , tkicedef=2.,  tkoceandef=5000.)
      parameter (pbl_eddydef=180000.)
      parameter (pbl_eddy_nlat_def=35.,pbl_eddy_slat_def=-35.)
      parameter (pbl_eddy_nwidth_def=15.,pbl_eddy_swidth_def=15.)
      parameter (pbl_eddy_xbox_def=40.,pbl_eddy_ybox_def=25.)
      parameter (hqdef=0.25, hfdef=2.0)
      parameter (c0_facdef=1.)

c   Common block for constants you might like to change

      common/amlice/ albedoocean, albedoice, tfreeze, cicemax, hsnow,
     +               sice, itermax, ssticemax, hicemin, tksnow, tkice,
     +               tkocean, hq, hf, albedof , pbl_eddy, slat,
     +               pbl_eddy_nlat, pbl_eddy_nwidth,
     +               pbl_eddy_slat, pbl_eddy_swidth,
     +               pbl_eddy_xbox, pbl_eddy_ybox, c0_fac
      

c  common block for test output 

      common /test / qao,qas,qio,qsi,qif,niter,albedo