anchor.diff:  c
         intrinsic sqrt
   c
  -c...  user-defined functions
  +c...  external routines
   c
         double precision sprod
         external sprod
  
positive.code: c*      write(6,*) "Hello from invar_f!"
  c
        sinv1=stn(1)+stn(2)+stn(3)
        smean=sinv1/3.0d0
        sdev(1)=stn(1)-smean
        sdev(2)=stn(2)-smean
        sdev(3)=stn(3)-smean
  
negative.code:       end
  c
  c version
  c $Id: invar.f,v 1.2 2004/07/07 15:31:08 willic3 Exp $
  c
  c Generated automatically by Fortran77Mill on Wed May 21 14:15:03 2003
  c
  