anchor.diff:    integer :: int_dir_name(len)
     character (len=len) :: dir
   
  -  integer :: i, un
  +  integer :: i, un, stat
   
     ! dir will be the string name of the checkpoint directory
     do i = 1, len
  
positive.code: 
    if (problem .eq. 3) then
  
       open (unit=un, file=trim(dir)//"/Rotation", status="old")
  
       read (un,100) rot_period
  
       close (un)
  
    endif
  
  
  
    open (unit=un, file=trim(dir)//"/Extrema", status="old")
  
    read (un,100) T_global_max
    read (un,100) rho_global_max
  
    close(un)
  
  end subroutine problem_restart
  
  
negative.code: 
    write (un,100) T_global_max
    write (un,100) rho_global_max
    write (un,100) ts_te_global_max
  
    close (un)
  
  