anchor.diff:  
           integer :: uid
   
  -        character*22 :: filename
  -        filename = './data/parameters.json'
  -        open(newunit=uid,file=filename)
  +        character(*), parameter :: filename = 'parameters.json'
  +        open(newunit=uid,file=write_dir//filename)
   
           write(uid,'(A)') '{'
   
  
positive.code:         write(uid,'(A,f11.8,A)') '"hy": ', hy, ','
  
          write(uid,'(A,f11.8,A)') '"uvel": ', uvel, ','
  
          write(uid,'(A,f11.8,A)') '"cfl": ', cfl, ','
          write(uid,'(A,f11.8,A)') '"dt": ', dt, ','
          write(uid,'(A,f11.8,A)') '"tf": ', tf, ','
          write(uid,'(A,I8,A)') '"nts": ', nts, ','
  
          write(uid,'(A,I4)') '"niter": ', niter
  
  
  
  
negative.code:         integer, parameter :: sizeofreal=4
          integer :: uid
          logical :: exist
          character*9 :: filename
  
          !! Save the full variables
          WRITE(filename,'(I0.5,A)') it,'.raw'
          open(newunit=uid,file=write_dir//filename, form='unformatted', &
               access='direct',recl=nnx*nny*sizeofreal*2)
          write(uid,rec=1) u,v
          close(uid)
  