Changeset 9684 for issm/trunk/test/NightlyRun/test1302.m
- Timestamp:
- 09/08/11 10:57:56 (14 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
TabularUnified issm/trunk/test/NightlyRun/test1302.m ¶
r9664 r9684 14 14 pos1=find(md.elementonbed); md.thermal.spctemperature(md.elements(pos1,1:3))=10; 15 15 pos2=find(md.elementonsurface); md.thermal.spctemperature(md.elements(pos2,4:6))=0; 16 md. vz=0.1*ones(md.numberofnodes,1);17 md. vel=sqrt( md.vx.^2+ md.vy.^2+ md.vz.^2);18 md. pressure=zeros(md.numberofnodes,1);16 md.initialization.vz=0.1*ones(md.numberofnodes,1); 17 md.initialization.vel=sqrt( md.initialization.vx.^2+ md.initialization.vy.^2+ md.initialization.vz.^2); 18 md.initialization.pressure=zeros(md.numberofnodes,1); 19 19 20 20 %analytical results … … 23 23 A=10/(exp(alpha*(-1000))-1); %A=T(bed)/(exp(alpha*bed)-1) with bed=-1000 T(bed)=10 24 24 B=-A; 25 md. temperature=A*exp(alpha*md.z)+B;25 md.initialization.temperature=A*exp(alpha*md.z)+B; 26 26 27 27 %modeled results … … 32 32 comp_temp=zeros(md.numberofnodes,1); 33 33 comp_temp(md.results.ThermalSolution.Temperature.index)=md.results.ThermalSolution.Temperature.value; 34 relative=abs((comp_temp-md. temperature)./md.temperature)*100;35 relative(find(comp_temp==md. temperature))=0;34 relative=abs((comp_temp-md.initialization.temperature)./md.initialization.temperature)*100; 35 relative(find(comp_temp==md.initialization.temperature))=0; 36 36 plotmodel(md,'data',comp_temp,'title','Modeled temperature [K]','data','temperature','view',3,... 37 'title','Analytical temperature [K]','view',3,'data',comp_temp-md. temperature,...37 'title','Analytical temperature [K]','view',3,'data',comp_temp-md.initialization.temperature,... 38 38 'title','Absolute error [K]','view',3,'data',relative,'title','Relative error [%]','view',3,... 39 39 'figposition','mathieu','FontSize#all',20)
Note:
See TracChangeset
for help on using the changeset viewer.