source: issm/oecreview/Archive/12678-13393/ISSM-12780-12781.diff@ 14312

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

Added 12678-13393
File size: 5.3 KB

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

     
    1111md.surfaceforcings.a_neg=-20000. - 0.00005*md.mesh.x + 0.0001*md.mesh.y;
    1212md.surfaceforcings.b_neg=250. + 0.000051*md.mesh.x - 0.00011*md.mesh.y;
    1313md.surfaceforcings.hc=(md.surfaceforcings.a_pos-md.surfaceforcings.a_neg)./(md.surfaceforcings.b_neg-md.surfaceforcings.b_pos);
     14md.transient.requested_outputs=TotalSmbEnum();
    1415md.cluster=generic('name',oshostname(),'np',3);
    1516md=solve(md,TransientSolutionEnum);
    1617
    1718%Fields and tolerances to track changes
    18 field_names     ={'Vx1','Vy1','Vz1','Vel1','Bed1','Surface1','Thickness1','Temperature1','SMB1','Vx2','Vy2','Vz2','Vel2','Bed2','Surface2','Thickness2','Temperature2','SMB2','Vx3','Vy3','Vz3','Vel3','Bed3','Surface3','Thickness3','Temperature3','SMB3'};
    19 field_tolerances={1e-09,1e-09,1e-09,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,...
    20                                                 1e-09,1e-09,1e-10,1e-09,1e-10,1e-10,1e-10,1e-10,1e-10,...
    21                                                 1e-09,5e-10,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10};
     19field_names     ={'Vx1','Vy1','Vz1','Vel1','Bed1','Surface1','Thickness1','Temperature1','SMB1','TotalSmb1','Vx2','Vy2','Vz2','Vel2','Bed2','Surface2','Thickness2','Temperature2','SMB2','TotalSmb1','Vx3','Vy3','Vz3','Vel3','Bed3','Surface3','Thickness3','Temperature3','SMB3','TotalSmb1'};
     20field_tolerances={1e-09,1e-09,1e-09,1e-09,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,...
     21                                                1e-09,1e-09,1e-10,1e-09,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,...
     22                                                1e-09,5e-09,1e-09,1e-09,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10};
    2223
    2324field_values={...
    2425        (md.results.TransientSolution(1).Vx),...
    … …  
    3031        (md.results.TransientSolution(1).Thickness),...
    3132        (md.results.TransientSolution(1).Temperature),...
    3233        (md.results.TransientSolution(1).SurfaceforcingsMassBalance),...
     34        (md.results.TransientSolution(1).TotalSmb),...
    3335        (md.results.TransientSolution(2).Vx),...
    3436        (md.results.TransientSolution(2).Vy),...
    3537        (md.results.TransientSolution(2).Vz),...
    … …  
    3941        (md.results.TransientSolution(2).Thickness),...
    4042        (md.results.TransientSolution(2).Temperature),...
    4143        (md.results.TransientSolution(2).SurfaceforcingsMassBalance),...
     44        (md.results.TransientSolution(2).TotalSmb),...
    4245        (md.results.TransientSolution(3).Vx),...
    4346        (md.results.TransientSolution(3).Vy),...
    4447        (md.results.TransientSolution(3).Vz),...
    … …  
    4851        (md.results.TransientSolution(3).Thickness),...
    4952        (md.results.TransientSolution(3).Temperature),...
    5053        (md.results.TransientSolution(3).SurfaceforcingsMassBalance),...
     54        (md.results.TransientSolution(3).TotalSmb),...
    5155        };

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

     
    1010md.surfaceforcings.a_neg=-20000. - 0.00005*md.mesh.x + 0.0001*md.mesh.y;
    1111md.surfaceforcings.b_neg=250. + 0.000051*md.mesh.x - 0.00011*md.mesh.y;
    1212md.surfaceforcings.hc=(md.surfaceforcings.a_pos-md.surfaceforcings.a_neg)./(md.surfaceforcings.b_neg-md.surfaceforcings.b_pos);
     13md.transient.requested_outputs=TotalSmbEnum();
    1314md.cluster=generic('name',oshostname(),'np',3);
    1415md=solve(md,TransientSolutionEnum);
    1516
    1617%Fields and tolerances to track changes
    17 field_names     ={'Vx1','Vy1','Vel1','Bed1','Surface1','Thickness1','SMB1','Vx2','Vy2','Vel2','Bed2','Surface2','Thickness2','SMB2','Vx3','Vy3','Vel3','Bed3','Surface3','Thickness3','SMB3'};
    18 field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13};
     18field_names     ={'Vx1','Vy1','Vel1','Bed1','Surface1','Thickness1','SMB1','TotalSmb1','Vx2','Vy2','Vel2','Bed2','Surface2','Thickness2','SMB2','TotalSmb2','Vx3','Vy3','Vel3','Bed3','Surface3','Thickness3','SMB3','TotalSmb3'};
     19field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13};
    1920field_values={...
    2021        (md.results.TransientSolution(1).Vx),...
    2122        (md.results.TransientSolution(1).Vy),...
    … …  
    2425        (md.results.TransientSolution(1).Surface),...
    2526        (md.results.TransientSolution(1).Thickness),...
    2627        (md.results.TransientSolution(1).SurfaceforcingsMassBalance),...
     28        (md.results.TransientSolution(1).TotalSmb),...
    2729        (md.results.TransientSolution(2).Vx),...
    2830        (md.results.TransientSolution(2).Vy),...
    2931        (md.results.TransientSolution(2).Vel),...
    3032        (md.results.TransientSolution(2).Bed),...
    3133        (md.results.TransientSolution(2).Surface),...
    3234        (md.results.TransientSolution(2).Thickness),...
     35        (md.results.TransientSolution(2).TotalSmb),...
    3336        (md.results.TransientSolution(2).SurfaceforcingsMassBalance),...
    3437        (md.results.TransientSolution(3).Vx),...
    3538        (md.results.TransientSolution(3).Vy),...
    … …  
    3841        (md.results.TransientSolution(3).Surface),...
    3942        (md.results.TransientSolution(3).Thickness),...
    4043        (md.results.TransientSolution(3).SurfaceforcingsMassBalance),...
     44        (md.results.TransientSolution(3).TotalSmb),...
    4145        };
Note: See TracBrowser for help on using the repository browser.