anchor.diff:  c
         include "implicit.inc"
   c
  -c...  subroutine arguments
  -c
  -      integer nstr,nstate,nddmat,nprop
  +c...  parameter definitions
  +c
  +      include "ndimens.inc"
  +      include "rconsts.inc"
  +c
  +c...  subroutine arguments
  +c
  +      integer nstate,nprop
         integer iddmat(nstr,nstr)
         double precision state(nstr,nstate),dmat(nddmat),prop(nprop)
         logical matchg
   c
  -c...  parameters
  -c
  -      integer inc
  -      double precision zero
  -      parameter(inc=1,zero=0.0d0)
  -c
   c...  local variables
   c
         integer i,j
  
positive.code:       call fill(dmat,zero,nddmat)
        e=prop(2)
        pr=prop(3)
        pr1=1.0d0-pr
        pr2=1.0d0+pr
        pr3=1.0d0-2.0d0*pr
        fac=e/(pr2*pr3)
        dd=pr1*fac
        od=pr*fac
        ss=0.5d0*pr3*fac
        do i=1,3
          dmat(iddmat(i,i))=dd
          dmat(iddmat(i+3,i+3))=ss
  
negative.code: c       
  
  c version
  c $Id: mat_1.f,v 1.2 2004/06/18 14:40:55 willic3 Exp $
  
  c Generated automatically by Fortran77Mill on Tue May 18 14:18:50 2004
  
  