104       real*8, allocatable :: mscale(:)  local
124       allocate (mscale(n))
129          mscale(i) = 1.0d0
168             mscale(i12(j,i)) = m2scale
171             mscale(i13(j,i)) = m3scale
174             mscale(i14(j,i)) = m4scale
177             mscale(i15(j,i)) = m5scale
232                   rr1 = f * mscale(k) / r
301             mscale(i12(j,i)) = 1.0d0
304             mscale(i13(j,i)) = 1.0d0
307             mscale(i14(j,i)) = 1.0d0
310             mscale(i15(j,i)) = 1.0d0
349                mscale(i12(j,i)) = m2scale
352                mscale(i13(j,i)) = m3scale
355                mscale(i14(j,i)) = m4scale
358                mscale(i15(j,i)) = m5scale
380      &                  mscale(k) = 1.0d0
415                         rr1 = f * mscale(k) / r
486                mscale(i12(j,i)) = 1.0d0
489                mscale(i13(j,i)) = 1.0d0
492                mscale(i14(j,i)) = 1.0d0
495                mscale(i15(j,i)) = 1.0d0
502       deallocate (mscale)
567       real*8, allocatable :: mscale(:)  local
587       allocate (mscale(n))
592          mscale(i) = 1.0d0
641             mscale(i12(j,i)) = m2scale
644             mscale(i13(j,i)) = m3scale
647             mscale(i14(j,i)) = m4scale
650             mscale(i15(j,i)) = m5scale
706                   rr1 = f * mscale(k) / r
775             mscale(i12(j,i)) = 1.0d0
778             mscale(i13(j,i)) = 1.0d0
781             mscale(i14(j,i)) = 1.0d0
784             mscale(i15(j,i)) = 1.0d0
795       deallocate (mscale)
973       real*8, allocatable :: mscale(:)  local
979       allocate (mscale(n))
984          mscale(i) = 1.0d0
1019             mscale(i12(j,i)) = m2scale
1022             mscale(i13(j,i)) = m3scale
1025             mscale(i14(j,i)) = m4scale
1028             mscale(i15(j,i)) = m5scale
1105                   scalek = mscale(k)
1106                   scalek = mscale(k)
1133                   scalek = 1.0d0 - mscale(k)
1152             mscale(i12(j,i)) = 1.0d0
1155             mscale(i13(j,i)) = 1.0d0
1158             mscale(i14(j,i)) = 1.0d0
1161             mscale(i15(j,i)) = 1.0d0
1196                mscale(i12(j,i)) = m2scale
1199                mscale(i13(j,i)) = m3scale
1202                mscale(i14(j,i)) = m4scale
1205                mscale(i15(j,i)) = m5scale
1219      &               mscale(k) = 1.0d0
1285                         scalek = mscale(k)
1314                         scalek = 1.0d0 - mscale(k)
1335                mscale(i12(j,i)) = 1.0d0
1338                mscale(i13(j,i)) = 1.0d0
1341                mscale(i14(j,i)) = 1.0d0
1344                mscale(i15(j,i)) = 1.0d0
1351       deallocate (mscale)
1528       real*8, allocatable :: mscale(:)  local
1534       allocate (mscale(n))
1539          mscale(i) = 1.0d0
1583             mscale(i12(j,i)) = m2scale
1586             mscale(i13(j,i)) = m3scale
1589             mscale(i14(j,i)) = m4scale
1592             mscale(i15(j,i)) = m5scale
1670                   scalek = mscale(k)
1697                   scalek = 1.0d0 - mscale(k)
1716             mscale(i12(j,i)) = 1.0d0
1719             mscale(i13(j,i)) = 1.0d0
1722             mscale(i14(j,i)) = 1.0d0
1725             mscale(i15(j,i)) = 1.0d0
1736       deallocate (mscale)