%
% file: v3dex_snr1987a.sl
% 
% Examples of V3D routines use applied to
%
% - - - - - SNR 1987A- - - - - 
%

% Set the scale (in x10^17 cm) and resolution:
v3d_setup(8.0,37);

% Orientation of the axis perp to equator of ring system
% The North (+Y) part of the ring is towards the observer.
variable ax_thphi = [260.0, -45.0];

% Make a dense inner ring
variable inring = v3d_torus(5.9,6.5, , ax_thphi);
% and some protrusions
variable prots = v3d_sphere_ring(5.7, 0.3, 27, , ax_thphi);
% the H II region
variable HIIreg = v3d_torus(4.5, 9.5, , ax_thphi);
% SN ejecta
variable snej = v3d_sphere_rlup(0.0,4.0, 4.5*[0.0,0.1,0.2,0.4,0.8,1.5],0.05*[10.0,5.0,3.0,1.0,0.7,0.5]);
% blastwave sphere - outside of the H II region
variable blast = v3d_merge_max( v3d_sphere(6.0,7.0) - HIIreg, 0.0*HIIreg );

% Sliced 3D view of these:
#ifexists volview
v3d_view((snej + HIIreg/10.0 + prots + inring + blast/10.0)*v3d_cube(10.0,10.0*[-1.,0.,0.]) , "v3dex_snr1987a_sliced.h5"); 
#endif

% View the components in projection with colors assigned
v3d_project( snej + 2.0*prots + inring, HIIreg/20.0, snej + blast/10.0);