auto_observation-VIC.f90~

  ke=1001

  iop = (/596/rate /)
  jop = (/634/rate ,636/rate, 638/rate/)


     WRITE(rdirname, '(a,I1.1,a,I1.1)') 'outfiles-VIC-',movex,'-',movey

     do j = 1, size(jop)
        do i = 1, size(iop)

           op1=iop(i)
           op2=jop(j)

          WRITE(filename, '(a,I3.3,a,I3.3,a,I3.3)') 'outE-',op1,'-',op2

          ke = ke+5
          
         open(unit=fu, &
     &           file=trim(trim(rdirname)//'/'//filename), &
     &           access='append', status='unknown')
 
         
          write(ke, *) ex(op1,op2), ey(op1,op2), 2*t

          close(unit=ke)

         
    end do
    end do