anchor.diff:       endif
     enddo
     !     printout theta and DF based on dzy
  +  ! The points on x=0 are passed to dzy, which does some output on the DF
  +  ! seen in the y-z plane, and the angle of inclination from the z-axis
     call dzy(zy,icount,klok, cap_i)
     DEALLOCATE (zy,XFP)
     return 
  
positive.code: !**********************************************************************
  subroutine wlindump(wlin,clock)
    implicit none
    integer, parameter :: lxng=$lngx$,lyng=$lngy$,lzng=$lngz$
    integer, parameter :: ngx=2**lxng,ngy=2**lyng,ngz=2**lzng
    integer i, clock
  
negative.code: !**********************************************************************
  subroutine padzeros(str)
    implicit none
    ! Takes a string and replaces any space characters in it with zeros.
    integer i
    character(len=*) str
    do i=1,len(str)
  