2.9 -0.65 1.15 : rr0m,zr0m,rpam ! rhot[m] = min(sqrt((r-rr0m)**2+(z-zr0m)**2), rpam); tor flux phi = pi*b0*rhot**2 -2.04 : b0 ! Bphi[T] @ rr0m[m] 1 3 2.881899 : q0,qa,alq ! q = q0 + (qa-q0)*sqrt(psin)**alq 60 : nlim ! number of points in first wall (rlim,zlim) polygon 1.70484 -1.77790 1.70484 1.87160 1.72500 1.89510 1.72500 2.09510 1.76960 2.17060 1.76960 2.55560 1.77112 2.56755 1.99463 2.88089 2.00394 2.88321 2.36973 3.00329 2.37846 2.99654 2.84841 2.87265 2.84793 2.77148 3.10321 2.61390 3.11286 2.61105 3.42780 2.36445 3.43529 2.35776 3.71576 2.07256 3.72232 2.06496 3.96360 1.74593 3.96913 1.73754 4.16719 1.39002 4.17159 1.38099 4.32311 1.01080 4.32631 1.00128 4.42876 0.61463 4.43069 0.60477 4.48235 0.20812 4.48300 0.19810 4.48300 -0.20190 4.48235 -0.21192 4.43069 -0.60857 4.42876 -0.61843 4.32631 -1.00508 4.32311 -1.01460 4.17159 -1.38479 4.16719 -1.39382 3.96913 -1.74134 3.96360 -1.74973 3.72232 -2.06876 3.71820 -2.07792 3.44162 -2.36689 3.42605 -2.36323 3.11119 -2.60994 3.10556 -2.61897 2.75674 -2.81473 2.74725 -2.81806 2.35254 -2.88291 2.36752 -3.03568 2.16964 -3.00662 2.16013 -3.00358 1.98211 -2.91242 1.97408 -2.90648 1.83484 -2.76291 1.82915 -2.75471 1.74347 -2.57399 1.74072 -2.56439 1.71772 -2.36572 1.72500 -2.35630 1.72500 -1.81130