changed output to separate raw and norm freqs
added faceting on frequency type
added a note
added a loop to accelerate calculation
updated with explicit job IDs
free y-scale for facets