81 xmonth2sec = xyear2day/12.*xday2sec
108 xalf1 = 0.55, xpow = 0.28
132 xmhofusn0 = 3.3355E+05
166 xrki1 = 2.11, xrki2 = -0.011, xrki3 = 0.09, xtempc = -1.8
237 dptfr1 = 0.8, dptfr2 = -1.E-4
277 hofusn0 = xmhofusn0 * rhoice0
280 hofusni0 = 1. / hofusn0
283 rhocpice0 = rhoice0*cpice0