anchor.diff:  
     close (un)
   
  +
  +
  +  open (unit=un, file=trim(dir)//"/Rotation", status="unknown")
  +
  +  write (un,100) rot_period
  +
  +  close (un)
  +
   end subroutine problem_checkpoint
   
   
  
positive.code:   read (un,100) vel_P(3), vel_S(3)
  
    close (un)
  
  end subroutine problem_restart
  
  
negative.code:   use bl_IO_module
    use probdata_module, only: com_P, com_S, vel_P, vel_S, mass_P, mass_S
    use prob_params_module, only: center
    use meth_params_module, only: rot_period
  
    implicit none
  
  