replace find with some index magic

Hopefully this will be marginally faster
OOPS Correct to use phase speed in calculating wave length

Reference the correct variable
add binary outputs, increase freq resulution
Make vertical slice of gravity wave velocity field
add outputs
this will work with any domain on w and k
make it so the output is not put in the temp directory
change source dag file
update spectra generation
make a 2D lie bracket of two functions
Merge all the NCEP files into one
change default souce files for plotspecVvar.m
Break out subfunctions into their own files,

figure ploting scripts are moved into subdirectory
more general use functions are moved to thier own file in the main
directory
Plot upper 1/2 of domain
use t0sim and tfsim as time axes limits instead of hard coding
Correct the sign of Shear Production
Move helper functions into subdirectory
calculate diff between model and observations
move palete and ticbtic setting into ..limits.plt
move fig6 and DAGtkeprofiles into subdirectory