added binary output for energy
streamlined energy.f and added to breaking wave case
Eb calculation works but overestimates
fixed bug in phi_d calc