Changeset 509

Timestamp:

08/04/13 19:19:35

Author:

cbc

Message:

Add Harvey's CTD QA code for Pelagia.

Files:

• gliderproc/trunk/gliderCTD_Generate_L1_Data.m (modified) (5 diffs)

gliderproc/trunk/gliderCTD_Generate_L1_Data.m

@@ -23 +23 @@
-%   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 +309 @@
-            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 +362 @@
-        avgDepthRate = avgDepthRate(i); glideAngle = glideAngle(i); 
-        ptime = ptime(i); temp = temp(i); cond = cond(i); pres = pres(i); 
-
+ hv = hv(i);
-
-        % scale up the pressure...
@@ -392 +402 @@
-        % top or bottom of profiles.  Use original velocity measure (hv) and pitch
-        % to identify points for exclusion
-
+
+
-
-        if(gliderIndex == 2)
@@ -399 +410 @@
-            ib = union(iv,ip);
-            salinCorrected(ib) = NaN;
-
+
+
+
+
+
+
+
-            % the step above likely removes many points, need to make sure that dataset
-            % is consistent, so use salinity to ID valid times going forward
-
-    
-    
-    
-end
+
+
+
+
-
-        % calculate density...should this use temp corrected?? HES - no, now have