anchor.diff:    integer i
     open(25,file=strfname,status='unknown')
     do i = 1,size(xfn,2)
  -     write(25,*) xfn(1,i),xfn(2,i),xfn(3,i)
  +     write(25,'(es24.17,3x,es24.17,3x,es24.17)') xfn(1,i),xfn(2,i),xfn(3,i)
     end do
     close(25)
   end subroutine saveallsolid
  
positive.code:   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))
          end do
       end do
    end do
  
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
  