From e2c73fd156fc3f2316e884e3541a9ff44ccebda8 Mon Sep 17 00:00:00 2001 From: Paul Romano Date: Mon, 15 Oct 2012 09:07:54 -0400 Subject: [PATCH] Surface sense calculation for cones. --- src/geometry.F90 | 30 ++++++++++++++++++++++++++++++ 1 file changed, 30 insertions(+) diff --git a/src/geometry.F90 b/src/geometry.F90 index 72e074d7e4..e74e356a60 100644 --- a/src/geometry.F90 +++ b/src/geometry.F90 @@ -1035,6 +1035,36 @@ contains z = z - z0 func = x*x + y*y + z*z - r*r + case (SURF_CONE_X) + x0 = surf % coeffs(1) + y0 = surf % coeffs(2) + z0 = surf % coeffs(3) + r = surf % coeffs(4) + x = x - x0 + y = y - y0 + z = z - z0 + func = y*y + z*z - r*x*x + + case (SURF_CONE_Y) + x0 = surf % coeffs(1) + y0 = surf % coeffs(2) + z0 = surf % coeffs(3) + r = surf % coeffs(4) + x = x - x0 + y = y - y0 + z = z - z0 + func = x*x + z*z - r*y*y + + case (SURF_CONE_Z) + x0 = surf % coeffs(1) + y0 = surf % coeffs(2) + z0 = surf % coeffs(3) + r = surf % coeffs(4) + x = x - x0 + y = y - y0 + z = z - z0 + func = x*x + y*y - r*z*z + case (SURF_BOX_X) y0 = surf % coeffs(1) z0 = surf % coeffs(2)