#	
#
  head 2 "Test of different wheel/rail couplings  cwr_coupl= $cwr_coupl"
  head 3 "axl_111 pure long.creep; axl_112 pure lat.creep; axl_113 pure spin"
  head 4 "Cylindric wheel profile  mu= $mu_"

##
## Analyse specific input data
##
  idebug= 1
  fresp_param= Fourier_CG1  1.e-4  fstart .2   fstop 20.  fstep -1.05
  modal_param= Schur_fact1  1.e-4
  quasi_param= Damped_Tens  1.e-4 0.1   20    2.
#
# tsim_param = heun  tstart 0.  tstop  500  tstep 1.00  tout 1.00
# tsim_param = heun_u .5  tstart 0.  tstop  2.50  tstep 0.001  tout .005
# tsim_param = heun       tstart 0.  tstop  2.50  tstep 0.005  tout .005
# tsim_param = heun       tstart 0.  tstop  5.00  tstep 0.005  tout .005
  tsim_param = heun_u .05 .2e-3  tstart 0.  tstop  3.00  tstep 0.0001 tout .005

##
##  Track irregularities and designed track geometry
##
  func	  intpl_r   lat_trac	  -200    0.   0    0.
  func	  intpl_r   spv_trac	  -200 1435.   0 1435.
  func	  intpl_r   vert_trac	  -200    0.   0    0.
  func	  intpl_r   fi_trac	  -200    0.   0    0.

# func	  intpl_r   krok_R  -200 0  10 0   20  .005  100  .005
  func	  intpl_r   krok_R  -200 0   0 0 #.001 .005  100  .005
  func	  intpl_r   krok_h  -200 0   0 0    
  func	  intpl_r   krok_z  -200 0   0 0    

##
##  Wheel-rail geometry
##
# insert file $genkpf/kona_0.000_uf.kpfr
# in_substruct kpf_kona_0.000 [ 1 ]    # Same wheel-rail geometry in all wheels vehicle #1
#
  insert file $genkpf/kpf_rkona.kpfr
  in_substruct kpf_rkona [ 1 0. ]    # Same wheel-rail geometry in all wheels vehicle #1

##
##  Define constants for the wheel-rail-coupling substructure wr_coupl
##
  func	const_block
   Vo = 10.   # Speed of the vehicles
   mu_= 0.5   # Friction coefficient, all wheels in the train
   ro_= 0.45  # Wheel radius, all wheels in the train
   bo = 0.75  # Semi distance between the wheel profile origins, all wheels in the train

   YMtrac= 1e-3        # From mm to m;   Scale factor 1.00; Lateral irregularities
   ZMtrac= 1e-3        # From mm to m;   Scale factor 1.00; Vertical irregularities
   GMtrac= 1e-3        # From mm to m;   Scale factor 1.00; Gauge irregularities
   CMtrac= 1e-3/(2*bo) # From mm to rad; Scale factor 1.00; Cant irregularities

   gauge_average= 1435.    # Average gauge of the track
   gauge_dev    = 0.000    # Deviation in gauge
  end_block

##
##  Create the wheelsets.
##  ---------------------
  lsys	 e_abs_bendrf  esys  Vo  0.  krok_R krok_h krok_z
                                     4 4    4 4    4 4

  substruct create_axles [
   lsys  l_local      lsa_$1  esys     0. 0. 0.
   mass  m_rigid_6    axl_$1  lsa_$1   0. 0. -ro_   1887 1887 1887  720  100  720
   mass  m_rigid_6    trc_$1  lsa_$1   0. 0. 0.    1.e3 1.e3 1.e3 1.e3 1.e3 1.e3
   mass  fixpoint_6   grd_$1  lsa_$1   0. 0. 0.
   body  w_set_lsys   axl_$1  lsa_$1   ro_ bo
   body  figure_lsys  trc_$1  lsa_$1   1 5
                                       .5 1 0  .5 -1 0   -.5 -1 0  -.5 1 0  .5 1 0
#
   constr fix_rigid_1  axl_$1  x  0.
   constr fix_rigid_1  axl_$1  y  0.010
#  constr fix_rigid_1  axl_$1  z  0.
#  constr fix_rigid_1  axl_$1  f  0.
   constr fix_free_1   axl_$1.vk  0.
   constr fix_rigid_1  axl_$1  p  0.
   constr fix_rigid_1  trc_$1  y  0.
#
   force rel_lsys1 vload_$1  axl_$1  0. 0. -ro_  0. 0. 200e3-1887*9.81  0. 0. 0.
#  coupl p_lin     czbap= 0.   50e3
#  coupl c         czba_$1r  axl_$1  0  bo -ro_ 
#                            grd_$1  0  bo -ro_ 
#                            czbap   esys z
#  coupl c         czba_$1l  axl_$1  0 -bo -ro_ 
#                            grd_$1  0 -bo -ro_ 
#                            czbap   esys z
  ]
#
  in_substruct create_axles [ 111 ]
  in_substruct create_axles [ 112 ]
  in_substruct create_axles [ 113 ]
#
  mass	fixpoint_6   car_1    lsa_111	0. 0. 0.
#
# coupl k_lin  ip_coupl axl_111 x 0. 0. -ro_
#                       trc_111 x 0. 0. -ro_   17

# coupl p_kfrkc  p_name   +-`F0 +-`ke +-`Ffm +-`x2 +-`kv +-`c  
# coupl p_kfrkc  p_rubb       0   100     10   .01    0    1  
# coupl kfrkc   c_rubb  axl_111     0      0   -ro_
#                       trc_111     0      0   -ro_
#                       p_rubb    esys     y
#
#coupl p_lin  stifx        0.   500e3  
#coupl p_lin  stify        0.   500e3  
#coupl p_lin  stifz    -100e3     1e6  
#coupl k      stiff   axl_111     0      0   oper_init( -ro_ - .1 )
#                     trc_111     0      0   oper_init( -ro_ + .1 )
#                     stifx     esys     y

#force rel_lsys1  ipforc  axl_111  0. 0. 0.  100 0 0   0 0 0

# coupl k_air3_exp  p_airip  axl_111     0      0   -ro_            
#                            trc_111     0      0   -ro_            
#                            esys   m
#                            stifx       0.
#                            stify       0.
#                            stifz                                 
#                            10   .01            10   .01
#                             0    1 
#                             0    1    1  10

##
##  Define pure creepage on the three different axles
##  -------------------------------------------------
# func operp axl_111.vx= -.05
# func operp axl_111.vx= time * -.25
  func operp axl_111.vk= time * -.25
# func operp axl_112.vy=        -.25
# func operp axl_112.vy= time * -.25
  func operp axl_112.p = time *  .1
# func operp axl_113.vp= time * -1.e+2
# func operp axl_113.vx= time * -.75e+2
# func operp axl_113.vp= time * -50
# func operp axl_113.vx= time * -37.5
  func operp axl_113.vp= time * -3.4
  func operp axl_113.vx= time * -2.55
# func operp axl_113.vp= 0.
# func operp axl_113.vx= 0.
  constr fix_rigid_1  axl_111 f  0.
  constr fix_rigid_1  axl_112 f  0.
  constr fix_rigid_1  axl_113 f  0.

# func operp axl_111.vx=  500 * .5e-3
# func operp axl_112.vy=  500 * .5e-3
# func operp axl_113.vp=  500 * 1.e-1
# func operp axl_113.vx=  500 * .75e-1

##
##	Calculate creepage and creep-forces
##	-----------------------------------
# func char cwr_coupl= wr_coupl_nl1
  func char cwr_coupl= wr_coupl_npol1
# func char cwr_coupl= wr_coupl_ne1
# func char cwr_coupl= wr_coupl_nr1
# func char cwr_coupl= wr_coupl_nra1
# func char cwr_coupl= wr_coupl_pe1
# func char cwr_coupl= wr_coupl_pr1
# func char cwr_coupl= wr_coupl_pra1
# func char cwr_coupl= wr_coupl_lt1

  if_then_char_init cwr_coupl .eq. wr_coupl_nl1 .or.
                    cwr_coupl .eq. wr_coupl_npol1 .or.
                    cwr_coupl .eq. wr_coupl_ne1 .or.
                    cwr_coupl .eq. wr_coupl_nr1 .or.
                    cwr_coupl .eq. wr_coupl_nra1  # {
#
   func const knwr.F0_= 100e3 # all wheels
#
   func  intpl_r knwr_    # Stiffness Tread-Rail all wheels
    -1.57     42e6
     0.      180e6
     1.57     42e6
   func  intpl_r cnwr_    # Damping Tread-Rail all wheels
    -1.57    60e3
     0.     250e3
     1.57    60e3
   func  intpl_r knfr_    # Stiffness Flange-Rail all wheels
    -1.57     42e6
     0.      180e6
     1.57     42e6
   func  intpl_r cnfr_    # Damping Flange-Rail all wheels
    -1.57    60e3
     0.     250e3
     1.57    60e3
  endif # }
#
  if_then_char_init cwr_coupl .eq. wr_coupl_nl1
   insert file $gensys/calc/insert_files/wr_coupl_nl1.ins
   in_substruct wr_coupl [ 1 111 ]
   in_substruct wr_coupl [ 1 112 ]
   in_substruct wr_coupl [ 1 113 ]
  endif
  if_then_char_init cwr_coupl .eq. wr_coupl_npol1
   insert file $gensys/calc/insert_files/wr_coupl_npol1.ins
   in_substruct wr_coupl [ 1 111 ]
   in_substruct wr_coupl [ 1 112 ]
   in_substruct wr_coupl [ 1 113 ]
  endif
  if_then_char_init cwr_coupl .eq. wr_coupl_ne1
   insert file $gensys/calc/insert_files/wr_coupl_ne1.ins
   in_substruct wr_coupl [ 1 111 ]
   in_substruct wr_coupl [ 1 112 ]
   in_substruct wr_coupl [ 1 113 ]
  endif
  if_then_char_init cwr_coupl .eq. wr_coupl_nr1
   insert file $gensys/calc/insert_files/wr_coupl_nr1.ins
   in_substruct wr_coupl [ 1 111 ]
   in_substruct wr_coupl [ 1 112 ]
   in_substruct wr_coupl [ 1 113 ]
  endif
  if_then_char_init cwr_coupl .eq. wr_coupl_nra1
   insert file $gensys/calc/insert_files/wr_coupl_nra1.ins
   in_substruct wr_coupl [ 1 111 ]
   in_substruct wr_coupl [ 1 112 ]
   in_substruct wr_coupl [ 1 113 ]
  endif
##
##
##
  if_then_char_init cwr_coupl .eq. wr_coupl_pe1  .or.
                    cwr_coupl .eq. wr_coupl_pr1  .or.
                    cwr_coupl .eq. wr_coupl_pra1 .or.
                    cwr_coupl .eq. wr_coupl_lt1         # {
#
   func const knwr.F0_= 100e3 # all wheels
#
   func  intpl_r knwr_      # Stiffness Tread-Rail all wheels
   -1.57     600e6
    0.       600e6
    1.57     600e6
#  func  intpl_r cnwr_      # Damping Tread-Rail all wheels
#  -1.57     600e3
#   0.       600e3
#   1.57     600e3
   func  intpl_r knfr_      # Stiffness Flange-Rail all wheels
   -1.57     600e6
    0.       600e6
    1.57     600e6
#  func  intpl_r cnfr_      # Damping Flange-Rail all wheels
#  -1.57     600e3
#   0.       600e3
#   1.57     600e3
   func const_block
    pknwr   = 1             # Linjär styvhet i kontaktpunkten
    kyrt    =  42e6*100.
    kzrt    = 180e6*100.
    kzrt.F0_= knwr.F0_
    cyrt    = 205e3*10#31.6
    czrt    = 250e3*10#31.6
   end_block
  endif # }
  if_then_char_init cwr_coupl .eq. wr_coupl_lt1         # {
   func  intpl_r cnwr_    # Damping Tread-Rail all wheels
    -1.57    60e3
     0.     250e3
     1.57    60e3
   func  intpl_r cnfr_    # Damping Flange-Rail all wheels
    -1.57    60e3
     0.     250e3
     1.57    60e3
  endif # }

##
  if_then_char_init cwr_coupl .eq. wr_coupl_pe1
   insert file $gensys/calc/insert_files/wr_coupl_pe1.ins
   in_substruct wr_coupl [ 1 111 ]
   in_substruct wr_coupl [ 1 112 ]
   in_substruct wr_coupl [ 1 113 ]
  endif
  if_then_char_init cwr_coupl .eq. wr_coupl_pr1
   insert file $gensys/calc/insert_files/wr_coupl_pr1.ins
   in_substruct wr_coupl [ 1 111 ]
   in_substruct wr_coupl [ 1 112 ]
   in_substruct wr_coupl [ 1 113 ]
  endif
  if_then_char_init cwr_coupl .eq. wr_coupl_pra1
   insert file $gensys/calc/insert_files/wr_coupl_pra1.ins
   in_substruct wr_coupl [ 1 111 ]
   in_substruct wr_coupl [ 1 112 ]
   in_substruct wr_coupl [ 1 113 ]
  endif

  substruct wr_coupl_lt1 [
   mass m_rigid_6f  ral_$1r lsa_$1  0.  bo  0.    0. 0. 0.  0. 0. 0.
   mass m_rigid_6f  ral_$1l lsa_$1  0. -bo  0.    0. 0. 0.  0. 0. 0.
   constr fix_rigid_1  ral_$1r x 0.  constr fix_rigid_1  ral_$1r y 0.  constr fix_rigid_1  ral_$1r z 0.
   constr fix_rigid_1  ral_$1r f 0.  constr fix_rigid_1  ral_$1r k 0.  constr fix_rigid_1  ral_$1r p 0.
   constr fix_rigid_1  ral_$1l x 0.  constr fix_rigid_1  ral_$1l y 0.  constr fix_rigid_1  ral_$1l z 0.
   constr fix_rigid_1  ral_$1l f 0.  constr fix_rigid_1  ral_$1l k 0.  constr fix_rigid_1  ral_$1l p 0.
#
   constr fix_free_1  axl_$1.k  0.
   constr fix_rigid_1 trc_$1 x  0.  constr fix_rigid_1 trc_$1 z  0.
   constr fix_rigid_1 trc_$1 f  0.  constr fix_rigid_1 trc_$1 k  0.
   constr fix_rigid_1 trc_$1 p  0.
#
   func const axl_$1.vk=  -Vo/ro_ 
#
   func  operp whe_$1r.x= axl_$1.x - bo * axl_$1.p
   func  operp whe_$1l.x= axl_$1.x + bo * axl_$1.p
#
   func  operp whe_$1r.vx= axl_$1.vx - bo * axl_$1.vp
   func  operp whe_$1l.vx= axl_$1.vx + bo * axl_$1.vp
#
   coupl creep_lookuptable_1 cpt_$1r
                axl_$1  0.0        bo  0.0
                ral_$1r whe_$1r.x  bo  0.0  esys m
                  0.    # tral$1r.y        # Lateral position track irregularity
                  0.    # tral$1r.z        # Vertical position track irregularity
                  0.    # tral$1r.vy       # Lateral velocity track irregularity
                  0.    # tral$1r.vz       # Vertical velocity track irregularity
                  1.    # mulfact_nux_$1r  # Longitudinal creep relaxation due to contaminated surfaces
                  1.    # mulfact_nuy_$1r  # Lateral      creep relaxation due to contaminated surfaces
                  1.    # mulfact_spin_$1r # Spin         creep relaxation due to contaminated surfaces
                  cpt_1.zfn            # Wheel lift geometric function
                  cpt_1.drfn           # Wheel radius geometric function
                  cpt_1.gamfn          # Contact angle geometric function
                  cpt_1.rofn           # Lateral curvature geometric function
                  cpt_1.poswfn         # Lateral position of contact point on wheel 
                  cpt_1.posrfn         # Lateral position of contact point on rail
                  knwr.F0_$1r          # Wheel/rail prestess force in normal direction of contact surface
                  knwr_$1r             # Wheel/rail stiffness in normal direction of contact surface
                  cnwr_$1r             # Wheel/rail damping   in normal direction of contact surface
                  2.05e11              # The combined modulus of elasticity in body1 and body2
                  .3                   # The combined Poisson's ratio in body1 and body2
                  mu_$1r               # The coefficient of friction in the contact surface
                  ro_$1r               # Wheel rolling radius, will vary for OOR-wheels
#
   coupl creep_lookuptable_1 cpt_$1l
                axl_$1  0.0       -bo  0.0
                ral_$1l whe_$1l.x -bo  0.0  esys m
                  0.    # tral$1l.y        # Lateral position track irregularity
                  0.    # tral$1l.z        # Vertical position track irregularity
                  0.    # tral$1l.vy       # Lateral velocity track irregularity
                  0.    # tral$1l.vz       # Vertical velocity track irregularity
                  1.    # mulfact_nux_$1l  # Longitudinal creep relaxation due to contaminated surfaces
                  1.    # mulfact_nuy_$1l  # Lateral      creep relaxation due to contaminated surfaces
                  1.    # mulfact_spin_$1l # Spin         creep relaxation due to contaminated surfaces
                  cpt_1.zfn            # Wheel lift geometric function
                  cpt_1.drfn           # Wheel radius geometric function
                  cpt_1.gamfn          # Contact angle geometric function
                  cpt_1.rofn           # Lateral curvature geometric function
                  cpt_1.poswfn         # Lateral position of contact point on wheel 
                  cpt_1.posrfn         # Lateral position of contact point on rail
                  knwr.F0_$1l          # Wheel/rail prestess force in normal direction of contact surface
                  knwr_$1l             # Wheel/rail stiffness in normal direction of contact surface
                  cnwr_$1l             # Wheel/rail damping   in normal direction of contact surface
                  2.05e11              # The combined modulus of elasticity in body1 and body2
                  .3                   # The combined Poisson's ratio in body1 and body2
                  mu_$1l               # The coefficient of friction in the contact surface
                  ro_$1l               # Wheel rolling radius, will vary for OOR-wheels
  ]
  if_then_char_init cwr_coupl .eq. wr_coupl_lt1
   in_substruct wr_coupl_lt1 [ 111 ]
   in_substruct wr_coupl_lt1 [ 112 ]
   in_substruct wr_coupl_lt1 [ 113 ]
  endif

##
##  Modal parameters of the wheelsets.
##  ----------------------------------
# mass  m_flex_1 axl_111  70   0.10               ## Torsionsform
# coupl m_flex_1 cpt_111r end_2  0.0636396 0. 0.  # 0. 0.14142 0.
# coupl m_flex_1 cpt_111l end_2  -.0636396 0. 0.  # 0. 0.14142 0.
# s_var sngl axl_111.c1
# s_var sngl axl_111.v1

##
##  Store variables for postprocessing
##  ----------------------------------
#
  func print06_init  cpt_111r.a/b  # Print the initial value of a variable to standard output
#
  s_var scalar_0 mu_
#
  s_var sngl ro_
  s_var sngl bo
  s_var sngl Vo
#
  s_var var_0 esys.x  
  s_var var_0 esys.y  
  s_var var_0 esys.z  
  s_var var_0 esys.f  
  s_var var_0 esys.k  
  s_var var_0 esys.p  
#
  s_var var_0 esys.vx 
  s_var var_0 esys.vy 
  s_var var_0 esys.vz 
  s_var var_0 esys.vf 
  s_var var_0 esys.vk 
  s_var var_0 esys.vp 
#
  s_var var_0 esys.ax 
  s_var var_0 esys.ay 
  s_var var_0 esys.az 
  s_var var_0 esys.af 
  s_var var_0 esys.ak 
  s_var var_0 esys.ap 
#
  substruct store_axle [
   s_var var_0 axl_$1.x
   s_var var_0 axl_$1.y
   s_var var_0 axl_$1.z
   s_var var_0 axl_$1.f
   s_var var_0 axl_$1.k
   s_var var_0 axl_$1.p
#
   s_var var_0 axl_$1.vx
   s_var var_0 axl_$1.vy
   s_var var_0 axl_$1.vz
   s_var var_0 axl_$1.vf
   s_var var_0 axl_$1.vk
   s_var var_0 axl_$1.vp
#
   s_var var_0 axl_$1.Fx
   s_var var_0 axl_$1.Fy
   s_var var_0 axl_$1.Fz
   s_var var_0 axl_$1.Mf
   s_var var_0 axl_$1.Mk
   s_var var_0 axl_$1.Mp
#
   s_var var_0 axl_$1.Ax
   s_var var_0 axl_$1.Ay
   s_var var_0 axl_$1.Az
   s_var var_0 axl_$1.Af
   s_var var_0 axl_$1.Ak
   s_var var_0 axl_$1.Ap
#
  s_var var_0 trc_$1.x     s_var var_0 trc_$1.y     s_var var_0 trc_$1.z
  s_var var_0 trc_$1.f     s_var var_0 trc_$1.k     s_var var_0 trc_$1.p
  s_var var_0 trc_$1.vx    s_var var_0 trc_$1.vy    s_var var_0 trc_$1.vz
  s_var var_0 trc_$1.vf    s_var var_0 trc_$1.vk    s_var var_0 trc_$1.vp
  ]
#
  substruct store_wheel [
   s_var var_0 ral_$1.y
   s_var var_0 ral_$1.z
   s_var var_0 ral_$1.vy
   s_var var_0 ral_$1.vz
#
   s_var var_0 cp_$1.eta   
#
   s_var var_0 tral$1.y    
   s_var var_0 tral$1.z    
   s_var var_0 tral$1.vy   
   s_var var_0 tral$1.vz   

   s_var var_0 cpt_$1.z    
   s_var var_0 cpt_$1.dr   
   s_var var_0 cpt_$1.gam  
   s_var var_0 cpt_$1.c    
   s_var var_0 cpt_$1.a/b  
   s_var var_0 cpt_$1.posw
   s_var var_0 cpt_$1.bo
#
   s_var var_0 cpt_$1.nux  
   s_var var_0 cpt_$1.nuy  
   s_var var_0 cpt_$1.spin 
#
   s_var var_0 cpt_$1.vux
   s_var var_0 cpt_$1.vuy
   s_var var_0 cpt_$1.vuz
   s_var var_0 cpt_$1.Vref
#
   s_var var_0 cpt_$1.Fn
   s_var var_0 cpt_$1.Fny
   s_var var_0 cpt_$1.Fx
   s_var var_0 cpt_$1.Fy
   s_var var_0 cpt_$1.Fz
   s_var var_0 cpt_$1.F1x
   s_var var_0 cpt_$1.F1y
   s_var var_0 cpt_$1.F1z
   s_var var_0 cpt_$1.F2x
   s_var var_0 cpt_$1.F2y
   s_var var_0 cpt_$1.F2z
   s_var var_0 cpt_$1.M1f
   s_var var_0 cpt_$1.M1k
   s_var var_0 cpt_$1.M1p
   s_var var_0 cpt_$1.M2f
   s_var var_0 cpt_$1.M2k
   s_var var_0 cpt_$1.M2p
#
   s_var var_0 cpt_$1.fwk
   s_var var_0 cpt_$1.wArchH
   s_var var_0 cpt_$1.wArchH_1
   s_var var_0 cpt_$1.wArchH_2
   s_var var_0 cpt_$1.wArchH_3
   s_var var_0 cpt_$1.wArchH_4
   s_var var_0 cpt_$1.wArchH_5
   s_var var_0 cpt_$1.wArchH_6
   s_var var_0 cpt_$1.wArchH_7
   s_var var_0 cpt_$1.wArchH_8
   s_var var_0 cpt_$1.wArchH_9
   s_var var_0 cpt_$1.wArchH_10
   s_var var_0 cpt_$1.wArchH_11
   s_var var_0 cpt_$1.wArchH_12
   s_var var_0 cpt_$1.wArchH_13
   s_var var_0 cpt_$1.wArchH_14
   s_var var_0 cpt_$1.wArchH_15
   s_var var_0 cpt_$1.wArchH_16
   s_var var_0 cpt_$1.wArchS
   s_var var_0 cpt_$1.wArchS_1
   s_var var_0 cpt_$1.wArchS_2
   s_var var_0 cpt_$1.wArchS_3
   s_var var_0 cpt_$1.wArchS_4
   s_var var_0 cpt_$1.wArchS_5
   s_var var_0 cpt_$1.wArchS_6
   s_var var_0 cpt_$1.wArchS_7
   s_var var_0 cpt_$1.wArchS_8
   s_var var_0 cpt_$1.wArchS_9
   s_var var_0 cpt_$1.wArchS_10
   s_var var_0 cpt_$1.wArchS_11
   s_var var_0 cpt_$1.wArchS_12
   s_var var_0 cpt_$1.wArchS_13
   s_var var_0 cpt_$1.wArchS_14
   s_var var_0 cpt_$1.wArchS_15
   s_var var_0 cpt_$1.wArchS_16
  ]
  in_substruct store_axle  [ 111  ]
  in_substruct store_axle  [ 112  ]
  in_substruct store_axle  [ 113  ]
  in_substruct store_wheel [ 111r ]
  in_substruct store_wheel [ 112r ]
  in_substruct store_wheel [ 113r ]
  in_substruct store_wheel [ 111l ]
  in_substruct store_wheel [ 112l ]
  in_substruct store_wheel [ 113l ]
#
# in_substruct store [ 111l ]
# in_substruct store [ 112l ]
# in_substruct store [ 113l ]
#
 eof
