; total = 172.668*p31
; d31=p31*1.599
;R
  (p31*2.000 ph0 d31):C
  (p31*1.556 ph0 3u):C
  (p31*3.556 ph2 d31):C
  (p31*1.000 ph2 3u):C
  (p31*3.000 ph0 d31):C
  (p31*0.222 ph0 3u):C
  (p31*2.222 ph2 d31):C
  (p31*0.944 ph2 3u):C
  (p31*0.333 ph0 3u):C
  (p31*1.389 ph2 d31):C
  (p31*1.333 ph2 3u):C
  (p31*3.333 ph0 d31):C
  (p31*0.833 ph0 3u):C
  (p31*2.833 ph2 d31):C
  (p31*0.111 ph2 3u):C
  (p31*2.111 ph0 d31):C
  (p31*2.000 ph0 d31):C
;R-bar
  (p31*2.000 ph2 d31):C
  (p31*1.556 ph2 3u):C
  (p31*3.556 ph0 d31):C
  (p31*1.000 ph0 3u):C
  (p31*3.000 ph2 d31):C
  (p31*0.222 ph2 3u):C
  (p31*2.222 ph0 d31):C
  (p31*0.944 ph0 3u):C
  (p31*0.333 ph2 3u):C
  (p31*1.389 ph0 d31):C
  (p31*1.333 ph0 3u):C
  (p31*3.333 ph2 d31):C
  (p31*0.833 ph2 3u):C
  (p31*2.833 ph0 d31):C
  (p31*0.111 ph0 3u):C
  (p31*2.111 ph2 d31):C
  (p31*2.000 ph2 d31):C
;R-bar
  (p31*2.000 ph2 d31):C
  (p31*1.556 ph2 3u):C
  (p31*3.556 ph0 d31):C
  (p31*1.000 ph0 3u):C
  (p31*3.000 ph2 d31):C
  (p31*0.222 ph2 3u):C
  (p31*2.222 ph0 d31):C
  (p31*0.944 ph0 3u):C
  (p31*0.333 ph2 3u):C
  (p31*1.389 ph0 d31):C
  (p31*1.333 ph0 3u):C
  (p31*3.333 ph2 d31):C
  (p31*0.833 ph2 3u):C
  (p31*2.833 ph0 d31):C
  (p31*0.111 ph0 3u):C
  (p31*2.111 ph2 d31):C
  (p31*2.000 ph2 d31):C
;R
  (p31*2.000 ph0 d31):C
  (p31*1.556 ph0 3u):C
  (p31*3.556 ph2 d31):C
  (p31*1.000 ph2 3u):C
  (p31*3.000 ph0 d31):C
  (p31*0.222 ph0 3u):C
  (p31*2.222 ph2 d31):C
  (p31*0.944 ph2 3u):C
  (p31*0.333 ph0 3u):C
  (p31*1.389 ph2 d31):C
  (p31*1.333 ph2 3u):C
  (p31*3.333 ph0 d31):C
  (p31*0.833 ph0 3u):C
  (p31*2.833 ph2 d31):C
  (p31*0.111 ph2 3u):C
  (p31*2.111 ph0 d31):C
  (p31*2.000 ph0 d31):C