anchor.diff:    do i =1,ngx
        do j = 1,ngy
           do k = 0,ngz-1
  -           write(25,*) i,j,k, dreal(wr(i,j,k)), dimag(wr(i,j,k))
  +           write(25,'(3(i4),2(es24.17))') i,j,k, dreal(wr(i,j,k)), &
  +                dimag(wr(i,j,k))
           end do
        end do
     end do
  
positive.code: 
    call makefilename('TaylorDF__', cap_i,'.txt',strfname)
    open(204, file=strfname, access='append')
    write(204,*) clock, DF
    close(204)
  
    return
  
negative.code:   do i =1,ngx
       do j = 1,ngy
          do k = 0,ngz-1
             write(25,'(3(i4),x,4(es24.17))') i,j,k, dreal(ur(i,j,k)), &
                  dreal(vr(i,j,k)), &
                  dreal(wr(i,j,k)), dreal(pr(i,j,k))
          end do
       end do
  