Move the part for ok_exlicitsnow from enerbil_fusion to the the end of thermosoil_main and thermosoilc_main. This does not make any change in the order of calculations.
temp_sol_new is now modified in thermosoil_main/thermosoilc_main (intent(inout)) but only for the case ok_explicitsnow
call enerbil_fusion after hydrol_main as done before changeset [3269]
protected variable fusion from writing for case ok_explicitsnow. This variable is calculated by enerbil_fusion