18 real temp(9), pi, tpi, zmu0
22 parameter(pi = 3.1415927)
23 parameter(tpi = 2.0*pi)
24 parameter(zmu0 = 4.0e-7*pi)
74 parms(7) = 1.0 / temp(2)
75 parms(8) = 1.0 / temp(3)
76 parms(9) = temp(3) - temp(2)
77 parms(10) = alog(temp(3)/temp(2))
78 parms(11) = temp(2) * tpi * temp(6)
79 parms(12) = tpi * temp(5)
80 parms(13) = tpi * temp(6)
81 parms(14) = tpi * temp(7)
82 parms(15) = 1.0 / parms(12)
83 parms(16) = 1.0 / parms(13)
84 parms(17) = 1.0 / parms(14)
86 zmu = (temp(8) - bh1) / (zmu0*(temp(6)-temp(5)))
87 parms(18) = 1.0 + temp(1)*(zmu-1.0)
88 parms(19) = parms(18) * parms(10)
89 parms(20) = parms(18) - 1.0
90 parms(21) = 1.0 - parms(18)
91 parms(22) = temp(4) * zmu0 / tpi
subroutine magparm(parms, nparms)