Changeset 27889 for issm/trunk/src/m/contrib/musselman/read_netCDF.m

Timestamp:

08/30/23 15:08:57 (19 months ago)

Author:

musselman

Message:

Removed global vars from all files. Updated read log statements to use "load" instead of "save"

File:

• issm/trunk/src/m/contrib/musselman/read_netCDF.m (modified) (18 diffs)

issm/trunk/src/m/contrib/musselman/read_netCDF.m

@@ -16 +16 @@
 
 
-function model_copy = read_netCDF(filename)
+function model_copy = read_netCDF(filename, varargin)
     if nargin > 1
         verbose = true;
@@ -27 +27 @@
     end
     % make a model framework to fill that is in the scope of this file
-    global model_copy;
     model_copy = model();
 
@@ -37 +36 @@
 
         % Open the given netCDF4 file
-        global NCData;
         NCData = netcdf.open(filename, 'NOWRITE');
         % Remove masks from netCDF data for easy conversion: NOT WORKING
@@ -44 +42 @@
         % see if results is in there, if it is we have to instantiate some classes
         try
-            results_group_id = netcdf.inqNcid(NCData, "results", verbose);
-            make_results_subclasses(verbose);
+            results_group_id = netcdf.inqNcid(NCData, "results");
+            model_copy = make_results_subclasses(model_copy, NCData, verbose);
         catch
         end % 'results' group doesn't exist 
@@ -52 +50 @@
         try
             inversion_group_id = netcdf.inqNcid(NCData, "inversion");
-            check_inversion_class(verbose);
+            model_copy = check_inversion_class(model_copy, NCData, verbose);
         catch
         end % 'inversion' group doesn't exist 
@@ -61 +59 @@
             %disp(netcdf.inqGrpNameFull(group_id))
             % hand off first level to recursive search
-            walk_nested_groups(group_id, verbose);
+            model_copy = walk_nested_groups(group_id, model_copy, NCData, verbose);
         end
         
@@ -75 +73 @@
 
 
-function make_results_subclasses(verbose)
-    global model_copy;
-    global NCData;
+function model_copy = make_results_subclasses(model_copy, NCData, verbose)
     resultsGroup = netcdf.inqNcid(NCData, "results");
     variables = netcdf.inqVarIDs(resultsGroup);
@@ -90 +86 @@
         %model_copy.results = setfield(model_copy.results, class_instance, class_instance_name);
     end
+    model_copy = model_copy;
     if verbose
         disp('Successfully recreated results structs:')
@@ -99 +96 @@
 
 
-function check_inversion_class(verbose)
+function model_copy = check_inversion_class(model_copy, NCData, verbose)
     % get the name of the inversion class: either inversion or m1qn3inversion or taoinversion
-    global model_copy;
-    global NCData;
     inversionGroup = netcdf.inqNcid(NCData, "inversion");
     varid = netcdf.inqVarID(inversionGroup, 'inversion_class_name');
@@ -121 +116 @@
         end
     end
-end
-
-
-function walk_nested_groups(group_location_in_file, verbose)
-    global model_copy;
-    global NCData;    
+    model_copy = model_copy;
+end
+
+
+function model_copy = walk_nested_groups(group_location_in_file, model_copy, NCData, verbose)  
     % we search the current group level for variables by getting this struct
     variables = netcdf.inqVarIDs(group_location_in_file); 
@@ -137 +131 @@
         if strcmp(varname, 'this_is_a_nested')
             is_nested = true;
-            copy_nested_struct(group_location_in_file, verbose)
+            model_copy = copy_nested_struct(group_location_in_file, model_copy, NCData, verbose);
         elseif strcmp(varname, 'solution')
             % band-aid pass..
         else
-            copy_variable_data_to_new_model(group_location_in_file, varname, xtype, verbose);
+            model_copy = copy_variable_data_to_new_model(group_location_in_file, varname, xtype, model_copy, NCData, verbose);
         end
     end
@@ -157 +151 @@
                 group_id = netcdf.inqNcid(group_location_in_file, netcdf.inqGrpName(group));
                 %disp(netcdf.inqGrpNameFull(group_id))
-                walk_nested_groups(group, verbose);
+                model_copy = walk_nested_groups(group, model_copy, NCData, verbose);
             end
         end
@@ -167 +161 @@
 
 
-function copy_nested_struct(group_location_in_file, verbose)
-    global model_copy;
-    global NCData;
+function model_copy = copy_nested_struct(group_location_in_file, model_copy, NCData, verbose)
     %{
         A common multidimensional struct array is the 1xn md.results.TransientSolution struct. 
@@ -235 +227 @@
         model_copy.(address_in_model).(name_of_struct)(current_layer);
         if verbose
-            fprintf("Successfully saved layer %s to multidimension struct array\n", num2str(current_layer))
-        end
-    end
+            fprintf("Successfully loaded layer %s to multidimension struct array\n", num2str(current_layer))
+        end
+    end
+    model_copy = model_copy;
     if verbose
         fprintf('Successfully recreated multidimensional structure array %s in md.%s\n', name_of_struct, address_in_model)
@@ -252 +245 @@
 %}
 
-function copy_variable_data_to_new_model(group_location_in_file, varname, xtype, verbose)
-    global model_copy;
-    global NCData;
+function model_copy = copy_variable_data_to_new_model(group_location_in_file, varname, xtype, model_copy, NCData, verbose)
     %disp(varname)
     % this is an inversion band-aid
@@ -305 +296 @@
                 arg_to_eval = ['model_copy', address_to_attr, '.', varname, ' = ' , 'double(data);'];
                 eval(arg_to_eval);
-                %disp('saved int64 as int16')
+                %disp('Loaded int64 as int16')
             else
                 arg_to_eval = ['model_copy', address_to_attr, '.', varname, ' = data;'];
@@ -315 +306 @@
                 %disp(xtype)
                 %class(data)
-                fprintf('Successfully saved %s to %s\n', varname, full_addy);
+                fprintf('Successfully loaded %s to %s\n', varname, full_addy);
             end
 
@@ -332 +323 @@
         end
     end
-end
+    model_copy = model_copy;
+end