minor phase separation cleanup

Author: evbauer

star/private/phase_separation.f90

@@ -54,7 +54,7 @@ subroutine do_phase_separation(s, dt, ierr)
          integer :: k, k_bound, kstart, net_ic12, net_io16
          logical :: save_Skye_use_ion_offsets
 
-         s% eps_phase_separation(:) = 0d0
+         s% eps_phase_separation(1:s%nz) = 0d0
 
          if(s% phase(s% nz) < eos_phase_boundary) then
             s% crystal_core_boundary_mass = 0d0

star/private/profile_getval.f90

@@ -1504,7 +1504,7 @@ subroutine getval_for_profile(s, c, k, val, int_flag, int_val)
                end if
 
             case (p_eps_phase_separation)
-               if (s% do_phase_separation_heating) val = s% eps_phase_separation(k)
+               if (s% do_phase_separation .and. s% do_phase_separation_heating) val = s% eps_phase_separation(k)
 
             case (p_eps_WD_sedimentation)
                if (s% do_element_diffusion) val = s% eps_WD_sedimentation(k)