Index: /issm/trunk/test/Exp/SquareHalfRight.exp =================================================================== --- /issm/trunk/test/Exp/SquareHalfRight.exp (revision 4960) +++ /issm/trunk/test/Exp/SquareHalfRight.exp (revision 4961) @@ -4,7 +4,7 @@ 5 1. # X pos Y pos -0 500000 -1000000 500000 +500000 0 +500000 1000000 1000000 1000000 -0 1000000 -0 500000 +1000000 0 +500000 0 Index: /issm/trunk/test/NightlyRun/Id2Name.m =================================================================== --- /issm/trunk/test/NightlyRun/Id2Name.m (revision 4960) +++ /issm/trunk/test/NightlyRun/Id2Name.m (revision 4961) @@ -6,4 +6,25 @@ if (id==201), name='SquareShelfDiagM2dSerial'; -elseif (id==202), name='SquareShelf...'; +elseif (id==202), name='SquareShelfDiagM2dParallel'; +elseif (id==203), name='SquareShelfDiagM3dSerial'; +elseif (id==204), name='SquareShelfDiagM3dParallel'; +elseif (id==205), name='SquareShelfDiagP3dSerial'; +elseif (id==206), name='SquareShelfDiagP3dParallel'; +elseif (id==207), name='SquareShelfDiagS3dSerial'; +elseif (id==208), name='SquareShelfDiagS3dParallel'; +elseif (id==209), name='SquareShelfDiagMP3dSerial'; +elseif (id==210), name='SquareShelfDiagMP3dParallel'; +elseif (id==211), name='SquareShelfTherSteaSerial'; +elseif (id==212), name='SquareShelfTherSteaParallel'; +elseif (id==213), name='SquareShelfTherTranSerial'; +elseif (id==214), name='SquareShelfTherTranParallel'; +elseif (id==215), name='SquareShelfTranM2dSerial'; +elseif (id==216), name='SquareShelfTranM2dParallel'; +elseif (id==217), name='SquareShelfTranM3dSerial'; +elseif (id==218), name='SquareShelfTranM3dParallel'; +elseif (id==219), name='SquareShelfTranP3dSerial'; +elseif (id==220), name='SquareShelfTranP3dParallel'; +elseif (id==221), name='SquareShelfTranS3dSerial'; +elseif (id==222), name='SquareShelfTranS3dParallel'; +elseif (id==702), name='SquareShelf...'; else name='N/A'; end Index: /issm/trunk/test/NightlyRun/test201.m =================================================================== --- /issm/trunk/test/NightlyRun/test201.m (revision 4960) +++ /issm/trunk/test/NightlyRun/test201.m (revision 4961) @@ -3,3 +3,4 @@ md=parameterize(md,'../Par/SquareShelf.par'); md=setelementstype(md,'macayeal','all'); +md.cluster='none'; md=solve(md,'analysis_type',DiagnosticSolutionEnum); Index: /issm/trunk/test/NightlyRun/test202.m =================================================================== --- /issm/trunk/test/NightlyRun/test202.m (revision 4961) +++ /issm/trunk/test/NightlyRun/test202.m (revision 4961) @@ -0,0 +1,6 @@ +md=mesh(model,'../Exp/Square.exp',150000); +md=geography(md,'all',''); +md=parameterize(md,'../Par/SquareShelf.par'); +md=setelementstype(md,'macayeal','all'); +md.cluster=oshostname; +md=solve(md,'analysis_type',DiagnosticSolutionEnum); Index: /issm/trunk/test/NightlyRun/test202_nightly.m =================================================================== --- /issm/trunk/test/NightlyRun/test202_nightly.m (revision 4961) +++ /issm/trunk/test/NightlyRun/test202_nightly.m (revision 4961) @@ -0,0 +1,8 @@ +field_names ={'Vx','Vy','Vel','Pressure'}; +field_tolerances={1e-13,1e-13,1e-13,1e-13}; +field_values={... + PatchToVec(md.results.DiagnosticSolution.Vx),... + PatchToVec(md.results.DiagnosticSolution.Vy),... + PatchToVec(md.results.DiagnosticSolution.Vel),... + PatchToVec(md.results.DiagnosticSolution.Pressure),... + }; Index: /issm/trunk/test/NightlyRun/test203.m =================================================================== --- /issm/trunk/test/NightlyRun/test203.m (revision 4961) +++ /issm/trunk/test/NightlyRun/test203.m (revision 4961) @@ -0,0 +1,7 @@ +md=mesh(model,'../Exp/Square.exp',180000); +md=geography(md,'all',''); +md=parameterize(md,'../Par/SquareShelf.par'); +md=extrude(md,3,2); +md=setelementstype(md,'macayeal','all'); +md.cluster='none'; +md=solve(md,'analysis_type',DiagnosticSolutionEnum); Index: /issm/trunk/test/NightlyRun/test203_nightly.m =================================================================== --- /issm/trunk/test/NightlyRun/test203_nightly.m (revision 4961) +++ /issm/trunk/test/NightlyRun/test203_nightly.m (revision 4961) @@ -0,0 +1,9 @@ +field_names ={'Vx','Vy','Vz','Vel','Pressure'}; +field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13}; +field_values={... + PatchToVec(md.results.DiagnosticSolution.Vx),... + PatchToVec(md.results.DiagnosticSolution.Vy),... + PatchToVec(md.results.DiagnosticSolution.Vz),... + PatchToVec(md.results.DiagnosticSolution.Vel),... + PatchToVec(md.results.DiagnosticSolution.Pressure),... + }; Index: /issm/trunk/test/NightlyRun/test204.m =================================================================== --- /issm/trunk/test/NightlyRun/test204.m (revision 4961) +++ /issm/trunk/test/NightlyRun/test204.m (revision 4961) @@ -0,0 +1,7 @@ +md=mesh(model,'../Exp/Square.exp',180000); +md=geography(md,'all',''); +md=parameterize(md,'../Par/SquareShelf.par'); +md=extrude(md,3,2); +md=setelementstype(md,'macayeal','all'); +md.cluster=oshostname; +md=solve(md,'analysis_type',DiagnosticSolutionEnum); Index: /issm/trunk/test/NightlyRun/test204_nightly.m =================================================================== --- /issm/trunk/test/NightlyRun/test204_nightly.m (revision 4961) +++ /issm/trunk/test/NightlyRun/test204_nightly.m (revision 4961) @@ -0,0 +1,9 @@ +field_names ={'Vx','Vy','Vz','Vel','Pressure'}; +field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13}; +field_values={... + PatchToVec(md.results.DiagnosticSolution.Vx),... + PatchToVec(md.results.DiagnosticSolution.Vy),... + PatchToVec(md.results.DiagnosticSolution.Vz),... + PatchToVec(md.results.DiagnosticSolution.Vel),... + PatchToVec(md.results.DiagnosticSolution.Pressure),... + }; Index: /issm/trunk/test/NightlyRun/test205.m =================================================================== --- /issm/trunk/test/NightlyRun/test205.m (revision 4961) +++ /issm/trunk/test/NightlyRun/test205.m (revision 4961) @@ -0,0 +1,8 @@ +md=mesh(model,'../Exp/Square.exp',180000); +md=geography(md,'all',''); +md=parameterize(md,'../Par/SquareShelf.par'); +md=extrude(md,3,2); +md=setelementstype(md,'pattyn','all'); +md.verbose=2; +md.cluster='none'; +md=solve(md,'analysis_type',DiagnosticSolutionEnum); Index: /issm/trunk/test/NightlyRun/test205_nightly.m =================================================================== --- /issm/trunk/test/NightlyRun/test205_nightly.m (revision 4961) +++ /issm/trunk/test/NightlyRun/test205_nightly.m (revision 4961) @@ -0,0 +1,9 @@ +field_names ={'Vx','Vy','Vz','Vel','Pressure'}; +field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13}; +field_values={... + PatchToVec(md.results.DiagnosticSolution.Vx),... + PatchToVec(md.results.DiagnosticSolution.Vy),... + PatchToVec(md.results.DiagnosticSolution.Vz),... + PatchToVec(md.results.DiagnosticSolution.Vel),... + PatchToVec(md.results.DiagnosticSolution.Pressure),... + }; Index: /issm/trunk/test/NightlyRun/test206.m =================================================================== --- /issm/trunk/test/NightlyRun/test206.m (revision 4961) +++ /issm/trunk/test/NightlyRun/test206.m (revision 4961) @@ -0,0 +1,7 @@ +md=mesh(model,'../Exp/Square.exp',180000); +md=geography(md,'all',''); +md=parameterize(md,'../Par/SquareShelf.par'); +md=extrude(md,3,2); +md=setelementstype(md,'pattyn','all'); +md.cluster=oshostname; +md=solve(md,'analysis_type',DiagnosticSolutionEnum); Index: /issm/trunk/test/NightlyRun/test206_nightly.m =================================================================== --- /issm/trunk/test/NightlyRun/test206_nightly.m (revision 4961) +++ /issm/trunk/test/NightlyRun/test206_nightly.m (revision 4961) @@ -0,0 +1,9 @@ +field_names ={'Vx','Vy','Vz','Vel','Pressure'}; +field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13}; +field_values={... + PatchToVec(md.results.DiagnosticSolution.Vx),... + PatchToVec(md.results.DiagnosticSolution.Vy),... + PatchToVec(md.results.DiagnosticSolution.Vz),... + PatchToVec(md.results.DiagnosticSolution.Vel),... + PatchToVec(md.results.DiagnosticSolution.Pressure),... + }; Index: /issm/trunk/test/NightlyRun/test207.m =================================================================== --- /issm/trunk/test/NightlyRun/test207.m (revision 4961) +++ /issm/trunk/test/NightlyRun/test207.m (revision 4961) @@ -0,0 +1,7 @@ +md=mesh(model,'../Exp/Square.exp',180000); +md=geography(md,'all',''); +md=parameterize(md,'../Par/SquareShelf.par'); +md=extrude(md,3,2); +md=setelementstype(md,'macayeal','all','stokes','all'); +md.cluster='none'; +md=solve(md,'analysis_type',DiagnosticSolutionEnum); Index: /issm/trunk/test/NightlyRun/test207_nightly.m =================================================================== --- /issm/trunk/test/NightlyRun/test207_nightly.m (revision 4961) +++ /issm/trunk/test/NightlyRun/test207_nightly.m (revision 4961) @@ -0,0 +1,9 @@ +field_names ={'Vx','Vy','Vz','Vel','Pressure'}; +field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13}; +field_values={... + PatchToVec(md.results.DiagnosticSolution.Vx),... + PatchToVec(md.results.DiagnosticSolution.Vy),... + PatchToVec(md.results.DiagnosticSolution.Vz),... + PatchToVec(md.results.DiagnosticSolution.Vel),... + PatchToVec(md.results.DiagnosticSolution.Pressure),... + }; Index: /issm/trunk/test/NightlyRun/test208.m =================================================================== --- /issm/trunk/test/NightlyRun/test208.m (revision 4961) +++ /issm/trunk/test/NightlyRun/test208.m (revision 4961) @@ -0,0 +1,7 @@ +md=mesh(model,'../Exp/Square.exp',180000); +md=geography(md,'all',''); +md=parameterize(md,'../Par/SquareShelf.par'); +md=extrude(md,3,2); +md=setelementstype(md,'macayeal','all','stokes','all'); +md.cluster=oshostname; +md=solve(md,'analysis_type',DiagnosticSolutionEnum); Index: /issm/trunk/test/NightlyRun/test208_nightly.m =================================================================== --- /issm/trunk/test/NightlyRun/test208_nightly.m (revision 4961) +++ /issm/trunk/test/NightlyRun/test208_nightly.m (revision 4961) @@ -0,0 +1,9 @@ +field_names ={'Vx','Vy','Vz','Vel','Pressure'}; +field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13}; +field_values={... + PatchToVec(md.results.DiagnosticSolution.Vx),... + PatchToVec(md.results.DiagnosticSolution.Vy),... + PatchToVec(md.results.DiagnosticSolution.Vz),... + PatchToVec(md.results.DiagnosticSolution.Vel),... + PatchToVec(md.results.DiagnosticSolution.Pressure),... + }; Index: /issm/trunk/test/NightlyRun/test209.m =================================================================== --- /issm/trunk/test/NightlyRun/test209.m (revision 4961) +++ /issm/trunk/test/NightlyRun/test209.m (revision 4961) @@ -0,0 +1,7 @@ +md=mesh(model,'../Exp/Square.exp',150000); +md=geography(md,'all',''); +md=parameterize(md,'../Par/SquareShelf.par'); +md=extrude(md,3,2); +md=setelementstype(md,'pattyn','../Exp/SquareHalfRight.exp','fill','macayeal'); +md.cluster='none'; +md=solve(md,'analysis_type',DiagnosticSolutionEnum); Index: /issm/trunk/test/NightlyRun/test209_nightly.m =================================================================== --- /issm/trunk/test/NightlyRun/test209_nightly.m (revision 4961) +++ /issm/trunk/test/NightlyRun/test209_nightly.m (revision 4961) @@ -0,0 +1,9 @@ +field_names ={'Vx','Vy','Vz','Vel','Pressure'}; +field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13}; +field_values={... + PatchToVec(md.results.DiagnosticSolution.Vx),... + PatchToVec(md.results.DiagnosticSolution.Vy),... + PatchToVec(md.results.DiagnosticSolution.Vz),... + PatchToVec(md.results.DiagnosticSolution.Vel),... + PatchToVec(md.results.DiagnosticSolution.Pressure),... + }; Index: /issm/trunk/test/NightlyRun/test210.m =================================================================== --- /issm/trunk/test/NightlyRun/test210.m (revision 4961) +++ /issm/trunk/test/NightlyRun/test210.m (revision 4961) @@ -0,0 +1,7 @@ +md=mesh(model,'../Exp/Square.exp',150000); +md=geography(md,'all',''); +md=parameterize(md,'../Par/SquareShelf.par'); +md=extrude(md,3,2); +md=setelementstype(md,'pattyn','../Exp/SquareHalfRight.exp','fill','macayeal'); +md.cluster=oshostname; +md=solve(md,'analysis_type',DiagnosticSolutionEnum); Index: /issm/trunk/test/NightlyRun/test210_nightly.m =================================================================== --- /issm/trunk/test/NightlyRun/test210_nightly.m (revision 4961) +++ /issm/trunk/test/NightlyRun/test210_nightly.m (revision 4961) @@ -0,0 +1,9 @@ +field_names ={'Vx','Vy','Vz','Vel','Pressure'}; +field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13}; +field_values={... + PatchToVec(md.results.DiagnosticSolution.Vx),... + PatchToVec(md.results.DiagnosticSolution.Vy),... + PatchToVec(md.results.DiagnosticSolution.Vz),... + PatchToVec(md.results.DiagnosticSolution.Vel),... + PatchToVec(md.results.DiagnosticSolution.Pressure),... + }; Index: /issm/trunk/test/NightlyRun/test211.m =================================================================== --- /issm/trunk/test/NightlyRun/test211.m (revision 4961) +++ /issm/trunk/test/NightlyRun/test211.m (revision 4961) @@ -0,0 +1,8 @@ +md=mesh(model,'../Exp/Square.exp',180000); +md=geography(md,'all',''); +md=parameterize(md,'../Par/SquareShelf.par'); +md=extrude(md,3,1); +md=setelementstype(md,'macayeal','all'); +md.dt=0; +md.cluster='none'; +md=solve(md,'analysis_type',ThermalSolutionEnum); Index: /issm/trunk/test/NightlyRun/test211_nighlty.m =================================================================== --- /issm/trunk/test/NightlyRun/test211_nighlty.m (revision 4961) +++ /issm/trunk/test/NightlyRun/test211_nighlty.m (revision 4961) @@ -0,0 +1,6 @@ +field_names ={'Temperature','MeltingRate'}; +field_tolerances={1e-13,1e-13}; +field_values={... + PatchToVec(md.results.ThermalSolution.Temperature),... + PatchToVec(md.results.ThermalSolution.MeltingRate),... + }; Index: /issm/trunk/test/NightlyRun/test212.m =================================================================== --- /issm/trunk/test/NightlyRun/test212.m (revision 4961) +++ /issm/trunk/test/NightlyRun/test212.m (revision 4961) @@ -0,0 +1,8 @@ +md=mesh(model,'../Exp/Square.exp',180000); +md=geography(md,'all',''); +md=parameterize(md,'../Par/SquareShelf.par'); +md=extrude(md,3,1); +md=setelementstype(md,'macayeal','all'); +md.dt=0; +md.cluster=oshostname; +md=solve(md,'analysis_type',ThermalSolutionEnum); Index: /issm/trunk/test/NightlyRun/test212_nightly.m =================================================================== --- /issm/trunk/test/NightlyRun/test212_nightly.m (revision 4961) +++ /issm/trunk/test/NightlyRun/test212_nightly.m (revision 4961) @@ -0,0 +1,6 @@ +field_names ={'Temperature','MeltingRate'}; +field_tolerances={1e-13,1e-13}; +field_values={... + PatchToVec(md.results.ThermalSolution.Temperature),... + PatchToVec(md.results.ThermalSolution.MeltingRate),... + }; Index: /issm/trunk/test/NightlyRun/test213.m =================================================================== --- /issm/trunk/test/NightlyRun/test213.m (revision 4961) +++ /issm/trunk/test/NightlyRun/test213.m (revision 4961) @@ -0,0 +1,7 @@ +md=mesh(model,'../Exp/Square.exp',180000); +md=geography(md,'all',''); +md=parameterize(md,'../Par/SquareShelf.par'); +md=extrude(md,3,1); +md=setelementstype(md,'macayeal','all'); +md.cluster='none'; +md=solve(md,'analysis_type',ThermalSolutionEnum); Index: /issm/trunk/test/NightlyRun/test213_nightly.m =================================================================== --- /issm/trunk/test/NightlyRun/test213_nightly.m (revision 4961) +++ /issm/trunk/test/NightlyRun/test213_nightly.m (revision 4961) @@ -0,0 +1,10 @@ +field_names ={'Temperature1','MeltingRate1','Temperature2','MeltingRate2','Temperature3','MeltingRate3'}; +field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13}; +field_values={... + PatchToVec(md.results.ThermalSolution(1).Temperature),... + PatchToVec(md.results.ThermalSolution(1).MeltingRate),... + PatchToVec(md.results.ThermalSolution(2).Temperature),... + PatchToVec(md.results.ThermalSolution(2).MeltingRate),... + PatchToVec(md.results.ThermalSolution(3).Temperature),... + PatchToVec(md.results.ThermalSolution(3).MeltingRate),... + }; Index: /issm/trunk/test/NightlyRun/test214.m =================================================================== --- /issm/trunk/test/NightlyRun/test214.m (revision 4961) +++ /issm/trunk/test/NightlyRun/test214.m (revision 4961) @@ -0,0 +1,7 @@ +md=mesh(model,'../Exp/Square.exp',180000); +md=geography(md,'all',''); +md=parameterize(md,'../Par/SquareShelf.par'); +md=extrude(md,3,1); +md=setelementstype(md,'macayeal','all'); +md.cluster=oshostname; +md=solve(md,'analysis_type',ThermalSolutionEnum); Index: /issm/trunk/test/NightlyRun/test214_nightly.m =================================================================== --- /issm/trunk/test/NightlyRun/test214_nightly.m (revision 4961) +++ /issm/trunk/test/NightlyRun/test214_nightly.m (revision 4961) @@ -0,0 +1,10 @@ +field_names ={'Temperature1','MeltingRate1','Temperature2','MeltingRate2','Temperature3','MeltingRate3'}; +field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13}; +field_values={... + PatchToVec(md.results.ThermalSolution(1).Temperature),... + PatchToVec(md.results.ThermalSolution(1).MeltingRate),... + PatchToVec(md.results.ThermalSolution(2).Temperature),... + PatchToVec(md.results.ThermalSolution(2).MeltingRate),... + PatchToVec(md.results.ThermalSolution(3).Temperature),... + PatchToVec(md.results.ThermalSolution(3).MeltingRate),... + }; Index: /issm/trunk/test/NightlyRun/test215.m =================================================================== --- /issm/trunk/test/NightlyRun/test215.m (revision 4961) +++ /issm/trunk/test/NightlyRun/test215.m (revision 4961) @@ -0,0 +1,6 @@ +md=mesh(model,'../Exp/Square.exp',150000); +md=geography(md,'all',''); +md=parameterize(md,'../Par/SquareShelf.par'); +md=setelementstype(md,'macayeal','all'); +md.cluster='none'; +md=solve(md,'analysis_type',Transient2DSolutionEnum); Index: /issm/trunk/test/NightlyRun/test215_nightly.m =================================================================== --- /issm/trunk/test/NightlyRun/test215_nightly.m (revision 4961) +++ /issm/trunk/test/NightlyRun/test215_nightly.m (revision 4961) @@ -0,0 +1,25 @@ +field_names ={'Vx1','Vy1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Vx2','Vy2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Vx3','Vy3','Vel3','Pressure3','Bed3','Surface3','Thickness3'}; +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}; +field_values={... + PatchToVec(md.results.Transient2DSolution(1).Vx),... + PatchToVec(md.results.Transient2DSolution(1).Vy),... + PatchToVec(md.results.Transient2DSolution(1).Vel),... + PatchToVec(md.results.Transient2DSolution(1).Pressure),... + PatchToVec(md.results.Transient2DSolution(1).Bed),... + PatchToVec(md.results.Transient2DSolution(1).Surface),... + PatchToVec(md.results.Transient2DSolution(1).Thickness),... + PatchToVec(md.results.Transient2DSolution(2).Vx),... + PatchToVec(md.results.Transient2DSolution(2).Vy),... + PatchToVec(md.results.Transient2DSolution(2).Vel),... + PatchToVec(md.results.Transient2DSolution(2).Pressure),... + PatchToVec(md.results.Transient2DSolution(2).Bed),... + PatchToVec(md.results.Transient2DSolution(2).Surface),... + PatchToVec(md.results.Transient2DSolution(2).Thickness),... + PatchToVec(md.results.Transient2DSolution(3).Vx),... + PatchToVec(md.results.Transient2DSolution(3).Vy),... + PatchToVec(md.results.Transient2DSolution(3).Vel),... + PatchToVec(md.results.Transient2DSolution(3).Pressure),... + PatchToVec(md.results.Transient2DSolution(3).Bed),... + PatchToVec(md.results.Transient2DSolution(3).Surface),... + PatchToVec(md.results.Transient2DSolution(3).Thickness),... + }; Index: /issm/trunk/test/NightlyRun/test216.m =================================================================== --- /issm/trunk/test/NightlyRun/test216.m (revision 4961) +++ /issm/trunk/test/NightlyRun/test216.m (revision 4961) @@ -0,0 +1,6 @@ +md=mesh(model,'../Exp/Square.exp',150000); +md=geography(md,'all',''); +md=parameterize(md,'../Par/SquareShelf.par'); +md=setelementstype(md,'macayeal','all'); +md.cluster=oshostname; +md=solve(md,'analysis_type',Transient2DSolutionEnum); Index: /issm/trunk/test/NightlyRun/test216_nightly.m =================================================================== --- /issm/trunk/test/NightlyRun/test216_nightly.m (revision 4961) +++ /issm/trunk/test/NightlyRun/test216_nightly.m (revision 4961) @@ -0,0 +1,25 @@ +field_names ={'Vx1','Vy1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Vx2','Vy2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Vx3','Vy3','Vel3','Pressure3','Bed3','Surface3','Thickness3'}; +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}; +field_values={... + PatchToVec(md.results.Transient2DSolution(1).Vx),... + PatchToVec(md.results.Transient2DSolution(1).Vy),... + PatchToVec(md.results.Transient2DSolution(1).Vel),... + PatchToVec(md.results.Transient2DSolution(1).Pressure),... + PatchToVec(md.results.Transient2DSolution(1).Bed),... + PatchToVec(md.results.Transient2DSolution(1).Surface),... + PatchToVec(md.results.Transient2DSolution(1).Thickness),... + PatchToVec(md.results.Transient2DSolution(2).Vx),... + PatchToVec(md.results.Transient2DSolution(2).Vy),... + PatchToVec(md.results.Transient2DSolution(2).Vel),... + PatchToVec(md.results.Transient2DSolution(2).Pressure),... + PatchToVec(md.results.Transient2DSolution(2).Bed),... + PatchToVec(md.results.Transient2DSolution(2).Surface),... + PatchToVec(md.results.Transient2DSolution(2).Thickness),... + PatchToVec(md.results.Transient2DSolution(3).Vx),... + PatchToVec(md.results.Transient2DSolution(3).Vy),... + PatchToVec(md.results.Transient2DSolution(3).Vel),... + PatchToVec(md.results.Transient2DSolution(3).Pressure),... + PatchToVec(md.results.Transient2DSolution(3).Bed),... + PatchToVec(md.results.Transient2DSolution(3).Surface),... + PatchToVec(md.results.Transient2DSolution(3).Thickness),... + }; Index: /issm/trunk/test/NightlyRun/test217.m =================================================================== --- /issm/trunk/test/NightlyRun/test217.m (revision 4961) +++ /issm/trunk/test/NightlyRun/test217.m (revision 4961) @@ -0,0 +1,7 @@ +md=mesh(model,'../Exp/Square.exp',180000); +md=geography(md,'all',''); +md=parameterize(md,'../Par/SquareShelf.par'); +md=extrude(md,3,1); +md=setelementstype(md,'macayeal','all'); +md.cluster='none'; +md=solve(md,'analysis_type',Transient3DSolutionEnum); Index: /issm/trunk/test/NightlyRun/test217_nightly.m =================================================================== --- /issm/trunk/test/NightlyRun/test217_nightly.m (revision 4961) +++ /issm/trunk/test/NightlyRun/test217_nightly.m (revision 4961) @@ -0,0 +1,33 @@ +field_names ={'Vx1','Vy1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1','MeltingRate1', ... + 'Vx2','Vy2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2','MeltingRate2', ... + 'Vx3','Vy3','Vel3','Pressure3','Bed3','Surface3','Thickness3','Temperature3','MeltingRate3'}; +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}; +field_values={... + PatchToVec(md.results.Transient2DSolution(1).Vx),... + PatchToVec(md.results.Transient2DSolution(1).Vy),... + PatchToVec(md.results.Transient2DSolution(1).Vel),... + PatchToVec(md.results.Transient2DSolution(1).Pressure),... + PatchToVec(md.results.Transient2DSolution(1).Bed),... + PatchToVec(md.results.Transient2DSolution(1).Surface),... + PatchToVec(md.results.Transient2DSolution(1).Thickness),... + PatchToVec(md.results.Transient2DSolution(1).Temperature),... + PatchToVec(md.results.Transient2DSolution(1).MeltingRate),... + PatchToVec(md.results.Transient2DSolution(2).Vx),... + PatchToVec(md.results.Transient2DSolution(2).Vy),... + PatchToVec(md.results.Transient2DSolution(2).Vel),... + PatchToVec(md.results.Transient2DSolution(2).Pressure),... + PatchToVec(md.results.Transient2DSolution(2).Bed),... + PatchToVec(md.results.Transient2DSolution(2).Surface),... + PatchToVec(md.results.Transient2DSolution(2).Thickness),... + PatchToVec(md.results.Transient2DSolution(2).Temperature),... + PatchToVec(md.results.Transient2DSolution(2).MeltingRate),... + PatchToVec(md.results.Transient2DSolution(3).Vx),... + PatchToVec(md.results.Transient2DSolution(3).Vy),... + PatchToVec(md.results.Transient2DSolution(3).Vel),... + PatchToVec(md.results.Transient2DSolution(3).Pressure),... + PatchToVec(md.results.Transient2DSolution(3).Bed),... + PatchToVec(md.results.Transient2DSolution(3).Surface),... + PatchToVec(md.results.Transient2DSolution(3).Thickness),... + PatchToVec(md.results.Transient2DSolution(3).Temperature),... + PatchToVec(md.results.Transient2DSolution(3).MeltingRate),... + }; Index: /issm/trunk/test/NightlyRun/test218.m =================================================================== --- /issm/trunk/test/NightlyRun/test218.m (revision 4961) +++ /issm/trunk/test/NightlyRun/test218.m (revision 4961) @@ -0,0 +1,7 @@ +md=mesh(model,'../Exp/Square.exp',180000); +md=geography(md,'all',''); +md=parameterize(md,'../Par/SquareShelf.par'); +md=extrude(md,3,1); +md=setelementstype(md,'macayeal','all'); +md.cluster=oshostname; +md=solve(md,'analysis_type',Transient3DSolutionEnum); Index: /issm/trunk/test/NightlyRun/test218_nightly.m =================================================================== --- /issm/trunk/test/NightlyRun/test218_nightly.m (revision 4961) +++ /issm/trunk/test/NightlyRun/test218_nightly.m (revision 4961) @@ -0,0 +1,33 @@ +field_names ={'Vx1','Vy1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1','MeltingRate1', ... + 'Vx2','Vy2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2','MeltingRate2', ... + 'Vx3','Vy3','Vel3','Pressure3','Bed3','Surface3','Thickness3','Temperature3','MeltingRate3'}; +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}; +field_values={... + PatchToVec(md.results.Transient2DSolution(1).Vx),... + PatchToVec(md.results.Transient2DSolution(1).Vy),... + PatchToVec(md.results.Transient2DSolution(1).Vel),... + PatchToVec(md.results.Transient2DSolution(1).Pressure),... + PatchToVec(md.results.Transient2DSolution(1).Bed),... + PatchToVec(md.results.Transient2DSolution(1).Surface),... + PatchToVec(md.results.Transient2DSolution(1).Thickness),... + PatchToVec(md.results.Transient2DSolution(1).Temperature),... + PatchToVec(md.results.Transient2DSolution(1).MeltingRate),... + PatchToVec(md.results.Transient2DSolution(2).Vx),... + PatchToVec(md.results.Transient2DSolution(2).Vy),... + PatchToVec(md.results.Transient2DSolution(2).Vel),... + PatchToVec(md.results.Transient2DSolution(2).Pressure),... + PatchToVec(md.results.Transient2DSolution(2).Bed),... + PatchToVec(md.results.Transient2DSolution(2).Surface),... + PatchToVec(md.results.Transient2DSolution(2).Thickness),... + PatchToVec(md.results.Transient2DSolution(2).Temperature),... + PatchToVec(md.results.Transient2DSolution(2).MeltingRate),... + PatchToVec(md.results.Transient2DSolution(3).Vx),... + PatchToVec(md.results.Transient2DSolution(3).Vy),... + PatchToVec(md.results.Transient2DSolution(3).Vel),... + PatchToVec(md.results.Transient2DSolution(3).Pressure),... + PatchToVec(md.results.Transient2DSolution(3).Bed),... + PatchToVec(md.results.Transient2DSolution(3).Surface),... + PatchToVec(md.results.Transient2DSolution(3).Thickness),... + PatchToVec(md.results.Transient2DSolution(3).Temperature),... + PatchToVec(md.results.Transient2DSolution(3).MeltingRate),... + }; Index: /issm/trunk/test/NightlyRun/test219.m =================================================================== --- /issm/trunk/test/NightlyRun/test219.m (revision 4961) +++ /issm/trunk/test/NightlyRun/test219.m (revision 4961) @@ -0,0 +1,7 @@ +md=mesh(model,'../Exp/Square.exp',200000); +md=geography(md,'all',''); +md=parameterize(md,'../Par/SquareShelf.par'); +md=extrude(md,3,1); +md=setelementstype(md,'pattyn','all'); +md.cluster='none'; +md=solve(md,'analysis_type',Transient3DSolutionEnum); Index: /issm/trunk/test/NightlyRun/test219_nightly.m =================================================================== --- /issm/trunk/test/NightlyRun/test219_nightly.m (revision 4961) +++ /issm/trunk/test/NightlyRun/test219_nightly.m (revision 4961) @@ -0,0 +1,33 @@ +field_names ={'Vx1','Vy1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1','MeltingRate1', ... + 'Vx2','Vy2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2','MeltingRate2', ... + 'Vx3','Vy3','Vel3','Pressure3','Bed3','Surface3','Thickness3','Temperature3','MeltingRate3'}; +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}; +field_values={... + PatchToVec(md.results.Transient2DSolution(1).Vx),... + PatchToVec(md.results.Transient2DSolution(1).Vy),... + PatchToVec(md.results.Transient2DSolution(1).Vel),... + PatchToVec(md.results.Transient2DSolution(1).Pressure),... + PatchToVec(md.results.Transient2DSolution(1).Bed),... + PatchToVec(md.results.Transient2DSolution(1).Surface),... + PatchToVec(md.results.Transient2DSolution(1).Thickness),... + PatchToVec(md.results.Transient2DSolution(1).Temperature),... + PatchToVec(md.results.Transient2DSolution(1).MeltingRate),... + PatchToVec(md.results.Transient2DSolution(2).Vx),... + PatchToVec(md.results.Transient2DSolution(2).Vy),... + PatchToVec(md.results.Transient2DSolution(2).Vel),... + PatchToVec(md.results.Transient2DSolution(2).Pressure),... + PatchToVec(md.results.Transient2DSolution(2).Bed),... + PatchToVec(md.results.Transient2DSolution(2).Surface),... + PatchToVec(md.results.Transient2DSolution(2).Thickness),... + PatchToVec(md.results.Transient2DSolution(2).Temperature),... + PatchToVec(md.results.Transient2DSolution(2).MeltingRate),... + PatchToVec(md.results.Transient2DSolution(3).Vx),... + PatchToVec(md.results.Transient2DSolution(3).Vy),... + PatchToVec(md.results.Transient2DSolution(3).Vel),... + PatchToVec(md.results.Transient2DSolution(3).Pressure),... + PatchToVec(md.results.Transient2DSolution(3).Bed),... + PatchToVec(md.results.Transient2DSolution(3).Surface),... + PatchToVec(md.results.Transient2DSolution(3).Thickness),... + PatchToVec(md.results.Transient2DSolution(3).Temperature),... + PatchToVec(md.results.Transient2DSolution(3).MeltingRate),... + }; Index: /issm/trunk/test/NightlyRun/test220.m =================================================================== --- /issm/trunk/test/NightlyRun/test220.m (revision 4961) +++ /issm/trunk/test/NightlyRun/test220.m (revision 4961) @@ -0,0 +1,7 @@ +md=mesh(model,'../Exp/Square.exp',200000); +md=geography(md,'all',''); +md=parameterize(md,'../Par/SquareShelf.par'); +md=extrude(md,3,1); +md=setelementstype(md,'pattyn','all'); +md.cluster=oshostname; +md=solve(md,'analysis_type',Transient3DSolutionEnum); Index: /issm/trunk/test/NightlyRun/test220_nightly.m =================================================================== --- /issm/trunk/test/NightlyRun/test220_nightly.m (revision 4961) +++ /issm/trunk/test/NightlyRun/test220_nightly.m (revision 4961) @@ -0,0 +1,33 @@ +field_names ={'Vx1','Vy1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1','MeltingRate1', ... + 'Vx2','Vy2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2','MeltingRate2', ... + 'Vx3','Vy3','Vel3','Pressure3','Bed3','Surface3','Thickness3','Temperature3','MeltingRate3'}; +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}; +field_values={... + PatchToVec(md.results.Transient2DSolution(1).Vx),... + PatchToVec(md.results.Transient2DSolution(1).Vy),... + PatchToVec(md.results.Transient2DSolution(1).Vel),... + PatchToVec(md.results.Transient2DSolution(1).Pressure),... + PatchToVec(md.results.Transient2DSolution(1).Bed),... + PatchToVec(md.results.Transient2DSolution(1).Surface),... + PatchToVec(md.results.Transient2DSolution(1).Thickness),... + PatchToVec(md.results.Transient2DSolution(1).Temperature),... + PatchToVec(md.results.Transient2DSolution(1).MeltingRate),... + PatchToVec(md.results.Transient2DSolution(2).Vx),... + PatchToVec(md.results.Transient2DSolution(2).Vy),... + PatchToVec(md.results.Transient2DSolution(2).Vel),... + PatchToVec(md.results.Transient2DSolution(2).Pressure),... + PatchToVec(md.results.Transient2DSolution(2).Bed),... + PatchToVec(md.results.Transient2DSolution(2).Surface),... + PatchToVec(md.results.Transient2DSolution(2).Thickness),... + PatchToVec(md.results.Transient2DSolution(2).Temperature),... + PatchToVec(md.results.Transient2DSolution(2).MeltingRate),... + PatchToVec(md.results.Transient2DSolution(3).Vx),... + PatchToVec(md.results.Transient2DSolution(3).Vy),... + PatchToVec(md.results.Transient2DSolution(3).Vel),... + PatchToVec(md.results.Transient2DSolution(3).Pressure),... + PatchToVec(md.results.Transient2DSolution(3).Bed),... + PatchToVec(md.results.Transient2DSolution(3).Surface),... + PatchToVec(md.results.Transient2DSolution(3).Thickness),... + PatchToVec(md.results.Transient2DSolution(3).Temperature),... + PatchToVec(md.results.Transient2DSolution(3).MeltingRate),... + }; Index: /issm/trunk/test/NightlyRun/test221.m =================================================================== --- /issm/trunk/test/NightlyRun/test221.m (revision 4961) +++ /issm/trunk/test/NightlyRun/test221.m (revision 4961) @@ -0,0 +1,8 @@ +md=mesh(model,'../Exp/Square.exp',200000); +md=geography(md,'all',''); +md=parameterize(md,'../Par/SquareShelf.par'); +md=extrude(md,3,1); +md=setelementstype(md,'macayeal','all','stokes','all'); +md.eps_rel=NaN; +md.cluster='none'; +md=solve(md,'analysis_type',Transient3DSolutionEnum); Index: /issm/trunk/test/NightlyRun/test221_nightly.m =================================================================== --- /issm/trunk/test/NightlyRun/test221_nightly.m (revision 4961) +++ /issm/trunk/test/NightlyRun/test221_nightly.m (revision 4961) @@ -0,0 +1,33 @@ +field_names ={'Vx1','Vy1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1','MeltingRate1', ... + 'Vx2','Vy2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2','MeltingRate2', ... + 'Vx3','Vy3','Vel3','Pressure3','Bed3','Surface3','Thickness3','Temperature3','MeltingRate3'}; +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}; +field_values={... + PatchToVec(md.results.Transient2DSolution(1).Vx),... + PatchToVec(md.results.Transient2DSolution(1).Vy),... + PatchToVec(md.results.Transient2DSolution(1).Vel),... + PatchToVec(md.results.Transient2DSolution(1).Pressure),... + PatchToVec(md.results.Transient2DSolution(1).Bed),... + PatchToVec(md.results.Transient2DSolution(1).Surface),... + PatchToVec(md.results.Transient2DSolution(1).Thickness),... + PatchToVec(md.results.Transient2DSolution(1).Temperature),... + PatchToVec(md.results.Transient2DSolution(1).MeltingRate),... + PatchToVec(md.results.Transient2DSolution(2).Vx),... + PatchToVec(md.results.Transient2DSolution(2).Vy),... + PatchToVec(md.results.Transient2DSolution(2).Vel),... + PatchToVec(md.results.Transient2DSolution(2).Pressure),... + PatchToVec(md.results.Transient2DSolution(2).Bed),... + PatchToVec(md.results.Transient2DSolution(2).Surface),... + PatchToVec(md.results.Transient2DSolution(2).Thickness),... + PatchToVec(md.results.Transient2DSolution(2).Temperature),... + PatchToVec(md.results.Transient2DSolution(2).MeltingRate),... + PatchToVec(md.results.Transient2DSolution(3).Vx),... + PatchToVec(md.results.Transient2DSolution(3).Vy),... + PatchToVec(md.results.Transient2DSolution(3).Vel),... + PatchToVec(md.results.Transient2DSolution(3).Pressure),... + PatchToVec(md.results.Transient2DSolution(3).Bed),... + PatchToVec(md.results.Transient2DSolution(3).Surface),... + PatchToVec(md.results.Transient2DSolution(3).Thickness),... + PatchToVec(md.results.Transient2DSolution(3).Temperature),... + PatchToVec(md.results.Transient2DSolution(3).MeltingRate),... + }; Index: /issm/trunk/test/NightlyRun/test222.m =================================================================== --- /issm/trunk/test/NightlyRun/test222.m (revision 4961) +++ /issm/trunk/test/NightlyRun/test222.m (revision 4961) @@ -0,0 +1,8 @@ +md=mesh(model,'../Exp/Square.exp',200000); +md=geography(md,'all',''); +md=parameterize(md,'../Par/SquareShelf.par'); +md=extrude(md,3,1); +md=setelementstype(md,'macayeal','all','stokes','all'); +md.eps_rel=NaN; +md.cluster=oshostname; +md=solve(md,'analysis_type',Transient3DSolutionEnum); Index: /issm/trunk/test/NightlyRun/test222_nightly.m =================================================================== --- /issm/trunk/test/NightlyRun/test222_nightly.m (revision 4961) +++ /issm/trunk/test/NightlyRun/test222_nightly.m (revision 4961) @@ -0,0 +1,33 @@ +field_names ={'Vx1','Vy1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1','MeltingRate1', ... + 'Vx2','Vy2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2','MeltingRate2', ... + 'Vx3','Vy3','Vel3','Pressure3','Bed3','Surface3','Thickness3','Temperature3','MeltingRate3'}; +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}; +field_values={... + PatchToVec(md.results.Transient2DSolution(1).Vx),... + PatchToVec(md.results.Transient2DSolution(1).Vy),... + PatchToVec(md.results.Transient2DSolution(1).Vel),... + PatchToVec(md.results.Transient2DSolution(1).Pressure),... + PatchToVec(md.results.Transient2DSolution(1).Bed),... + PatchToVec(md.results.Transient2DSolution(1).Surface),... + PatchToVec(md.results.Transient2DSolution(1).Thickness),... + PatchToVec(md.results.Transient2DSolution(1).Temperature),... + PatchToVec(md.results.Transient2DSolution(1).MeltingRate),... + PatchToVec(md.results.Transient2DSolution(2).Vx),... + PatchToVec(md.results.Transient2DSolution(2).Vy),... + PatchToVec(md.results.Transient2DSolution(2).Vel),... + PatchToVec(md.results.Transient2DSolution(2).Pressure),... + PatchToVec(md.results.Transient2DSolution(2).Bed),... + PatchToVec(md.results.Transient2DSolution(2).Surface),... + PatchToVec(md.results.Transient2DSolution(2).Thickness),... + PatchToVec(md.results.Transient2DSolution(2).Temperature),... + PatchToVec(md.results.Transient2DSolution(2).MeltingRate),... + PatchToVec(md.results.Transient2DSolution(3).Vx),... + PatchToVec(md.results.Transient2DSolution(3).Vy),... + PatchToVec(md.results.Transient2DSolution(3).Vel),... + PatchToVec(md.results.Transient2DSolution(3).Pressure),... + PatchToVec(md.results.Transient2DSolution(3).Bed),... + PatchToVec(md.results.Transient2DSolution(3).Surface),... + PatchToVec(md.results.Transient2DSolution(3).Thickness),... + PatchToVec(md.results.Transient2DSolution(3).Temperature),... + PatchToVec(md.results.Transient2DSolution(3).MeltingRate),... + }; Index: /issm/trunk/test/Par/SquareSheetConstrained.par =================================================================== --- /issm/trunk/test/Par/SquareSheetConstrained.par (revision 4960) +++ /issm/trunk/test/Par/SquareSheetConstrained.par (revision 4961) @@ -15,4 +15,6 @@ md.vy=InterpFromMeshToMesh2d(index,x,y,vy,md.x,md.y,0); clear vx vy x y index; +md.vz=zeros(md.numberofgrids,1); +md.pressure=zeros(md.numberofgrids,1); %Materials @@ -20,4 +22,5 @@ md.rheology_B=paterson(md.observed_temperature); md.rheology_n=3*ones(md.numberofelements,1); +md.temperature=md.observed_temperature; %Friction @@ -38,4 +41,6 @@ md.eps_rel=0.005; md.eps_abs=NaN; +md.dt=1; +md.ndt=3; %Boundary conditions: Index: /issm/trunk/test/Par/SquareShelf.par =================================================================== --- /issm/trunk/test/Par/SquareShelf.par (revision 4960) +++ /issm/trunk/test/Par/SquareShelf.par (revision 4961) @@ -10,9 +10,11 @@ md.surface=md.bed+md.thickness; -%Initial velocity +%Initial velocity and pressure load('./../Data/SquareShelf.data','-mat'); md.vx=InterpFromMeshToMesh2d(index,x,y,vx,md.x,md.y,0); md.vy=InterpFromMeshToMesh2d(index,x,y,vy,md.x,md.y,0); clear vx vy x y index; +md.vz=zeros(md.numberofgrids,1); +md.pressure=zeros(md.numberofgrids,1); %Materials @@ -20,4 +22,5 @@ md.rheology_B=paterson(md.observed_temperature); md.rheology_n=3*ones(md.numberofelements,1); +md.temperature=md.observed_temperature; %Friction @@ -38,4 +41,6 @@ md.eps_rel=0.02; md.eps_abs=NaN; +md.dt=1; +md.ndt=3; %Boundary conditions: Index: /issm/trunk/test/Par/SquareShelfConstrained.par =================================================================== --- /issm/trunk/test/Par/SquareShelfConstrained.par (revision 4960) +++ /issm/trunk/test/Par/SquareShelfConstrained.par (revision 4961) @@ -15,4 +15,6 @@ md.vy=InterpFromMeshToMesh2d(index,x,y,vy,md.x,md.y,0); clear vx vy x y index; +md.vz=zeros(md.numberofgrids,1); +md.pressure=zeros(md.numberofgrids,1); %Materials @@ -20,4 +22,5 @@ md.rheology_B=paterson(md.observed_temperature); md.rheology_n=3*ones(md.numberofelements,1); +md.temperature=md.observed_temperature; %Friction @@ -38,4 +41,6 @@ md.eps_rel=0.005; md.eps_abs=NaN; +md.dt=1; +md.ndt=3; %Deal with boundary conditions: