Frequency-dependent internal wave drag

Implementation of frequency-dependent drag for M2 and K1 velocities,
utilizing the streaming banspass filter. Currently, it only works
for scalar form of internal wave drag (e.g., JSL). Implementation
for directional internal wave drag will follow in another commit.