Changeset 4993

Timestamp:

08/04/10 19:42:13 (15 years ago)

Author:

Eric.Larour

Message:

Added circles to test for 0 mass flux with constant velocity

Location:

issm/trunk/test/Verification/test29_MassFlux

Files:

• Archive_Diag_None_1_0_nan_0.mat (modified) ( previous)
• Archive_Diag_None_1_0_nan_1.mat (modified) ( previous)
• Lrel29.m (modified) (1 diff)
• MassFlux4.exp (added)
• MassFlux5.exp (added)
• MassFlux6.exp (added)
• README (modified) (1 diff)
• length_profiles.m (modified) (1 diff)
• testpostsolve.m (modified) (1 diff)

issm/trunk/test/Verification/test29_MassFlux/Lrel29.m

@@ -9 +9 @@
 md.responses.rf (end+1)=response_function('MassFlux2',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
 md.responses.rf (end+1)=response_function('MassFlux3',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
+md.responses.rf (end+1)=response_function('MassFlux4',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
+md.responses.rf (end+1)=response_function('MassFlux5',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
+md.responses.rf (end+1)=response_function('MassFlux6',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
 
 %mass flux profiles
-md.qmu_mass_flux_profiles={'MassFlux1.exp','MassFlux2.exp','MassFlux3.exp'};
+md.qmu_mass_flux_profiles={'MassFlux1.exp','MassFlux2.exp','MassFlux3.exp','MassFlux4.exp','MassFlux5.exp','MassFlux6.exp'};
 md.qmu_mass_flux_profile_directory=pwd;
 

issm/trunk/test/Verification/test29_MassFlux/README

@@ -1 +1 @@
 Qmu sensitivity analysis on MassFlux response. 
+MassFlux1, Massflux2 and Massflux3 are 3 sides of a triangle. MassFlux1 + MassFlux2 should be equal to MassFlux3
+MassFlux4 to 6 are circles, 10,100 and 1000 grids, .4,.3 and .2 radius. Mass flux for thoses circlyes should be 0.

issm/trunk/test/Verification/test29_MassFlux/length_profiles.m

@@ -1 @@
-for i=1:3,
-        a=md.qmu_mass_flux_segments{i};
-        %sum(sqrt((a(:,3)-a(:,1)).^2+ (a(:,4)-a(:,2)).^2))
-        [sum(a(:,3)-a(:,1)) sum(a(:,4)-a(:,2))]
+num=2;
+
+
+if num==1,
+        for i=1:length(md.qmu_mass_flux_segments),
+                a=md.qmu_mass_flux_segments{i};
+                %sum(sqrt((a(:,3)-a(:,1)).^2+ (a(:,4)-a(:,2)).^2))
+                [sum(a(:,3)-a(:,1)) sum(a(:,4)-a(:,2))]
+        end
 end
+
+
+if num==2,
+        hold on;
+        %for i=1:length(md.qmu_mass_flux_segments),
+        for i=4:length(md.qmu_mass_flux_segments),
+                a=md.qmu_mass_flux_segments{i};
+                for j=1:length(a),
+                        plot([a(j,1) a(j,3)],[a(j,2) a(j,4)]);
+                end
+        end
+end

issm/trunk/test/Verification/test29_MassFlux/testpostsolve.m

@@ -6 +6 @@
 %also, check that the stddev are 0.
 
-md.results.dakota.importancefactors=[md.results.dakota.dresp_out(1).mean md.results.dakota.dresp_out(2).mean md.results.dakota.dresp_out(3).mean md.results.dakota.dresp_out(4).mean ...
-                                     md.results.dakota.dresp_out(1).stddev md.results.dakota.dresp_out(2).stddev md.results.dakota.dresp_out(3).stddev md.results.dakota.dresp_out(4).stddev];
+md.results.dakota.importancefactors=[];
+for i=1:7,
+        md.results.dakota.importancefactors=[md.results.dakota.importancefactors md.results.dakota.dresp_out(i).mean];
+end
+for i=1:7,
+        md.results.dakota.importancefactors=[md.results.dakota.importancefactors md.results.dakota.dresp_out(i).stddev];
+end