Processing HI cube - II

A) We now have a cube V_los(x,y,z)

1. Compute V_min and V_max (extrema) of V_los matrix,
2. Create new array V_obs, where velocity values (V_los) range from V_min to V_max, in steps of 1 km/s,

B) We want an output image matrix, OUT(x,y,z), where x_i corresponds to V_obs(i)

1. Scan the V_los(x,y,z) matrix; say, V = V_los(x,y,z)
2. If V1 < V < V2, then OUT(v,y,z) += rho(x,y,z)
   
3. Repeat a & b for energy and pressure.