Merge branch 'main' of github.com:NOAA-PMEL/BGC_Argo_Mat_Toolbox into full_argo

Author: hfrenzel-noaa

auxil/depth_interp.m

@@ -87,6 +87,11 @@
 try
     assert(strcmp(raw, 'no'));
     X = Data.PRES_ADJUSTED;
+    good_vals = sum(isfinite(X));
+    % a somewhat arbitrary criterion: at least half of the profiles
+    % must have valid adjusted pressure values, or switch to
+    % raw pressure
+    assert(sum(good_vals) > 0.5 * length(good_vals));
 catch
     X = Data.PRES;
 end

auxil/plot_profiles.m

@@ -280,6 +280,11 @@
         try
             assert(strcmp(raw, 'no')) % use PRES if raw values are used
             PRES = Data.(floats{f}).PRES_ADJUSTED;
+            good_vals = sum(isfinite(PRES));
+            % a somewhat arbitrary criterion: at least half of the profiles
+            % must have valid adjusted pressure values, or switch to
+            % raw pressure
+            assert(sum(good_vals) > 0.5 * length(good_vals));
         catch
             PRES = Data.(floats{f}).PRES;
         end

auxil/plot_sections.m

@@ -125,9 +125,9 @@ function plot_sections(Data, Mdata, variables, nvars, plot_isopyc, ...
 
 % unless 'raw' is specified, plot adjusted data
 if strncmpi(raw,'y',1)
-    title_add = ' [raw values]';
+    [title_added{1:nvars}] = deal(' [raw values]');
 else
-    title_add = '';
+    [title_added{1:nvars}] = deal(''); % default; may change later
     for v = 1:nvars
         % if all floats have adjusted values available for a variable,
         % they will be used instead of raw values
@@ -142,7 +142,7 @@ function plot_sections(Data, Mdata, variables, nvars, plot_isopyc, ...
                 warning(['adjusted values for %s for float %s are not available,',...
                     ' showing raw value profiles instead'], ...
                     variables{v}, floats{f});
-                title_add = ' [raw values]';
+                title_added{v} = ' [raw values]';
             end
         end
         if has_adj == nfloats
@@ -208,8 +208,8 @@ function plot_sections(Data, Mdata, variables, nvars, plot_isopyc, ...
             ylim(depth);
         end
         title(sprintf('Float %d: %s %s%s', ...
-            Mdata.(float_ids{f}).WMO_NUMBER, long_name, units, title_add),...
-            'FontSize',14);
+            Mdata.(float_ids{f}).WMO_NUMBER, long_name, units, ...
+            title_added{v}), 'FontSize', 14);
         ylabel('Pressure (dbar)');
         set(gca,'Ydir','reverse');
         colorbar;

auxil/try_download.m

@@ -51,9 +51,9 @@
     catch
         if Settings.verbose
             fprintf('failure!\n');
-            if exist([dest_path, '.html'], 'file')
-                delete([dest_path, '.html']);
-            end
+        end
+        if exist([dest_path, '.html'], 'file')
+            delete([dest_path, '.html']);
         end
         if h == length(Settings.hosts)
             if exist(dest_path, 'file') ~= 2