anchor.diff:    call padzeros(strfname)
     open(25, file=strfname, access='append')
     DO i=1, nnode
  -     IF((XFN(1,i) >fngx/2.d0-1.d0).AND.(XFN(1,i) &
  -          <fngx/2.d0+1.d0)) THEN
  -        write(25,'(es24.17,x,es24.17)') XFN(2,i),XFN(3,i)
  +     IF((XFN(3,i) > 0.5d0*(fngz-1.d0) - 0.5d0) .AND.(XFN(3,i) &
  +          < 0.5d0*(fngz-1.d0) + 0.5d0)) THEN
  +        write(25,'(es24.17,x,es24.17)') XFN(1,i),XFN(2,i)
        ENDIF
     ENDDO
     close(25)
  
positive.code:         zy(2,icount) = yc
       endif
    enddo
    !     printout profile in the zy plane and value of Dxy
    call dzy(zy,icount,klok, cap_i)
    DEALLOCATE (zy,XFP)
    return 
  
negative.code: ! that plots zy in shape should be used.
  subroutine profile(cap_i, xfn, clock)
    implicit none
    integer, parameter :: lngz = $lngx$
    integer, parameter :: ngz = 2**lngz - 1
    double precision, parameter :: fngz = ngz
    character*21 strfname
    integer cap_i, clock, i
    double precision :: xfn(:, :)
  