Index: gliderproc/trunk/gliderCTD_Generate_L1_Data.m
===================================================================
--- gliderproc/trunk/gliderCTD_Generate_L1_Data.m (revision 508)
+++ gliderproc/trunk/gliderCTD_Generate_L1_Data.m (revision 509)
@@ -23,4 +23,5 @@
 %   June 2013, implemented some QC procedures for ramses, HES
 %   20130703 - Iterate throuth gliders and deployments. - CBC
+%   20130725 - Implemented some QC procedures for pelagia. - HES
 %
 %//////////////////////////////////////////////////////////////////////////
@@ -308,4 +309,13 @@
             ptime = ptime(i);  temp = temp(i);  pres = pres(i);  cond = cond(i); 
             ctd_time = ctd_time(i);
+        elseif(gliderIndex == 1)
+            ib=find(abs(diff(temp))>1.5);
+            ib2=find(abs(diff(cond))>0.1);
+            ibb=union(ib,ib2);
+            temp(ibb+1)=NaN;
+            cond(ibb+1)=NaN;
+            i=find(~isnan(temp));
+            ptime = ptime(i);  temp = temp(i);  pres = pres(i);  cond = cond(i); 
+            ctd_time = ctd_time(i);
         end
 
@@ -352,5 +362,5 @@
         avgDepthRate = avgDepthRate(i); glideAngle = glideAngle(i); 
         ptime = ptime(i); temp = temp(i); cond = cond(i); pres = pres(i); 
-        ctd_time = ctd_time(i);
+        ctd_time = ctd_time(i); hv = hv(i);
 
         % scale up the pressure...
@@ -392,5 +402,6 @@
         % top or bottom of profiles.  Use original velocity measure (hv) and pitch
         % to identify points for exclusion
-        % these values set by look at ramses; may need alternate set for pelagia
+        % these values set by look at ramses; may need alternate set for
+        % pelagia (now set based on deployment 1)
 
         if(gliderIndex == 2)
@@ -399,12 +410,18 @@
             ib = union(iv,ip);
             salinCorrected(ib) = NaN;
-
+        elseif(gliderIndex == 1)
+            iv = find(hv<0.1);
+            ip = find(pitch>5 & pitch < 15);
+            ib = union(iv,ip);
+            salinCorrected(ib) = NaN;
+        end
+            
             % the step above likely removes many points, need to make sure that dataset
             % is consistent, so use salinity to ID valid times going forward
 
-            i = find(~isnan(salinCorrected));
-            ptime=ptime(i);  temp=temp(i);  tempCorrected=tempCorrected(i); salin=salin(i);
-            salinCorrected=salinCorrected(i); pres=pres(i); dens=dens(i);
-        end
+        i = find(~isnan(salinCorrected));
+        ptime=ptime(i);  temp=temp(i);  tempCorrected=tempCorrected(i); salin=salin(i);
+        salinCorrected=salinCorrected(i); pres=pres(i); dens=dens(i);
+        
 
         % calculate density...should this use temp corrected?? HES - no, now have