Changeset 519


Ignore:
Timestamp:
05/20/09 09:51:48 (16 years ago)
Author:
Mathieu Morlighem
Message:

use nightly run utilities to compute analysis type, sub analysis type package and fields to be checked

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

Legend:

Unmodified
Added
Removed

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

    r487 r519  
    4242                md=parameterize(md,'Square.par');
    4343                md=setelementstype(md,'macayeal','all');
    44                 if strcmpi(package,'cielo_parallel'), md.cluster='wilkes'; end
    4544                if md.numberofgrids==388
    4645                        load Velocities; md.vx=0.5*vx; md.vy=0.5*vy;
    … …  
    4847
    4948                %compute solution
    50                 if (strcmpi(package,'cielo_serial') | strcmpi(package,'cielo_parallel')),
    51                         md=solve(md,'analysis_type',solution,'package','cielo');
    52                 else
    53                         md=solve(md,'analysis_type',solution,'package',package);
    54                 end
     49                [analysis_type sub_analysis_type]=testsgetanalysis(solution);
     50                [md packagefinal]=testsgetpackage(md,package);
     51                md=solve(md,'analysis_type',analysis_type,'sub_analysis_type',sub_analysis_type,'package',package);
    5552
    56                 %compute field to be checked
    57                 if strcmpi(solution,'diagnostic'),
    58                         fields={'vel'};
    59                 elseif strcmpi(solution,'prognostic'),
    60                         fields={'new_thickness'};
    61                 end
     53                %compute fields to be checked
     54                fields=testsgetfields(md.type,solution);
    6255
    6356                %load archive

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

    r487 r519  
    3535                md=parameterize(md,'Square.par');
    3636                md=setelementstype(md,'macayeal','all');
    37                 if strcmpi(package,'cielo_parallel'), md.cluster='wilkes'; end
    3837                if md.numberofgrids==388
    3938                        load Velocities; md.vx=0.5*vx; md.vy=0.5*vy;
    … …  
    4140
    4241                %compute solution
    43                 if (strcmpi(package,'cielo_serial') | strcmpi(package,'cielo_parallel')),
    44                         md=solve(md,'analysis_type',solution,'package','cielo');
    45                 else
    46                         md=solve(md,'analysis_type',solution,'package',package);
    47                 end
     42                [analysis_type sub_analysis_type]=testsgetanalysis(solution);
     43                [md packagefinal]=testsgetpackage(md,package);
     44                md=solve(md,'analysis_type',analysis_type,'sub_analysis_type',sub_analysis_type,'package',package);
    4845
    49                 %field to be saved
    50                 if strcmpi(solution,'diagnostic'),
    51                         fields={'vel'};
    52                 elseif strcmpi(solution,'prognostic'),
    53                         fields={'new_thickness'};
    54                 end
     46                %compute fields to be checked
     47                fields=testsgetfields(md.type,solution);
    5548
    5649                %save new archive

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

    r489 r519  
    4545                md=extrude(md,8,4);
    4646                md=setelementstype(md,'macayeal','all');
    47                 if strcmpi(package,'cielo_parallel'), md.cluster='wilkes'; end
    4847                if md.numberofgrids==832
    4948                        load Velocities; md.vx=0.8*vx; md.vy=0.8*vy; md.vz=0.8*vz; md.temperature=temperature-1; md.pressure=pressure;
    … …  
    5150
    5251                %compute solution
    53                 if (strcmpi(package,'cielo_serial') | strcmpi(package,'cielo_parallel')),
    54                         md=solve(md,'analysis_type',solution,'package','cielo');
    55                 else
    56                         md=solve(md,'analysis_type',solution,'package',package);
    57                 end
     52                [analysis_type sub_analysis_type]=testsgetanalysis(solution);
     53                [md packagefinal]=testsgetpackage(md,package);
     54                md=solve(md,'analysis_type',analysis_type,'sub_analysis_type',sub_analysis_type,'package',package);
    5855
    59                 %compute field to be checked
    60                 if strcmpi(solution,'diagnostic'),
    61                         fields={'vy','vz'};
    62                 elseif strcmpi(solution,'thermalsteady'),
    63                         fields={'temperature','melting'};
    64                 elseif strcmpi(solution,'thermaltransient'),
    65                         fields={'thermaltransient_results(end).temperature','thermaltransient_results(end).melting'};
    66                 elseif strcmpi(solution,'prognostic'),
    67                         fields={'new_thickness'};
    68                 elseif strcmpi(solution,'transient'),
    69                         fields={'transient_results(end).vy','transient_results(end).vz','transient_results(end).temperature','transient_results(end).melting','transient_results(end).thickness'};
    70                 end
     56                %compute fields to be checked
     57                fields=testsgetfields(md.type,solution);
    7158
    7259                %load archive

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

    r489 r519  
    3838                md=extrude(md,8,4);
    3939                md=setelementstype(md,'macayeal','all');
    40                 if strcmpi(package,'cielo_parallel'), md.cluster='wilkes'; end
    4140                if md.numberofgrids==832
    4241                        load Velocities; md.vx=0.8*vx; md.vy=0.8*vy; md.vz=0.8*vz; md.temperature=temperature-1; md.pressure=pressure;
    … …  
    4443
    4544                %compute solution
    46                 if (strcmpi(package,'cielo_serial') | strcmpi(package,'cielo_parallel')),
    47                         md=solve(md,'analysis_type',solution,'package','cielo');
    48                 else
    49                         md=solve(md,'analysis_type',solution,'package',package);
    50                 end
     45                [analysis_type sub_analysis_type]=testsgetanalysis(solution);
     46                [md packagefinal]=testsgetpackage(md,package);
     47                md=solve(md,'analysis_type',analysis_type,'sub_analysis_type',sub_analysis_type,'package',package);
    5148
    52                 %compute field to be checked
    53                 if strcmpi(solution,'diagnostic'),
    54                         fields={'vy','vz'};
    55                 elseif strcmpi(solution,'thermalsteady'),
    56                         fields={'temperature','melting'};
    57                 elseif strcmpi(solution,'thermaltransient'),
    58                         fields={'thermaltransient_results(end).temperature','thermaltransient_results(end).melting'};
    59                 elseif strcmpi(solution,'prognostic'),
    60                         fields={'new_thickness'};
    61                 elseif strcmpi(solution,'transient'),
    62                         fields={'transient_results(end).vy','transient_results(end).vz','transient_results(end).temperature','transient_results(end).melting','transient_results(end).thickness'};
    63                 end
     49                %compute fields to be checked
     50                fields=testsgetfields(md.type,solution);
    6451
    6552                %save new archive

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

    r489 r519  
    4545                md=extrude(md,8,4);
    4646                md=setelementstype(md,'pattyn','Pattyn.exp','fill','macayeal');
    47                 if strcmpi(package,'cielo_parallel'), md.cluster='wilkes'; end
    4847                if md.numberofgrids==832
    4948                        load Velocities; md.vx=0.8*vx; md.vy=0.8*vy; md.vz=0.8*vz; md.temperature=temperature-1; md.pressure=pressure;
    … …  
    5150
    5251                %compute solution
    53                 if (strcmpi(package,'cielo_serial') | strcmpi(package,'cielo_parallel')),
    54                         md=solve(md,'analysis_type',solution,'package','cielo');
    55                 else
    56                         md=solve(md,'analysis_type',solution,'package',package);
    57                 end
     52                [analysis_type sub_analysis_type]=testsgetanalysis(solution);
     53                [md packagefinal]=testsgetpackage(md,package);
     54                md=solve(md,'analysis_type',analysis_type,'sub_analysis_type',sub_analysis_type,'package',package);
    5855
    59                 %compute field to be checked
    60                 if strcmpi(solution,'diagnostic'),
    61                         fields={'vy','vz'};
    62                 elseif strcmpi(solution,'thermalsteady'),
    63                         fields={'temperature','melting'};
    64                 elseif strcmpi(solution,'thermaltransient'),
    65                         fields={'thermaltransient_results(end).temperature','thermaltransient_results(end).melting'};
    66                 elseif strcmpi(solution,'prognostic'),
    67                         fields={'new_thickness'};
    68                 elseif strcmpi(solution,'transient'),
    69                         fields={'transient_results(end).vy','transient_results(end).vz','transient_results(end).temperature','transient_results(end).melting','transient_results(end).thickness'};
    70                 end
     56                %compute fields to be checked
     57                fields=testsgetfields(md.type,solution);
    7158
    7259                %load archive

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

    r489 r519  
    3838                md=extrude(md,8,4);
    3939                md=setelementstype(md,'pattyn','Pattyn.exp','fill','macayeal');
    40                 if strcmpi(package,'cielo_parallel'), md.cluster='wilkes'; end
    4140                if md.numberofgrids==832
    4241                        load Velocities; md.vx=0.8*vx; md.vy=0.8*vy; md.vz=0.8*vz; md.temperature=temperature-1; md.pressure=pressure;
    … …  
    4443
    4544                %compute solution
    46                 if (strcmpi(package,'cielo_serial') | strcmpi(package,'cielo_parallel')),
    47                         md=solve(md,'analysis_type',solution,'package','cielo');
    48                 else
    49                         md=solve(md,'analysis_type',solution,'package',package);
    50                 end
     45                [analysis_type sub_analysis_type]=testsgetanalysis(solution);
     46                [md packagefinal]=testsgetpackage(md,package);
     47                md=solve(md,'analysis_type',analysis_type,'sub_analysis_type',sub_analysis_type,'package',package);
    5148
    52                 %compute field to be checked
    53                 if strcmpi(solution,'diagnostic'),
    54                         fields={'vy','vz'};
    55                 elseif strcmpi(solution,'thermalsteady'),
    56                         fields={'temperature','melting'};
    57                 elseif strcmpi(solution,'thermaltransient'),
    58                         fields={'thermaltransient_results(end).temperature','thermaltransient_results(end).melting'};
    59                 elseif strcmpi(solution,'prognostic'),
    60                         fields={'new_thickness'};
    61                 elseif strcmpi(solution,'transient'),
    62                         fields={'transient_results(end).vy','transient_results(end).vz','transient_results(end).temperature','transient_results(end).melting','transient_results(end).thickness'};
    63                 end
     49                %compute fields to be checked
     50                fields=testsgetfields(md.type,solution);
    6451
    6552                %save new archive

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

    r489 r519  
    4545                md=extrude(md,8,4);
    4646                md=setelementstype(md,'pattyn','all');
    47                 if strcmpi(package,'cielo_parallel'), md.cluster='wilkes'; end
    4847                if md.numberofgrids==832
    4948                        load Velocities; md.vx=0.8*vx; md.vy=0.8*vy; md.vz=0.8*vy; md.temperature=temperature-1; md.pressure=pressure;
    … …  
    5150
    5251                %compute solution
    53                 if (strcmpi(package,'cielo_serial') | strcmpi(package,'cielo_parallel')),
    54                         md=solve(md,'analysis_type',solution,'package','cielo');
    55                 else
    56                         md=solve(md,'analysis_type',solution,'package',package);
    57                 end
     52                [analysis_type sub_analysis_type]=testsgetanalysis(solution);
     53                [md packagefinal]=testsgetpackage(md,package);
     54                md=solve(md,'analysis_type',analysis_type,'sub_analysis_type',sub_analysis_type,'package',package);
    5855
    59                 %compute field to be checked
    60                 if strcmpi(solution,'diagnostic'),
    61                         fields={'vy','vz'};
    62                 elseif strcmpi(solution,'thermalsteady'),
    63                         fields={'temperature','melting'};
    64                 elseif strcmpi(solution,'thermaltransient'),
    65                         fields={'thermaltransient_results(end).temperature','thermaltransient_results(end).melting'};
    66                 elseif strcmpi(solution,'prognostic'),
    67                         fields={'new_thickness'};
    68                 elseif strcmpi(solution,'transient'),
    69                         fields={'transient_results(end).vy','transient_results(end).vz','transient_results(end).temperature','transient_results(end).melting','transient_results(end).thickness'};
    70                 end
     56                %compute fields to be checked
     57                fields=testsgetfields(md.type,solution);
    7158
    7259                %load archive

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

    r489 r519  
    3838                md=extrude(md,8,4);
    3939                md=setelementstype(md,'pattyn','all');
    40                 if strcmpi(package,'cielo_parallel'), md.cluster='wilkes'; end
    4140                if md.numberofgrids==832
    4241                        load Velocities; md.vx=0.8*vx; md.vy=0.8*vy; md.vz=0.8*vy; md.temperature=temperature-1; md.pressure=pressure;
    … …  
    4443
    4544                %compute solution
    46                 if (strcmpi(package,'cielo_serial') | strcmpi(package,'cielo_parallel')),
    47                         md=solve(md,'analysis_type',solution,'package','cielo');
    48                 else
    49                         md=solve(md,'analysis_type',solution,'package',package);
    50                 end
     45                [analysis_type sub_analysis_type]=testsgetanalysis(solution);
     46                [md packagefinal]=testsgetpackage(md,package);
     47                md=solve(md,'analysis_type',analysis_type,'sub_analysis_type',sub_analysis_type,'package',package);
    5148
    52                 %compute field to be checked
    53                 if strcmpi(solution,'diagnostic'),
    54                         fields={'vy','vz'};
    55                 elseif strcmpi(solution,'thermalsteady'),
    56                         fields={'temperature','melting'};
    57                 elseif strcmpi(solution,'thermaltransient'),
    58                         fields={'thermaltransient_results(end).temperature','thermaltransient_results(end).melting'};
    59                 elseif strcmpi(solution,'prognostic'),
    60                         fields={'new_thickness'};
    61                 elseif strcmpi(solution,'transient'),
    62                         fields={'transient_results(end).vy','transient_results(end).vz','transient_results(end).temperature','transient_results(end).melting','transient_results(end).thickness'};
    63                 end
     49                %compute fields to be checked
     50                fields=testsgetfields(md.type,solution);
    6451
    6552                %save new archive

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

    r487 r519  
    4343                md=extrude(md,8,4);
    4444                md=setelementstype(md,'macayeal','all','stokes','all');
    45                 if strcmpi(package,'cielo_parallel'), md.cluster='wilkes'; end
    4645                if md.numberofgrids==832
    4746                        load Velocities; md.vx=0.8*vx; md.vy=0.8*vy; md.vz=0.8*vz; md.temperature=temperature-1; md.pressure=pressure;
    … …  
    4948
    5049                %compute solution
    51                 if (strcmpi(package,'cielo_serial') | strcmpi(package,'cielo_parallel')),
    52                         md=solve(md,'analysis_type',solution,'package','cielo');
    53                 else
    54                         md=solve(md,'analysis_type',solution,'package',package);
    55                 end
     50                [analysis_type sub_analysis_type]=testsgetanalysis(solution);
     51                [md packagefinal]=testsgetpackage(md,package);
     52                md=solve(md,'analysis_type',analysis_type,'sub_analysis_type',sub_analysis_type,'package',package);
    5653
    57                 %compute field to be checked
    58                 if strcmpi(solution,'diagnostic'),
    59                         fields={'vy','vz'};
    60                 elseif strcmpi(solution,'thermalsteady'),
    61                         fields={'temperature','melting'};
    62                 elseif strcmpi(solution,'thermaltransient'),
    63                         fields={'thermaltransient_results(end).temperature','thermaltransient_results(end).melting'};
    64                 elseif strcmpi(solution,'prognostic'),
    65                         fields={'new_thickness'};
    66                 elseif strcmpi(solution,'transient'),
    67                         fields={'transient_results(end).vy','transient_results(end).vz','transient_results(end).temperature','transient_results(end).melting','transient_results(end).thickness'};
    68                 end
     54                %compute fields to be checked
     55                fields=testsgetfields(md.type,solution);
    6956
    7057                %load archive

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

    r487 r519  
    3636                md=extrude(md,8,4);
    3737                md=setelementstype(md,'macayeal','all','stokes','all');
    38                 if strcmpi(package,'cielo_parallel'), md.cluster='wilkes'; end
    3938                if md.numberofgrids==832
    4039                        load Velocities; md.vx=0.8*vx; md.vy=0.8*vy; md.vz=0.8*vz; md.temperature=temperature-1; md.pressure=pressure;
    … …  
    4241
    4342                %compute solution
    44                 if (strcmpi(package,'cielo_serial') | strcmpi(package,'cielo_parallel')),
    45                         md=solve(md,'analysis_type',solution,'package','cielo');
    46                 else
    47                         md=solve(md,'analysis_type',solution,'package',package);
    48                 end
     43                [analysis_type sub_analysis_type]=testsgetanalysis(solution);
     44                [md packagefinal]=testsgetpackage(md,package);
     45                md=solve(md,'analysis_type',analysis_type,'sub_analysis_type',sub_analysis_type,'package',package);
    4946
    50                 %compute field to be checked
    51                 if strcmpi(solution,'diagnostic'),
    52                         fields={'vy','vz'};
    53                 elseif strcmpi(solution,'thermalsteady'),
    54                         fields={'temperature','melting'};
    55                 elseif strcmpi(solution,'thermaltransient'),
    56                         fields={'thermaltransient_results(end).temperature','thermaltransient_results(end).melting'};
    57                 elseif strcmpi(solution,'prognostic'),
    58                         fields={'new_thickness'};
    59                 elseif strcmpi(solution,'transient'),
    60                         fields={'transient_results(end).vy','transient_results(end).vz','transient_results(end).temperature','transient_results(end).melting','transient_results(end).thickness'};
    61                 end
     47                %compute fields to be checked
     48                fields=testsgetfields(md.type,solution);
    6249
    6350                %save new archive

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

    r497 r519  
    4242                md=parameterize(md,'Square.par');
    4343                md=setelementstype(md,'hutter','all');
    44                 if strcmpi(package,'cielo_parallel'), md.cluster='wilkes'; end
    4544                if md.numberofgrids==388
    4645                        load Velocities; md.vx=0.5*vx; md.vy=0.5*vy;
    … …  
    4847
    4948                %compute solution
    50                 if (strcmpi(package,'cielo_serial') | strcmpi(package,'cielo_parallel')),
    51                         md=solve(md,'analysis_type',solution,'package','cielo');
    52                 else
    53                         md=solve(md,'analysis_type',solution,'package',package);
    54                 end
     49                [analysis_type sub_analysis_type]=testsgetanalysis(solution);
     50                [md packagefinal]=testsgetpackage(md,package);
     51                md=solve(md,'analysis_type',analysis_type,'sub_analysis_type',sub_analysis_type,'package',package);
    5552
    56                 %compute field to be checked
    57                 if strcmpi(solution,'diagnostic'),
    58                         fields={'vel'};
    59                 elseif strcmpi(solution,'prognostic'),
    60                         fields={'new_thickness'};
    61                 end
     53                %compute fields to be checked
     54                fields=testsgetfields(md.type,solution);
    6255
    6356                %load archive

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

    r497 r519  
    3535                md=parameterize(md,'Square.par');
    3636                md=setelementstype(md,'hutter','all');
    37                 if strcmpi(package,'cielo_parallel'), md.cluster='wilkes'; end
    3837                if md.numberofgrids==388
    3938                        load Velocities; md.vx=0.5*vx; md.vy=0.5*vy;
    … …  
    4140
    4241                %compute solution
    43                 if (strcmpi(package,'cielo_serial') | strcmpi(package,'cielo_parallel')),
    44                         md=solve(md,'analysis_type',solution,'package','cielo');
    45                 else
    46                         md=solve(md,'analysis_type',solution,'package',package);
    47                 end
     42                [analysis_type sub_analysis_type]=testsgetanalysis(solution);
     43                [md packagefinal]=testsgetpackage(md,package);
     44                md=solve(md,'analysis_type',analysis_type,'sub_analysis_type',sub_analysis_type,'package',package);
    4845
    49                 %field to be saved
    50                 if strcmpi(solution,'diagnostic'),
    51                         fields={'vel'};
    52                 elseif strcmpi(solution,'prognostic'),
    53                         fields={'new_thickness'};
    54                 end
     46                %compute fields to be checked
     47fields=testsgetfields(md.type,solution);
    5548
    5649                %save new archive

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

    r489 r519  
    4545                md=extrude(md,8,4);
    4646                md=setelementstype(md,'hutter','all');
    47                 if strcmpi(package,'cielo_parallel'), md.cluster='wilkes'; end
    4847                if md.numberofgrids==832
    4948                        load Velocities; md.vx=0.8*vx; md.vy=0.8*vy; md.vz=0.8*vz; md.temperature=temperature-1; md.pressure=pressure;
    … …  
    5150
    5251                %compute solution
    53                 if (strcmpi(package,'cielo_serial') | strcmpi(package,'cielo_parallel')),
    54                         md=solve(md,'analysis_type',solution,'package','cielo');
    55                 else
    56                         md=solve(md,'analysis_type',solution,'package',package);
    57                 end
     52                [analysis_type sub_analysis_type]=testsgetanalysis(solution);
     53                [md packagefinal]=testsgetpackage(md,package);
     54                md=solve(md,'analysis_type',analysis_type,'sub_analysis_type',sub_analysis_type,'package',package);
    5855
    59                 %compute field to be checked
    60                 if strcmpi(solution,'diagnostic'),
    61                         fields={'vy','vz'};
    62                 elseif strcmpi(solution,'thermalsteady'),
    63                         fields={'temperature','melting'};
    64                 elseif strcmpi(solution,'thermaltransient'),
    65                         fields={'thermaltransient_results(end).temperature','thermaltransient_results(end).melting'};
    66                 elseif strcmpi(solution,'prognostic'),
    67                         fields={'new_thickness'};
    68                 elseif strcmpi(solution,'transient'),
    69                         fields={'transient_results(end).vy','transient_results(end).vz','transient_results(end).temperature','transient_results(end).melting','transient_results(end).thickness'};
    70                 end
     56                %compute fields to be checked
     57                fields=testsgetfields(md.type,solution);
    7158
    7259                %load archive

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

    r489 r519  
    3838                md=extrude(md,8,4);
    3939                md=setelementstype(md,'hutter','all');
    40                 if strcmpi(package,'cielo_parallel'), md.cluster='wilkes'; end
    4140                if md.numberofgrids==832
    4241                        load Velocities; md.vx=0.8*vx; md.vy=0.8*vy; md.vz=0.8*vz; md.temperature=temperature-1; md.pressure=pressure;
    … …  
    4443
    4544                %compute solution
    46                 if (strcmpi(package,'cielo_serial') | strcmpi(package,'cielo_parallel')),
    47                         md=solve(md,'analysis_type',solution,'package','cielo');
    48                 else
    49                         md=solve(md,'analysis_type',solution,'package',package);
    50                 end
     45                [analysis_type sub_analysis_type]=testsgetanalysis(solution);
     46                [md packagefinal]=testsgetpackage(md,package);
     47                md=solve(md,'analysis_type',analysis_type,'sub_analysis_type',sub_analysis_type,'package',package);
    5148
    52                 %compute field to be checked
    53                 if strcmpi(solution,'diagnostic'),
    54                         fields={'vy','vz'};
    55                 elseif strcmpi(solution,'thermalsteady'),
    56                         fields={'temperature','melting'};
    57                 elseif strcmpi(solution,'thermaltransient'),
    58                         fields={'thermaltransient_results(end).temperature','thermaltransient_results(end).melting'};
    59                 elseif strcmpi(solution,'prognostic'),
    60                         fields={'new_thickness'};
    61                 elseif strcmpi(solution,'transient'),
    62                         fields={'transient_results(end).vy','transient_results(end).vz','transient_results(end).temperature','transient_results(end).melting','transient_results(end).thickness'};
    63                 end
     49                %compute fields to be checked
     50                fields=testsgetfields(md.type,solution);
    6451
    6552                %save new archive

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

    r487 r519  
    4242                md=parameterize(md,'Square.par');
    4343                md=setelementstype(md,'hutter','Hutter.exp','fill','MacAyeal');
    44                 if strcmpi(package,'cielo_parallel'), md.cluster='wilkes'; end
    4544                if md.numberofgrids==388
    4645                        load Velocities; md.vx=0.5*vx; md.vy=0.5*vy;
    … …  
    4847
    4948                %compute solution
    50                 if (strcmpi(package,'cielo_serial') | strcmpi(package,'cielo_parallel')),
    51                         md=solve(md,'analysis_type',solution,'package','cielo');
    52                 else
    53                         md=solve(md,'analysis_type',solution,'package',package);
    54                 end
     49                [analysis_type sub_analysis_type]=testsgetanalysis(solution);
     50                [md packagefinal]=testsgetpackage(md,package);
     51                md=solve(md,'analysis_type',analysis_type,'sub_analysis_type',sub_analysis_type,'package',package);
    5552
    56                 %compute field to be checked
    57                 if strcmpi(solution,'diagnostic'),
    58                         fields={'vel'};
    59                 elseif strcmpi(solution,'prognostic'),
    60                         fields={'new_thickness'};
    61                 end
     53                %compute fields to be checked
     54                fields=testsgetfields(md.type,solution);
    6255
    6356                %load archive

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

    r487 r519  
    3535                md=parameterize(md,'Square.par');
    3636                md=setelementstype(md,'hutter','Hutter.exp','fill','MacAyeal');
    37                 if strcmpi(package,'cielo_parallel'), md.cluster='wilkes'; end
    3837                if md.numberofgrids==388
    3938                        load Velocities; md.vx=0.5*vx; md.vy=0.5*vy;
    … …  
    4140
    4241                %compute solution
    43                 if (strcmpi(package,'cielo_serial') | strcmpi(package,'cielo_parallel')),
    44                         md=solve(md,'analysis_type',solution,'package','cielo');
    45                 else
    46                         md=solve(md,'analysis_type',solution,'package',package);
    47                 end
     42                [analysis_type sub_analysis_type]=testsgetanalysis(solution);
     43                [md packagefinal]=testsgetpackage(md,package);
     44                md=solve(md,'analysis_type',analysis_type,'sub_analysis_type',sub_analysis_type,'package',package);
    4845
    49                 %field to be saved
    50                 if strcmpi(solution,'diagnostic'),
    51                         fields={'vel'};
    52                 elseif strcmpi(solution,'prognostic'),
    53                         fields={'new_thickness'};
    54                 end
     46                %compute fields to be checked
     47fields=testsgetfields(md.type,solution);
    5548
    5649                %save new archive

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

    r487 r519  
    4343                md=extrude(md,8,4);
    4444                md=setelementstype(md,'hutter','Hutter.exp','fill','macayeal');
    45                 if strcmpi(package,'cielo_parallel'), md.cluster='wilkes'; end
    4645                if md.numberofgrids==832
    4746                        load Velocities; md.vx=0.8*vx; md.vy=0.8*vy; md.vz=0.8*vz; md.temperature=temperature-1; md.pressure=pressure;
    … …  
    4948
    5049                %compute solution
    51                 if (strcmpi(package,'cielo_serial') | strcmpi(package,'cielo_parallel')),
    52                         md=solve(md,'analysis_type',solution,'package','cielo');
    53                 else
    54                         md=solve(md,'analysis_type',solution,'package',package);
    55                 end
     50                [analysis_type sub_analysis_type]=testsgetanalysis(solution);
     51                [md packagefinal]=testsgetpackage(md,package);
     52                md=solve(md,'analysis_type',analysis_type,'sub_analysis_type',sub_analysis_type,'package',package);
    5653
    57                 %compute field to be checked
    58                 if strcmpi(solution,'diagnostic'),
    59                         fields={'vy','vz'};
    60                 elseif strcmpi(solution,'thermalsteady'),
    61                         fields={'temperature','melting'};
    62                 elseif strcmpi(solution,'thermaltransient'),
    63                         fields={'thermaltransient_results(end).temperature','thermaltransient_results(end).melting'};
    64                 elseif strcmpi(solution,'prognostic'),
    65                         fields={'new_thickness'};
    66                 elseif strcmpi(solution,'transient'),
    67                         fields={'transient_results(end).vy','transient_results(end).vz','transient_results(end).temperature','transient_results(end).melting','transient_results(end).thickness'};
    68                 end
     54                %compute fields to be checked
     55                fields=testsgetfields(md.type,solution);
    6956
    7057                %load archive

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

    r487 r519  
    3636                md=extrude(md,8,4);
    3737                md=setelementstype(md,'hutter','Hutter.exp','fill','macayeal');
    38                 if strcmpi(package,'cielo_parallel'), md.cluster='wilkes'; end
    3938                if md.numberofgrids==832
    4039                        load Velocities; md.vx=0.8*vx; md.vy=0.8*vy; md.vz=0.8*vz; md.temperature=temperature-1; md.pressure=pressure;
    … …  
    4241
    4342                %compute solution
    44                 if (strcmpi(package,'cielo_serial') | strcmpi(package,'cielo_parallel')),
    45                         md=solve(md,'analysis_type',solution,'package','cielo');
    46                 else
    47                         md=solve(md,'analysis_type',solution,'package',package);
    48                 end
     43                [analysis_type sub_analysis_type]=testsgetanalysis(solution);
     44                [md packagefinal]=testsgetpackage(md,package);
     45                md=solve(md,'analysis_type',analysis_type,'sub_analysis_type',sub_analysis_type,'package',package);
    4946
    50                 %compute field to be checked
    51                 if strcmpi(solution,'diagnostic'),
    52                         fields={'vy','vz'};
    53                 elseif strcmpi(solution,'thermalsteady'),
    54                         fields={'temperature','melting'};
    55                 elseif strcmpi(solution,'thermaltransient'),
    56                         fields={'thermaltransient_results(end).temperature','thermaltransient_results(end).melting'};
    57                 elseif strcmpi(solution,'prognostic'),
    58                         fields={'new_thickness'};
    59                 elseif strcmpi(solution,'transient'),
    60                         fields={'transient_results(end).vy','transient_results(end).vz','transient_results(end).temperature','transient_results(end).melting','transient_results(end).thickness'};
    61                 end
     47                %compute fields to be checked
     48                fields=testsgetfields(md.type,solution);
    6249
    6350                %save new archive

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

    r487 r519  
    4242                md=parameterize(md,'Square.par');
    4343                md=setelementstype(md,'macayeal','all');
    44                 if strcmpi(package,'cielo_parallel'), md.cluster='wilkes'; end
    4544                if md.numberofgrids==388
    4645                        load Velocities; md.vx=0.5*vx; md.vy=0.5*vy;
    … …  
    4847
    4948                %compute solution
    50                 if (strcmpi(package,'cielo_serial') | strcmpi(package,'cielo_parallel')),
    51                         md=solve(md,'analysis_type',solution,'package','cielo');
    52                 else
    53                         md=solve(md,'analysis_type',solution,'package',package);
    54                 end
     49                [analysis_type sub_analysis_type]=testsgetanalysis(solution);
     50                [md packagefinal]=testsgetpackage(md,package);
     51                md=solve(md,'analysis_type',analysis_type,'sub_analysis_type',sub_analysis_type,'package',package);
    5552
    56                 %compute field to be checked
    57                 if strcmpi(solution,'diagnostic'),
    58                         fields={'vel'};
    59                 elseif strcmpi(solution,'prognostic'),
    60                         fields={'new_thickness'};
    61                 end
     53                %compute fields to be checked
     54                fields=testsgetfields(md.type,solution);
    6255
    6356                %load archive

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

    r487 r519  
    3535                md=parameterize(md,'Square.par');
    3636                md=setelementstype(md,'macayeal','all');
    37                 if strcmpi(package,'cielo_parallel'), md.cluster='wilkes'; end
    3837                if md.numberofgrids==388
    3938                        load Velocities; md.vx=0.5*vx; md.vy=0.5*vy;
    … …  
    4140
    4241                %compute solution
    43                 if (strcmpi(package,'cielo_serial') | strcmpi(package,'cielo_parallel')),
    44                         md=solve(md,'analysis_type',solution,'package','cielo');
    45                 else
    46                         md=solve(md,'analysis_type',solution,'package',package);
    47                 end
     42                [analysis_type sub_analysis_type]=testsgetanalysis(solution);
     43                [md packagefinal]=testsgetpackage(md,package);
     44                md=solve(md,'analysis_type',analysis_type,'sub_analysis_type',sub_analysis_type,'package',package);
    4845
    49                 %field to be saved
    50                 if strcmpi(solution,'diagnostic'),
    51                         fields={'vel'};
    52                 elseif strcmpi(solution,'prognostic'),
    53                         fields={'new_thickness'};
    54                 end
     46                %compute fields to be checked
     47fields=testsgetfields(md.type,solution);
    5548
    5649                %save new archive

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

    r489 r519  
    4545                md=extrude(md,8,4);
    4646                md=setelementstype(md,'macayeal','all');
    47                 if strcmpi(package,'cielo_parallel'), md.cluster='wilkes'; end
    4847                if md.numberofgrids==832
    4948                        load Velocities; md.vx=0.8*vx; md.vy=0.8*vy; md.vz=0.8*vz; md.temperature=temperature-1; md.pressure=pressure;
    … …  
    5150
    5251                %compute solution
    53                 if (strcmpi(package,'cielo_serial') | strcmpi(package,'cielo_parallel')),
    54                         md=solve(md,'analysis_type',solution,'package','cielo');
    55                 else
    56                         md=solve(md,'analysis_type',solution,'package',package);
    57                 end
     52                [analysis_type sub_analysis_type]=testsgetanalysis(solution);
     53                [md packagefinal]=testsgetpackage(md,package);
     54                md=solve(md,'analysis_type',analysis_type,'sub_analysis_type',sub_analysis_type,'package',package);
    5855
    59                 %compute field to be checked
    60                 if strcmpi(solution,'diagnostic'),
    61                         fields={'vy','vz'};
    62                 elseif strcmpi(solution,'thermalsteady'),
    63                         fields={'temperature','melting'};
    64                 elseif strcmpi(solution,'thermaltransient'),
    65                         fields={'thermaltransient_results(end).temperature','thermaltransient_results(end).melting'};
    66                 elseif strcmpi(solution,'prognostic'),
    67                         fields={'new_thickness'};
    68                 elseif strcmpi(solution,'transient'),
    69                         fields={'transient_results(end).vy','transient_results(end).vz','transient_results(end).temperature','transient_results(end).melting','transient_results(end).thickness'};
    70                 end
     56                %compute fields to be checked
     57                fields=testsgetfields(md.type,solution);
    7158
    7259                %load archive

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

    r489 r519  
    3838                md=extrude(md,8,4);
    3939                md=setelementstype(md,'macayeal','all');
    40                 if strcmpi(package,'cielo_parallel'), md.cluster='wilkes'; end
    4140                if md.numberofgrids==832
    4241                        load Velocities; md.vx=0.8*vx; md.vy=0.8*vy; md.vz=0.8*vz; md.temperature=temperature-1; md.pressure=pressure;
    … …  
    4443
    4544                %compute solution
    46                 if (strcmpi(package,'cielo_serial') | strcmpi(package,'cielo_parallel')),
    47                         md=solve(md,'analysis_type',solution,'package','cielo');
    48                 else
    49                         md=solve(md,'analysis_type',solution,'package',package);
    50                 end
     45                [analysis_type sub_analysis_type]=testsgetanalysis(solution);
     46                [md packagefinal]=testsgetpackage(md,package);
     47                md=solve(md,'analysis_type',analysis_type,'sub_analysis_type',sub_analysis_type,'package',package);
    5148
    52                 %compute field to be checked
    53                 if strcmpi(solution,'diagnostic'),
    54                         fields={'vy','vz'};
    55                 elseif strcmpi(solution,'thermalsteady'),
    56                         fields={'temperature','melting'};
    57                 elseif strcmpi(solution,'thermaltransient'),
    58                         fields={'thermaltransient_results(end).temperature','thermaltransient_results(end).melting'};
    59                 elseif strcmpi(solution,'prognostic'),
    60                         fields={'new_thickness'};
    61                 elseif strcmpi(solution,'transient'),
    62                         fields={'transient_results(end).vy','transient_results(end).vz','transient_results(end).temperature','transient_results(end).melting','transient_results(end).thickness'};
    63                 end
     49                %compute fields to be checked
     50                fields=testsgetfields(md.type,solution);
    6451
    6552                %save new archive

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

    r489 r519  
    4545                md=extrude(md,8,4);
    4646                md=setelementstype(md,'pattyn','all');
    47                 if strcmpi(package,'cielo_parallel'), md.cluster='wilkes'; end
    4847                if md.numberofgrids==832
    4948                        load Velocities; md.vx=0.8*vx; md.vy=0.8*vy; md.vz=0.8*vz; md.temperature=temperature-1; md.pressure=pressure;
    … …  
    5150
    5251                %compute solution
    53                 if (strcmpi(package,'cielo_serial') | strcmpi(package,'cielo_parallel')),
    54                         md=solve(md,'analysis_type',solution,'package','cielo');
    55                 else
    56                         md=solve(md,'analysis_type',solution,'package',package);
    57                 end
     52                [analysis_type sub_analysis_type]=testsgetanalysis(solution);
     53                [md packagefinal]=testsgetpackage(md,package);
     54                md=solve(md,'analysis_type',analysis_type,'sub_analysis_type',sub_analysis_type,'package',package);
    5855
    59                 %compute field to be checked
    60                 if strcmpi(solution,'diagnostic'),
    61                         fields={'vy','vz'};
    62                 elseif strcmpi(solution,'thermalsteady'),
    63                         fields={'temperature','melting'};
    64                 elseif strcmpi(solution,'thermaltransient'),
    65                         fields={'thermaltransient_results(end).temperature','thermaltransient_results(end).melting'};
    66                 elseif strcmpi(solution,'prognostic'),
    67                         fields={'new_thickness'};
    68                 elseif strcmpi(solution,'transient'),
    69                         fields={'transient_results(end).vy','transient_results(end).vz','transient_results(end).temperature','transient_results(end).melting','transient_results(end).thickness'};
    70                 end
     56                %compute fields to be checked
     57                fields=testsgetfields(md.type,solution);
    7158
    7259                %load archive

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

    r489 r519  
    3838                md=extrude(md,8,4);
    3939                md=setelementstype(md,'pattyn','all');
    40                 if strcmpi(package,'cielo_parallel'), md.cluster='wilkes'; end
    4140                if md.numberofgrids==832
    4241                        load Velocities; md.vx=0.8*vx; md.vy=0.8*vy; md.vz=0.8*vz; md.temperature=temperature-1; md.pressure=pressure;
    … …  
    4443
    4544                %compute solution
    46                 if (strcmpi(package,'cielo_serial') | strcmpi(package,'cielo_parallel')),
    47                         md=solve(md,'analysis_type',solution,'package','cielo');
    48                 else
    49                         md=solve(md,'analysis_type',solution,'package',package);
    50                 end
     45                [analysis_type sub_analysis_type]=testsgetanalysis(solution);
     46                [md packagefinal]=testsgetpackage(md,package);
     47                md=solve(md,'analysis_type',analysis_type,'sub_analysis_type',sub_analysis_type,'package',package);
    5148
    52                 %compute field to be checked
    53                 if strcmpi(solution,'diagnostic'),
    54                         fields={'vy','vz'};
    55                 elseif strcmpi(solution,'thermalsteady'),
    56                         fields={'temperature','melting'};
    57                 elseif strcmpi(solution,'thermaltransient'),
    58                         fields={'thermaltransient_results(end).temperature','thermaltransient_results(end).melting'};
    59                 elseif strcmpi(solution,'prognostic'),
    60                         fields={'new_thickness'};
    61                 elseif strcmpi(solution,'transient'),
    62                         fields={'transient_results(end).vy','transient_results(end).vz','transient_results(end).temperature','transient_results(end).melting','transient_results(end).thickness'};
    63                 end
     49                %compute fields to be checked
     50                fields=testsgetfields(md.type,solution);
    6451
    6552                %save new archive

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

    r487 r519  
    4343                md=extrude(md,8,4);
    4444                md=setelementstype(md,'pattyn','all','stokes','all');
    45                 if strcmpi(package,'cielo_parallel'), md.cluster='wilkes'; end
    4645                if md.numberofgrids==832
    4746                        load Velocities; md.vx=0.8*vx; md.vy=0.8*vy; md.vz=0.8*vz; md.temperature=temperature-1; md.pressure=pressure;
    … …  
    4948
    5049                %compute solution
    51                 if (strcmpi(package,'cielo_serial') | strcmpi(package,'cielo_parallel')),
    52                         md=solve(md,'analysis_type',solution,'package','cielo');
    53                 else
    54                         md=solve(md,'analysis_type',solution,'package',package);
    55                 end
     50                [analysis_type sub_analysis_type]=testsgetanalysis(solution);
     51                [md packagefinal]=testsgetpackage(md,package);
     52                md=solve(md,'analysis_type',analysis_type,'sub_analysis_type',sub_analysis_type,'package',package);
    5653
    57                 %compute field to be checked
    58                 if strcmpi(solution,'diagnostic'),
    59                         fields={'vy','vz'};
    60                 elseif strcmpi(solution,'thermalsteady'),
    61                         fields={'temperature','melting'};
    62                 elseif strcmpi(solution,'thermaltransient'),
    63                         fields={'thermaltransient_results(end).temperature','thermaltransient_results(end).melting'};
    64                 elseif strcmpi(solution,'prognostic'),
    65                         fields={'new_thickness'};
    66                 elseif strcmpi(solution,'transient'),
    67                         fields={'transient_results(end).vy','transient_results(end).vz','transient_results(end).temperature','transient_results(end).melting','transient_results(end).thickness'};
    68                 end
     54                %compute fields to be checked
     55                fields=testsgetfields(md.type,solution);
    6956
    7057                %load archive

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

    r487 r519  
    3636                md=extrude(md,8,4);
    3737                md=setelementstype(md,'pattyn','all','stokes','all');
    38                 if strcmpi(package,'cielo_parallel'), md.cluster='wilkes'; end
    3938                if md.numberofgrids==832
    4039                        load Velocities; md.vx=0.8*vx; md.vy=0.8*vy; md.vz=0.8*vz; md.temperature=temperature-1; md.pressure=pressure;
    … …  
    4241
    4342                %compute solution
    44                 if (strcmpi(package,'cielo_serial') | strcmpi(package,'cielo_parallel')),
    45                         md=solve(md,'analysis_type',solution,'package','cielo');
    46                 else
    47                         md=solve(md,'analysis_type',solution,'package',package);
    48                 end
     43                [analysis_type sub_analysis_type]=testsgetanalysis(solution);
     44                [md packagefinal]=testsgetpackage(md,package);
     45                md=solve(md,'analysis_type',analysis_type,'sub_analysis_type',sub_analysis_type,'package',package);
    4946
    50                 %compute field to be checked
    51                 if strcmpi(solution,'diagnostic'),
    52                         fields={'vy','vz'};
    53                 elseif strcmpi(solution,'thermalsteady'),
    54                         fields={'temperature','melting'};
    55                 elseif strcmpi(solution,'thermaltransient'),
    56                         fields={'thermaltransient_results(end).temperature','thermaltransient_results(end).melting'};
    57                 elseif strcmpi(solution,'prognostic'),
    58                         fields={'new_thickness'};
    59                 elseif strcmpi(solution,'transient'),
    60                         fields={'transient_results(end).vy','transient_results(end).vz','transient_results(end).temperature','transient_results(end).melting','transient_results(end).thickness'};
    61                 end
     47                %compute fields to be checked
     48                fields=testsgetfields(md.type,solution);
    6249
    6350                %save new archive

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

    r487 r519  
    4242                md=parameterize(md,'Square.par');
    4343                md=setelementstype(md,'macayeal','all');
    44                 if strcmpi(package,'cielo_parallel'), md.cluster='wilkes'; end
    4544                if md.numberofgrids==388
    4645                        load Velocities; md.vx=0.5*vx; md.vy=0.5*vy;
    … …  
    4847
    4948                %compute solution
    50                 if (strcmpi(package,'cielo_serial') | strcmpi(package,'cielo_parallel')),
    51                         md=solve(md,'analysis_type',solution,'package','cielo');
    52                 else
    53                         md=solve(md,'analysis_type',solution,'package',package);
    54                 end
     49                [analysis_type sub_analysis_type]=testsgetanalysis(solution);
     50                [md packagefinal]=testsgetpackage(md,package);
     51                md=solve(md,'analysis_type',analysis_type,'sub_analysis_type',sub_analysis_type,'package',package);
    5552
    56                 %compute field to be checked
    57                 if strcmpi(solution,'diagnostic'),
    58                         fields={'vel'};
    59                 elseif strcmpi(solution,'prognostic'),
    60                         fields={'new_thickness'};
    61                 end
     53                %compute fields to be checked
     54                fields=testsgetfields(md.type,solution);
    6255
    6356                %load archive

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

    r487 r519  
    3535                md=parameterize(md,'Square.par');
    3636                md=setelementstype(md,'macayeal','all');
    37                 if strcmpi(package,'cielo_parallel'), md.cluster='wilkes'; end
    3837                if md.numberofgrids==388
    3938                        load Velocities; md.vx=0.5*vx; md.vy=0.5*vy;
    … …  
    4140
    4241                %compute solution
    43                 if (strcmpi(package,'cielo_serial') | strcmpi(package,'cielo_parallel')),
    44                         md=solve(md,'analysis_type',solution,'package','cielo');
    45                 else
    46                         md=solve(md,'analysis_type',solution,'package',package);
    47                 end
     42                [analysis_type sub_analysis_type]=testsgetanalysis(solution);
     43                [md packagefinal]=testsgetpackage(md,package);
     44                md=solve(md,'analysis_type',analysis_type,'sub_analysis_type',sub_analysis_type,'package',package);
    4845
    49                 %field to be saved
    50                 if strcmpi(solution,'diagnostic'),
    51                         fields={'vel'};
    52                 elseif strcmpi(solution,'prognostic'),
    53                         fields={'new_thickness'};
    54                 end
     46                %compute fields to be checked
     47fields=testsgetfields(md.type,solution);
    5548
    5649                %save new archive

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

    r489 r519  
    4545                md=extrude(md,8,4);
    4646                md=setelementstype(md,'macayeal','all');
    47                 if strcmpi(package,'cielo_parallel'), md.cluster='wilkes'; end
    4847                if md.numberofgrids==832
    4948                        load Velocities; md.vx=0.8*vx; md.vy=0.8*vy; md.vz=0.8*vz; md.temperature=temperature-1; md.pressure=pressure;
    … …  
    5150
    5251                %compute solution
    53                 if (strcmpi(package,'cielo_serial') | strcmpi(package,'cielo_parallel')),
    54                         md=solve(md,'analysis_type',solution,'package','cielo');
    55                 else
    56                         md=solve(md,'analysis_type',solution,'package',package);
    57                 end
     52                [analysis_type sub_analysis_type]=testsgetanalysis(solution);
     53                [md packagefinal]=testsgetpackage(md,package);
     54                md=solve(md,'analysis_type',analysis_type,'sub_analysis_type',sub_analysis_type,'package',package);
    5855
    59                 %compute field to be checked
    60                 if strcmpi(solution,'diagnostic'),
    61                         fields={'vy','vz'};
    62                 elseif strcmpi(solution,'thermalsteady'),
    63                         fields={'temperature','melting'};
    64                 elseif strcmpi(solution,'thermaltransient'),
    65                         fields={'thermaltransient_results(end).temperature','thermaltransient_results(end).melting'};
    66                 elseif strcmpi(solution,'prognostic'),
    67                         fields={'new_thickness'};
    68                 elseif strcmpi(solution,'transient'),
    69                         fields={'transient_results(end).vy','transient_results(end).vz','transient_results(end).temperature','transient_results(end).melting','transient_results(end).thickness'};
    70                 end
     56                %compute fields to be checked
     57                fields=testsgetfields(md.type,solution);
    7158
    7259                %load archive

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

    r489 r519  
    3838                md=extrude(md,8,4);
    3939                md=setelementstype(md,'macayeal','all');
    40                 if strcmpi(package,'cielo_parallel'), md.cluster='wilkes'; end
    4140                if md.numberofgrids==832
    4241                        load Velocities; md.vx=0.8*vx; md.vy=0.8*vy; md.vz=0.8*vz; md.temperature=temperature-1; md.pressure=pressure;
    … …  
    4443
    4544                %compute solution
    46                 if (strcmpi(package,'cielo_serial') | strcmpi(package,'cielo_parallel')),
    47                         md=solve(md,'analysis_type',solution,'package','cielo');
    48                 else
    49                         md=solve(md,'analysis_type',solution,'package',package);
    50                 end
     45                [analysis_type sub_analysis_type]=testsgetanalysis(solution);
     46                [md packagefinal]=testsgetpackage(md,package);
     47                md=solve(md,'analysis_type',analysis_type,'sub_analysis_type',sub_analysis_type,'package',package);
    5148
    52                 %compute field to be checked
    53                 if strcmpi(solution,'diagnostic'),
    54                         fields={'vy','vz'};
    55                 elseif strcmpi(solution,'thermalsteady'),
    56                         fields={'temperature','melting'};
    57                 elseif strcmpi(solution,'thermaltransient'),
    58                         fields={'thermaltransient_results(end).temperature','thermaltransient_results(end).melting'};
    59                 elseif strcmpi(solution,'prognostic'),
    60                         fields={'new_thickness'};
    61                 elseif strcmpi(solution,'transient'),
    62                         fields={'transient_results(end).vy','transient_results(end).vz','transient_results(end).temperature','transient_results(end).melting','transient_results(end).thickness'};
    63                 end
     49                %compute fields to be checked
     50                fields=testsgetfields(md.type,solution);
    6451
    6552                %save new archive

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

    r489 r519  
    4545                md=extrude(md,8,4);
    4646                md=setelementstype(md,'pattyn','Pattyn.exp','fill','macayeal');
    47                 if strcmpi(package,'cielo_parallel'), md.cluster='wilkes'; end
    4847                if md.numberofgrids==832
    4948                        load Velocities; md.vx=0.8*vx; md.vy=0.8*vy; md.vz=0.8*vz; md.temperature=temperature-1; md.pressure=pressure;
    … …  
    5150
    5251                %compute solution
    53                 if (strcmpi(package,'cielo_serial') | strcmpi(package,'cielo_parallel')),
    54                         md=solve(md,'analysis_type',solution,'package','cielo');
    55                 else
    56                         md=solve(md,'analysis_type',solution,'package',package);
    57                 end
     52                [analysis_type sub_analysis_type]=testsgetanalysis(solution);
     53                [md packagefinal]=testsgetpackage(md,package);
     54                md=solve(md,'analysis_type',analysis_type,'sub_analysis_type',sub_analysis_type,'package',package);
    5855
    59                 %compute field to be checked
    60                 if strcmpi(solution,'diagnostic'),
    61                         fields={'vy','vz'};
    62                 elseif strcmpi(solution,'thermalsteady'),
    63                         fields={'temperature','melting'};
    64                 elseif strcmpi(solution,'thermaltransient'),
    65                         fields={'thermaltransient_results(end).temperature','thermaltransient_results(end).melting'};
    66                 elseif strcmpi(solution,'prognostic'),
    67                         fields={'new_thickness'};
    68                 elseif strcmpi(solution,'transient'),
    69                         fields={'transient_results(end).vy','transient_results(end).vz','transient_results(end).temperature','transient_results(end).melting','transient_results(end).thickness'};
    70                 end
     56                %compute fields to be checked
     57                fields=testsgetfields(md.type,solution);
    7158
    7259                %load archive

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

    r489 r519  
    3838                md=extrude(md,8,4);
    3939                md=setelementstype(md,'pattyn','Pattyn.exp','fill','macayeal');
    40                 if strcmpi(package,'cielo_parallel'), md.cluster='wilkes'; end
    4140                if md.numberofgrids==832
    4241                        load Velocities; md.vx=0.8*vx; md.vy=0.8*vy; md.vz=0.8*vz; md.temperature=temperature-1; md.pressure=pressure;
    … …  
    4443
    4544                %compute solution
    46                 if (strcmpi(package,'cielo_serial') | strcmpi(package,'cielo_parallel')),
    47                         md=solve(md,'analysis_type',solution,'package','cielo');
    48                 else
    49                         md=solve(md,'analysis_type',solution,'package',package);
    50                 end
     45                [analysis_type sub_analysis_type]=testsgetanalysis(solution);
     46                [md packagefinal]=testsgetpackage(md,package);
     47                md=solve(md,'analysis_type',analysis_type,'sub_analysis_type',sub_analysis_type,'package',package);
    5148
    52                 %compute field to be checked
    53                 if strcmpi(solution,'diagnostic'),
    54                         fields={'vy','vz'};
    55                 elseif strcmpi(solution,'thermalsteady'),
    56                         fields={'temperature','melting'};
    57                 elseif strcmpi(solution,'thermaltransient'),
    58                         fields={'thermaltransient_results(end).temperature','thermaltransient_results(end).melting'};
    59                 elseif strcmpi(solution,'prognostic'),
    60                         fields={'new_thickness'};
    61                 elseif strcmpi(solution,'transient'),
    62                         fields={'transient_results(end).vy','transient_results(end).vz','transient_results(end).temperature','transient_results(end).melting','transient_results(end).thickness'};
    63                 end
     49                %compute fields to be checked
     50                fields=testsgetfields(md.type,solution);
    6451
    6552                %save new archive

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

    r487 r519  
    4343                md=extrude(md,8,4);
    4444                md=setelementstype(md,'pattyn','Pattyn.exp','fill','macayeal','stokes','Stokes.exp');
    45                 if strcmpi(package,'cielo_parallel'), md.cluster='wilkes'; end
    4645                if md.numberofgrids==832
    4746                        load Velocities; md.vx=0.8*vx; md.vy=0.8*vy; md.vz=0.8*vz; md.temperature=temperature-1; md.pressure=pressure;
    … …  
    4948
    5049                %compute solution
    51                 if (strcmpi(package,'cielo_serial') | strcmpi(package,'cielo_parallel')),
    52                         md=solve(md,'analysis_type',solution,'package','cielo');
    53                 else
    54                         md=solve(md,'analysis_type',solution,'package',package);
    55                 end
     50                [analysis_type sub_analysis_type]=testsgetanalysis(solution);
     51                [md packagefinal]=testsgetpackage(md,package);
     52                md=solve(md,'analysis_type',analysis_type,'sub_analysis_type',sub_analysis_type,'package',package);
    5653
    57                 %compute field to be checked
    58                 if strcmpi(solution,'diagnostic'),
    59                         fields={'vy','vz'};
    60                 elseif strcmpi(solution,'thermalsteady'),
    61                         fields={'temperature','melting'};
    62                 elseif strcmpi(solution,'thermaltransient'),
    63                         fields={'thermaltransient_results(end).temperature','thermaltransient_results(end).melting'};
    64                 elseif strcmpi(solution,'prognostic'),
    65                         fields={'new_thickness'};
    66                 elseif strcmpi(solution,'transient'),
    67                         fields={'transient_results(end).vy','transient_results(end).vz','transient_results(end).temperature','transient_results(end).melting','transient_results(end).thickness'};
    68                 end
     54                %compute fields to be checked
     55                fields=testsgetfields(md.type,solution);
    6956
    7057                %load archive

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

    r487 r519  
    3636                md=extrude(md,8,4);
    3737                md=setelementstype(md,'pattyn','Pattyn.exp','fill','macayeal','stokes','Stokes.exp');
    38                 if strcmpi(package,'cielo_parallel'), md.cluster='wilkes'; end
    3938                if md.numberofgrids==832
    4039                        load Velocities; md.vx=0.8*vx; md.vy=0.8*vy; md.vz=0.8*vz; md.temperature=temperature-1; md.pressure=pressure;
    … …  
    4241
    4342                %compute solution
    44                 if (strcmpi(package,'cielo_serial') | strcmpi(package,'cielo_parallel')),
    45                         md=solve(md,'analysis_type',solution,'package','cielo');
    46                 else
    47                         md=solve(md,'analysis_type',solution,'package',package);
    48                 end
     43                [analysis_type sub_analysis_type]=testsgetanalysis(solution);
     44                [md packagefinal]=testsgetpackage(md,package);
     45                md=solve(md,'analysis_type',analysis_type,'sub_analysis_type',sub_analysis_type,'package',package);
    4946
    50                 %compute field to be checked
    51                 if strcmpi(solution,'diagnostic'),
    52                         fields={'vy','vz'};
    53                 elseif strcmpi(solution,'thermalsteady'),
    54                         fields={'temperature','melting'};
    55                 elseif strcmpi(solution,'thermaltransient'),
    56                         fields={'thermaltransient_results(end).temperature','thermaltransient_results(end).melting'};
    57                 elseif strcmpi(solution,'prognostic'),
    58                         fields={'new_thickness'};
    59                 elseif strcmpi(solution,'transient'),
    60                         fields={'transient_results(end).vy','transient_results(end).vz','transient_results(end).temperature','transient_results(end).melting','transient_results(end).thickness'};
    61                 end
     47                %compute fields to be checked
     48                fields=testsgetfields(md.type,solution);
    6249
    6350                %save new archive

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

    r489 r519  
    4545                md=extrude(md,8,4);
    4646                md=setelementstype(md,'pattyn','all');
    47                 if strcmpi(package,'cielo_parallel'), md.cluster='wilkes'; end
    4847                if md.numberofgrids==832
    4948                        load Velocities; md.vx=0.8*vx; md.vy=0.8*vy; md.vz=0.8*vz; md.temperature=temperature-1; md.pressure=pressure;
    … …  
    5150
    5251                %compute solution
    53                 if (strcmpi(package,'cielo_serial') | strcmpi(package,'cielo_parallel')),
    54                         md=solve(md,'analysis_type',solution,'package','cielo');
    55                 else
    56                         md=solve(md,'analysis_type',solution,'package',package);
    57                 end
     52                [analysis_type sub_analysis_type]=testsgetanalysis(solution);
     53                [md packagefinal]=testsgetpackage(md,package);
     54                md=solve(md,'analysis_type',analysis_type,'sub_analysis_type',sub_analysis_type,'package',package);
    5855
    59                 %compute field to be checked
    60                 if strcmpi(solution,'diagnostic'),
    61                         fields={'vy','vz'};
    62                 elseif strcmpi(solution,'thermalsteady'),
    63                         fields={'temperature','melting'};
    64                 elseif strcmpi(solution,'thermaltransient'),
    65                         fields={'thermaltransient_results(end).temperature','thermaltransient_results(end).melting'};
    66                 elseif strcmpi(solution,'prognostic'),
    67                         fields={'new_thickness'};
    68                 elseif strcmpi(solution,'transient'),
    69                         fields={'transient_results(end).vy','transient_results(end).vz','transient_results(end).temperature','transient_results(end).melting','transient_results(end).thickness'};
    70                 end
     56                %compute fields to be checked
     57                fields=testsgetfields(md.type,solution);
    7158
    7259                %load archive

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

    r489 r519  
    3838                md=extrude(md,8,4);
    3939                md=setelementstype(md,'pattyn','all');
    40                 if strcmpi(package,'cielo_parallel'), md.cluster='wilkes'; end
    4140                if md.numberofgrids==832
    4241                        load Velocities; md.vx=0.8*vx; md.vy=0.8*vy; md.vz=0.8*vz; md.temperature=temperature-1; md.pressure=pressure;
    … …  
    4443
    4544                %compute solution
    46                 if (strcmpi(package,'cielo_serial') | strcmpi(package,'cielo_parallel')),
    47                         md=solve(md,'analysis_type',solution,'package','cielo');
    48                 else
    49                         md=solve(md,'analysis_type',solution,'package',package);
    50                 end
     45                [analysis_type sub_analysis_type]=testsgetanalysis(solution);
     46                [md packagefinal]=testsgetpackage(md,package);
     47                md=solve(md,'analysis_type',analysis_type,'sub_analysis_type',sub_analysis_type,'package',package);
    5148
    52                 %compute field to be checked
    53                 if strcmpi(solution,'diagnostic'),
    54                         fields={'vy','vz'};
    55                 elseif strcmpi(solution,'thermalsteady'),
    56                         fields={'temperature','melting'};
    57                 elseif strcmpi(solution,'thermaltransient'),
    58                         fields={'thermaltransient_results(end).temperature','thermaltransient_results(end).melting'};
    59                 elseif strcmpi(solution,'prognostic'),
    60                         fields={'new_thickness'};
    61                 elseif strcmpi(solution,'transient'),
    62                         fields={'transient_results(end).vy','transient_results(end).vz','transient_results(end).temperature','transient_results(end).melting','transient_results(end).thickness'};
    63                 end
     49                %compute fields to be checked
     50                fields=testsgetfields(md.type,solution);
    6451
    6552                %save new archive

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

    r487 r519  
    4343                md=extrude(md,8,4);
    4444                md=setelementstype(md,'pattyn','all','stokes','all');
    45                 if strcmpi(package,'cielo_parallel'), md.cluster='wilkes'; end
    4645                if md.numberofgrids==832
    4746                        load Velocities; md.vx=0.8*vx; md.vy=0.8*vy; md.vz=0.8*vz; md.temperature=temperature-1; md.pressure=pressure;
    … …  
    4948
    5049                %compute solution
    51                 if (strcmpi(package,'cielo_serial') | strcmpi(package,'cielo_parallel')),
    52                         md=solve(md,'analysis_type',solution,'package','cielo');
    53                 else
    54                         md=solve(md,'analysis_type',solution,'package',package);
    55                 end
     50                [analysis_type sub_analysis_type]=testsgetanalysis(solution);
     51                [md packagefinal]=testsgetpackage(md,package);
     52                md=solve(md,'analysis_type',analysis_type,'sub_analysis_type',sub_analysis_type,'package',package);
    5653
    57                 %compute field to be checked
    58                 if strcmpi(solution,'diagnostic'),
    59                         fields={'vy','vz'};
    60                 elseif strcmpi(solution,'thermalsteady'),
    61                         fields={'temperature','melting'};
    62                 elseif strcmpi(solution,'thermaltransient'),
    63                         fields={'thermaltransient_results(end).temperature','thermaltransient_results(end).melting'};
    64                 elseif strcmpi(solution,'prognostic'),
    65                         fields={'new_thickness'};
    66                 elseif strcmpi(solution,'transient'),
    67                         fields={'transient_results(end).vy','transient_results(end).vz','transient_results(end).temperature','transient_results(end).melting','transient_results(end).thickness'};
    68                 end
     54                %compute fields to be checked
     55                fields=testsgetfields(md.type,solution);
    6956
    7057                %load archive

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

    r487 r519  
    3636                md=extrude(md,8,4);
    3737                md=setelementstype(md,'pattyn','all','stokes','all');
    38                 if strcmpi(package,'cielo_parallel'), md.cluster='wilkes'; end
    3938                if md.numberofgrids==832
    4039                        load Velocities; md.vx=0.8*vx; md.vy=0.8*vy; md.vz=0.8*vz; md.temperature=temperature-1; md.pressure=pressure;
    … …  
    4241
    4342                %compute solution
    44                 if (strcmpi(package,'cielo_serial') | strcmpi(package,'cielo_parallel')),
    45                         md=solve(md,'analysis_type',solution,'package','cielo');
    46                 else
    47                         md=solve(md,'analysis_type',solution,'package',package);
    48                 end
     43                [analysis_type sub_analysis_type]=testsgetanalysis(solution);
     44                [md packagefinal]=testsgetpackage(md,package);
     45                md=solve(md,'analysis_type',analysis_type,'sub_analysis_type',sub_analysis_type,'package',package);
    4946
    50                 %compute field to be checked
    51                 if strcmpi(solution,'diagnostic'),
    52                         fields={'vy','vz'};
    53                 elseif strcmpi(solution,'thermalsteady'),
    54                         fields={'temperature','melting'};
    55                 elseif strcmpi(solution,'thermaltransient'),
    56                         fields={'thermaltransient_results(end).temperature','thermaltransient_results(end).melting'};
    57                 elseif strcmpi(solution,'prognostic'),
    58                         fields={'new_thickness'};
    59                 elseif strcmpi(solution,'transient'),
    60                         fields={'transient_results(end).vy','transient_results(end).vz','transient_results(end).temperature','transient_results(end).melting','transient_results(end).thickness'};
    61                 end
     47                %compute fields to be checked
     48                fields=testsgetfields(md.type,solution);
    6249
    6350                %save new archive

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

    r487 r519  
    5050                end
    5151                md=modelextract(md,md.elementoniceshelf);
    52                 if strcmpi(package,'cielo_parallel'), md.cluster='wilkes'; end
    5352
    5453                %launch control method
    … …  
    5958                end
    6059
    61                 %compute field to be checked
     60                %compute fields to be checked
     61fields=testsgetfields(md.type,solution);
    6262                fields={'cont_vel','cont_J','cont_parameter'};
    6363

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

    r487 r519  
    4545                end
    4646                md=modelextract(md,md.elementoniceshelf);
    47                 if strcmpi(package,'cielo_parallel'), md.cluster='wilkes'; end
    4847
    4948                %launch control method
    … …  
    5453                end
    5554
    56                 %compute field to be checked
     55                %compute fields to be checked
     56fields=testsgetfields(md.type,solution);
    5757                fields={'cont_vel','cont_J','cont_parameter'};
    5858

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

    r487 r519  
    4949                        md.fit(:)=2;
    5050                end
    51                 if strcmpi(package,'cielo_parallel'), md.cluster='wilkes'; end
    5251
    5352                %launch control method
    … …  
    5857                end
    5958
    60                 %compute field to be checked
     59                %compute fields to be checked
     60fields=testsgetfields(md.type,solution);
    6161                fields={'cont_vel','cont_J','cont_parameter'};
    6262

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

    r487 r519  
    4444                        md.fit(:)=2;
    4545                end
    46                 if strcmpi(package,'cielo_parallel'), md.cluster='wilkes'; end
    4746
    4847                %launch control method
    … …  
    5352                end
    5453
    55                 %compute field to be checked
     54                %compute fields to be checked
     55fields=testsgetfields(md.type,solution);
    5656                fields={'cont_vel','cont_J','cont_parameter'};
    5757
Note: See TracChangeset for help on using the changeset viewer.