Changeset 4281
- Timestamp:
- 06/28/10 16:19:30 (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk/src/m/utils/Nightly/testsgetfields.m
r4279 r4281 23 23 if (analysis_type==DiagnosticSolutionEnum), 24 24 if control, 25 fields={,'DiagnosticAnalysis.Vel.value','DiagnosticAnalysis.J.value','DiagnosticAnalysis.parameter.value'}; 25 field_names={'Vel','J','parameter'}; 26 vel=zeros(md.numberofgrids,1); 27 vel(md.results.DiagnosticAnalysis.Vel.index)=md.results.DiagnosticAnalysis.Vel.value; 28 j=zeros(md.numberofgrids,1); 29 j(md.results.DiagnosticAnalysis.J.index)=md.results.DiagnosticAnalysis.J.value; 30 parameter=zeros(md.numberofgrids,1); 31 parameter(md.results.DiagnosticAnalysis.parameter.index)=md.results.DiagnosticAnalysis.parameter.value; 32 field_values={vel,j,parameter}; 26 33 elseif qmu, 27 34 fields={'dakota.importancefactors'}; 28 35 else 29 36 if md.dim==3 30 fields={'DiagnosticAnalysis.Vy.value','DiagnosticAnalysis.Vz.value'}; 37 field_names={'Vy','Vz'}; 38 vy=zeros(md.numberofgrids,1); 39 vy(md.results.DiagnosticAnalysis.Vy.index)=md.results.DiagnosticAnalysis.Vy.value; 40 vz=zeros(md.numberofgrids,1); 41 vz(md.results.DiagnosticAnalysis.Vz.index)=md.results.DiagnosticAnalysis.Vz.value; 42 field_values={vy,vz}; 31 43 else 32 44 field_names={'Vel'}; … … 38 50 39 51 elseif (analysis_type==ThermalSolutionEnum), 40 fields={'ThermalSolution(end).temperature','ThermalSolution(end).melting'}; 52 field_names={'Temperature','Melting'}; 53 temperature=zeros(md.numberofgrids,1); 54 temperature(md.results.ThermalAnalysis(end).Temperature.index)=md.results.ThermalAnalysis(end).Temperature.value; 55 melting=zeros(md.numberofgrids,1); 56 melting(md.results.ThermalAnalysis(end).Melting.index)=md.results.ThermalAnalysis(end).Melting.value; 57 field_values={temperature,melting}; 41 58 42 59 elseif (analysis_type==PrognosticSolutionEnum), … … 47 64 48 65 elseif (analysis_type==Prognostic2SolutionEnum), 49 fields={'Prognostic2Solution.thickness'}; 66 field_names={'Thickness'}; 67 thickness=zeros(md.numberofgrids,1); 68 thickness(md.results.Prognostic2Analysis.Thickness.index)=md.results.Prognostic2Analysis.Thickness.value; 69 field_values={thickness}; 50 70 51 71 elseif (analysis_type==Transient2DSolutionEnum), 52 fields={'TransientSolution(end).vel','TransientSolution(end).thickness','TransientSolution(end).surface','TransientSolution(end).bed'};53 72 field_names={'Vel','Thickness','Surface','Bed'}; 54 73 vel=zeros(md.numberofgrids,1); … … 63 82 64 83 elseif (analysis_type==Transient3DSolutionEnum), 65 fields={'TransientSolution(end).vel','TransientSolution(end).pressure','TransientSolution(end).temperature','TransientSolution(end).melting','TransientSolution(end).thickness','TransientSolution(end).surface','TransientSolution(end).bed'}; 84 field_names={'Vel','Pressure','Temperature','Melting','Thickness','Surface','Bed'}; 85 vel=zeros(md.numberofgrids,1); 86 vel(md.results.Transient3DAnalysis(end).Vel.index)=md.results.Transient3DAnalysis(end).Vel.value; 87 pressure=zeros(md.numberofgrids,1); 88 pressure(md.results.Transient3DAnalysis(end).Pressure.index)=md.results.Transient3DAnalysis(end).Pressure.value; 89 temperature=zeros(md.numberofgrids,1); 90 temperature(md.results.Transient3DAnalysis(end).Temperature.index)=md.results.Transient3DAnalysis(end).Temperature.value; 91 melting=zeros(md.numberofgrids,1); 92 melting(md.results.Transient3DAnalysis(end).Melting.index)=md.results.Transient3DAnalysis(end).Melting.value; 93 thickness=zeros(md.numberofgrids,1); 94 thickness(md.results.Transient3DAnalysis(end).Thickness.index)=md.results.Transient3DAnalysis(end).Thickness.value; 95 surface=zeros(md.numberofgrids,1); 96 surface(md.results.Transient3DAnalysis(end).Surface.index)=md.results.Transient3DAnalysis(end).Surface.value; 97 bed=zeros(md.numberofgrids,1); 98 bed(md.results.Transient3DAnalysis(end).Bed.index)=md.results.Transient3DAnalysis(end).Bed.value; 99 field_values={vel,pressure,temperature,melting,thickness,surface,bed}; 66 100 67 101 elseif (analysis_type==SteadystateSolutionEnum), 68 102 if control, 69 fields={'SteadystateSolution.vel','SteadystateSolution.J','SteadystateSolution.parameter'}; 103 field_names={'Vel','J','parameter'}; 104 vel=zeros(md.numberofgrids,1); 105 vel(md.results.SteadystateAnalysis.Vel.index)=md.results.SteadystateAnalysis.Vel.value; 106 j=zeros(md.numberofgrids,1); 107 j(md.results.SteadystateAnalysis.J.index)=md.results.SteadystateAnalysis.J.value; 108 parameter=zeros(md.numberofgrids,1); 109 parameter(md.results.SteadystateAnalysis.parameter.index)=md.results.SteadystateAnalysis.parameter.value; 110 field_values={vel,j,parameter}; 70 111 else 71 fields={'SteadystateSolution.vel','SteadystateSolution.pressure','SteadystateSolution.temperature'}; 112 field_names={'Vel','Pressure','Temperature'}; 113 vel=zeros(md.numberofgrids,1); 114 vel(md.results.SteadystateAnalysis.Vel.index)=md.results.SteadystateAnalysis.Vel.value; 115 presure=zeros(md.numberofgrids,1); 116 presure(md.results.SteadystateAnalysis.Pressure.index)=md.results.SteadystateAnalysis.Pressure.value; 117 temperature=zeros(md.numberofgrids,1); 118 temperature(md.results.SteadystateAnalysis.Temperature.index)=md.results.SteadystateAnalysis.Temperature.value; 119 field_values={vel,pressure,temperature}; 72 120 end 73 121 … … 79 127 80 128 elseif (analysis_type==SlopeSolutionEnum), 81 fields={'SlopeSolution.slopex','SlopeSolution.slopey'}; 129 field_names={'Slopex','Slopey'}; 130 slopex=zeros(md.numberofgrids,1); 131 slopex(md.results.DiagnosticAnalysis.Slopex.index)=md.results.DiagnosticAnalysis.Slopex.value; 132 slopey=zeros(md.numberofgrids,1); 133 slopey(md.results.DiagnosticAnalysis.Slopey.index)=md.results.DiagnosticAnalysis.Slopey.value; 134 field_values={slopex,slopey}; 82 135 83 136 else
Note:
See TracChangeset
for help on using the changeset viewer.