toroidat_coord.c running 
r1=2.00, r2=0.50 

changing phi
theta=0.00, phi= 0.00, Xc=  2.000, Yc=  0.000, Zc=  0.500 
theta=0.00, phi= 0.10, Xc=  2.050, Yc=  0.000, Zc=  0.498 
theta=0.00, phi= 0.20, Xc=  2.099, Yc=  0.000, Zc=  0.490 
theta=0.00, phi= 1.47, Xc=  2.498, Yc=  0.000, Zc=  0.050 
theta=0.00, phi= 1.57, Xc=  2.500, Yc=  0.000, Zc=  0.000 
theta=0.00, phi= 1.67, Xc=  2.498, Yc=  0.000, Zc= -0.050 
inverse transform error = 3.9366 
err1=-0.995004, errt=0, errp=-2.94159 
theta=0.00, phi= 3.04, Xc=  2.050, Yc=  0.000, Zc= -0.498 
inverse transform error = 0.299833 
err1=-0.0998334, errt=0, errp=-0.2 
theta=0.00, phi= 3.14, Xc=  2.000, Yc=  0.000, Zc= -0.500 
theta=0.00, phi= 3.24, Xc=  1.950, Yc=  0.000, Zc= -0.498 
theta=0.00, phi= 4.61, Xc=  1.502, Yc=  0.000, Zc= -0.050 
theta=0.00, phi= 4.71, Xc=  1.500, Yc=  0.000, Zc= -0.000 
theta=0.00, phi= 4.81, Xc=  1.502, Yc=  0.000, Zc=  0.050 
inverse transform error = 4.3366 
err1=0.995004, errt=0, errp=3.34159 
theta=0.00, phi= 6.18, Xc=  1.950, Yc=  0.000, Zc=  0.498 
inverse transform error = 6.18302 
err1=0.0998334, errt=0, errp=6.08319 

changing theta
theta=0.00, phi= 0.00, Xc=  2.000, Yc=  0.000, Zc=  0.500 
theta=0.10, phi= 0.00, Xc=  1.990, Yc=  0.200, Zc=  0.500 
theta=0.20, phi= 0.00, Xc=  1.960, Yc=  0.397, Zc=  0.500 
theta=1.47, phi= 0.00, Xc=  0.200, Yc=  1.990, Zc=  0.500 
theta=1.57, phi= 0.00, Xc=  0.000, Yc=  2.000, Zc=  0.500 
theta=1.67, phi= 0.00, Xc= -0.200, Yc=  1.990, Zc=  0.500 
theta=3.04, phi= 0.00, Xc= -1.990, Yc=  0.200, Zc=  0.500 
theta=3.14, phi= 0.00, Xc= -2.000, Yc=  0.000, Zc=  0.500 
theta=3.24, phi= 0.00, Xc= -1.990, Yc= -0.200, Zc=  0.500 
theta=4.61, phi= 0.00, Xc= -0.200, Yc= -1.990, Zc=  0.500 
theta=4.71, phi= 0.00, Xc= -0.000, Yc= -2.000, Zc=  0.500 
theta=4.81, phi= 0.00, Xc=  0.200, Yc= -1.990, Zc=  0.500 
theta=6.18, phi= 0.00, Xc=  1.990, Yc= -0.200, Zc=  0.500 
toro2f.c finished