add in extra eddington limit calculation at formation of CO (#569)

Author: katiebreivik

cosmic/src/evolv2.f

@@ -1318,7 +1318,7 @@ SUBROUTINE evolv2(kstar,mass,tb,ecc,z,tphysf,
      &                      bacc(1),bacc(2),tacc(1),tacc(2),epoch(1),
      &                      epoch(2),bhspin(1),bhspin(2))
                CALL kick(kw,mass(k),mt,0.d0,0.d0,-1.d0,0.d0,vk,k,
-     &                  0.d0,fallback,sigmahold,kick_info,disrupt,bkick)
+     &                  0.d0,fallback,sigmahold,kick_info,disrupt,bkick)     
 
                sigma = sigmahold !reset sigma after possible ECSN kick dist. Remove this if u want some kick link to the intial pulsar values...
 * set kick values for the bcm array
@@ -1398,11 +1398,14 @@ SUBROUTINE evolv2(kstar,mass,tb,ecc,z,tphysf,
             mass0(k) = m0
             epoch(k) = tphys - age
          endif
+
 *
 *
 * Force new NS or BH to have a birth spin peirod and magnetic field.
-*
+* 
          if(kstar(k).eq.13.or.kstar(k).eq.14)then
+* re-calculate the Eddington limit of the newly formed CO
+            dme = 2.08d-03*eddfac*(1.d0/(1.d0 + zpars(11)))*rad(k)  
             if(tphys-epoch(k).lt.tiny)then
                if(kstar(k).eq.13.and.pulsar.gt.0)then
 *                  write(93,*)'birth start: ',tphys,k,B_0(k),ospin(k)
@@ -1903,6 +1906,7 @@ SUBROUTINE evolv2(kstar,mass,tb,ecc,z,tphysf,
 * Eddington limit for accretion on to the secondary in one orbit.
 *
  8    dme = 2.08d-03*eddfac*(1.d0/(1.d0 + zpars(11)))*rad(j2)*tb
+      
       supedd = .false.
       novae = .false.
       disk = .false.
@@ -2909,6 +2913,7 @@ SUBROUTINE evolv2(kstar,mass,tb,ecc,z,tphysf,
          dm2 = km*dm2
          dm22 = km*dm22
          dme = km*dme
+
 *
 * Calculate orbital angular momentum change due to system mass loss.
 *
@@ -3647,6 +3652,7 @@ SUBROUTINE evolv2(kstar,mass,tb,ecc,z,tphysf,
      &                    bacc(1),bacc(2),tacc(1),tacc(2),epoch(1),
      &                    epoch(2),bhspin(1),bhspin(2))
       endif
+
 *
 * Test whether the primary still fills its Roche lobe.
 *