Changeset 4875
- Timestamp:
- 07/29/10 10:09:21 (15 years ago)
- Location:
- issm/trunk/test/Validation/PattynStokes
- Files:
-
- 5 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk/test/Validation/PattynStokes/test1_iceshelf/runme.m
r899 r4875 11 11 md=setelementstype(md,'Pattyn','all'); 12 12 13 md=solve(md,'analysis_type','diagnostic','package','ice'); 14 vel_pattyn=md.results.diagnostic.vel; 15 vx_pattyn=md.results.diagnostic.vx; 16 vy_pattyn=md.results.diagnostic.vy; 17 vz_pattyn=md.results.diagnostic.vz; 18 pressure_pattyn=md.results.diagnostic.pressure; 13 md=solve(md,'analysis_type',DiagnosticSolutionEnum); 14 vel_pattyn=zeros(md.numberofgrids,1); 15 vel_pattyn(md.results.DiagnosticSolution.Vel.index)=md.results.DiagnosticSolution.Vel.value; 16 vx_pattyn=zeros(md.numberofgrids,1); 17 vx_pattyn(md.results.DiagnosticSolution.Vx.index)=md.results.DiagnosticSolution.Vx.value; 18 vy_pattyn=zeros(md.numberofgrids,1); 19 vy_pattyn(md.results.DiagnosticSolution.Vy.index)=md.results.DiagnosticSolution.Vy.value; 20 vz_pattyn=zeros(md.numberofgrids,1); 21 vz_pattyn(md.results.DiagnosticSolution.Vz.index)=md.results.DiagnosticSolution.Vz.value; 22 pressure_pattyn=zeros(md.numberofgrids,1); 23 pressure_pattyn(md.results.DiagnosticSolution.Pressure.index)=md.results.DiagnosticSolution.Pressure.value; 19 24 20 25 mdpattyn=md; … … 29 34 md=setelementstype(md,'Stokes','all','Pattyn','all'); 30 35 31 md=solve(md,'analysis_type','diagnostic','package','ice'); 32 vel_stokes=md.vel; 33 vx_stokes=md.vx; 34 vy_stokes=md.vy; 35 vz_stokes=md.vz; 36 pressure_stokes=md.pressure; 36 md=solve(md,'analysis_type',DiagnosticSolutionEnum); 37 vel_stokes=zeros(md.numberofgrids,1); 38 vel_stokes(md.results.DiagnosticSolution.Vel.index)=md.results.DiagnosticSolution.Vel.value; 39 vx_stokes=zeros(md.numberofgrids,1); 40 vx_stokes(md.results.DiagnosticSolution.Vx.index)=md.results.DiagnosticSolution.Vx.value; 41 vy_stokes=zeros(md.numberofgrids,1); 42 vy_stokes(md.results.DiagnosticSolution.Vy.index)=md.results.DiagnosticSolution.Vy.value; 43 vz_stokes=zeros(md.numberofgrids,1); 44 vz_stokes(md.results.DiagnosticSolution.Vz.index)=md.results.DiagnosticSolution.Vz.value; 45 pressure_stokes=zeros(md.numberofgrids,1); 46 pressure_stokes(md.results.DiagnosticSolution.Pressure.index)=md.results.DiagnosticSolution.Pressure.value; 37 47 38 48 mdstokes=md; 39 save modelstokes mdstokes; 49 save modelstokes mdstokes 50 51 %plot the two velocities 52 plotmodel(md,'data',vel_pattyn,'data',vel_stokes,'data',abs((vel_pattyn-vel_stokes)./vel_pattyn)*100,'data',vel_pattyn-vel_stokes,'title','Pattyn velocity','title','Stokes velocity','title','Relative difference','title','Absolute difference','view#all',2) -
issm/trunk/test/Validation/PattynStokes/test2_icesheet/Square.par
r3768 r4875 34 34 md.np=3; 35 35 md.time=50; 36 md.waitonlock=1 ;36 md.waitonlock=10; -
issm/trunk/test/Validation/PattynStokes/test2_icesheet/runme.m
r911 r4875 11 11 md=setelementstype(md,'Pattyn','all'); 12 12 13 md=solve(md,'analysis_type','diagnostic','package','ice'); 14 vel_pattyn=md.results.diagnostic.vel; 15 vx_pattyn=md.results.diagnostic.vx; 16 vy_pattyn=md.results.diagnostic.vy; 17 vz_pattyn=md.results.diagnostic.vz; 18 pressure_pattyn=md.results.diagnostic.pressure; 13 md=solve(md,'analysis_type',DiagnosticSolutionEnum); 14 vel_pattyn=zeros(md.numberofgrids,1); 15 vel_pattyn(md.results.DiagnosticSolution.Vel.index)=md.results.DiagnosticSolution.Vel.value; 16 vx_pattyn=zeros(md.numberofgrids,1); 17 vx_pattyn(md.results.DiagnosticSolution.Vx.index)=md.results.DiagnosticSolution.Vx.value; 18 vy_pattyn=zeros(md.numberofgrids,1); 19 vy_pattyn(md.results.DiagnosticSolution.Vy.index)=md.results.DiagnosticSolution.Vy.value; 20 vz_pattyn=zeros(md.numberofgrids,1); 21 vz_pattyn(md.results.DiagnosticSolution.Vz.index)=md.results.DiagnosticSolution.Vz.value; 22 pressure_pattyn=zeros(md.numberofgrids,1); 23 pressure_pattyn(md.results.DiagnosticSolution.Pressure.index)=md.results.DiagnosticSolution.Pressure.value; 19 24 20 25 mdpattyn=md; … … 29 34 md=setelementstype(md,'Pattyn','all','Stokes','all'); 30 35 31 md=solve(md,'analysis_type','diagnostic','package','ice'); 32 vel_stokes=md.results.diagnostic.vel; 33 vx_stokes=md.results.diagnostic.vx; 34 vy_stokes=md.results.diagnostic.vy; 35 vz_stokes=md.results.diagnostic.vz; 36 pressure_stokes=md.results.diagnostic.pressure; 36 md=solve(md,'analysis_type',DiagnosticSolutionEnum); 37 vel_stokes=zeros(md.numberofgrids,1); 38 vel_stokes(md.results.DiagnosticSolution.Vel.index)=md.results.DiagnosticSolution.Vel.value; 39 vx_stokes=zeros(md.numberofgrids,1); 40 vx_stokes(md.results.DiagnosticSolution.Vx.index)=md.results.DiagnosticSolution.Vx.value; 41 vy_stokes=zeros(md.numberofgrids,1); 42 vy_stokes(md.results.DiagnosticSolution.Vy.index)=md.results.DiagnosticSolution.Vy.value; 43 vz_stokes=zeros(md.numberofgrids,1); 44 vz_stokes(md.results.DiagnosticSolution.Vz.index)=md.results.DiagnosticSolution.Vz.value; 45 pressure_stokes=zeros(md.numberofgrids,1); 46 pressure_stokes(md.results.DiagnosticSolution.Pressure.index)=md.results.DiagnosticSolution.Pressure.value; 37 47 38 48 mdstokes=md; 39 49 save modelstokes mdstokes; 40 50 51 %plot the two velocities 52 plotmodel(md,'data',vel_pattyn,'data',vel_stokes,'data',abs((vel_pattyn-vel_stokes)./vel_pattyn)*100,'data',vel_pattyn-vel_stokes,'title','Pattyn velocity','title','Stokes velocity','title','Relative difference','title','Absolute difference','view#all',2) -
issm/trunk/test/Validation/PattynStokes/test3_icesheet_iceshelf/Square.par
r3768 r4875 34 34 md.np=3; 35 35 md.time=50; 36 md.waitonlock=1 ;36 md.waitonlock=10; -
issm/trunk/test/Validation/PattynStokes/test3_icesheet_iceshelf/runme.m
r899 r4875 11 11 md=setelementstype(md,'Pattyn','all'); 12 12 13 md=solve(md,'analysis_type','diagnostic','package','ice'); 14 vel_pattyn=md.results.diagnostic.vel; 15 vx_pattyn=md.results.diagnostic.vx; 16 vy_pattyn=md.results.diagnostic.vy; 17 vz_pattyn=md.results.diagnostic.vz; 18 pressure_pattyn=md.results.diagnostic.pressure; 13 md=solve(md,'analysis_type',DiagnosticSolutionEnum); 14 vel_pattyn=zeros(md.numberofgrids,1); 15 vel_pattyn(md.results.DiagnosticSolution.Vel.index)=md.results.DiagnosticSolution.Vel.value; 16 vx_pattyn=zeros(md.numberofgrids,1); 17 vx_pattyn(md.results.DiagnosticSolution.Vx.index)=md.results.DiagnosticSolution.Vx.value; 18 vy_pattyn=zeros(md.numberofgrids,1); 19 vy_pattyn(md.results.DiagnosticSolution.Vy.index)=md.results.DiagnosticSolution.Vy.value; 20 vz_pattyn=zeros(md.numberofgrids,1); 21 vz_pattyn(md.results.DiagnosticSolution.Vz.index)=md.results.DiagnosticSolution.Vz.value; 22 pressure_pattyn=zeros(md.numberofgrids,1); 23 pressure_pattyn(md.results.DiagnosticSolution.Pressure.index)=md.results.DiagnosticSolution.Pressure.value; 19 24 20 25 mdpattyn=md; … … 29 34 md=setelementstype(md,'Pattyn','all','Stokes','all'); 30 35 31 md=solve(md,'analysis_type','diagnostic','package','ice'); 32 vel_stokes=md.results.diagnostic.vel; 33 vx_stokes=md.results.diagnostic.vx; 34 vy_stokes=md.results.diagnostic.vy; 35 vz_stokes=md.results.diagnostic.vz; 36 pressure_stokes=md.results.diagnostic.pressure; 36 md=solve(md,'analysis_type',DiagnosticSolutionEnum); 37 vel_stokes=zeros(md.numberofgrids,1); 38 vel_stokes(md.results.DiagnosticSolution.Vel.index)=md.results.DiagnosticSolution.Vel.value; 39 vx_stokes=zeros(md.numberofgrids,1); 40 vx_stokes(md.results.DiagnosticSolution.Vx.index)=md.results.DiagnosticSolution.Vx.value; 41 vy_stokes=zeros(md.numberofgrids,1); 42 vy_stokes(md.results.DiagnosticSolution.Vy.index)=md.results.DiagnosticSolution.Vy.value; 43 vz_stokes=zeros(md.numberofgrids,1); 44 vz_stokes(md.results.DiagnosticSolution.Vz.index)=md.results.DiagnosticSolution.Vz.value; 45 pressure_stokes=zeros(md.numberofgrids,1); 46 pressure_stokes(md.results.DiagnosticSolution.Pressure.index)=md.results.DiagnosticSolution.Pressure.value; 37 47 38 48 mdstokes=md; 39 49 save modelstokes mdstokes; 50 51 %plot the two velocities 52 plotmodel(md,'data',vel_pattyn,'data',vel_stokes,'data',abs((vel_pattyn-vel_stokes)./vel_pattyn)*100,'data',vel_pattyn-vel_stokes,'title','Pattyn velocity','title','Stokes velocity','title','Relative difference','title','Absolute difference','view#all',2)
Note:
See TracChangeset
for help on using the changeset viewer.