anchor.diff:          deallocate(a_prestress, nvec_i)
           deallocate(fineness, nnode, nelm)
           deallocate(cap_n_start, cap_n_end, cap_e_start, cap_e_end)
  -        deallocate(rad, xcenter, ycenter, zcenter, &
  +        deallocate(xcenter, ycenter, zcenter, &
                xcenterold, ycenterold, zcenterold)
  +     deallocate(ellipn1, ellipn2, ellipn3)
  +     deallocate(ellipa, ellipb, ellipc)
  +     if (nsph > 0) deallocate(rad) 
        end if
      END PROGRAM cell
  -       !**********************************************************************
  +!**********************************************************************
  
positive.code:   INTEGER KLOK,KLOK1,KLOKEND,NSTEP
    double precision :: T,H,h64,TD,VSC,TIME,RHO,PI,RADX,FOSTAR
    double precision, allocatable :: rad(:), a_prestress(:)
    double precision, allocatable :: xcenter(:), ycenter(:), &
         zcenter(:)
    double precision, allocatable :: xcenterold(:), ycenterold(:), & 
         zcenterold(:)
    integer, allocatable :: nvec_i(:,:)
    integer ncap
    double precision ::  LCUBE,NU,MU,MASS,LENGTH
    !     Velocities are always expressed in program units
    double COMPLEX,ALLOCATABLE :: UR(:,:,:),VR(:,:,:),WR(:,:,:)
  
negative.code:         xcenterold(i) = xcenter(i)
          ycenterold(i) = ycenter(i)
          zcenterold(i) = zcenter(i)
          ! Save the n vectors, the six vectors from the center of the capsule
          ! to the six tracer points.
          call dumpnvec(klok, xfn(:,cap_n_start(i):cap_n_end(i)), xcenter(i), &
               ycenter(i), zcenter(i), nvec_i(:, i), i)
          call shape(h64,klok, i, xfn(1:3, cap_n_start(i): &
  