anchor.diff:    do i =1,ngx
        do j = 1,ngy
           do k = 0,ngz-1
  -           write(25,*) i,j,k, dreal(ur(i,j,k)), dreal(vr(i,j,k)), &
  +           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
  
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 j = 1, ngy
       k = ngz/2-1
       !            do k = 0, ngz-1
       write(25,'(3(i5x),x,e24.17,x,e24.17)') i,j,k, dble(wr(i,j,k)), dimag(wr(i,j,k))
       !            end do
    end do
    !      end do
  