Index: /issm/trunk-jpl/test/NightlyRun/IdToName.m =================================================================== --- /issm/trunk-jpl/test/NightlyRun/IdToName.m (revision 18001) +++ /issm/trunk-jpl/test/NightlyRun/IdToName.m (revision 18002) @@ -67,4 +67,6 @@ case 236, name='SquareShelfTranIspddIsdeltaSSA2d'; case 237, name='SquareShelfTranIspddIsdeltaSSA3d'; + case 238, name='SquareShelfTranForceNoInterp2d'; + case 239, name='SquareShelfTranForceNoInterp3d'; case 270, name='SquareShelfStressSSA2dDamage'; case 271, name='SquareShelfStressSSA2dDamageTran'; @@ -109,4 +111,8 @@ case 334, name='SquareSheetConstrainedExtrudedHydrologyDC'; case 335, name='SquareSheetExtrudedHydrologyDCTwoLayers'; + case 336, name='SquareSheetConstrainedSmbComponents2d'; + case 337, name='SquareSheetConstrainedSmbComponents3d'; + case 338, name='SquareSheetConstrainedSmbMeltComponents2d'; + case 339, name='SquareSheetConstrainedSmbMeltComponents3d'; case 401, name='SquareSheetShelfStressSSA2d'; case 402, name='SquareSheetShelfStressSSA3d'; Index: /issm/trunk-jpl/test/NightlyRun/test238.m =================================================================== --- /issm/trunk-jpl/test/NightlyRun/test238.m (revision 18002) +++ /issm/trunk-jpl/test/NightlyRun/test238.m (revision 18002) @@ -0,0 +1,64 @@ +md=triangle(model(),'../Exp/Square.exp',150000.); +md=setmask(md,'all',''); +md=parameterize(md,'../Par/SquareShelf.par'); +md=setflowequation(md,'SSA','all'); +md.cluster=generic('name',oshostname(),'np',3); + +md.timestepping.time_step=1.; +md.settings.output_frequency=1; +md.timestepping.final_time=4.; +md.timestepping.interp_forcings=0; + +%Set up transient +smb=ones(md.mesh.numberofvertices,1)*3.6; +smb=[ smb smb*-1. ]; + +md.surfaceforcings.mass_balance=smb; +md.surfaceforcings.mass_balance(end+1,:)=[1.5 3.]; +md.transient.isthermal=0; + +md=solve(md,TransientSolutionEnum()); + +%Fields and tolerances to track changes +field_names={'Vx1','Vy1','Vel1','Pressure1','Bed1','Surface1','Thickness1','SurfaceforcingsMassBalance1', ... + 'Vx2','Vy2','Vel2','Pressure2','Bed2','Surface2','Thickness2','SurfaceforcingsMassBalance2', ... + 'Vx3','Vy3','Vel3','Pressure3','Bed3','Surface3','Thickness3','SurfaceforcingsMassBalance3', ... + 'Vx4','Vy4','Vel4','Pressure4','Bed4','Surface4','Thickness4','SurfaceforcingsMassBalance4'}; +field_tolerances={1e-09,1e-09,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,... + 1e-09,1e-09,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,... + 1e-09,1e-09,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,... + 1e-09,1e-09,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10}; +field_values={... + (md.results.TransientSolution(1).Vx),... + (md.results.TransientSolution(1).Vy),... + (md.results.TransientSolution(1).Vel),... + (md.results.TransientSolution(1).Pressure),... + (md.results.TransientSolution(1).Base),... + (md.results.TransientSolution(1).Surface),... + (md.results.TransientSolution(1).Thickness),... + (md.results.TransientSolution(1).SurfaceforcingsMassBalance),... + (md.results.TransientSolution(2).Vx),... + (md.results.TransientSolution(2).Vy),... + (md.results.TransientSolution(2).Vel),... + (md.results.TransientSolution(2).Pressure),... + (md.results.TransientSolution(2).Base),... + (md.results.TransientSolution(2).Surface),... + (md.results.TransientSolution(2).Thickness),... + (md.results.TransientSolution(2).SurfaceforcingsMassBalance),... + (md.results.TransientSolution(3).Vx),... + (md.results.TransientSolution(3).Vy),... + (md.results.TransientSolution(3).Vel),... + (md.results.TransientSolution(3).Pressure),... + (md.results.TransientSolution(3).Base),... + (md.results.TransientSolution(3).Surface),... + (md.results.TransientSolution(3).Thickness),... + (md.results.TransientSolution(3).SurfaceforcingsMassBalance),... + (md.results.TransientSolution(4).Vx),... + (md.results.TransientSolution(4).Vy),... + (md.results.TransientSolution(4).Vel),... + (md.results.TransientSolution(4).Pressure),... + (md.results.TransientSolution(4).Base),... + (md.results.TransientSolution(4).Surface),... + (md.results.TransientSolution(4).Thickness),... + (md.results.TransientSolution(4).SurfaceforcingsMassBalance),... + }; Index: /issm/trunk-jpl/test/NightlyRun/test238.py =================================================================== --- /issm/trunk-jpl/test/NightlyRun/test238.py (revision 18002) +++ /issm/trunk-jpl/test/NightlyRun/test238.py (revision 18002) @@ -0,0 +1,73 @@ +import numpy +from model import * +from triangle import * +from setmask import * +from parameterize import * +from setflowequation import * +from EnumDefinitions import * +from solve import * +from MatlabFuncs import * + +md=triangle(model(),'../Exp/Square.exp',150000.) +md=setmask(md,'all','') +md=parameterize(md,'../Par/SquareShelf.py') +md=setflowequation(md,'SSA','all') +md.cluster=generic('name',oshostname(),'np',3) + +md.timestepping.time_step=1. +md.settings.output_frequency=1 +md.timestepping.final_time=4. +md.timestepping.interp_forcings=False + +#Set up transient +smb=numpy.ones((md.mesh.numberofvertices,1))*3.6 +smb=numpy.hstack((smb,smb*-1.)) + +md.surfaceforcings.mass_balance=numpy.vstack((smb,[1.5,3.])) +md.transient.isthermal=False + +md=solve(md,TransientSolutionEnum()) + +#Fields and tolerances to track changes +field_names=['Vx1','Vy1','Vel1','Pressure1','Bed1','Surface1','Thickness1','SurfaceforcingsMassBalance1', \ + 'Vx2','Vy2','Vel2','Pressure2','Bed2','Surface2','Thickness2','SurfaceforcingsMassBalance2', \ + 'Vx3','Vy3','Vel3','Pressure3','Bed3','Surface3','Thickness3','SurfaceforcingsMassBalance3', \ + 'Vx4','Vy4','Vel4','Pressure4','Bed4','Surface4','Thickness4','SurfaceforcingsMassBalance4'] +field_tolerances=[1e-09,1e-09,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,\ + 1e-09,1e-09,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,\ + 1e-09,1e-09,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,\ + 1e-09,1e-09,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10] +field_values=[\ + md.results.TransientSolution[0].Vx,\ + md.results.TransientSolution[0].Vy,\ + md.results.TransientSolution[0].Vel,\ + md.results.TransientSolution[0].Pressure,\ + md.results.TransientSolution[0].Base,\ + md.results.TransientSolution[0].Surface,\ + md.results.TransientSolution[0].Thickness,\ + md.results.TransientSolution[0].SurfaceforcingsMassBalance,\ + md.results.TransientSolution[1].Vx,\ + md.results.TransientSolution[1].Vy,\ + md.results.TransientSolution[1].Vel,\ + md.results.TransientSolution[1].Pressure,\ + md.results.TransientSolution[1].Base,\ + md.results.TransientSolution[1].Surface,\ + md.results.TransientSolution[1].Thickness,\ + md.results.TransientSolution[1].SurfaceforcingsMassBalance,\ + md.results.TransientSolution[2].Vx,\ + md.results.TransientSolution[2].Vy,\ + md.results.TransientSolution[2].Vel,\ + md.results.TransientSolution[2].Pressure,\ + md.results.TransientSolution[2].Base,\ + md.results.TransientSolution[2].Surface,\ + md.results.TransientSolution[2].Thickness,\ + md.results.TransientSolution[2].SurfaceforcingsMassBalance,\ + md.results.TransientSolution[3].Vx,\ + md.results.TransientSolution[3].Vy,\ + md.results.TransientSolution[3].Vel,\ + md.results.TransientSolution[3].Pressure,\ + md.results.TransientSolution[3].Base,\ + md.results.TransientSolution[3].Surface,\ + md.results.TransientSolution[3].Thickness,\ + md.results.TransientSolution[3].SurfaceforcingsMassBalance,\ + ] Index: /issm/trunk-jpl/test/NightlyRun/test239.m =================================================================== --- /issm/trunk-jpl/test/NightlyRun/test239.m (revision 18002) +++ /issm/trunk-jpl/test/NightlyRun/test239.m (revision 18002) @@ -0,0 +1,69 @@ +md=triangle(model(),'../Exp/Square.exp',350000.); +md=setmask(md,'all',''); +md=parameterize(md,'../Par/SquareShelf.par'); +md=extrude(md,3,1.); +md=setflowequation(md,'HO','all'); +md.cluster=generic('name',oshostname(),'np',3); + +md.timestepping.time_step=1; +md.settings.output_frequency=1; +md.timestepping.final_time=4; +md.timestepping.interp_forcings=0; + +%Set up transient +smb=ones(md.mesh.numberofvertices,1)*3.6; +smb=[ smb smb*-1. ]; + +md.surfaceforcings.mass_balance=smb; +md.surfaceforcings.mass_balance(end+1,:)=[1.5 3.]; +md.transient.isthermal=0; + +md=solve(md,TransientSolutionEnum()); + +%Fields and tolerances to track changes +field_names ={'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','SurfaceforcingsMassBalance1', ... + 'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','SurfaceforcingsMassBalance2', ... + 'Vx3','Vy3','Vz3','Vel3','Pressure3','Bed3','Surface3','Thickness3','SurfaceforcingsMassBalance3', ... + 'Vx4','Vy4','Vz4','Vel4','Pressure4','Bed4','Surface4','Thickness4','SurfaceforcingsMassbalance4'}; +field_tolerances={1e-09,1e-09,1e-09,1e-09,1e-10,1e-10,1e-10,1e-10,1e-13,... + 1e-09,1e-09,1e-09,1e-09,1e-10,1e-10,1e-10,1e-10,1e-13,... + 1e-09,1e-09,1e-09,1e-09,1e-10,1e-10,1e-10,1e-10,1e-13,... + 1e-09,1e-09,1e-09,1e-09,1e-10,1e-10,1e-10,1e-10,1e-13}; +field_values={... + (md.results.TransientSolution(1).Vx),... + (md.results.TransientSolution(1).Vy),... + (md.results.TransientSolution(1).Vz),... + (md.results.TransientSolution(1).Vel),... + (md.results.TransientSolution(1).Pressure),... + (md.results.TransientSolution(1).Base),... + (md.results.TransientSolution(1).Surface),... + (md.results.TransientSolution(1).Thickness),... + (md.results.TransientSolution(1).SurfaceforcingsMassBalance),... + (md.results.TransientSolution(2).Vx),... + (md.results.TransientSolution(2).Vy),... + (md.results.TransientSolution(2).Vz),... + (md.results.TransientSolution(2).Vel),... + (md.results.TransientSolution(2).Pressure),... + (md.results.TransientSolution(2).Base),... + (md.results.TransientSolution(2).Surface),... + (md.results.TransientSolution(2).Thickness),... + (md.results.TransientSolution(2).SurfaceforcingsMassBalance),... + (md.results.TransientSolution(3).Vx),... + (md.results.TransientSolution(3).Vy),... + (md.results.TransientSolution(3).Vz),... + (md.results.TransientSolution(3).Vel),... + (md.results.TransientSolution(3).Pressure),... + (md.results.TransientSolution(3).Base),... + (md.results.TransientSolution(3).Surface),... + (md.results.TransientSolution(3).Thickness),... + (md.results.TransientSolution(3).SurfaceforcingsMassBalance),... + (md.results.TransientSolution(4).Vx),... + (md.results.TransientSolution(4).Vy),... + (md.results.TransientSolution(4).Vz),... + (md.results.TransientSolution(4).Vel),... + (md.results.TransientSolution(4).Pressure),... + (md.results.TransientSolution(4).Base),... + (md.results.TransientSolution(4).Surface),... + (md.results.TransientSolution(4).Thickness),... + (md.results.TransientSolution(4).SurfaceforcingsMassBalance),... + }; Index: /issm/trunk-jpl/test/NightlyRun/test239.py =================================================================== --- /issm/trunk-jpl/test/NightlyRun/test239.py (revision 18002) +++ /issm/trunk-jpl/test/NightlyRun/test239.py (revision 18002) @@ -0,0 +1,78 @@ +import numpy +from model import * +from triangle import * +from setmask import * +from parameterize import * +from setflowequation import * +from EnumDefinitions import * +from solve import * +from MatlabFuncs import * + +md=triangle(model(),'../Exp/Square.exp',350000.) +md=setmask(md,'all','') +md=parameterize(md,'../Par/SquareShelf.py') +md.extrude(3,1.) +md=setflowequation(md,'HO','all') +md.cluster=generic('name',oshostname(),'np',3) + +md.timestepping.time_step=1. +md.settings.output_frequency=1 +md.timestepping.final_time=4. +md.timestepping.interp_forcings=False + +#Set up transient +smb=numpy.ones((md.mesh.numberofvertices,1))*3.6 +smb=numpy.hstack((smb,smb*-1.)) + +md.surfaceforcings.mass_balance=numpy.vstack((smb,[1.5,3.])) +md.transient.isthermal=False + +md=solve(md,TransientSolutionEnum()) + +#Fields and tolerances to track changes +field_names =['Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','SurfaceforcingsMassBalance1', \ + 'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','SurfaceforcingsMassBalance2', \ + 'Vx3','Vy3','Vz3','Vel3','Pressure3','Bed3','Surface3','Thickness3','SurfaceforcingsMassBalance3', \ + 'Vx4','Vy4','Vz4','Vel4','Pressure4','Bed4','Surface4','Thickness4','SurfaceforcingsMassbalance4'] +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,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] +field_values=[\ + md.results.TransientSolution[0].Vx,\ + md.results.TransientSolution[0].Vy,\ + md.results.TransientSolution[0].Vz,\ + md.results.TransientSolution[0].Vel,\ + md.results.TransientSolution[0].Pressure,\ + md.results.TransientSolution[0].Base,\ + md.results.TransientSolution[0].Surface,\ + md.results.TransientSolution[0].Thickness,\ + md.results.TransientSolution[0].SurfaceforcingsMassBalance,\ + md.results.TransientSolution[1].Vx,\ + md.results.TransientSolution[1].Vy,\ + md.results.TransientSolution[1].Vz,\ + md.results.TransientSolution[1].Vel,\ + md.results.TransientSolution[1].Pressure,\ + md.results.TransientSolution[1].Base,\ + md.results.TransientSolution[1].Surface,\ + md.results.TransientSolution[1].Thickness,\ + md.results.TransientSolution[1].SurfaceforcingsMassBalance,\ + md.results.TransientSolution[2].Vx,\ + md.results.TransientSolution[2].Vy,\ + md.results.TransientSolution[2].Vz,\ + md.results.TransientSolution[2].Vel,\ + md.results.TransientSolution[2].Pressure,\ + md.results.TransientSolution[2].Base,\ + md.results.TransientSolution[2].Surface,\ + md.results.TransientSolution[2].Thickness,\ + md.results.TransientSolution[2].SurfaceforcingsMassBalance,\ + md.results.TransientSolution[3].Vx,\ + md.results.TransientSolution[3].Vy,\ + md.results.TransientSolution[3].Vz,\ + md.results.TransientSolution[3].Vel,\ + md.results.TransientSolution[3].Pressure,\ + md.results.TransientSolution[3].Base,\ + md.results.TransientSolution[3].Surface,\ + md.results.TransientSolution[3].Thickness,\ + md.results.TransientSolution[3].SurfaceforcingsMassBalance,\ + ] Index: /issm/trunk-jpl/test/NightlyRun/test336.m =================================================================== --- /issm/trunk-jpl/test/NightlyRun/test336.m (revision 18002) +++ /issm/trunk-jpl/test/NightlyRun/test336.m (revision 18002) @@ -0,0 +1,65 @@ +md=triangle(model(),'../Exp/Square.exp',150000.); +md=setmask(md,'all',''); +md=parameterize(md,'../Par/SquareShelf.par'); +md=setflowequation(md,'SSA','all'); +md.cluster=generic('name',oshostname(),'np',3); + +md.timestepping.time_step=1.; +md.settings.output_frequency=1; +md.timestepping.final_time=4.; + +%Set up transient +smb=ones(md.mesh.numberofvertices,1)*3.6; +smb=[ smb smb*-1. ]; + +md.surfaceforcings=SMBcomponents(); +md.surfaceforcings.accumulation=[smb*2; [1.5 3.]]; +md.surfaceforcings.runoff=[smb/2; [1.5 3.]]; +md.surfaceforcings.evaporation=[smb/2; [1.5 3.]]; +md.transient.isthermal=0; + +md=solve(md,TransientSolutionEnum()); + +%Fields and tolerances to track changes +field_names={'Vx1','Vy1','Vel1','Pressure1','Bed1','Surface1','Thickness1','SurfaceforcingsMassBalance1', ... + 'Vx2','Vy2','Vel2','Pressure2','Bed2','Surface2','Thickness2','SurfaceforcingsMassBalance2', ... + 'Vx3','Vy3','Vel3','Pressure3','Bed3','Surface3','Thickness3','SurfaceforcingsMassBalance3', ... + 'Vx4','Vy4','Vel4','Pressure4','Bed4','Surface4','Thickness4','SurfaceforcingsMassBalance4'}; +field_tolerances={1e-09,1e-09,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,... + 1e-09,1e-09,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,... + 1e-09,1e-09,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,... + 1e-09,1e-09,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10}; +field_values={... + (md.results.TransientSolution(1).Vx),... + (md.results.TransientSolution(1).Vy),... + (md.results.TransientSolution(1).Vel),... + (md.results.TransientSolution(1).Pressure),... + (md.results.TransientSolution(1).Base),... + (md.results.TransientSolution(1).Surface),... + (md.results.TransientSolution(1).Thickness),... + (md.results.TransientSolution(1).SurfaceforcingsMassBalance),... + (md.results.TransientSolution(2).Vx),... + (md.results.TransientSolution(2).Vy),... + (md.results.TransientSolution(2).Vel),... + (md.results.TransientSolution(2).Pressure),... + (md.results.TransientSolution(2).Base),... + (md.results.TransientSolution(2).Surface),... + (md.results.TransientSolution(2).Thickness),... + (md.results.TransientSolution(2).SurfaceforcingsMassBalance),... + (md.results.TransientSolution(3).Vx),... + (md.results.TransientSolution(3).Vy),... + (md.results.TransientSolution(3).Vel),... + (md.results.TransientSolution(3).Pressure),... + (md.results.TransientSolution(3).Base),... + (md.results.TransientSolution(3).Surface),... + (md.results.TransientSolution(3).Thickness),... + (md.results.TransientSolution(3).SurfaceforcingsMassBalance),... + (md.results.TransientSolution(4).Vx),... + (md.results.TransientSolution(4).Vy),... + (md.results.TransientSolution(4).Vel),... + (md.results.TransientSolution(4).Pressure),... + (md.results.TransientSolution(4).Base),... + (md.results.TransientSolution(4).Surface),... + (md.results.TransientSolution(4).Thickness),... + (md.results.TransientSolution(4).SurfaceforcingsMassBalance),... + }; Index: /issm/trunk-jpl/test/NightlyRun/test336.py =================================================================== --- /issm/trunk-jpl/test/NightlyRun/test336.py (revision 18002) +++ /issm/trunk-jpl/test/NightlyRun/test336.py (revision 18002) @@ -0,0 +1,75 @@ +import numpy +from model import * +from triangle import * +from setmask import * +from parameterize import * +from setflowequation import * +from EnumDefinitions import * +from solve import * +from MatlabFuncs import * + +md=triangle(model(),'../Exp/Square.exp',150000.) +md=setmask(md,'all','') +md=parameterize(md,'../Par/SquareShelf.py') +md=setflowequation(md,'SSA','all') +md.cluster=generic('name',oshostname(),'np',3) + +md.timestepping.time_step=1. +md.settings.output_frequency=1 +md.timestepping.final_time=4. + +#Set up transient +smb=numpy.ones((md.mesh.numberofvertices,1))*3.6 +smb=numpy.hstack((smb,smb*-1.)) + +md.surfaceforcings=SMBcomponents(); +md.surfaceforcings.accumulation=numpy.vstack((smb*2, [1.5,3.])); +md.surfaceforcings.runoff=numpy.vstack((smb/2, [1.5,3.])); +md.surfaceforcings.evaporation=numpy.vstack((smb/2, [1.5,3.])); +md.transient.isthermal=False + +md=solve(md,TransientSolutionEnum()) + +#Fields and tolerances to track changes +field_names=['Vx1','Vy1','Vel1','Pressure1','Bed1','Surface1','Thickness1','SurfaceforcingsMassBalance1', \ + 'Vx2','Vy2','Vel2','Pressure2','Bed2','Surface2','Thickness2','SurfaceforcingsMassBalance2', \ + 'Vx3','Vy3','Vel3','Pressure3','Bed3','Surface3','Thickness3','SurfaceforcingsMassBalance3', \ + 'Vx4','Vy4','Vel4','Pressure4','Bed4','Surface4','Thickness4','SurfaceforcingsMassBalance4'] +field_tolerances=[1e-09,1e-09,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,\ + 1e-09,1e-09,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,\ + 1e-09,1e-09,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,\ + 1e-09,1e-09,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10] +field_values=[\ + md.results.TransientSolution[0].Vx,\ + md.results.TransientSolution[0].Vy,\ + md.results.TransientSolution[0].Vel,\ + md.results.TransientSolution[0].Pressure,\ + md.results.TransientSolution[0].Base,\ + md.results.TransientSolution[0].Surface,\ + md.results.TransientSolution[0].Thickness,\ + md.results.TransientSolution[0].SurfaceforcingsMassBalance,\ + md.results.TransientSolution[1].Vx,\ + md.results.TransientSolution[1].Vy,\ + md.results.TransientSolution[1].Vel,\ + md.results.TransientSolution[1].Pressure,\ + md.results.TransientSolution[1].Base,\ + md.results.TransientSolution[1].Surface,\ + md.results.TransientSolution[1].Thickness,\ + md.results.TransientSolution[1].SurfaceforcingsMassBalance,\ + md.results.TransientSolution[2].Vx,\ + md.results.TransientSolution[2].Vy,\ + md.results.TransientSolution[2].Vel,\ + md.results.TransientSolution[2].Pressure,\ + md.results.TransientSolution[2].Base,\ + md.results.TransientSolution[2].Surface,\ + md.results.TransientSolution[2].Thickness,\ + md.results.TransientSolution[2].SurfaceforcingsMassBalance,\ + md.results.TransientSolution[3].Vx,\ + md.results.TransientSolution[3].Vy,\ + md.results.TransientSolution[3].Vel,\ + md.results.TransientSolution[3].Pressure,\ + md.results.TransientSolution[3].Base,\ + md.results.TransientSolution[3].Surface,\ + md.results.TransientSolution[3].Thickness,\ + md.results.TransientSolution[3].SurfaceforcingsMassBalance,\ + ] Index: /issm/trunk-jpl/test/NightlyRun/test337.m =================================================================== --- /issm/trunk-jpl/test/NightlyRun/test337.m (revision 18002) +++ /issm/trunk-jpl/test/NightlyRun/test337.m (revision 18002) @@ -0,0 +1,66 @@ +md=triangle(model(),'../Exp/Square.exp',150000.); +md=setmask(md,'all',''); +md=parameterize(md,'../Par/SquareShelf.par'); +md=extrude(md,3,1.); +md=setflowequation(md,'HO','all'); +md.cluster=generic('name',oshostname(),'np',3); + +md.timestepping.time_step=1.; +md.settings.output_frequency=1; +md.timestepping.final_time=4.; + +%Set up transient +smb=ones(md.mesh.numberofvertices,1)*3.6; +smb=[ smb smb*-1. ]; + +md.surfaceforcings=SMBcomponents(); +md.surfaceforcings.accumulation=[smb*2; [1.5 3.]]; +md.surfaceforcings.runoff=[smb/2; [1.5 3.]]; +md.surfaceforcings.evaporation=[smb/2; [1.5 3.]]; +md.transient.isthermal=0; + +md=solve(md,TransientSolutionEnum()); + +%Fields and tolerances to track changes +field_names={'Vx1','Vy1','Vel1','Pressure1','Bed1','Surface1','Thickness1','SurfaceforcingsMassBalance1', ... + 'Vx2','Vy2','Vel2','Pressure2','Bed2','Surface2','Thickness2','SurfaceforcingsMassBalance2', ... + 'Vx3','Vy3','Vel3','Pressure3','Bed3','Surface3','Thickness3','SurfaceforcingsMassBalance3', ... + 'Vx4','Vy4','Vel4','Pressure4','Bed4','Surface4','Thickness4','SurfaceforcingsMassBalance4'}; +field_tolerances={1e-09,1e-09,1e-09,1e-10,1e-10,1e-10,1e-10,1e-10,... + 1e-09,1e-09,1e-09,1e-10,1e-10,1e-10,1e-10,1e-10,... + 1e-09,1e-09,1e-09,1e-10,1e-10,1e-10,1e-10,1e-10,... + 1e-09,1e-09,1e-09,1e-10,1e-10,1e-10,1e-10,1e-10}; +field_values={... + (md.results.TransientSolution(1).Vx),... + (md.results.TransientSolution(1).Vy),... + (md.results.TransientSolution(1).Vel),... + (md.results.TransientSolution(1).Pressure),... + (md.results.TransientSolution(1).Base),... + (md.results.TransientSolution(1).Surface),... + (md.results.TransientSolution(1).Thickness),... + (md.results.TransientSolution(1).SurfaceforcingsMassBalance),... + (md.results.TransientSolution(2).Vx),... + (md.results.TransientSolution(2).Vy),... + (md.results.TransientSolution(2).Vel),... + (md.results.TransientSolution(2).Pressure),... + (md.results.TransientSolution(2).Base),... + (md.results.TransientSolution(2).Surface),... + (md.results.TransientSolution(2).Thickness),... + (md.results.TransientSolution(2).SurfaceforcingsMassBalance),... + (md.results.TransientSolution(3).Vx),... + (md.results.TransientSolution(3).Vy),... + (md.results.TransientSolution(3).Vel),... + (md.results.TransientSolution(3).Pressure),... + (md.results.TransientSolution(3).Base),... + (md.results.TransientSolution(3).Surface),... + (md.results.TransientSolution(3).Thickness),... + (md.results.TransientSolution(3).SurfaceforcingsMassBalance),... + (md.results.TransientSolution(4).Vx),... + (md.results.TransientSolution(4).Vy),... + (md.results.TransientSolution(4).Vel),... + (md.results.TransientSolution(4).Pressure),... + (md.results.TransientSolution(4).Base),... + (md.results.TransientSolution(4).Surface),... + (md.results.TransientSolution(4).Thickness),... + (md.results.TransientSolution(4).SurfaceforcingsMassBalance),... + }; Index: /issm/trunk-jpl/test/NightlyRun/test337.py =================================================================== --- /issm/trunk-jpl/test/NightlyRun/test337.py (revision 18002) +++ /issm/trunk-jpl/test/NightlyRun/test337.py (revision 18002) @@ -0,0 +1,76 @@ +import numpy +from model import * +from triangle import * +from setmask import * +from parameterize import * +from setflowequation import * +from EnumDefinitions import * +from solve import * +from MatlabFuncs import * + +md=triangle(model(),'../Exp/Square.exp',150000.) +md=setmask(md,'all','') +md=parameterize(md,'../Par/SquareShelf.py') +md.extrude(3,1.) +md=setflowequation(md,'HO','all') +md.cluster=generic('name',oshostname(),'np',3) + +md.timestepping.time_step=1. +md.settings.output_frequency=1 +md.timestepping.final_time=4. + +#Set up transient +smb=numpy.ones((md.mesh.numberofvertices,1))*3.6 +smb=numpy.hstack((smb,smb*-1.)) + +md.surfaceforcings=SMBcomponents(); +md.surfaceforcings.accumulation=numpy.vstack((smb*2, [1.5,3.])); +md.surfaceforcings.runoff=numpy.vstack((smb/2, [1.5,3.])); +md.surfaceforcings.evaporation=numpy.vstack((smb/2, [1.5,3.])); +md.transient.isthermal=False + +md=solve(md,TransientSolutionEnum()) + +#Fields and tolerances to track changes +field_names=['Vx1','Vy1','Vel1','Pressure1','Bed1','Surface1','Thickness1','SurfaceforcingsMassBalance1', \ + 'Vx2','Vy2','Vel2','Pressure2','Bed2','Surface2','Thickness2','SurfaceforcingsMassBalance2', \ + 'Vx3','Vy3','Vel3','Pressure3','Bed3','Surface3','Thickness3','SurfaceforcingsMassBalance3', \ + 'Vx4','Vy4','Vel4','Pressure4','Bed4','Surface4','Thickness4','SurfaceforcingsMassBalance4'] +field_tolerances=[1e-09,1e-09,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,\ + 1e-09,1e-09,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,\ + 1e-09,1e-09,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,\ + 1e-09,1e-09,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10] +field_values=[\ + md.results.TransientSolution[0].Vx,\ + md.results.TransientSolution[0].Vy,\ + md.results.TransientSolution[0].Vel,\ + md.results.TransientSolution[0].Pressure,\ + md.results.TransientSolution[0].Base,\ + md.results.TransientSolution[0].Surface,\ + md.results.TransientSolution[0].Thickness,\ + md.results.TransientSolution[0].SurfaceforcingsMassBalance,\ + md.results.TransientSolution[1].Vx,\ + md.results.TransientSolution[1].Vy,\ + md.results.TransientSolution[1].Vel,\ + md.results.TransientSolution[1].Pressure,\ + md.results.TransientSolution[1].Base,\ + md.results.TransientSolution[1].Surface,\ + md.results.TransientSolution[1].Thickness,\ + md.results.TransientSolution[1].SurfaceforcingsMassBalance,\ + md.results.TransientSolution[2].Vx,\ + md.results.TransientSolution[2].Vy,\ + md.results.TransientSolution[2].Vel,\ + md.results.TransientSolution[2].Pressure,\ + md.results.TransientSolution[2].Base,\ + md.results.TransientSolution[2].Surface,\ + md.results.TransientSolution[2].Thickness,\ + md.results.TransientSolution[2].SurfaceforcingsMassBalance,\ + md.results.TransientSolution[3].Vx,\ + md.results.TransientSolution[3].Vy,\ + md.results.TransientSolution[3].Vel,\ + md.results.TransientSolution[3].Pressure,\ + md.results.TransientSolution[3].Base,\ + md.results.TransientSolution[3].Surface,\ + md.results.TransientSolution[3].Thickness,\ + md.results.TransientSolution[3].SurfaceforcingsMassBalance,\ + ] Index: /issm/trunk-jpl/test/NightlyRun/test338.m =================================================================== --- /issm/trunk-jpl/test/NightlyRun/test338.m (revision 18002) +++ /issm/trunk-jpl/test/NightlyRun/test338.m (revision 18002) @@ -0,0 +1,66 @@ +md=triangle(model(),'../Exp/Square.exp',150000.); +md=setmask(md,'all',''); +md=parameterize(md,'../Par/SquareShelf.par'); +md=setflowequation(md,'SSA','all'); +md.cluster=generic('name',oshostname(),'np',3); + +md.timestepping.time_step=1.; +md.settings.output_frequency=1; +md.timestepping.final_time=4.; + +%Set up transient +smb=ones(md.mesh.numberofvertices,1)*3.6; +smb=[ smb smb*-1. ]; + +md.surfaceforcings=SMBmeltcomponents(); +md.surfaceforcings.accumulation=[smb; [1.5 3.]]; +md.surfaceforcings.melt=[smb/2; [1.5 3.]]; +md.surfaceforcings.refreeze=[smb; [1.5 3.]]; +md.surfaceforcings.evaporation=[smb/2; [1.5 3.]]; +md.transient.isthermal=0; + +md=solve(md,TransientSolutionEnum()); + +%Fields and tolerances to track changes +field_names={'Vx1','Vy1','Vel1','Pressure1','Bed1','Surface1','Thickness1','SurfaceforcingsMassBalance1', ... + 'Vx2','Vy2','Vel2','Pressure2','Bed2','Surface2','Thickness2','SurfaceforcingsMassBalance2', ... + 'Vx3','Vy3','Vel3','Pressure3','Bed3','Surface3','Thickness3','SurfaceforcingsMassBalance3', ... + 'Vx4','Vy4','Vel4','Pressure4','Bed4','Surface4','Thickness4','SurfaceforcingsMassBalance4'}; +field_tolerances={1e-09,1e-09,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,... + 1e-09,1e-09,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,... + 1e-09,1e-09,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,... + 1e-09,1e-09,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10}; +field_values={... + (md.results.TransientSolution(1).Vx),... + (md.results.TransientSolution(1).Vy),... + (md.results.TransientSolution(1).Vel),... + (md.results.TransientSolution(1).Pressure),... + (md.results.TransientSolution(1).Base),... + (md.results.TransientSolution(1).Surface),... + (md.results.TransientSolution(1).Thickness),... + (md.results.TransientSolution(1).SurfaceforcingsMassBalance),... + (md.results.TransientSolution(2).Vx),... + (md.results.TransientSolution(2).Vy),... + (md.results.TransientSolution(2).Vel),... + (md.results.TransientSolution(2).Pressure),... + (md.results.TransientSolution(2).Base),... + (md.results.TransientSolution(2).Surface),... + (md.results.TransientSolution(2).Thickness),... + (md.results.TransientSolution(2).SurfaceforcingsMassBalance),... + (md.results.TransientSolution(3).Vx),... + (md.results.TransientSolution(3).Vy),... + (md.results.TransientSolution(3).Vel),... + (md.results.TransientSolution(3).Pressure),... + (md.results.TransientSolution(3).Base),... + (md.results.TransientSolution(3).Surface),... + (md.results.TransientSolution(3).Thickness),... + (md.results.TransientSolution(3).SurfaceforcingsMassBalance),... + (md.results.TransientSolution(4).Vx),... + (md.results.TransientSolution(4).Vy),... + (md.results.TransientSolution(4).Vel),... + (md.results.TransientSolution(4).Pressure),... + (md.results.TransientSolution(4).Base),... + (md.results.TransientSolution(4).Surface),... + (md.results.TransientSolution(4).Thickness),... + (md.results.TransientSolution(4).SurfaceforcingsMassBalance),... + }; Index: /issm/trunk-jpl/test/NightlyRun/test338.py =================================================================== --- /issm/trunk-jpl/test/NightlyRun/test338.py (revision 18002) +++ /issm/trunk-jpl/test/NightlyRun/test338.py (revision 18002) @@ -0,0 +1,76 @@ +import numpy +from model import * +from triangle import * +from setmask import * +from parameterize import * +from setflowequation import * +from EnumDefinitions import * +from solve import * +from MatlabFuncs import * + +md=triangle(model(),'../Exp/Square.exp',150000.) +md=setmask(md,'all','') +md=parameterize(md,'../Par/SquareShelf.py') +md=setflowequation(md,'SSA','all') +md.cluster=generic('name',oshostname(),'np',3) + +md.timestepping.time_step=1. +md.settings.output_frequency=1 +md.timestepping.final_time=4. + +#Set up transient +smb=numpy.ones((md.mesh.numberofvertices,1))*3.6 +smb=numpy.hstack((smb,smb*-1.)) + +md.surfaceforcings=SMBmeltcomponents(); +md.surfaceforcings.accumulation=numpy.vstack((smb, [1.5,3.])); +md.surfaceforcings.evaporation=numpy.vstack((smb/2, [1.5,3.])); +md.surfaceforcings.melt=numpy.vstack((smb/2, [1.5,3.])); +md.surfaceforcings.refreeze=numpy.vstack((smb, [1.5,3.])); +md.transient.isthermal=False + +md=solve(md,TransientSolutionEnum()) + +#Fields and tolerances to track changes +field_names=['Vx1','Vy1','Vel1','Pressure1','Bed1','Surface1','Thickness1','SurfaceforcingsMassBalance1', \ + 'Vx2','Vy2','Vel2','Pressure2','Bed2','Surface2','Thickness2','SurfaceforcingsMassBalance2', \ + 'Vx3','Vy3','Vel3','Pressure3','Bed3','Surface3','Thickness3','SurfaceforcingsMassBalance3', \ + 'Vx4','Vy4','Vel4','Pressure4','Bed4','Surface4','Thickness4','SurfaceforcingsMassBalance4'] +field_tolerances=[1e-09,1e-09,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,\ + 1e-09,1e-09,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,\ + 1e-09,1e-09,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,\ + 1e-09,1e-09,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10] +field_values=[\ + md.results.TransientSolution[0].Vx,\ + md.results.TransientSolution[0].Vy,\ + md.results.TransientSolution[0].Vel,\ + md.results.TransientSolution[0].Pressure,\ + md.results.TransientSolution[0].Base,\ + md.results.TransientSolution[0].Surface,\ + md.results.TransientSolution[0].Thickness,\ + md.results.TransientSolution[0].SurfaceforcingsMassBalance,\ + md.results.TransientSolution[1].Vx,\ + md.results.TransientSolution[1].Vy,\ + md.results.TransientSolution[1].Vel,\ + md.results.TransientSolution[1].Pressure,\ + md.results.TransientSolution[1].Base,\ + md.results.TransientSolution[1].Surface,\ + md.results.TransientSolution[1].Thickness,\ + md.results.TransientSolution[1].SurfaceforcingsMassBalance,\ + md.results.TransientSolution[2].Vx,\ + md.results.TransientSolution[2].Vy,\ + md.results.TransientSolution[2].Vel,\ + md.results.TransientSolution[2].Pressure,\ + md.results.TransientSolution[2].Base,\ + md.results.TransientSolution[2].Surface,\ + md.results.TransientSolution[2].Thickness,\ + md.results.TransientSolution[2].SurfaceforcingsMassBalance,\ + md.results.TransientSolution[3].Vx,\ + md.results.TransientSolution[3].Vy,\ + md.results.TransientSolution[3].Vel,\ + md.results.TransientSolution[3].Pressure,\ + md.results.TransientSolution[3].Base,\ + md.results.TransientSolution[3].Surface,\ + md.results.TransientSolution[3].Thickness,\ + md.results.TransientSolution[3].SurfaceforcingsMassBalance,\ + ] Index: /issm/trunk-jpl/test/NightlyRun/test339.m =================================================================== --- /issm/trunk-jpl/test/NightlyRun/test339.m (revision 18002) +++ /issm/trunk-jpl/test/NightlyRun/test339.m (revision 18002) @@ -0,0 +1,67 @@ +md=triangle(model(),'../Exp/Square.exp',150000.); +md=setmask(md,'all',''); +md=parameterize(md,'../Par/SquareShelf.par'); +md=extrude(md,3,1.); +md=setflowequation(md,'HO','all'); +md.cluster=generic('name',oshostname(),'np',3); + +md.timestepping.time_step=1.; +md.settings.output_frequency=1; +md.timestepping.final_time=4.; + +%Set up transient +smb=ones(md.mesh.numberofvertices,1)*3.6; +smb=[ smb smb*-1. ]; + +md.surfaceforcings=SMBmeltcomponents(); +md.surfaceforcings.accumulation=[smb; [1.5 3.]]; +md.surfaceforcings.melt=[smb/2; [1.5 3.]]; +md.surfaceforcings.refreeze=[smb; [1.5 3.]]; +md.surfaceforcings.evaporation=[smb/2; [1.5 3.]]; +md.transient.isthermal=0; + +md=solve(md,TransientSolutionEnum()); + +%Fields and tolerances to track changes +field_names={'Vx1','Vy1','Vel1','Pressure1','Bed1','Surface1','Thickness1','SurfaceforcingsMassBalance1', ... + 'Vx2','Vy2','Vel2','Pressure2','Bed2','Surface2','Thickness2','SurfaceforcingsMassBalance2', ... + 'Vx3','Vy3','Vel3','Pressure3','Bed3','Surface3','Thickness3','SurfaceforcingsMassBalance3', ... + 'Vx4','Vy4','Vel4','Pressure4','Bed4','Surface4','Thickness4','SurfaceforcingsMassBalance4'}; +field_tolerances={1e-09,1e-09,1e-09,1e-10,1e-10,1e-10,1e-10,1e-10,... + 1e-09,1e-09,1e-09,1e-10,1e-10,1e-10,1e-10,1e-10,... + 1e-09,1e-09,1e-09,1e-10,1e-10,1e-10,1e-10,1e-10,... + 1e-09,1e-09,1e-09,1e-10,1e-10,1e-10,1e-10,1e-10}; +field_values={... + (md.results.TransientSolution(1).Vx),... + (md.results.TransientSolution(1).Vy),... + (md.results.TransientSolution(1).Vel),... + (md.results.TransientSolution(1).Pressure),... + (md.results.TransientSolution(1).Base),... + (md.results.TransientSolution(1).Surface),... + (md.results.TransientSolution(1).Thickness),... + (md.results.TransientSolution(1).SurfaceforcingsMassBalance),... + (md.results.TransientSolution(2).Vx),... + (md.results.TransientSolution(2).Vy),... + (md.results.TransientSolution(2).Vel),... + (md.results.TransientSolution(2).Pressure),... + (md.results.TransientSolution(2).Base),... + (md.results.TransientSolution(2).Surface),... + (md.results.TransientSolution(2).Thickness),... + (md.results.TransientSolution(2).SurfaceforcingsMassBalance),... + (md.results.TransientSolution(3).Vx),... + (md.results.TransientSolution(3).Vy),... + (md.results.TransientSolution(3).Vel),... + (md.results.TransientSolution(3).Pressure),... + (md.results.TransientSolution(3).Base),... + (md.results.TransientSolution(3).Surface),... + (md.results.TransientSolution(3).Thickness),... + (md.results.TransientSolution(3).SurfaceforcingsMassBalance),... + (md.results.TransientSolution(4).Vx),... + (md.results.TransientSolution(4).Vy),... + (md.results.TransientSolution(4).Vel),... + (md.results.TransientSolution(4).Pressure),... + (md.results.TransientSolution(4).Base),... + (md.results.TransientSolution(4).Surface),... + (md.results.TransientSolution(4).Thickness),... + (md.results.TransientSolution(4).SurfaceforcingsMassBalance),... + }; Index: /issm/trunk-jpl/test/NightlyRun/test339.py =================================================================== --- /issm/trunk-jpl/test/NightlyRun/test339.py (revision 18002) +++ /issm/trunk-jpl/test/NightlyRun/test339.py (revision 18002) @@ -0,0 +1,77 @@ +import numpy +from model import * +from triangle import * +from setmask import * +from parameterize import * +from setflowequation import * +from EnumDefinitions import * +from solve import * +from MatlabFuncs import * + +md=triangle(model(),'../Exp/Square.exp',150000.) +md=setmask(md,'all','') +md=parameterize(md,'../Par/SquareShelf.py') +md.extrude(3,1.) +md=setflowequation(md,'HO','all') +md.cluster=generic('name',oshostname(),'np',3) + +md.timestepping.time_step=1. +md.settings.output_frequency=1 +md.timestepping.final_time=4. + +#Set up transient +smb=numpy.ones((md.mesh.numberofvertices,1))*3.6 +smb=numpy.hstack((smb,smb*-1.)) + +md.surfaceforcings=SMBmeltcomponents(); +md.surfaceforcings.accumulation=numpy.vstack((smb, [1.5,3.])); +md.surfaceforcings.evaporation=numpy.vstack((smb/2, [1.5,3.])); +md.surfaceforcings.melt=numpy.vstack((smb/2, [1.5,3.])); +md.surfaceforcings.refreeze=numpy.vstack((smb, [1.5,3.])); +md.transient.isthermal=False + +md=solve(md,TransientSolutionEnum()) + +#Fields and tolerances to track changes +field_names=['Vx1','Vy1','Vel1','Pressure1','Bed1','Surface1','Thickness1','SurfaceforcingsMassBalance1', \ + 'Vx2','Vy2','Vel2','Pressure2','Bed2','Surface2','Thickness2','SurfaceforcingsMassBalance2', \ + 'Vx3','Vy3','Vel3','Pressure3','Bed3','Surface3','Thickness3','SurfaceforcingsMassBalance3', \ + 'Vx4','Vy4','Vel4','Pressure4','Bed4','Surface4','Thickness4','SurfaceforcingsMassBalance4'] +field_tolerances=[1e-09,1e-09,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,\ + 1e-09,1e-09,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,\ + 1e-09,1e-09,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,\ + 1e-09,1e-09,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10] +field_values=[\ + md.results.TransientSolution[0].Vx,\ + md.results.TransientSolution[0].Vy,\ + md.results.TransientSolution[0].Vel,\ + md.results.TransientSolution[0].Pressure,\ + md.results.TransientSolution[0].Base,\ + md.results.TransientSolution[0].Surface,\ + md.results.TransientSolution[0].Thickness,\ + md.results.TransientSolution[0].SurfaceforcingsMassBalance,\ + md.results.TransientSolution[1].Vx,\ + md.results.TransientSolution[1].Vy,\ + md.results.TransientSolution[1].Vel,\ + md.results.TransientSolution[1].Pressure,\ + md.results.TransientSolution[1].Base,\ + md.results.TransientSolution[1].Surface,\ + md.results.TransientSolution[1].Thickness,\ + md.results.TransientSolution[1].SurfaceforcingsMassBalance,\ + md.results.TransientSolution[2].Vx,\ + md.results.TransientSolution[2].Vy,\ + md.results.TransientSolution[2].Vel,\ + md.results.TransientSolution[2].Pressure,\ + md.results.TransientSolution[2].Base,\ + md.results.TransientSolution[2].Surface,\ + md.results.TransientSolution[2].Thickness,\ + md.results.TransientSolution[2].SurfaceforcingsMassBalance,\ + md.results.TransientSolution[3].Vx,\ + md.results.TransientSolution[3].Vy,\ + md.results.TransientSolution[3].Vel,\ + md.results.TransientSolution[3].Pressure,\ + md.results.TransientSolution[3].Base,\ + md.results.TransientSolution[3].Surface,\ + md.results.TransientSolution[3].Thickness,\ + md.results.TransientSolution[3].SurfaceforcingsMassBalance,\ + ]