anchor.diff:        ieleven=11
         itmp = int(log10(float(myrank+1)))+1
   
  +      iotime = TMRC()
         if( input_mode .eq. -1 ) then
           call streamio_setup_read(fhandle, geomRestartStream)
         else if( input_mode .eq. 0 ) then
  
positive.code:       endif
  
        call phio_closefile(fhandle);
  c.... Read restart files
        if( input_mode .eq. -1 ) then
          call streamio_setup_read(fhandle, geomRestartStream)
        else if( input_mode .eq. 0 ) then
  
negative.code:       real*8, target, allocatable :: xread(:,:), qread(:,:), acread(:,:)
        real*8, target, allocatable :: uread(:,:)
        real*8, target, allocatable :: BCinpread(:,:)
        real*8 :: iotime
        integer, target, allocatable :: iperread(:), iBCtmpread(:)
        integer, target, allocatable :: ilworkread(:), nBCread(:)
        character*10 cname2
  