anchor.diff:          write(uid,'(A,f11.8,A)') '"hy": ', hy, ','
   
           write(uid,'(A,f11.8,A)') '"uvel": ', uvel, ','
  +        write(uid,'(A,f11.9,A)') '"dpg": ', dpg, ','
   
           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
  +        write(uid,'(A,I4,A)') '"niter": ', niter, ','
  +
  +        write(uid,'(A,I4,A)') '"write_freq": ', write_freq, ','
  +        write(uid,'(A,A,A)') '"write_dir": "', write_dir, '",'
  +
  +        write(uid,'(A,I4,A)') '"ndisks": ', ndisks, ','
  +        write(uid,'(A,A,A)') '"adisk_fname": "', adisk_fname, '",'
  +        write(uid,'(A,I4,A)') '"r": ', r, ','
  +        write(uid,'(A,f11.8,A)') '"alpha": ', alpha
   
   
   
  
positive.code: 
          integer :: uid
  
          character*22 :: filename
          filename = './data/parameters.json'
          open(newunit=uid,file=filename)
  
          write(uid,'(A)') '{'
  
  
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)
  