anchor.diff:            zo = real(coordinate_Zcylindr(wd_index) * METERS_IN_PARSEC)
   C         TODO: find out the meaning of ecini
             ecini = 0.5 * wo * wo
  -          call get_potential_energy(zo, epoti)
             time_increment = (galacticDiskAge - starBirthTime(wd_index)) 
        &                     / float(NJUMPS)
   C         Time in seconds
  
positive.code:      &        htry,
       &        initial_time,
       &        ecinf,
       &        epotf,
       &        epoti,
       &        f
        real :: uu(MAX_STARS_COUNT),
       &        vv(MAX_STARS_COUNT),
  
negative.code:      &                HMIN, NOK, NBAD, DERIVS, RKQC, yscal, y, dydx)
  C         TODO: find out the meaning of ecinf    
            ecinf = 0.5 * y(2) * y(2)
            coordinate_Zcylindr(wd_index) = y(1) / METERS_IN_PARSEC
            ww(wd_index) = y(2) + WOSUN
        end do
        end subroutine
  
  
  