Fast empirical models of the broad emission line region (BLR) are a powerful tool to interpret velocity-resolved reverberation mapping (RM) data, estimate the mass of the supermassive black holes, and gain insight into its geometry and kinematics. Much of the effort so far has been devoted to describing the emissivity of one emission line at a time. We present here an alternative approach aimed at describing the underlying BLR gas distribution, by exploiting simple numerical recipes to connect it with emissivity. This approach is a step toward describing multiple emission lines originating from the same gas and allows us to clarify some issues related to the interpretation of RM data. We illustrate this approach—implemented in the code CARAMEL-gas—using three data sets covering the Hβ emission line (Mrk 50, Mrk 1511, Arp 151) that have been modeled using the emissivity-based version of the code. As expected, we find differences in the parameters describing the BLR gas and emissivity distribution, but the emissivity-weighted lag measurements and all other model parameters including black hole mass and overall BLR morphology and kinematics are consistent with the previous measurements. We also model the Hα emission line for Arp 151 using both the gas- and emissivity-based BLR models. We find ionization stratification in the BLR with Hα arising at larger radii than Hβ, while all other model parameters are consistent within the uncertainties.