Changeset 44


Ignore:
Timestamp:
04/24/09 15:18:31 (16 years ago)
Author:
seroussi
Message:

removed accelerated triaelem model in tests 2d

Location:
issm/trunk/test/Verification
Files:
16 edited

Legend:

Unmodified
Added
Removed

  • issm/trunk/test/Verification/test1_icesheet_noicefront/runme.m

    r16 r44  
    1 
    2 % This file can be run to check that the current version of macayeal and ice models are giving
     1% This file can be run to check that the current version of macayeal and issm models are giving
    32% coherent results. This test deals with an icesheet without icefront for a 2d model. The geometry
    4 % is square. Just run this file in Matlab, with a properly setup Ice code.
     3% is square. Just run this file in Matlab, with a properly setup ISSM code.
    54% The results of this test will indicate if there is a difference between current computations
    65% and archived results.
    … …  
    3332md=solve(md,'diagnostic','ice');
    3433vel_ice=md.vel;
    35 
    36 %Compute solution with accelerated Ice model
    37 md.acceleration=1;
    38 md=solve(md,'diagnostic','ice');
    39 vel_acc=md.vel;
    4034
    4135%Compute solution with Hutter elements in Ice model
    … …  
    7367end
    7468
    75 %Accelerated Ice
    76 load archive_acc
    77 pos=find(archive_acc-vel_acc);
    78 error_vel=abs(norm((archive_acc(pos)-vel_acc(pos))./archive_acc(pos),2));
    79 
    80 disp(sprintf('numerical difference between old and new version of Accelerated Ice is : %d', error_vel));
    81 
    82 if (error_vel>tolerance);
    83         logoutput(logstring,sprintf('%s\n','    ERROR. Results from Accelerated Ice model differ from the archive version'));
    84 else
    85         logoutput(logstring,sprintf('%s\n','    SUCCESS. Results from Accelerated Ice model are comform with the archive version'));
    86 end
    87 
    8869%Ice with Hutter elements
    8970load archive_hutter

  • issm/trunk/test/Verification/test1_icesheet_noicefront/updatearchive.m

    r16 r44  
    2020
    2121%Compute solution with Ice model
    22 md.acceleration=0;
    2322md=solve(md,'diagnostic','ice');
    2423archive_ice=md.vel;
    … …  
    2625%Save the solution in the directory
    2726save archive_ice archive_ice
    28 
    29 %Compute solution with accelerated Ice model
    30 md.acceleration=1;
    31 md=solve(md,'diagnostic','ice');
    32 archive_acc=md.vel;
    33 
    34 %Save the solution in the directory
    35 save archive_acc archive_acc
    3627
    3728%Compute solution with Hutter elements in Ice model

  • issm/trunk/test/Verification/test2_icesheet_icefront/runme.m

    r16 r44  
    1616tolerance=10^-13;
    1717
    18 
    1918md=model;
    2019md=mesh(md,'DomainOutline.exp',50000);
    … …  
    2827
    2928%Compute solution with Ice model
    30 md.acceleration=0;
    3129md=solve(md,'diagnostic','ice');
    3230vel_ice=md.vel;
    33 
    34 %Compute solution with accelerated Ice model
    35 md.acceleration=1;
    36 md=solve(md,'diagnostic','ice');
    37 vel_acc=md.vel;
    3831
    3932%Compute prognostic with Ice model
    … …  
    7669end
    7770
    78 %Accelerated Ice
    79 load archive_acc
    80 pos=find(archive_acc-vel_acc);
    81 error_vel=abs(norm((archive_acc(pos)-vel_acc(pos))./archive_ice(pos),2));
    82 
    83 disp(sprintf('numerical difference between old and new version of Accelerated Ice is : %d', error_vel));
    84 
    85 if (error_vel>tolerance);
    86         logoutput(logstring,sprintf('%s\n','    ERROR. Results from Accelerated Ice model differ from the archive version'));
    87 else
    88         logoutput(logstring,sprintf('%s\n','    SUCCESS. Results from Accelerated Ice model are comform with the archive version'));
    89 end
    90 
    9171%Ice prognostic
    9272load archive_thickness

  • issm/trunk/test/Verification/test2_icesheet_icefront/updatearchive.m

    r16 r44  
    2020
    2121%Compute solution with Ice model
    22 md.acceleration=0;
    2322md=solve(md,'diagnostic','ice');
    2423archive_ice=md.vel;
    … …  
    2625%Save the solution in the directory
    2726save archive_ice archive_ice
    28 
    29 %Compute solution with accelerated Ice model
    30 md.acceleration=1;
    31 md=solve(md,'diagnostic','ice');
    32 archive_acc=md.vel;
    33 
    34 %Save the solution in the directory
    35 save archive_acc archive_acc
    3627
    3728%Compute prognostic with Ice model

  • issm/trunk/test/Verification/test3_iceshelf_icefront/runme.m

    r16 r44  
    2929
    3030%Compute solution with Ice model
    31 md.acceleration=0;
    3231md=solve(md,'diagnostic','ice');
    3332vel_ice=md.vel;
    34 
    35 %Compute solution with accelerated Ice model
    36 md.acceleration=1;
    37 md=solve(md,'diagnostic','ice');
    38 vel_acc=md.vel;
    3933
    4034%Compute prognostic with Ice model
    … …  
    7266end
    7367
    74 %Accelerated Ice
    75 load archive_acc
    76 pos=find(archive_acc-vel_acc);
    77 error_vel=abs(norm((archive_acc(pos)-vel_acc(pos))./archive_ice(pos),2));
    78 
    79 disp(sprintf('numerical difference between old and new version of Accelerated Ice is : %d', error_vel));
    80 
    81 if (error_vel>tolerance);
    82         logoutput(logstring,sprintf('%s\n','    ERROR. Results from Accelerated Ice model differ from the archive version'));
    83 else
    84         logoutput(logstring,sprintf('%s\n','    SUCCESS. Results from Accelerated Ice model are comform with the archive version'));
    85 end
    86 
    8768%Ice prognostic
    8869load archive_thickness

  • issm/trunk/test/Verification/test3_iceshelf_icefront/updatearchive.m

    r16 r44  
    2020
    2121%Compute solution with Ice model
    22 md.acceleration=0;
    2322md=solve(md,'diagnostic','ice');
    2423archive_ice=md.vel;
    … …  
    2625%Save the solution in the directory
    2726save archive_ice archive_ice
    28 
    29 %Compute solution with accelerated Ice model
    30 md.acceleration=1;
    31 md=solve(md,'diagnostic','ice');
    32 archive_acc=md.vel;
    33 
    34 %Save the solution in the directory
    35 save archive_acc archive_acc
    3627
    3728%Compute prognostic with Ice model
Note: See TracChangeset for help on using the changeset viewer.