# Contains all information for calculation of surface flux dependence on variables

[iimelt]
summer_formula1=SW_down * (l_snow.0.01 * (ac.alpha_cs - ac.alpha_s) + (l_snow.0.01 * -1. + 1.) * (ac.alpha_ci - ac.alpha_i))
winter_formula1=SW_down * (l_snow.0.01 * (ac.alpha_cs - ac.alpha_s) + (l_snow.0.01 * -1. + 1.) * (ac.alpha_ci - ac.alpha_i)) * complex_scale_factor
cvnames=iimelt,SW_down,l_snow.0.01,complex_scale_factor
obs_cvnames=iimelt

[Ice_area]
summer_formula1=SW_down * (ac.alpha_o - ac.alpha_i)
winter_formula1=(f_atmos_ice + B_lw_dep * (TTsfc*-1. + ac.Tbot)) * complex_scale_factor - f_atmos_ocean
cvnames=Ice_area,SW_down,complex_scale_factor,f_atmos_ice,B_lw_dep,TTsfc,f_atmos_ocean
obs_cvnames=Ice_area

[SW_down]
winter_formula1=Ice_area * complex_scale_factor * (ice_albedo * -1. + 1.) + (Ice_area * -1. + 1.) * (1. - ac.alpha_o)
summer_formula1=Ice_area * (ice_albedo * -1. + 1.) + (Ice_area * -1. + 1.) * (1. - ac.alpha_o)
cvnames=SW_down,Ice_area,complex_scale_factor,l_snow.0.01,ice_albedo
obs_cvnames=SW_down

[LW_down]
winter_formula1=complex_scale_factor * Ice_area + (Ice_area * -1. + 1.)
summer_formula1=1.
cvnames=LW_down,complex_scale_factor,Ice_area
obs_cvnames=LW_down

[Ice_thickness.cat.1]
winter_formula1=Ice_area.cat.1*complex_scale_factor.cat.1 ** 2. * (LW_down - C_lw_ref * 1. + B_lw_dep * (TTsfc*-1. + ac.Tbot) + Sensible_heat_flux_over_ice) * B_lw_dep / ac.kice
summer_formula1=0.
winter_formula2=Ice_area.cat.1*complex_scale_factor.cat.1 ** 3. * (LW_down - C_lw_ref * 1. + B_lw_dep * (TTsfc*-1. + ac.Tbot) + Sensible_heat_flux_over_ice) * (B_lw_dep / ac.kice) ** 2.
summer_formula2=0.
cvnames=LW_down,Ice_thickness.cat.1,Sensible_heat_flux_over_ice,B_lw_dep,C_lw_ref,TTsfc,complex_scale_factor.cat.1,Ice_area.cat.1
obs_cvnames=LW_down,Ice_thickness.cat.1

[Ice_thickness.cat.2]
winter_formula1=Ice_area.cat.2*complex_scale_factor.cat.2 ** 2. * (LW_down - C_lw_ref * 1. + B_lw_dep * (TTsfc*-1. + ac.Tbot) + Sensible_heat_flux_over_ice) * B_lw_dep / ac.kice
summer_formula1=0.
winter_formula2=Ice_area.cat.2*complex_scale_factor.cat.2 ** 3. * (LW_down - C_lw_ref * 1. + B_lw_dep * (TTsfc*-1. + ac.Tbot) + Sensible_heat_flux_over_ice) * (B_lw_dep / ac.kice) ** 2.
summer_formula2=0.
cvnames=LW_down,Ice_thickness.cat.2,Sensible_heat_flux_over_ice,B_lw_dep,C_lw_ref,TTsfc,complex_scale_factor.cat.2,Ice_area.cat.2
obs_cvnames=LW_down,Ice_thickness.cat.2

[Ice_thickness.cat.3]
winter_formula1=Ice_area.cat.3*complex_scale_factor.cat.3 ** 2. * (LW_down - C_lw_ref * 1. + B_lw_dep * (TTsfc*-1. + ac.Tbot) + Sensible_heat_flux_over_ice) * B_lw_dep / ac.kice
summer_formula1=0.
winter_formula2=Ice_area.cat.3*complex_scale_factor.cat.3 ** 3. * (LW_down - C_lw_ref * 1. + B_lw_dep * (TTsfc*-1. + ac.Tbot) + Sensible_heat_flux_over_ice) * (B_lw_dep / ac.kice) ** 2.
summer_formula2=0.
cvnames=LW_down,Ice_thickness.cat.3,Sensible_heat_flux_over_ice,B_lw_dep,C_lw_ref,TTsfc,complex_scale_factor.cat.3,Ice_area.cat.3
obs_cvnames=LW_down,Ice_thickness.cat.3

[Ice_thickness.cat.4]
winter_formula1=Ice_area.cat.4*complex_scale_factor.cat.4 ** 2. * (LW_down - C_lw_ref * 1. + B_lw_dep * (TTsfc*-1. + ac.Tbot) + Sensible_heat_flux_over_ice) * B_lw_dep / ac.kice
summer_formula1=0.
winter_formula2=Ice_area.cat.4*complex_scale_factor.cat.4 ** 3. * (LW_down - C_lw_ref * 1. + B_lw_dep * (TTsfc*-1. + ac.Tbot) + Sensible_heat_flux_over_ice) * (B_lw_dep / ac.kice) ** 2.
summer_formula2=0.
cvnames=LW_down,Ice_thickness.cat.4,Sensible_heat_flux_over_ice,B_lw_dep,C_lw_ref,TTsfc,complex_scale_factor.cat.4,Ice_area.cat.4
obs_cvnames=LW_down,Ice_thickness.cat.4

[Ice_thickness.cat.5]
winter_formula1=Ice_area.cat.5*complex_scale_factor.cat.5 ** 2. * (LW_down - C_lw_ref * 1. + B_lw_dep * (TTsfc*-1. + ac.Tbot) + Sensible_heat_flux_over_ice) * B_lw_dep / ac.kice
summer_formula1=0.
winter_formula2=Ice_area.cat.5*complex_scale_factor.cat.5 ** 3. * (LW_down - C_lw_ref * 1. + B_lw_dep * (TTsfc*-1. + ac.Tbot) + Sensible_heat_flux_over_ice) * (B_lw_dep / ac.kice) ** 2.
summer_formula2=0.
cvnames=LW_down,Ice_thickness.cat.5,Sensible_heat_flux_over_ice,B_lw_dep,C_lw_ref,TTsfc,complex_scale_factor.cat.5,Ice_area.cat.5
obs_cvnames=LW_down,Ice_thickness.cat.5
