 filename = 'WF-CANONCruiseSept2010_PARModifiedTime.xlsx';
 latValue = xlsread(filename,'K4:K464287');
 lonValue = xlsread(filename,'J4:J464287');
 time = xlsread(filename, 'H4:H464287');

%Define a second set of lat/lon vectors that we can edit the bad data out
%of
lonValue2 = lonValue;
latValue2 = latValue;

%Find the bad lons
badLons = find( -121 < lonValue | lonValue < - 126 );
%Get rid of any lat/lon with a bad lon
lonValue2(badLons) = NaN;
latValue2(badLons) = NaN;

%Find the bad lats
badLats = find( 36.5 < latValue | latValue < 35.0 );
%Get rid of any lat/lon with a bad lat
latValue2(badLats) = NaN;
lonValue2(badLons) = NaN;

latValue3 = latValue2(~isnan(latValue2));
lonValue3 = lonValue2(~isnan(lonValue2));

%Convert lat/lon to UTM Zone 10N, WGS84 northings and eastings
% a = 6378137;
% f = 1/298.257223563;
% e = sqrt(2*f - f^2);
% ellipsoid = [a e];
% utmstruct= defaultm('utm');
% utmstruct.zone = '10N';
% utmstruct.geoid = ellipsoid;
% utmstruct = defaultm(utm(utmstruct));
% 
% [easting, northing] = mfwdtran(utmstruct, latValue2, lonValue2);
