Changeset 4773 for issm/trunk/src/m/qmu/process_qmu_response_data.m

Timestamp:

07/22/10 18:57:52 (15 years ago)

Author:

Eric.Larour

Message:

New MeshProfileIntersection module, with local routines.
Needed a new object, Segment, because this is highly dynamic in nature.

More qmu debugging, especially dealing with segments for Mass flux computation.

File:

• issm/trunk/src/m/qmu/process_qmu_response_data.m (modified) (1 diff)

issm/trunk/src/m/qmu/process_qmu_response_data.m

@@ -7 +7 @@
 
 
+%preliminary data
+process_mass_flux_profiles=0;
 
+
+%loop through response descriptors, and act accordingly
 for i=1:numel(md.responsedescriptors),
-        if strncmpi(md.responsedescriptors{i},'mass_flux',9),
-                indx=str2int(md.responsedescriptors{i}(10:end));
-                if isempty(indx) || ~indx
-                        indx=1;
-                end
 
-                %we need a profile of points on which to compute the mass_flux, is it here? 
-                if isnans(md.qmu_mass_flux_profile),
-                        error('process_qmu_response_data error message: could not find a mass_flux exp profile!');
-                end
-                if ~(ischar(md.qmu_mass_flux_profile) && (indx == 1)) && ...
-                   ~(iscell(md.qmu_mass_flux_profile) && length(md.qmu_mass_flux_profile) >= indx && ischar(md.qmu_mass_flux_profile{indx})),
-                        error('process_qmu_response_data error message: mass_flux exp profile should be a domain outline name');
-                end
-
-                %ok, process the qmu_mass_flux_profile to build a list of segments: 
-                if ischar(md.qmu_mass_flux_profile),
-                        md.qmu_mass_flux_segments=MassFluxProcessProfile(md);
-                        md.qmu_mass_flux_num_segments=size(md.qmu_mass_flux_segments,1);
-                else
-                        if isempty(md.qmu_mass_flux_segments) || ~iscell(md.qmu_mass_flux_segments)
-                                md.qmu_mass_flux_segments=cell(numel(md.qmu_mass_flux_profile),1);
-                        end
-                        if isempty(md.qmu_mass_flux_num_segments) || ...
-                           (numel(md.qmu_mass_flux_num_segments) == 1 && isnan(md.qmu_mass_flux_num_segments))
-                                md.qmu_mass_flux_num_segments=zeros(numel(md.qmu_mass_flux_profile),1);
-                        end
-                        save=md.qmu_mass_flux_profile;
-                        md.qmu_mass_flux_profile=md.qmu_mass_flux_profile{indx};
-                        md.qmu_mass_flux_segments{indx,1}=MassFluxProcessProfile(md);
-                        md.qmu_mass_flux_num_segments(indx,1)=size(md.qmu_mass_flux_segments{indx},1);
-                        md.qmu_mass_flux_profile=save;
-                        clear save
-                end
+        %Do we have to process  mass flux profiles?
+        if strncmpi(md.responsedescriptors{i},'MassFlux',8),
+                process_mass_flux_profiles=1;
         end
 end
 
-% for now, concatenate all segments from cells to make a double array for parallel
-if iscell(md.qmu_mass_flux_segments)
-        segments=zeros(0,5);
-        for i=1:length(md.qmu_mass_flux_num_segments)
-                if md.qmu_mass_flux_num_segments(i)
-                        segments=[segments; md.qmu_mass_flux_segments{i}];
-                end
+
+%deal with mass flux profiles
+if process_mass_flux_profiles,
+
+        %we need a profile of points on which to compute the mass_flux, is it here? 
+        if isnans(md.qmu_mass_flux_profiles),
+                error('process_qmu_response_data error message: could not find a mass_flux exp profile!');
         end
-        md.qmu_mass_flux_segments=segments;
-        clear segments
+        
+        if ~iscell(md.qmu_mass_flux_profiles),
+                error('process_qmu_response_data error message: qmu_mass_flux_profiles field should be a cell array of domain outline names');
+        end
+
+        if isempty(md.qmu_mass_flux_profiles),
+                error('process_qmu_response_data error message: qmu_mass_flux_profiles cannot be empty!');
+        end
+
+        %ok, process the domains named in qmu_mass_flux_profiles,  to build a list of segments
+        md.qmu_mass_flux_segments=cell(numel(md.qmu_mass_flux_profiles),1);
+
+        for i=1:numel(md.qmu_mass_flux_profiles),
+                md.qmu_mass_flux_segments{i}=MassFluxProcessProfile(md,md.qmu_mass_flux_profile_directory,md.qmu_mass_flux_profiles{i});
+        end
 end
 
+error;
+