%
% file: v3dex_HydraA.sl
%
% Create a 3D model of the supercavities in Hydra A
% based on Wise et al 2006 (astro-ph/0612100 )
%

% Set "radius" of the 3D cube space
% units are kpc
v3d_setup(350.0, 37);   
v3d_user_unit = "kpc";

% define the locations of the cavities
% in R_sky, PA_sky, and z coordinates:

% kpc from radio core (Table 1 R)
variable rskys = [24.9, 25.6, 100.8, 59.3, 225.6, 104.3];
% angle CCW from North (guessing)
variable pas = [30.0, 215.0, 35.0, 220.0, 15.0, 185.0];
%
% convert to X (plus toward W), Y (plus toward N)
variable DTOR = PI/180.0;
variable xskys = -1.0 * rskys * sin(DTOR*pas);
variable yskys = rskys * cos(DTOR*pas);
%
% the cavity z's in kpc (Table 2 z) 
% not sure what the sign of z is
% here, +z is toward the observer...
% maybe the table/analysis only gives |z|
% put E (5th one) at a negative z to be more opposite F
variable zskys = [45.0, 19.0, 26.0, 50.0, -101.0, 80.0];

% and the cavity effective radii (Table 2 r_eff)
variable reffs = [15.97, 16.06, 38.58, 24.62, 102.34, 58.26];

% Now fill a v3d array with these 6 spheres:
%
% Blank array = 0 times r3d-squared values
variable r3dvals = v3d_r3dsq([0.,0.,0.]);
variable harr = 0.0*r3dvals;
variable this_cavity;

variable ic;
_for ic (0,length(rskys)-1,1)
{
  this_cavity = v3d_sphere(0.0, reffs[ic], [xskys[ic], yskys[ic], zskys[ic]] );
  % scale by 1/r so that the bigger cavities are a bit dimmer (locally)
  this_cavity = this_cavity/(reffs[ic]) ;
  harr = v3d_merge_max(harr, this_cavity);
}

% Make the AGN at 0,0,0
variable agn = v3d_sphere(0.0, 10.0);

% Make a big sphere of residual X-ray emission
variable resid = v3d_sphere(0.0, 375.0);
% and remove the cavities
resid = resid - 1.e4*harr;
% and keep just positive values
resid = v3d_merge_max(resid, 0.0*r3dvals);

% and view it 
%
% in projection
%  the cavities in green, resid in red and AGN blue:
v3d_project(resid/1000.0, harr, agn/10.0);

% and 3D arrangment of the cavities with volview:
v3d_view(harr);