# Integrate the line area
line_area_obs = pm.Deterministic('linearea',
tt.sum(emission_obs[mask_loc]) /
(obs_wave[mask_loc[-1]] - obs_wave[mask_loc[0]]) / 1e15)
The dividing by 1e15 part.
# Integrate the line area
line_area_obs = pm.Deterministic('linearea',
tt.sum(emission_obs[mask_loc]) /
(obs_wave[mask_loc[-1]] - obs_wave[mask_loc[0]]) / 1e15)
The dividing by 1e15 part.