anchor.diff:          enddo ! end collect total force
   
           acc = totalForce(1) / objMass
  +        if (myrank .eq. master) then
  +          write(*,*) "projectile acceleration is:", acc
  +        endif
   
           do i = 1,numnp
             if ( (ibits(iBC(i),14,3) .eq. 7) .and.
  
positive.code:       integer   answer
  c
        if (elasFDC .gt. 0) then
          write(*,*) "use Force-driven case:", elasFDC
          casenumber = elasFDC
        else
          write(*,*) "Please use Force-driven as Mesh Elas Model"
        endif
  c
  c.... Update BC value based on geom and iBC
  
negative.code:         xtsl     = 2.0000000000000e-4
          dyn_org  = 2.6250000000000e-3 + DBLE(lstep) * xtsl
          dyn_lnt  = 2.0000000000000e-3 * 0.9**lstep
          if (myrank .eq. master) then
            write(*,*) "current lstep:", lstep, "dyn_org:", dyn_org
          endif
          rotf     = 90.0
          shrkfactor = 2.0
          do i = 1,numnp
  