anchor.diff:    call random_seed(put=seed)
   
     write(argv,'(A,I0.3,A,I0.3,A,I0.3)') odir//"/EE",realization,"_n",beami,"_v",beamj
  -  open(newunit=uEE,file=trim(argv), status='unknown',action='write')
  -  print *,'writing to',trim(argv)
  +  open(newunit=uEE,file=trim(argv), status='new',action='write',access='stream')
  +  print *,'writing to ',trim(argv)
   
     write(argv,'(A,I0.3,A,I0.3,A,I0.3)') odir//"/nn",realization, "_n" ,beami, "_v",beamj
  -  open(newunit=uNN,file=trim(argv), status='unknown',action='write')
  -  print *,'writing to',trim(argv)
  +  open(newunit=uNN,file=trim(argv), status='new',action='write',access='stream')
  +  print *,'writing to ',trim(argv)
   
   !   main loops
   
  
positive.code:       vv(c2,pp,2)=vv(c1,pp,2)+(kv1(2)+2*kv2(2)+2*kv3(2)+kv4(2))/6+SSn(2)*Tstep;
        nn(c2,pp,1)=nn(c1,pp,1)+(kn1(1)+2*kn2(1)+2*kn3(1)+kn4(1))/6;
        nn(c2,pp,2)=nn(c1,pp,2)+(kn1(2)+2*kn2(2)+2*kn3(2)+kn4(2))/6;
        nn(c2,N/2,1)=0.0;
        nn(c2,N/2,2)=0.0;
  
        if (pp>=1) then
          nn(c2,N-pp,1)=nn(c2,pp,1);
  
negative.code:   write(argv,'(A,I0.3,A,I0.3)')  odir//"/parameters_n" , beami, "_v" , beamj
    ofn = trim(argv)
  
    open(newunit=u, file=ofn, status='new',action='write', access='stream')
    write(u) parameters, SEED
    close(u)
    print *, "Wrote parameters to ",ofn
  
    do ii=1,N
      p(ii)=ii-N/2
  