Steps for glider diel cycle code:

1. Load in spray data.
2. Index the experimental portion of the data.
	2. Using ODV to find the date range of the data that we want.
3. QC if necessary.
4. Calculate o2 saturation concentraion from t,s,p and reshape array to match the rest of the data.
	4a. Yui uses like satconc(t,s,p)
5. Calculate o2 anomaly. 
	5a. anom = dissolved oxygen - o2 saturation concentration.
	5b. this normalizes o2 data because it is temp dependent.
6. Calculate mean of o2 saturation percentage in the upper 5m depth. (This is indexed as 1:3 or 1:2 in barone code).
7. Calculate max of PAR in upper 5m.
8. Parse_to_diel the mean of o2 saturation percentage and max of Ypar. This function basically averages data into hourly sections. (Can improve this function potentially).
9. Calculate the amplitude of diel o2 saturation percentage (max - min), the mean of diel PAR. Integrate diel PAR.
10. Fit a line to the PAR and Amplitude.
	10a. Barone fits 3 lines and uses all of them to estimate the diel cycle.
11. Calculate diff O2 anom and quality control.
12. Calculate PAR.
13. Calculate hour in local time.
14. Find indices of morning hours.
15. Define bins for PAR and bin O2 anom accordingly.