$
$****************************************
$
&dimen -dimen feet kips
&device -auxin lmu.dat
$
$****************************************
$
inmodel
$
$****************************************
$
medit
  ~frict  gspr         x 1e3 1e4 y 1  1 z 1 2 -friction .1
  ~x_damp gspr         x 1e3 1e4 y 1  1 z 1 2 -x_damp 1e5
  ~y_damp gspr         x 1e3 1e4 y 1  1 z 1 2 -y_damp 1e5
  ~z_damp gspr         x 1e3 1e4 y 1  1 z 1 2 -z_damp 1e5
$
  connector frict     ~frict  *1
  connector x_damp    ~x_damp  *1
  connector y_damp    ~y_damp  *1
  connector z_damp    ~z_damp  *1
$
  &set gbs_hor     = 1e10
  &set gbs_hor_py  =  0.0     0.000  \
                      250     0.250  \
                     1000     0.500  \
                     2000     0.750  \
                     3000     1.000
  &set gbs_ver_py  =  0.0     0.000  \
                     1000     0.200  \
                     2070     0.400  \
                     2700     0.500  \
                     3400     0.600  \
                    83400     1.600
  ~gspr   gspr comp 1000 1000 0 x 1e9 y %gbs_hor z %gbs_hor \
                         -x_py   %gbs_ver_py          \
                         -y_py   %gbs_hor_py          \
                         -z_py   %gbs_hor_py 
$
  connector gspr ~gspr -euler -z  *1
end
$
$****************************************
$
&connect @damp frict -inactive
&type x_py &class(x_py ~gspr) 
&type y_py &class(y_py ~gspr) 
&cutype x = 1 z = +1
&instat -loc 1 0 +1 
&stat f_c
&loop x ( -1 -.75 -.5 -.25 0 .25 .5  .75 1. )
      &cutype x = %x z = 0
      &instat -loc %x 0 -.01
      &stat f_c
&endloop
$
$****************************************
$
&connect @damp frict -active gspr -inactive
&cutype Move -2 in X Vel 10 10 10
&instat -move -2 -vel 10 10 10
&stat f_c
$
$****************************************
$
&fini