anchor.diff:          end if
   
           ! Check for case in which the model iteration has passed a write time
  -       if(tt/tseas >= dble(Ni)/dble(N)) then
  +        if(tt/tseas >= dble(Ni)/dble(N)) then
   
               ! Time interpolation weights
               ! Time at the output time minus the time for x0 over the time window, weight for x0, etc
  
positive.code: ! controls the iterative stepping between the two input model outputs
  dstep = 1.d0/dble(iter)
  dtmin = dstep*tseas 
  
  flag = 0
  
  
negative.code:                                                            x0, y0, z0, x1, y1, z1, &
                                                             tt, tss, ts, rwn, rwp
  integer                                                 :: ntrac, niter, ia, ja, ka, &
                                                             iam, ib, jb, kb, errCode, &
                                                             Ni, wrapx, wrapy
  real*8,     parameter                                   :: UNDEF=1.d20
  
  real*8, dimension(6,2)                                    :: upr
  