A script that starts from raw DN data and produces the toa_rad, wl_rad and refl data
Subsets and converts to TOARAD
Full wang correction mostly there. Needs testing. Interpolation takes forever with current setup
Logic works, but results are weird. Could be a problem with DN_2_rad codes...
Got rid of previous do_wang.py script
Update do_wang.py to soemthing that works