anchor.diff:  
     ! The following defaults are already set in global_varaibles.f90.
     ! The user is free to change here too.
  -  dT         = 0.01_dp
  +  dT         = 0.001_dp
     start_time = 0.0_dp
     end_time   = 1.0_dp
   
  
positive.code: 
       end do newton
  
       ! print the summary of the newton iteration
       write(*,'(f15.2,e15.6,e15.6,e15.6,e15.6,e15.6,xi4,xi8)') &
            &time, update_norm, res_norm, body1%KE, body1%PE, &
            & body1%KE + body1%PE, newton_cnt, fcnt
  
       
    end do time_march
  
    call print_body(body1)
  
negative.code:        &  '             KE', '          PE','               TE', &
         &'       Niter',' FCNT'
  
    !-------------------------------------------------------------------!
    ! Time integration loop
    !-------------------------------------------------------------------!
  
    time = start_time
  
    time_march: do while ( time .le. end_time)  ! loop for time-marching
  
       time = time + dT ! update the time
  
       !----------------------------------------------------------------!
       ! Newton iteration loop
       !----------------------------------------------------------------!
  
       newton_cnt = 0
  
       newton: do k = 1, MAX_NEWTON_ITER
  