Changeset 5271


Ignore:
Timestamp:
08/16/10 10:32:42 (15 years ago)
Author:
seroussi
Message:

thermal validation continued

Location:
issm/trunk/test/NightlyRun
Files:
5 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk/test/NightlyRun/runme.m

    r5104 r5271  
    1313%      'exclude'       ids to be excluded from the test
    1414%      'benchmark'     'nightly' (nightly run/ daily run)
     15%                      'ismip'  : validation of ismip-hom tests
     16%                      'eismint': validation of eismint tests
     17%                      'thermal': validation of thermal tests
     18%                      'mesh'   : validation of mesh tests
    1519%                      ...
    1620%      'procedure'     'check' : run the test (default)
     
    4347%GET benchmark {{{1
    4448benchmark=getfieldvalue(options,'benchmark','nightly');
    45 if ~ismember(benchmark,{'all','nightly'})
     49if ~ismember(benchmark,{'all','nightly','ismip','eismint','thermal','mesh'})
    4650        disp('runme warnig: benchmark not supported, defaulting to test ''nightly''')
    4751        benchmark='nighlty';
     
    9195if strcmpi(benchmark,'nightly'),
    9296        test_ids=intersect(test_ids,[1:999]);
     97else if strcmpi(benchmark,'ismip'),
     98        test_ids=intersect(test_ids,[1101:1199]);
     99else if strcmpi(benchmark,'eismint'),
     100        test_ids=intersect(test_ids,[1201:1299]);
     101else if strcmpi(benchmark,'thermal'),
     102        test_ids=intersect(test_ids,[1301:1399]);
     103else if strcmpi(benchmark,'mesh'),
     104        test_ids=intersect(test_ids,[1401:1499]);
    93105end
    94106% }}}
  • issm/trunk/test/NightlyRun/test1301.m

    r5098 r5271  
    3232relative=abs((comp_melting-melting)./melting)*100;
    3333relative(find(comp_melting==melting))=0;
    34 plotmodel(md,'data',comp_melting,'title','modeled melting','data',melting,'title','analytical melting','data',comp_melting-melting,'title','absolute error','data',relative,'title','relative error [%]','layer#all',1,'caxis#2',[1.02964 1.02966]*10^-4)
     34plotmodel(md,'data',comp_melting,'title','Modeled melting','data',melting,'title','Analytical melting',...
     35        'data',comp_melting-melting,'title','Absolute error','data',relative,'title','Relative error [%]',...
     36        'layer#all',1,'caxis#2',[1.02964 1.02966]*10^-4,'FontSize#all',20,'figposition','mathieu','visible','off')
     37set(gcf,'Color','w')
     38printmodel('thermalmelting','png','margin','on','marginsize',25,'resolution',0.7,'hardcopy','off');
     39
    3540
    3641%Fields and tolerances to track changes
  • issm/trunk/test/NightlyRun/test1302.m

    r5098 r5271  
    3333relative=abs((comp_temp-md.observed_temperature)./md.observed_temperature)*100;
    3434relative(find(comp_temp==md.observed_temperature))=0;
    35 plotmodel(md,'data',comp_temp,'title','modeled temperature','data','observed_temperature','view',3,'title','analytical temperature','view',3,'data',comp_temp-md.observed_temperature,'title','absolute error','view',3,'data',relative,'title','relative error [%]','view',3)
     35plotmodel(md,'data',comp_temp,'title','Modeled temperature [K]','data','observed_temperature','view',3,...
     36        'title','Analytical temperature [K]','view',3,'data',comp_temp-md.observed_temperature,...
     37        'title','Absolute error [K]','view',3,'data',relative,'title','Relative error [%]','view',3,...
     38        'visible','off','figposition','mathieu','FontSize#all',20)
     39set(gcf,'Color','w')
     40printmodel('thermaladvection','png','margin','on','marginsize',25,'resolution',0.7,'hardcopy','off');
    3641
    3742%Fields and tolerances to track changes
  • issm/trunk/test/NightlyRun/test1303.m

    r5098 r5271  
    2828relative=abs((comp_temp-md.observed_temperature)./md.observed_temperature)*100;
    2929relative(find(comp_temp==md.observed_temperature))=0;
    30 plotmodel(md,'data',comp_temp,'title','modeled temperature','data','observed_temperature','view',3,'title','analytical temperature','view',3,'data',comp_temp-md.observed_temperature,'title','absolute error','view',3,'data',relative,'title','relative error [%]','view',3)
     30plotmodel(md,'data',comp_temp,'title','Modeled temperature [K]','data','observed_temperature','view',3,...
     31        'title','Analytical temperature [K]','view',3,'data',comp_temp-md.observed_temperature,...
     32        'title','Absolute error [K]','view',3,'data',relative,'title','Relative error [%]','view',3,...
     33        'visible','off','figposition','mathieu','FontSize#all',20)
     34set(gcf,'Color','w')
     35printmodel('thermalconduction','png','margin','on','marginsize',25,'resolution',0.7,'hardcopy','off');
    3136
    3237%Fields and tolerances to track changes
  • issm/trunk/test/NightlyRun/test1304.m

    r5098 r5271  
    2828relative=abs((comp_temp-md.observed_temperature)./md.observed_temperature)*100;
    2929relative(find(comp_temp==md.observed_temperature))=0;
    30 plotmodel(md,'data',comp_temp,'title','modeled temperature','data','observed_temperature','view',3,'title','analytical temperature','view',3,'data',comp_temp-md.observed_temperature,'title','absolute error','view',3,'data',relative,'title','relative error [%]','view',3)
     30plotmodel(md,'data',comp_temp,'title','Modeled temperature [K]','data','observed_temperature','view',3,...
     31        'title','Analytical temperature','view',3,'data',comp_temp-md.observed_temperature,...
     32        'title','Absolute error [K]','view',3,'data',relative,'title','Relative error [%]','view',3,...
     33        'visible','off','figposition','mathieu','FontSize#all',20)
     34set(gcf,'Color','w')
     35printmodel('thermalgeothermalflux','png','margin','on','marginsize',25,'resolution',0.7,'hardcopy','off');
    3136
    3237%Fields and tolerances to track changes
Note: See TracChangeset for help on using the changeset viewer.