Home    Search    Reference Manuals    Return   


#
#
#       Substructure file t_irr_p1.ins   rev.971113
#
#       ----------------------------------------------------------
#       Stiff axle.
#
#       Track irregularities will be interpolated by cubic splines
#       read from fields generated by command "func intpl_r",
#	"func fifo_mem" or "fifo_mem2".
#
#       Individual wheel/rail geometry functions for each wheel, will be
#	read if the geometry functions contains the wheel number.
#
#	All operations are made in one big function, in order to speed up
#	the calculations.
#
#       ----------------------------------------------------------
#
#       This substructure do the following:
#         1) Set constraints to the axle in pitch motion, and in all
#            directions for the track except in lateral direction.
#         2) Interpolate the track irregularities by cubic splines
#         3) Calculate cp_$2r.eta and cp_$2l.eta as the lateral
#	     displacement between wheel and rail
#         4) Interpolate the wheel-rail geometry functions linearly
#         5) Calculate the deformations and normal forces in the
#	     contact points.
#         6) Calculate creep and creep forces in the contact points.
#
#       ----------------------------------------------------------
#
#
#       Variables which must be defined before this substructure file can
#       be used:
#
#      Vo       = Nominal speed of the vehicle in m/s
#
#      Ro$2h    = Nominal rolling radius, axle number $2, right side.
#      ro$2r    = If Ro$2h not can be found
#      Ro$2v    = Nominal rolling radius, axle number $2, left side.
#      ro$2l    = If Ro$2v not can be found
#      Ro$1     = If not the variables above can be found
#      ro$1     = If Ro$1 not can be found
#      Ro       = If ro$1 not can be found
#      ro       = If Ro not can be found
#
#      Boh      = Lateral distance between track center line and nominal 
#                 rolling circle of right wheel
#      Bo       = If Boh not can be found
#      bo       = If Bo  not can be found
#      0.75     = If bo  not can be found
#
#      Bov      = Lateral distance between track center line and nominal
#                 rolling circle of left wheel
#      -Bo      = If Bov not can be found
#      -bo      = If Bo  not can be found
#      -.75     = If bo  not can be found
#
#      kph$2.my  = Friction coefficient between wheel and rail, tread right side.
#      cpt_$2r.mu= If kph$2.my   not can be found
#      my        = If cpt_$2r.mu not can be found
#      mu        = If my         not can be found
#
#      kpv$2.my  = Friction coefficient between wheel and rail, tread left side.
#      cpt_$2l.mu= If kpv$2.my   not can be found
#      my        = If cpt_$2l.mu not can be found
#      mu        = If my         not can be found
#
#      kfh$2.my  = Friction coefficient between wheel and rail, flange left side.
#      cpf_$2r.mu= If kfh$2.my   not can be found
#      myf       = If cpf_$2r.mu not can be found
#      muf       = If myf        not can be found
#      my        = If muf        not can be found
#      mu        = If my         not can be found
#
#      kfv$2.my  = Friction coefficient between wheel and rail, flange left side.
#      cpf_$2l.mu= If kfv$2.my   not can be found
#      myf       = If cpf_$2l.mu not can be found
#      muf       = If myf        not can be found
#      my        = If muf        not can be found
#      mu        = If my         not can be found
#
#
#      mulfact_nux_tread  = Longitudinal creepage reduction factor on tread
#      mulfact_nux        = If mulfact_nux_tread   not can be found
#      mulfact_nuy_tread  = Lateral      creepage reduction factor on tread
#      mulfact_nux        = If mulfact_nuy_tread   not can be found
#      mulfact_spin_tread = Spin         creepage reduction factor on tread
#      mulfact_nux        = If mulfact_spin_tread  not can be found
#      mulfact_nux_flange = Longitudinal creepage reduction factor on flange
#      mulfact_nux        = If mulfact_nux_flange  not can be found
#      mulfact_nuy_flange = Lateral      creepage reduction factor on flange
#      mulfact_nux        = If mulfact_nuy_flange  not can be found
#      mulfact_spin_flange= Spin         creepage reduction factor on flange
#      mulfact_nux        = If mulfact_spin_flange not can be found
#      If one or more of the above multiplication factors are missing,
#      they will be replaced with the value 1.
#
#
#	Surface stiffness between rail and wheel
#	========================================
#      knwr$2h_r= Stiffness perpendicular to the contact surface
#		  on tread right wheel, defined by "func intpl_r".
#                Ex:  func intpl_r knwr$2h_r -1.57 Horizontal stiffness
#					      0.   Vertical stiffness
#                                             1.57 Horizontal stiffness
#      knwr$2r_r= If knwr$2h_r  not can be found
#      knwr$1_r = If knwr$2r_r  not can be found
#      knwr_r   = If knwr$1_r   not can be found
#
#      knwr$2h_r.F0= Pre-stress force in spring knwr$2h_r.
#      knwr$2r_r.F0= If knwr$2h_r.F0 not can be found
#      knwr$1_r.F0 = If knwr$2r_r.F0 not can be found
#      knwr_r.F0   = If knwr$1_r.F0  not can be found
#
#
#      knwr$2v_r= Stiffness perpendicular to the contact surface
#                 on tread left wheel, defined by "func intpl_r".
#      knwr$2l_r= If knwr$2v_r  not can be found
#      knwr$1_r = If knwr$2l_r  not can be found
#      knwr_r   = If knwr$1_r   not can be found
#
#      knwr$2v_r.F0= Pre-stress force in spring knwr$2v_r.
#      knwr$2l_r.F0= If knwr$2v_r.F0 not can be found
#      knwr$1_r.F0 = If knwr$2l_r.F0 not can be found
#      knwr_r.F0   = If knwr$1_r.F0  not can be found
#
#      pknwr =	The value to which the rail-wheel deformation shall be
#		raised to when calculating the rail-wheel contact force.
#		In Hertz theory this value should be set to 1.5.
#		If not pknwr is defined pknwr will be set to 1. i.e.
#		a linear stiffness between rail-wheels are assumed.
#		Normally the rail-wheel stiffness is so high compared to
#		other stiffnesses in the model that it only work as a
#		method of calculating how much vertical load is resting
#		on the tread and how much vertical load is resting
#		on the flange.
#
#
#       Surface stiffness between rail and wheel
#       ======================================== 
#      kyrt$2h= Lateral stiffness between right rail and track
#      kyrt$2r= If kyrt$2h  not can be found
#      kyrt$1 = If kyrt$2r  not can be found
#      kyrt   = If kyrt$1   not can be found
#
#      kzrt$2h= Vertical stiffness between right rail and track
#      kzrt$2r= If kzrt$2h  not can be found
#      kzrt$1 = If kzrt$2r  not can be found
#      kzrt   = If kzrt$1   not can be found
#
#      kzrt$2h.F0= Pre-stress force in spring kzrt$2h
#      kzrt$2r.F0= If kzrt$2h.F0  not can be found
#      kzrt$1.F0 = If kzrt$2r.F0  not can be found
#      kzrt.F0   = If kzrt$1.F0   not can be found
#      0.        = If kzrt.F0     not can be found
#
#      cyrt$2h= Lateral damping between right rail and track
#      cyrt$2r= If cyrt$2h  not can be found
#      cyrt$1 = If cyrt$2r  not can be found
#      cyrt   = If cyrt$1   not can be found
#
#      czrt$2h= Vertical damping between right rail and track
#      czrt$2r= If czrt$2h  not can be found 
#      czrt$1 = If czrt$2r  not can be found
#      czrt   = If czrt$1   not can be found
#
#      kyrt$2v= Lateral stiffness between left rail and track
#      kyrt$2l= If kyrt$2v  not can be found
#      kyrt$1 = If kyrt$2l  not can be found
#      kyrt   = If kyrt$1   not can be found
#
#      kzrt$2v= Vertical stiffness between left rail and track
#      kzrt$2l= If kzrt$2v  not can be found
#      kzrt$1 = If kzrt$2l  not can be found
#      kzrt   = If kzrt$1   not can be found
#
#      kzrt$2v.F0= Pre-stress force in spring kzrt$2v
#      kzrt$2l.F0= If kzrt$2v.F0  not can be found
#      kzrt$1.F0 = If kzrt$2l.F0  not can be found
#      kzrt.F0   = If kzrt$1.F0   not can be found
#      0.        = If kzrt.F0     not can be found
#
#      cyrt$2v= Lateral damping between left rail and track
#      cyrt$2l= If cyrt$2v  not can be found
#      cyrt$1 = If cyrt$2l  not can be found
#      cyrt   = If cyrt$1   not can be found
#
#      czrt$2v= Vertical damping between left rail and track
#      czrt$2l= If czrt$2v  not can be found 
#      czrt$1 = If czrt$2l  not can be found
#      czrt   = If czrt$1   not can be found
#
#      sv_trac        = Average gauge of the track defined in the data field spv_trac
#      gauge_average  = If sv_trac       not can be found
#      sv_v$1         = Deviation in gauge for vehicle number $1
#			A positive value in sv_v$1 entails a wide gauge
#			A negative value in sv_v$1 entails a narrow gauge
#      gauge_dev_v$1  = If sv_v$1        not can be found
#      sv_v           = If gauge_dev_v$1 not can be found
#      gauge_dev      = If sv_v          not can be found
#
#      YMtrac   = Multiplying factor for lateral track irregularities
#      ZMtrac   = Multiplying factor for vertical track irregularities
#      FItrac   = Multiplying factor for cant irregularities
#      CMtrac   = If FItrac not can be found
#      SVtrac   = Multiplying factor for gauge irregularities
#      GMtrac   = If SVtrac not can be found
#
#      lsa_$2.pn= Position along the track of the linear local
#		  coordinate system lsa_$2
#      axl_$2.? = Variables containing the motions of the axle (wheelset)
#      trc_$2.? = Variables containing the motions of the track
#
#     Track irregularities will be interpolated from data fields with the
#     following names:
#
#      lat_trac  = Lateral irregularities (positive direction = right)
#      vert_trac = Vertical irregularities (positive direction = down)
#      fi_trac   = Cant irregularities (positive direction = positive
#		   rotation round the x-axle)
#      spv_trac  = Gauge irregularities (positive direction = wide gauge)
#      gauge_trac= If spv_trac not can be found
#
#
#     Wheel-rail geometry functions describing the tread of the right
#     wheel and rail shall be defined by command `func intpl_r` and must
#     have one of the following names:
#
#      drkp_r$2h    = Change in rolling radius, tread, right wheel
#      drkp_r$2r    = If drkp_r$2h    not can be found
#      cpt_$2r.drfn = If drkp_r$2r    not can be found
#      drkp_r$1     = If cpt_$2r.drfn not can be found
#      cpt_$1.drfn  = If drkp_r$1     not can be found
#      drkp_r       = If cpt_$1.drfn  not can be found
#      cpt_.drfn    = If drkp_r       not can be found
#
#      gamkph_r$2h  = Angle of contact surface, tread, right wheel
#      gamkph_r$2r  = If gamkph_r$2h  not can be found
#      cpt_$2r.gamfn= If gamkph_r$2r   not can be found
#      gamkph_r$1   = If cpt_$2r.gamfn not can be found
#      cpt_$1.gamfn = If gamkph_r$1    not can be found
#      gamkph_r     = If cpt_$1.gamfn  not can be found
#      cpt_.gamfn   = If gamkph_r      not can be found
#
#      zkp_r$2h     = Wheel lift, tread, right wheel
#      zkp_r$2r     = If zkp_r$2h     not can be found
#      cpt_$2r.zfn  = If zkp_r$2r      not can be found
#      zkp_r$1      = If cpt_$2r.zfn   not can be found
#      cpt_$1.zfn   = If zkp_r$1       not can be found
#      zkp_r        = If cpt_$1.zfn    not can be found
#      cpt_.zfn     = If zkp_r         not can be found
#
#      rohr_r$2h    = Wheel-rail curvature difference, tread, right wheel
#      rohr_r$2r    = If rohr_r$2h    not can be found
#      cpt_$2r.rofn = If rohr_r$2r     not can be found
#      rohr_r$1     = If cpt_$2r.rofn  not can be found
#      cpt_$1.rofn  = If rohr_r$1      not can be found
#      rohr_r       = If cpt_$1.rofn   not can be found
#      cpt_.rofn    = If rohr_r        not can be found
#
#      poskph_r$2h   = Position of contact surface of wheel, tread, right wheel
#      poskph_r$2r   = If poskph_r$2h    not can be found
#      cpt_$2r.poswfn= If poskph_r$2r    not can be found
#      poskph_r$1    = If cpt_$2r.poswfn not can be found
#      cpt_$1.poswfn = If poskph_r$1     not can be found
#      poskph_r      = If cpt_$1.poswfn  not can be found
#      cpt_.poswfn   = If poskph_r       not can be found
#      0.            = If cpt_.poswfn    not can be found
#
#
#     Wheel-rail geometry functions describing the tread of the left
#     wheel and rail shall be defined by command `func intpl_r` and must
#     have one of the following names:
#
#      drkp_r$2v    = Change in rolling radius, tread, left wheel
#      drkp_r$2l    = If drkp_r$2v    not can be found
#      cpt_$2l.drfn = If drkp_r$2l    not can be found
#      drkp_r$1     = If cpt_$2l.drfn not can be found
#      cpt_$1.drfn  = If drkp_r$1     not can be found
#      drkp_r       = If cpt_$1.drfn  not can be found
#      cpt_.drfn    = If drkp_r       not can be found
#
#      gamkpv_r$2v  = Angle of contact surface, tread, left wheel
#      gamkpv_r$2l  = If gamkph_r$2v   not can be found
#      cpt_$2l.gamfn= If gamkph_r$2l   not can be found
#      gamkpv_r$1   = If cpt_$2l.gamfn not can be found
#      cpt_$1.gamfn = If gamkph_r$1    not can be found
#      gamkpv_r     = If cpt_$1.gamfn  not can be found
#      cpt_.gamfn   = If gamkph_r      not can be found
#
#      zkp_r$2v     = Wheel lift, tread, left wheel
#      zkp_r$2l     = If zkp_r$2v      not can be found
#      cpt_$2l.zfn  = If zkp_r$2l      not can be found
#      zkp_r$1      = If cpt_$2l.zfn   not can be found
#      cpt_$1.zfn   = If zkp_r$1       not can be found
#      zkp_r        = If cpt_$1.zfn    not can be found
#      cpt_.zfn     = If zkp_r         not can be found
#
#      rohr_r$2v    = Wheel-rail curvature difference, tread, left wheel
#      rohr_r$2l    = If rohr_r$2v     not can be found
#      cpt_$2l.rofn = If rohr_r$2l     not can be found
#      rohr_r$1     = If cpt_$2l.rofn  not can be found
#      cpt_$1.rofn  = If rohr_r$1      not can be found
#      rohr_r       = If cpt_$1.rofn   not can be found
#      cpt_.rofn    = If rohr_r        not can be found
#
#      poskph_r$2v   = Position of contact surface of wheel, tread, left wheel
#      poskph_r$2l   = If poskph_r$2v    not can be found
#      cpt_$2l.poswfn= If poskph_r$2l    not can be found
#      poskph_r$1    = If cpt_$2l.poswfn not can be found
#      cpt_$1.poswfn = If poskph_r$1     not can be found
#      poskph_r      = If cpt_$1.poswfn  not can be found
#      cpt_.poswfn   = If poskph_r       not can be found
#      0.            = If cpt_.poswfn    not can be found
#
#     Wheel-rail geometry functions describing the flange of the right
#     wheel and rail shall be defined by command `func intpl_r` and must
#     have one of the following names:
#
#      drkp_f$2h    = Change in rolling radius, flange right wheel
#      drkp_f$2r    = If drkp_f$2h    not can be found
#      cpf_$2r.drfn = If drkp_f$2r    not can be found
#      drkp_f$1     = If cpf_$2r.drfn not can be found
#      cpf_$1.drfn  = If drkp_f$1     not can be found
#      drkp_f       = If cpf_$1.drfn  not can be found
#      cpf_.drfn    = If drkp_f       not can be found
#
#      gamkph_f$2h  = Angle of contact surface, flange right wheel
#      gamkph_f$2r  = If gamkph_f$2h  not can be found
#      cpf_$2r.gamfn= If gamkph_f$2r   not can be found
#      gamkph_f$1   = If cpf_$2r.gamfn not can be found
#      cpf_$1.gamfn = If gamkph_f$1    not can be found
#      gamkph_f     = If cpf_$1.gamfn  not can be found
#      cpf_.gamfn   = If gamkph_f      not can be found
#
#      zkp_f$2h     = Wheel lift, flange right wheel
#      zkp_f$2r     = If zkp_f$2h     not can be found
#      cpf_$2r.zfn  = If zkp_f$2r      not can be found
#      zkp_f$1      = If cpf_$2r.zfn   not can be found
#      cpf_$1.zfn   = If zkp_f$1       not can be found
#      zkp_f        = If cpf_$1.zfn    not can be found
#      cpf_.zfn     = If zkp_f         not can be found
#
#      rohr_f$2h    = Wheel-rail curvature difference, flange right wheel
#      rohr_f$2r    = If rohr_f$2h    not can be found
#      cpf_$2r.rofn = If rohr_f$2r     not can be found
#      rohr_f$1     = If cpf_$2r.rofn  not can be found
#      cpf_$1.rofn  = If rohr_f$1      not can be found
#      rohr_f       = If cpf_$1.rofn   not can be found
#      cpf_.rofn    = If rohr_f        not can be found
#
#      poskph_f$2h   = Position of contact surface of wheel, flange right wheel
#      poskph_f$2r   = If poskph_f$2h    not can be found
#      cpf_$2r.poswfn= If poskph_f$2r    not can be found
#      poskph_f$1    = If cpf_$2r.poswfn not can be found
#      cpf_$1.poswfn = If poskph_f$1     not can be found
#      poskph_f      = If cpf_$1.poswfn  not can be found
#      cpf_.poswfn   = If poskph_f       not can be found
#      0.            = If cpf_.poswfn    not can be found
#
#     Wheel-rail geometry functions describing the flange of the left
#     wheel and rail shall be defined by command `func intpl_r` and must
#     have one of the following names:
#
#      drkp_f$2h    = Change in rolling radius, flange left wheel
#      drkp_f$2r    = If drkp_f$2h    not can be found
#      cpf_$2r.drfn = If drkp_f$2r    not can be found
#      drkp_f$1     = If cpf_$2r.drfn not can be found
#      cpf_$1.drfn  = If drkp_f$1     not can be found
#      drkp_f       = If cpf_$1.drfn  not can be found
#      cpf_.drfn    = If drkp_f       not can be found
#
#      gamkph_f$2h  = Angle of contact surface, flange left wheel
#      gamkph_f$2r  = If gamkph_f$2h   not can be found
#      cpf_$2r.gamfn= If gamkph_f$2r   not can be found
#      gamkph_f$1   = If cpf_$2r.gamfn not can be found
#      cpf_$1.gamfn = If gamkph_f$1    not can be found
#      gamkph_f     = If cpf_$1.gamfn  not can be found
#      cpf_.gamfn   = If gamkph_f      not can be found
#
#      zkp_f$2h     = Wheel lift, flange left wheel
#      zkp_f$2r     = If zkp_f$2h      not can be found
#      cpf_$2r.zfn  = If zkp_f$2r      not can be found
#      zkp_f$1      = If cpf_$2r.zfn   not can be found
#      cpf_$1.zfn   = If zkp_f$1       not can be found
#      zkp_f        = If cpf_$1.zfn    not can be found
#      cpf_.zfn     = If zkp_f         not can be found
#
#      rohr_f$2h    = Wheel-rail curvature difference, flange left wheel
#      rohr_f$2r    = If rohr_f$2h     not can be found
#      cpf_$2r.rofn = If rohr_f$2r     not can be found
#      rohr_f$1     = If cpf_$2r.rofn  not can be found
#      cpf_$1.rofn  = If rohr_f$1      not can be found
#      rohr_f       = If cpf_$1.rofn   not can be found
#      cpf_.rofn    = If rohr_f        not can be found
#
#      poskph_f$2v   = Position of contact surface of wheel, flange left wheel
#      poskph_f$2l   = If poskph_f$2v    not can be found
#      cpf_$2l.poswfn= If poskph_f$2l    not can be found
#      poskph_f$1    = If cpf_$2l.poswfn not can be found
#      cpf_$1.poswfn = If poskph_f$1     not can be found
#      poskph_f      = If cpf_$1.poswfn  not can be found
#      cpf_.poswfn   = If poskph_f       not can be found
#      0.            = If cpf_.poswfn    not can be found
#
#
#       Output:
#       =======
#       This substructure file will generate the following variables:
#
#       tral$2.y   = Lateral track irregularities, track center line
#       tral$2r.y  = Lateral track irregularities, right rail
#       tral$2l.y  = Lateral track irregularities, left rail
#       tral$2r.z  = Vertical track irregularities, right rail
#       tral$2l.z  = Vertical track irregularities, left rail
#       tral$2r.vy = First order derivative of lateral track irregularities, right rail
#       tral$2l.vy = First order derivative of lateral track irregularities, left rail
#       tral$2r.vz = First order derivative of vertical track irregularities, right rail
#       tral$2l.vz = First order derivative of vertical track irregularities, left rail
#       tral$2.f   = Track irregularities in roll direction
#       tral$2r.k  = Pitch irregularities right rail
#       tral$2l.k  = Pitch irregularities left rail
#       tral$2r.p  = Yaw irregularities right rail
#       tral$2l.p  = Yaw irregularities left rail
#
#       cp_$2r.eta = Lateral shift between wheel and rail right wheel
#       cpt_$2r.ksi= Longitudinal position of the contact point, tread right wheel
#       cpt_$2r.dr = Change in wheel radius due to wheel-rail geometry, tread right wheel
#       cpt_$2r.gam= Angle of contact point, tread right wheel
#       cpt_$2r.z  = Vertical motion of wheel due to wheel-rail geometry, tread right wheel
#       cpt_$2r.irx= Longitudinal wheel-rail curvature difference, tread right wheel
#       cpt_$2r.iry= Lateral wheel-rail curvature difference, tread right wheel
#       cpt_$2r.bo = Lateral distance track center line to contact point, tread right wheel
#
#       cpt_$2r.nux = Longitudinal creepage, the division has been made with regard to the 
#			longitudinal velocity of the car-body, tread right wheel
#       cpt_$2r.nuy = Lateral creepage, the division has been made with regard to the
#                       longitudinal velocity of the car-body, tread right wheel
#       cpt_$2r.spin= Spin creepage, the division has been made with regard to the
#                       longitudinal velocity of the car-body, tread right wheel
#       cpt_$2r.nuxm= Modified long. creepage, due to greasy contact surface, tread right wheel
#       cpt_$2r.nuym= Modified lat. creepage, due to greasy contact surface, tread right wheel
#       cpt_$2r.spim= Modified spin creepage, due to greasy contact surface, tread right wheel
#
#       cpt_$2r.a/b = The a/b-ratio of the contact ellipse, tread right wheel
#       cpt_$2r.c   = Geom. average radius of the contact ellipse c=sqrt(a*b), tread right wheel
#
#       cpt_$2r.Fn  = Contact force perpendicular to the contact surface, tread right wheel
#       cpt_$2r.Fny = Lateral creep force in plane with the contact surface, tread right wheel
#       cpt_$2r.Fx  = Longitudinal force between wheel and track, tread right wheel
#       cpt_$2r.Fy  = Lateral force between wheel and track, tread right wheel
#       cpt_$2r.Fz  = Vertical force between wheel and track, tread right wheel
#
#	ral_$2r.y  = Lateral  position of the massless rail-head, tread right wheel
#	ral_$2r.vy = Lateral  velocity of the massless rail-head, tread right wheel
#	ral_$2r.z  = Vertical position of the massless rail-head, tread right wheel
#	ral_$2r.vz = Vertical velocity of the massless rail-head, tread right wheel
#
#	The same naming convention has also been used on the wheels on the left hand side
#	of the vehicle. The names of the variables on the left side can be created by
#	changing the name ral_$2r.* into ral_$2l.*.
#
#	If the wheel-rail geometry functions contain a separate description of the flange
#	(two-point contact), an equal amount of variables as for the tread will be generated.
#	The names of the variables will be the same as for the tread, just the introduction
#	cpt_* will be changed into cpf_*.
#       ---------------------------------------------------------------------------------------
#
#
#
#                       $1=car #, $2=axle #
   substruct trac_irr [
note
note
note       Start of substructure file t_irr_p1.ins     Rev.980212
note       -------------------------------------------------------
#
#
#       Set constraints to axle and track
#       ---------------------------------
  constr fix_free_1  axl_$2.k  0.
#
  constr fix_rigid_1 trc_$2  x 0.  constr fix_rigid_1 trc_$2  z 0.
  constr fix_rigid_1 trc_$2  f 0.  constr fix_rigid_1 trc_$2  k 0.
  constr fix_rigid_1 trc_$2  p 0.
#
  func t_irr_p1  $2  $1
#
#
       cpt_$2r     trc_$2  cpt_$2r.ksi cpt_$2r.bo  0.
                   axl_$2  cpt_$2r.ksi cpt_$2r.bo  0.
#      cpt_$2r     trc_$2  cpt_$2r.ksi  { Boh | Bo | bo }  0.
#                  axl_$2  cpt_$2r.ksi  { Boh | Bo | bo }  0.
       0.          fzwrnr$2 cpt_$2r.gam cpt_$2r.nuxm cpt_$2r.nuym cpt_$2r.spim
#      1/Rh+1/Rr   1/Ryh+1/Ryr    my
       cpt_$2r.irx cpt_$2r.iry { kph$2.my | cpt_$2r.mu | myh | mutr | my | mu }
#      E        nu  jvkt.rikt
       2.05e11 .25    z
#
       cpt_$2l     trc_$2  cpt_$2l.ksi cpt_$2l.bo 0.
                   axl_$2  cpt_$2l.ksi cpt_$2l.bo 0.
#      cpt_$2l     trc_$2  cpt_$2l.ksi { Bov | -Bo | -bo } 0.
#                  axl_$2  cpt_$2l.ksi { Bov | -Bo | -bo } 0.
       0.          fzwrnl$2 cpt_$2l.gam cpt_$2l.nuxm cpt_$2l.nuym cpt_$2l.spim
#      1/Rh+1/Rr   1/Ryh+1/Ryr    my
       cpt_$2l.irx cpt_$2l.iry { kpv$2.my | cpt_$2l.mu | myv | mutl | my | mu }
#      E        nu  jvkt.rikt
       2.05e11 .25    z
#
#
       cpf_$2r     trc_$2  cpf_$2r.ksi cpf_$2r.bo 0.
                   axl_$2  cpf_$2r.ksi cpf_$2r.bo 0.
#      cpf_$2r     trc_$2  cpf_$2r.ksi { Boh | Bo | bo } 0.
#                  axl_$2  cpf_$2r.ksi { Boh | Bo | bo } 0.
       0.          fzfrnr$2 cpf_$2r.gam cpf_$2r.nuxm cpf_$2r.nuym cpf_$2r.spim
#      1/Rh+1/Rr   1/Ryh+1/Ryr
       cpf_$2r.irx cpf_$2r.iry
       { kfh$2.my | cpf_$2r.mu | myfh | mufr | myf | muf | my | mu }
#      E        nu  jvkt.rikt
       2.05e11 .25    z
#
       cpf_$2l     trc_$2  cpf_$2l.ksi cpf_$2l.bo 0.
                   axl_$2  cpf_$2l.ksi cpf_$2l.bo 0.
#      cpf_$2l     trc_$2  cpf_$2l.ksi { Bov | -Bo | -bo } 0.
#                  axl_$2  cpf_$2l.ksi { Bov | -Bo | -bo } 0.
       0.          fzfrnl$2 cpf_$2l.gam cpf_$2l.nuxm cpf_$2l.nuym cpf_$2l.spim
#      1/Rh+1/Rr   1/Ryh+1/Ryr
       cpf_$2l.irx cpf_$2l.iry
       { kfv$2.my | cpf_$2l.mu | myfv | mufl | myf | muf | my | mu }
#      E        nu  jvkt.rikt
       2.05e11 .25    z
#
#  func print06_all cpt_$2r.nux
#  func print06_all cpt_$2r.nuy
#  func print06_all cpt_$2r.spi
#  func print06_all cpt_$2r.Fnx
#  func print06_all cpt_$2l.Fn
#  func print06_all cpt_$2l.Fnx
#  func print06_all cpt_$2l.Fny
#  func print06_all cpt_$2r.Fn
#  func print06_all cpt_$2r.Fny
#
note
note       -------------------------------------------
note       End   of substructure file t_irr_p1.ins
note       -------------------------------------------
  ]