anchor.diff:        ntimdat(1)=nstep
         maxitp=maxit(ii)
         ntimdat(2)=maxitp
  -      maxitcp=maxitc(ii)
  -      ntimdat(3)=maxitcp
  +      ntdinitp=ntdinit(ii)
  +      ntimdat(3)=ntdinitp
         lgdefp=lgdef(ii)
         ntimdat(4)=lgdefp
         ibbarp=ibbar(ii)
  
positive.code: c~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  c
  c
        subroutine const(maxstp,delt,alfa,maxit,maxitc,lgdef,ibbar,utol,
       & ftol,etol,itmax,nintg,i,naxstp,nfirst,rtimdat,deltp,alfap,
       & ntimdat,nstep,maxitp,maxitcp,lgdefp,ibbarp,itmaxp,gtol)
  c
  c...this subroutine defines the constants for each time step group.
  c
  
negative.code:       end
  c
  c version
  c $Id: const.f,v 1.2 2004/06/23 20:52:09 willic3 Exp $
  c
  c Generated automatically by Fortran77Mill on Wed May 21 14:15:03 2003
  c
  