source: issm/oecreview/Archive/12678-13393/ISSM-12807-12808.diff@ 13394

Last change on this file since 13394 was 13394, checked in by Mathieu Morlighem, 13 years ago

Added 12678-13393
File size: 5.7 KB

  • u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/test/NightlyRun/test110.m

     
    33md=parameterize(md,'../Par/SquareShelfConstrained.par');
    44md=setflowequation(md,'macayeal','all');
    55md.cluster=generic('name',oshostname(),'np',3);
     6md.transient.requested_outputs=IceVolumeEnum();
    67
    78md=solve(md,TransientSolutionEnum);
    89
    910%Fields and tolerances to track changes
    10 field_names     ={'Vx1','Vy1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Vx2','Vy2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Vx3','Vy3','Vel3','Pressure3','Bed3','Surface3','Thickness3'};
    11 field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,...
    12                                                 1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,...
    13                                                 1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13};
     11field_names     ={'Vx1','Vy1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Volume1','Vx2','Vy2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Volume2','Vx3','Vy3','Vel3','Pressure3','Bed3','Surface3','Thickness3','Volume3'};
     12field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,...
     13                                                1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,...
     14                                                1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13};
    1415field_values={...
    1516        (md.results.TransientSolution(1).Vx),...
    1617        (md.results.TransientSolution(1).Vy),...
    … …  
    1920        (md.results.TransientSolution(1).Bed),...
    2021        (md.results.TransientSolution(1).Surface),...
    2122        (md.results.TransientSolution(1).Thickness),...
     23        (md.results.TransientSolution(1).IceVolume),...
    2224        (md.results.TransientSolution(2).Vx),...
    2325        (md.results.TransientSolution(2).Vy),...
    2426        (md.results.TransientSolution(2).Vel),...
    … …  
    2628        (md.results.TransientSolution(2).Bed),...
    2729        (md.results.TransientSolution(2).Surface),...
    2830        (md.results.TransientSolution(2).Thickness),...
     31        (md.results.TransientSolution(2).IceVolume),...
    2932        (md.results.TransientSolution(3).Vx),...
    3033        (md.results.TransientSolution(3).Vy),...
    3134        (md.results.TransientSolution(3).Vel),...
    … …  
    3336        (md.results.TransientSolution(3).Bed),...
    3437        (md.results.TransientSolution(3).Surface),...
    3538        (md.results.TransientSolution(3).Thickness),...
     39        (md.results.TransientSolution(3).IceVolume),...
    3640        };

  • u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/test/NightlyRun/test111.m

     
    44md=extrude(md,3,1);
    55md=setflowequation(md,'pattyn','all');
    66md.cluster=generic('name',oshostname(),'np',3);
     7md.transient.requested_outputs=IceVolumeEnum();
    78md=solve(md,TransientSolutionEnum);
    89
    910%Fields and tolerances to track changes
    1011field_names={...
    11         'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1','BasalforcingsMeltingRate1', ...
    12         'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2','BasalforcingsMeltingRate2', ...
    13         'Vx3','Vy3','Vz3','Vel3','Pressure3','Bed3','Surface3','Thickness3','Temperature3','BasalforcingsMeltingRate3'};
     12        'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface2','Thickness1','Temperature1','BasalforcingsMeltingRate1','Volume1', ...
     13        'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2','BasalforcingsMeltingRate2','Volume2', ...
     14        'Vx3','Vy3','Vz3','Vel3','Pressure3','Bed3','Surface3','Thickness3','Temperature3','BasalforcingsMeltingRate3','Volume3'};
    1415field_tolerances={...
    15         1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-05,...
    16         1e-09,1e-09,1e-08,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-05,...
    17         1e-09,1e-09,1e-08,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-05};
     16        1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-05,1e-13,...
     17        1e-09,1e-09,1e-08,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-05,1e-13,...
     18        1e-09,1e-09,1e-08,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-05,1e-13};
    1819field_values={...
    1920        (md.results.TransientSolution(1).Vx),...
    2021        (md.results.TransientSolution(1).Vy),...
    … …  
    2627        (md.results.TransientSolution(1).Thickness),...
    2728        (md.results.TransientSolution(1).Temperature),...
    2829        (md.results.TransientSolution(1).BasalforcingsMeltingRate),...
     30        (md.results.TransientSolution(1).IceVolume),...
    2931        (md.results.TransientSolution(2).Vx),...
    3032        (md.results.TransientSolution(2).Vy),...
    3133        (md.results.TransientSolution(2).Vz),...
    … …  
    3638        (md.results.TransientSolution(2).Thickness),...
    3739        (md.results.TransientSolution(2).Temperature),...
    3840        (md.results.TransientSolution(2).BasalforcingsMeltingRate),...
     41        (md.results.TransientSolution(2).IceVolume),...
    3942        (md.results.TransientSolution(3).Vx),...
    4043        (md.results.TransientSolution(3).Vy),...
    4144        (md.results.TransientSolution(3).Vz),...
    … …  
    4649        (md.results.TransientSolution(3).Thickness),...
    4750        (md.results.TransientSolution(3).Temperature),...
    4851        (md.results.TransientSolution(3).BasalforcingsMeltingRate),...
     52        (md.results.TransientSolution(3).IceVolume),...
    4953        };
Note: See TracBrowser for help on using the repository browser.