use abs cause non-even power used down the line

Author: Matthias Fabry

star/private/hydro_rotation.f90

@@ -642,7 +642,7 @@ subroutine eval_fp_ft(id, nz, xm, r, rho, aw, ft, fp, r_polar, r_equatorial, rep
          do j=1, s% nz
             ! Compute fp, ft, re and rp using fits to the Roche geometry of a single star.
             ! by this point in the code, w_div_w_crit_roche is set
-            w = s% w_div_w_crit_roche(j)
+            w = abs(s% w_div_w_crit_roche(j))
             w% d1val1 = 1d0
 
             w2 = pow2(w)
@@ -701,7 +701,7 @@ subroutine eval_i_rot(s, k, r00, w_div_w_crit_roche, i_rot)
             i_rot% d1Array(i_lnR_00) = 2d0 * i_rot% val
             i_rot% d1Array(i_w_div_wc_00) = 0d0
          else
-            w = w_div_w_crit_roche
+            w = abs(w_div_w_crit_roche)
             w% d1val1 = 1d0
 
             r = r00