Index: /issm/trunk-jpl/jenkins/javascript/karma/scripts/m2js.sh =================================================================== --- /issm/trunk-jpl/jenkins/javascript/karma/scripts/m2js.sh (revision 21203) +++ /issm/trunk-jpl/jenkins/javascript/karma/scripts/m2js.sh (revision 21204) @@ -9,8 +9,4 @@ echo "Exiting..." exit - fi - - if [[ $OUTFILE == test*.js ]]; then - sed -i '2s/^/var md = new model();\n/' $OUTFILE # initialize the model fi @@ -62,7 +58,7 @@ - sed -i "s/\([^']\)NaN\([^']\)/\1null\2/g" $OUTFILE # NaN translates to null + #sed -i "s/\([^']\)NaN\([^']\)/\1null\2/g" $OUTFILE # NaN translates to null - sed -i 's/function\s*\(.*\)\s*=\s*\(.*\)(\(.*\)/this.\2 = function(\3)/' $OUTFILE + sed -i 's/function\s\(.*\)\s*=\s*\(.*\)(\(.*\)/this.\2 = function(\3)/' $OUTFILE sed -i 's/(self,*\(.*\))/(\1)/' $OUTFILE @@ -87,6 +83,14 @@ sed -i 's/self/this/g' $OUTFILE # self -> this + sed -i 's/md\.transient/md.trans/g' $OUTFILE # self -> this sed -i 's/function(\([^ ]*\))/function(\1) {/' $OUTFILE + + #sed -i 's/\(\d+\)\^\(\d+\)/Math.pow(\1,\2)/' $OUTFILE + sed -i 's/\([0-9]*\)\^\([0-9]*\)/Math.pow(\1,\2)/g' $OUTFILE + + #if [[ $OUTFILE == test*.js ]]; then + sed -i '2s/^/var md = new model();\n/' $OUTFILE # initialize the model + #fi #sed -i 'N;s/properties.*\n\s*\(.*\)\(\s*\)=\s*\(.*\)/this/' $OUTFILE Index: /issm/trunk-jpl/jenkins/javascript/karma/scripts/specs/issm.spec.js =================================================================== --- /issm/trunk-jpl/jenkins/javascript/karma/scripts/specs/issm.spec.js (revision 21203) +++ /issm/trunk-jpl/jenkins/javascript/karma/scripts/specs/issm.spec.js (revision 21204) @@ -1,2 +1,11 @@ +//var json:any; +it("should load a fixture", function () { + jasmine.getFixtures().fixturesPath = "base/Archives/" + var f = readFixtures("Archive101.json"); + var json = JSON.parse(f); + expect(json).toBeDefined(); + //console.log(json); +}); + describe("test101", function() { it("contains test101", function() { @@ -13,4 +22,7 @@ md.verbose.solution=2; md.verbose.convergence=0; + console.log("MESH"); + console.log(md.mesh.domaintype()); + md=solve(md,StressbalanceSolutionEnum(),'checkconsistency','no'); @@ -19,187 +31,187 @@ }) -describe("test105", function() { - it("contains test105", function() { - //Test Name: SquareShelfConstrainedMasstransp2d - var md = new model(); - triangle(md,square[0],150000.); - setmask(md,'all',''); - parameterize(md); - setflowequation(md,'SSA','all'); - //md.cluster=generic('name',oshostname(),'np',3); - md=solve(md,MasstransportSolutionEnum()); - - //Fields and tolerances to track changes - field_names =['Thickness']; - field_tolerances=[1e-13]; - field_values=[ - (md.results.MasstransportSolution[0].Thickness), - ]; - }); -}) - -describe("test110", function() { - it("contains test110", function() { - //Test Name: SquareShelfConstrainedTranSSA2d - var md = new model(); - triangle(md,square[0],150000.); - setmask(md,'all',''); - parameterize(md); - setflowequation(md,'SSA','all'); - //md.cluster=generic('name',oshostname(),'np',3); - md.trans.requested_outputs=['IceVolume']; - - md=solve(md,TransientSolutionEnum()); - - //Fields and tolerances to track changes - field_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']; - 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=[ - (md.results.TransientSolution[0](1).Vx), - (md.results.TransientSolution[0](1).Vy), - (md.results.TransientSolution[0](1).Vel), - (md.results.TransientSolution[0](1).Pressure), - (md.results.TransientSolution[0](1).Base), - (md.results.TransientSolution[0](1).Surface), - (md.results.TransientSolution[0](1).Thickness), - (md.results.TransientSolution[0](1).IceVolume), - (md.results.TransientSolution[0](2).Vx), - (md.results.TransientSolution[0](2).Vy), - (md.results.TransientSolution[0](2).Vel), - (md.results.TransientSolution[0](2).Pressure), - (md.results.TransientSolution[0](2).Base), - (md.results.TransientSolution[0](2).Surface), - (md.results.TransientSolution[0](2).Thickness), - (md.results.TransientSolution[0](2).IceVolume), - (md.results.TransientSolution[0](3).Vx), - (md.results.TransientSolution[0](3).Vy), - (md.results.TransientSolution[0](3).Vel), - (md.results.TransientSolution[0](3).Pressure), - (md.results.TransientSolution[0](3).Base), - (md.results.TransientSolution[0](3).Surface), - (md.results.TransientSolution[0](3).Thickness), - (md.results.TransientSolution[0](3).IceVolume), - ]; - }); -}) - -describe("test112", function() { - it("contains test112", function() { - //Test Name: SquareShelfConstrainedSurfSlop2d - var md = new model(); - triangle(md,square[0],150000.); - setmask(md,'all',''); - parameterize(md); - setflowequation(md,'SSA','all'); - //md.cluster=generic('name',oshostname(),'np',3); - md=solve(md,SurfaceSlopeSolutionEnum()); - - //Fields and tolerances to track changes - field_names =['SurfaceSlopeX','SurfaceSlopeY']; - field_tolerances=[1e-13,1e-13]; - field_values=[ - (md.results.SurfaceSlopeSolution[0].SurfaceSlopeX), - (md.results.SurfaceSlopeSolution[0].SurfaceSlopeY), - ]; - }); -}) - -describe("test114", function() { - it("contains test114", function() { - //Test Name: SquareShelfConstrainedBedSlop2d - var md = new model(); - triangle(md,square[0],150000.); - setmask(md,'all',''); - parameterize(md); - setflowequation(md,'SSA','all'); - //md.cluster=generic('name',oshostname(),'np',3); - md=solve(md,BedSlopeSolutionEnum()); - - //Fields and tolerances to track changes - field_names =['BedSlopeX','BedSlopeY']; - field_tolerances=[1e-13,1e-13]; - field_values=[ - (md.results.BedSlopeSolution[0].BedSlopeX), - (md.results.BedSlopeSolution[0].BedSlopeY), - ]; - }); -}) - -describe("test201", function() { - it("contains test201", function() { - //Test Name: SquareShelfStressSSA2d - var md = new model(); - triangle(md,square[0],150000.); - setmask(md,'all',''); - parameterize(md); - setflowequation(md,'SSA','all'); - //md.cluster=generic('name',oshostname(),'np',3); - md=solve(md,StressbalanceSolutionEnum()); - - //Fields and tolerances to track changes - field_names =['Vx','Vy','Vel','Pressure']; - field_tolerances=[1e-13,1e-13,1e-13,1e-13]; - field_values=[ - (md.results.StressbalanceSolution[0].Vx), - (md.results.StressbalanceSolution[0].Vy), - (md.results.StressbalanceSolution[0].Vel), - (md.results.StressbalanceSolution[0].Pressure), - ]; - }); -}) - -describe("test208", function() { - it("contains test208", function() { - //Test Name: SquareShelfTranSSA2d - var md = new model(); - triangle(md,square[0],150000.); - setmask(md,'all',''); - parameterize(md); - setflowequation(md,'SSA','all'); - //md.cluster=generic('name',oshostname(),'np',3); - md.trans.requested_outputs=['default','FloatingArea','GroundedArea','TotalGroundedBmb','TotalFloatingBmb']; - for (var i = 0; i < md.basalforcings.floatingice_melting_rate.length; ++i) { - md.basalforcings.floatingice_melting_rate[i] = 1; - } - md=solve(md,TransientSolutionEnum()); - - //Fields and tolerances to track changes - field_names =['Vx1','Vy1','Vel1','Pressure1','Bed1','Surface1','Thickness1','TotalGroundedBmb1','TotalFloatingBmb1', - 'Vx2','Vy2','Vel2','Pressure2','Bed2','Surface2','Thickness2','TotalGroundedBmb2','TotalFloatingBmb2', - 'Vx3','Vy3','Vel3','Pressure3','Bed3','Surface3','Thickness3','TotalGroundedBmb3','TotalFloatingBmb3']; - 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]; - field_values=[ - (md.results.TransientSolution[0](1).Vx), - (md.results.TransientSolution[0](1).Vy), - (md.results.TransientSolution[0](1).Vel), - (md.results.TransientSolution[0](1).Pressure), - (md.results.TransientSolution[0](1).Base), - (md.results.TransientSolution[0](1).Surface), - (md.results.TransientSolution[0](1).Thickness), - (md.results.TransientSolution[0](1).TotalGroundedBmb), - (md.results.TransientSolution[0](1).TotalFloatingBmb), - (md.results.TransientSolution[0](2).Vx), - (md.results.TransientSolution[0](2).Vy), - (md.results.TransientSolution[0](2).Vel), - (md.results.TransientSolution[0](2).Pressure), - (md.results.TransientSolution[0](2).Base), - (md.results.TransientSolution[0](2).Surface), - (md.results.TransientSolution[0](2).Thickness), - (md.results.TransientSolution[0](2).TotalGroundedBmb), - (md.results.TransientSolution[0](2).TotalFloatingBmb), - (md.results.TransientSolution[0](3).Vx), - (md.results.TransientSolution[0](3).Vy), - (md.results.TransientSolution[0](3).Vel), - (md.results.TransientSolution[0](3).Pressure), - (md.results.TransientSolution[0](3).Base), - (md.results.TransientSolution[0](3).Surface), - (md.results.TransientSolution[0](3).Thickness), - (md.results.TransientSolution[0](3).TotalGroundedBmb), - (md.results.TransientSolution[0](3).TotalFloatingBmb), - ]; - }); -}) +//describe("test105", function() { + //it("contains test105", function() { + ////Test Name: SquareShelfConstrainedMasstransp2d + //var md = new model(); + //triangle(md,square[0],150000.); + //setmask(md,'all',''); + //parameterize(md); + //setflowequation(md,'SSA','all'); + ////md.cluster=generic('name',oshostname(),'np',3); + //md=solve(md,MasstransportSolutionEnum()); + + ////Fields and tolerances to track changes + //field_names =['Thickness']; + //field_tolerances=[1e-13]; + //field_values=[ + //(md.results.MasstransportSolution[0].Thickness), + //]; + //}); +//}) + +//describe("test110", function() { + //it("contains test110", function() { + ////Test Name: SquareShelfConstrainedTranSSA2d + //var md = new model(); + //triangle(md,square[0],150000.); + //setmask(md,'all',''); + //parameterize(md); + //setflowequation(md,'SSA','all'); + ////md.cluster=generic('name',oshostname(),'np',3); + //md.trans.requested_outputs=['IceVolume']; + + //md=solve(md,TransientSolutionEnum()); + + ////Fields and tolerances to track changes + //field_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']; + //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=[ + //(md.results.TransientSolution[0](1).Vx), + //(md.results.TransientSolution[0](1).Vy), + //(md.results.TransientSolution[0](1).Vel), + //(md.results.TransientSolution[0](1).Pressure), + //(md.results.TransientSolution[0](1).Base), + //(md.results.TransientSolution[0](1).Surface), + //(md.results.TransientSolution[0](1).Thickness), + //(md.results.TransientSolution[0](1).IceVolume), + //(md.results.TransientSolution[0](2).Vx), + //(md.results.TransientSolution[0](2).Vy), + //(md.results.TransientSolution[0](2).Vel), + //(md.results.TransientSolution[0](2).Pressure), + //(md.results.TransientSolution[0](2).Base), + //(md.results.TransientSolution[0](2).Surface), + //(md.results.TransientSolution[0](2).Thickness), + //(md.results.TransientSolution[0](2).IceVolume), + //(md.results.TransientSolution[0](3).Vx), + //(md.results.TransientSolution[0](3).Vy), + //(md.results.TransientSolution[0](3).Vel), + //(md.results.TransientSolution[0](3).Pressure), + //(md.results.TransientSolution[0](3).Base), + //(md.results.TransientSolution[0](3).Surface), + //(md.results.TransientSolution[0](3).Thickness), + //(md.results.TransientSolution[0](3).IceVolume), + //]; + //}); +//}) + +//describe("test112", function() { + //it("contains test112", function() { + ////Test Name: SquareShelfConstrainedSurfSlop2d + //var md = new model(); + //triangle(md,square[0],150000.); + //setmask(md,'all',''); + //parameterize(md); + //setflowequation(md,'SSA','all'); + ////md.cluster=generic('name',oshostname(),'np',3); + //md=solve(md,SurfaceSlopeSolutionEnum()); + + ////Fields and tolerances to track changes + //field_names =['SurfaceSlopeX','SurfaceSlopeY']; + //field_tolerances=[1e-13,1e-13]; + //field_values=[ + //(md.results.SurfaceSlopeSolution[0].SurfaceSlopeX), + //(md.results.SurfaceSlopeSolution[0].SurfaceSlopeY), + //]; + //}); +//}) + +//describe("test114", function() { + //it("contains test114", function() { + ////Test Name: SquareShelfConstrainedBedSlop2d + //var md = new model(); + //triangle(md,square[0],150000.); + //setmask(md,'all',''); + //parameterize(md); + //setflowequation(md,'SSA','all'); + ////md.cluster=generic('name',oshostname(),'np',3); + //md=solve(md,BedSlopeSolutionEnum()); + + ////Fields and tolerances to track changes + //field_names =['BedSlopeX','BedSlopeY']; + //field_tolerances=[1e-13,1e-13]; + //field_values=[ + //(md.results.BedSlopeSolution[0].BedSlopeX), + //(md.results.BedSlopeSolution[0].BedSlopeY), + //]; + //}); +//}) + +//describe("test201", function() { + //it("contains test201", function() { + ////Test Name: SquareShelfStressSSA2d + //var md = new model(); + //triangle(md,square[0],150000.); + //setmask(md,'all',''); + //parameterize(md); + //setflowequation(md,'SSA','all'); + ////md.cluster=generic('name',oshostname(),'np',3); + //md=solve(md,StressbalanceSolutionEnum()); + + ////Fields and tolerances to track changes + //field_names =['Vx','Vy','Vel','Pressure']; + //field_tolerances=[1e-13,1e-13,1e-13,1e-13]; + //field_values=[ + //(md.results.StressbalanceSolution[0].Vx), + //(md.results.StressbalanceSolution[0].Vy), + //(md.results.StressbalanceSolution[0].Vel), + //(md.results.StressbalanceSolution[0].Pressure), + //]; + //}); +//}) + +//describe("test208", function() { + //it("contains test208", function() { + ////Test Name: SquareShelfTranSSA2d + //var md = new model(); + //triangle(md,square[0],150000.); + //setmask(md,'all',''); + //parameterize(md); + //setflowequation(md,'SSA','all'); + ////md.cluster=generic('name',oshostname(),'np',3); + //md.trans.requested_outputs=['default','FloatingArea','GroundedArea','TotalGroundedBmb','TotalFloatingBmb']; + //for (var i = 0; i < md.basalforcings.floatingice_melting_rate.length; ++i) { + //md.basalforcings.floatingice_melting_rate[i] = 1; + //} + //md=solve(md,TransientSolutionEnum()); + + ////Fields and tolerances to track changes + //field_names =['Vx1','Vy1','Vel1','Pressure1','Bed1','Surface1','Thickness1','TotalGroundedBmb1','TotalFloatingBmb1', + //'Vx2','Vy2','Vel2','Pressure2','Bed2','Surface2','Thickness2','TotalGroundedBmb2','TotalFloatingBmb2', + //'Vx3','Vy3','Vel3','Pressure3','Bed3','Surface3','Thickness3','TotalGroundedBmb3','TotalFloatingBmb3']; + //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]; + //field_values=[ + //(md.results.TransientSolution[0](1).Vx), + //(md.results.TransientSolution[0](1).Vy), + //(md.results.TransientSolution[0](1).Vel), + //(md.results.TransientSolution[0](1).Pressure), + //(md.results.TransientSolution[0](1).Base), + //(md.results.TransientSolution[0](1).Surface), + //(md.results.TransientSolution[0](1).Thickness), + //(md.results.TransientSolution[0](1).TotalGroundedBmb), + //(md.results.TransientSolution[0](1).TotalFloatingBmb), + //(md.results.TransientSolution[0](2).Vx), + //(md.results.TransientSolution[0](2).Vy), + //(md.results.TransientSolution[0](2).Vel), + //(md.results.TransientSolution[0](2).Pressure), + //(md.results.TransientSolution[0](2).Base), + //(md.results.TransientSolution[0](2).Surface), + //(md.results.TransientSolution[0](2).Thickness), + //(md.results.TransientSolution[0](2).TotalGroundedBmb), + //(md.results.TransientSolution[0](2).TotalFloatingBmb), + //(md.results.TransientSolution[0](3).Vx), + //(md.results.TransientSolution[0](3).Vy), + //(md.results.TransientSolution[0](3).Vel), + //(md.results.TransientSolution[0](3).Pressure), + //(md.results.TransientSolution[0](3).Base), + //(md.results.TransientSolution[0](3).Surface), + //(md.results.TransientSolution[0](3).Thickness), + //(md.results.TransientSolution[0](3).TotalGroundedBmb), + //(md.results.TransientSolution[0](3).TotalFloatingBmb), + //]; + //}); +//}) Index: sm/trunk-jpl/jenkins/javascript/karma/scripts/test101.js =================================================================== --- /issm/trunk-jpl/jenkins/javascript/karma/scripts/test101.js (revision 21203) +++ (revision ) @@ -1,15 +1,0 @@ -console.log('creating model'); -var md = new model(); - -console.log('meshing'); -triangle(md,square[0],40000); - -console.log('parameterization'); -setmask(md,'all',''); -parameterize(md); -setflowequation(md,'SSA','all'); -md.verbose.solution=2; md.verbose.convergence=0; - -md=solve(md,StressbalanceSolutionEnum(),'checkconsistency','no'); - -console.log(md.results['StressbalanceSolution'][0]['Vel']); Index: sm/trunk-jpl/jenkins/javascript/karma/scripts/test102.js =================================================================== --- /issm/trunk-jpl/jenkins/javascript/karma/scripts/test102.js (revision 21203) +++ (revision ) @@ -1,23 +1,0 @@ -//Test Name: SquareShelfConstrainedStressSSA3d -var md = new model(); - -triangle(md,square[0],180000.); - -setmask(md,'all',''); -parameterize(md); - -md.extrude(md,3,2.); -setflowequation(md,'SSA','all'); -//md.cluster=generic('name',oshostname(),'np',3); -md=solve(md,StressbalanceSolutionEnum()); - -//Fields and tolerances to track changes -field_names =['Vx','Vy','Vz','Vel','Pressure']; -field_tolerances=[1e-13,1e-13,1e-13,1e-13,1e-13]; -field_values=[ - (md.results.StressbalanceSolution[0].Vx), - (md.results.StressbalanceSolution[0].Vy), - (md.results.StressbalanceSolution[0].Vz), - (md.results.StressbalanceSolution[0].Vel), - (md.results.StressbalanceSolution[0].Pressure), - ]; Index: sm/trunk-jpl/jenkins/javascript/karma/scripts/test103.js =================================================================== --- /issm/trunk-jpl/jenkins/javascript/karma/scripts/test103.js (revision 21203) +++ (revision ) @@ -1,29 +1,0 @@ -//Test Name: SquareShelfConstrainedStressHO -var md = new model(); -triangle(md,square[0],180000.); -setmask(md,'all',''); -parameterize(md); -md.extrude(md,3,2.); -setflowequation(md,'HO','all'); -//md.cluster=generic('name',oshostname(),'np',3); -md.stressbalance.requested_outputs=['default','StressTensorxx','StressTensoryy','StressTensorzz','StressTensorxy','StressTensorxz','StressTensoryz']; -md=solve(md,StressbalanceSolutionEnum()); - -//Fields and tolerances to track changes -field_names =['Vx','Vy','Vz','Vel','Pressure', - 'StressTensorxx','StressTensoryy','StressTensorzz','StressTensorxy','StressTensorxz','StressTensoryz']; -field_tolerances=[1e-09,1e-09,1e-09,1e-09,1e-09, - 1e-09,1e-09,1e-09,1e-09,1e-09,1e-08]; -field_values=[ - (md.results.StressbalanceSolution[0].Vx), - (md.results.StressbalanceSolution[0].Vy), - (md.results.StressbalanceSolution[0].Vz), - (md.results.StressbalanceSolution[0].Vel), - (md.results.StressbalanceSolution[0].Pressure), - (md.results.StressbalanceSolution[0].StressTensorxx), - (md.results.StressbalanceSolution[0].StressTensoryy), - (md.results.StressbalanceSolution[0].StressTensorzz), - (md.results.StressbalanceSolution[0].StressTensorxy), - (md.results.StressbalanceSolution[0].StressTensorxz), - (md.results.StressbalanceSolution[0].StressTensoryz), - ]; Index: sm/trunk-jpl/jenkins/javascript/karma/scripts/test104.js =================================================================== --- /issm/trunk-jpl/jenkins/javascript/karma/scripts/test104.js (revision 21203) +++ (revision ) @@ -1,20 +1,0 @@ -//Test Name: SquareShelfConstrainedStressFS -var md = new model(); -triangle(md,square[0],180000.); -setmask(md,'all',''); -parameterize(md); -md.extrude(md,3,2.); -setflowequation(md,'FS','all'); -//md.cluster=generic('name',oshostname(),'np',3); -md=solve(md,StressbalanceSolutionEnum(),'checkconsistency','no'); - -//Fields and tolerances to track changes -field_names =['Vx','Vy','Vz','Vel','Pressure']; -field_tolerances=[1e-08,1e-08,1e-06,1e-08,1e-08]; -field_values=[ - (md.results.StressbalanceSolution[0].Vx), - (md.results.StressbalanceSolution[0].Vy), - (md.results.StressbalanceSolution[0].Vz), - (md.results.StressbalanceSolution[0].Vel), - (md.results.StressbalanceSolution[0].Pressure), - ]; Index: sm/trunk-jpl/jenkins/javascript/karma/scripts/test105.js =================================================================== --- /issm/trunk-jpl/jenkins/javascript/karma/scripts/test105.js (revision 21203) +++ (revision ) @@ -1,15 +1,0 @@ -//Test Name: SquareShelfConstrainedMasstransp2d -var md = new model(); -triangle(md,square[0],150000.); -setmask(md,'all',''); -parameterize(md); -setflowequation(md,'SSA','all'); -//md.cluster=generic('name',oshostname(),'np',3); -md=solve(md,MasstransportSolutionEnum()); - -//Fields and tolerances to track changes -field_names =['Thickness']; -field_tolerances=[1e-13]; -field_values=[ - (md.results.MasstransportSolution[0].Thickness), - ]; Index: sm/trunk-jpl/jenkins/javascript/karma/scripts/test106.js =================================================================== --- /issm/trunk-jpl/jenkins/javascript/karma/scripts/test106.js (revision 21203) +++ (revision ) @@ -1,18 +1,0 @@ -//Test Name: SquareShelfConstrainedMasstransp2dDG -var md = new model(); -triangle(md,square[0],150000.); -meshconvert(md); -setmask(md,'all',''); -parameterize(md); -setflowequation(md,'SSA','all'); -//md.cluster=generic('name',oshostname(),'np',3); -md.masstransport.stabilization=3; -md.masstransport.spcthickness=md.geometry.thickness; -md=solve(md,MasstransportSolutionEnum()); - -//Fields and tolerances to track changes -field_names =['Thickness']; -field_tolerances=[1e-13]; -field_values=[ - (md.results.MasstransportSolution[0].Thickness), - ]; Index: sm/trunk-jpl/jenkins/javascript/karma/scripts/test107.js =================================================================== --- /issm/trunk-jpl/jenkins/javascript/karma/scripts/test107.js (revision 21203) +++ (revision ) @@ -1,16 +1,0 @@ -//Test Name: SquareShelfConstrainedMasstransp3d -var md = new model(); -triangle(md,square[0],150000.); -setmask(md,'all',''); -parameterize(md); -setflowequation(md,'SSA','all'); -md.extrude(md,5,3.); -//md.cluster=generic('name',oshostname(),'np',3); -md=solve(md,MasstransportSolutionEnum()); - -//Fields and tolerances to track changes -field_names =['Thickness']; -field_tolerances=[1e-13]; -field_values=[ - (md.results.MasstransportSolution[0].Thickness), - ]; Index: sm/trunk-jpl/jenkins/javascript/karma/scripts/test108.js =================================================================== --- /issm/trunk-jpl/jenkins/javascript/karma/scripts/test108.js (revision 21203) +++ (revision ) @@ -1,18 +1,0 @@ -//Test Name: SquareShelfConstrainedTherStea -var md = new model(); -triangle(md,square[0],180000.); -setmask(md,'all',''); -parameterize(md); -md.extrude(md,3,1.); -setflowequation(md,'SSA','all'); -md.timestepping.time_step=0; -//md.cluster=generic('name',oshostname(),'np',3); -md=solve(md,ThermalSolutionEnum()); - -//Fields and tolerances to track changes -field_names =['Temperature','BasalforcingsGroundediceMeltingRate']; -field_tolerances=[1e-13,1e-5]; -field_values=[ - (md.results.ThermalSolution[0].Temperature), - (md.results.ThermalSolution[0].BasalforcingsGroundediceMeltingRate), - ]; Index: sm/trunk-jpl/jenkins/javascript/karma/scripts/test109.js =================================================================== --- /issm/trunk-jpl/jenkins/javascript/karma/scripts/test109.js (revision 21203) +++ (revision ) @@ -1,22 +1,0 @@ -//Test Name: SquareShelfConstrainedTherTran -var md = new model(); -triangle(md,square[0],180000.); -setmask(md,'all',''); -parameterize(md); -md.extrude(md,3,1.); -setflowequation(md,'SSA','all'); -//md.cluster=generic('name',oshostname(),'np',3); -md.transient.isstressbalance=0; -md.transient.ismasstransport=0; -md.transient.issmb=1; -md.transient.isthermal=1; -md.transient.isgroundingline=0; -md=solve(md,TransientSolutionEnum()); - -//Fields and tolerances to track changes -field_names =['Temperature','BasalforcingsGroundediceMeltingRate']; -field_tolerances=[1e-13,1e-13]; -field_values=[ - (md.results.TransientSolution[0](1).Temperature), - (md.results.TransientSolution[0](1).BasalforcingsGroundediceMeltingRate), - ]; Index: sm/trunk-jpl/jenkins/javascript/karma/scripts/test110.js =================================================================== --- /issm/trunk-jpl/jenkins/javascript/karma/scripts/test110.js (revision 21203) +++ (revision ) @@ -1,42 +1,0 @@ -//Test Name: SquareShelfConstrainedTranSSA2d -var md = new model(); -triangle(md,square[0],150000.); -setmask(md,'all',''); -parameterize(md); -setflowequation(md,'SSA','all'); -//md.cluster=generic('name',oshostname(),'np',3); -md.transient.requested_outputs=['IceVolume']; - -md=solve(md,TransientSolutionEnum()); - -//Fields and tolerances to track changes -field_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']; -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=[ - (md.results.TransientSolution[0](1).Vx), - (md.results.TransientSolution[0](1).Vy), - (md.results.TransientSolution[0](1).Vel), - (md.results.TransientSolution[0](1).Pressure), - (md.results.TransientSolution[0](1).Base), - (md.results.TransientSolution[0](1).Surface), - (md.results.TransientSolution[0](1).Thickness), - (md.results.TransientSolution[0](1).IceVolume), - (md.results.TransientSolution[0](2).Vx), - (md.results.TransientSolution[0](2).Vy), - (md.results.TransientSolution[0](2).Vel), - (md.results.TransientSolution[0](2).Pressure), - (md.results.TransientSolution[0](2).Base), - (md.results.TransientSolution[0](2).Surface), - (md.results.TransientSolution[0](2).Thickness), - (md.results.TransientSolution[0](2).IceVolume), - (md.results.TransientSolution[0](3).Vx), - (md.results.TransientSolution[0](3).Vy), - (md.results.TransientSolution[0](3).Vel), - (md.results.TransientSolution[0](3).Pressure), - (md.results.TransientSolution[0](3).Base), - (md.results.TransientSolution[0](3).Surface), - (md.results.TransientSolution[0](3).Thickness), - (md.results.TransientSolution[0](3).IceVolume), - ]; Index: sm/trunk-jpl/jenkins/javascript/karma/scripts/test111.js =================================================================== --- /issm/trunk-jpl/jenkins/javascript/karma/scripts/test111.js (revision 21203) +++ (revision ) @@ -1,55 +1,0 @@ -//Test Name: SquareShelfConstrainedTranHO -var md = new model(); -triangle(md,square[0],200000.); -setmask(md,'all',''); -parameterize(md); -md.extrude(md,3,1.); -setflowequation(md,'HO','all'); -//md.cluster=generic('name',oshostname(),'np',3); -md.transient.requested_outputs=['IceVolume']; -md=solve(md,TransientSolutionEnum()); - -//Fields and tolerances to track changes -field_names=[ - 'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface2','Thickness1','Temperature1','BasalforcingsGroundediceMeltingRate1','Volume1', - 'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2','BasalforcingsGroundediceMeltingRate2','Volume2', - 'Vx3','Vy3','Vz3','Vel3','Pressure3','Bed3','Surface3','Thickness3','Temperature3','BasalforcingsGroundediceMeltingRate3','Volume3']; -field_tolerances=[ - 1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-05,1e-13, - 1e-09,1e-08,1e-08,1e-08,1e-09,1e-09,1e-09,1e-09,1e-09,1e-05,2e-11, - 1e-09,1e-08,1e-08,1e-08,1e-09,1e-09,1e-09,1e-09,1e-09,1e-05,1e-11]; -field_values=[ - (md.results.TransientSolution[0](1).Vx), - (md.results.TransientSolution[0](1).Vy), - (md.results.TransientSolution[0](1).Vz), - (md.results.TransientSolution[0](1).Vel), - (md.results.TransientSolution[0](1).Pressure), - (md.results.TransientSolution[0](1).Base), - (md.results.TransientSolution[0](1).Surface), - (md.results.TransientSolution[0](1).Thickness), - (md.results.TransientSolution[0](1).Temperature), - (md.results.TransientSolution[0](1).BasalforcingsGroundediceMeltingRate), - (md.results.TransientSolution[0](1).IceVolume), - (md.results.TransientSolution[0](2).Vx), - (md.results.TransientSolution[0](2).Vy), - (md.results.TransientSolution[0](2).Vz), - (md.results.TransientSolution[0](2).Vel), - (md.results.TransientSolution[0](2).Pressure), - (md.results.TransientSolution[0](2).Base), - (md.results.TransientSolution[0](2).Surface), - (md.results.TransientSolution[0](2).Thickness), - (md.results.TransientSolution[0](2).Temperature), - (md.results.TransientSolution[0](2).BasalforcingsGroundediceMeltingRate), - (md.results.TransientSolution[0](2).IceVolume), - (md.results.TransientSolution[0](3).Vx), - (md.results.TransientSolution[0](3).Vy), - (md.results.TransientSolution[0](3).Vz), - (md.results.TransientSolution[0](3).Vel), - (md.results.TransientSolution[0](3).Pressure), - (md.results.TransientSolution[0](3).Base), - (md.results.TransientSolution[0](3).Surface), - (md.results.TransientSolution[0](3).Thickness), - (md.results.TransientSolution[0](3).Temperature), - (md.results.TransientSolution[0](3).BasalforcingsGroundediceMeltingRate), - (md.results.TransientSolution[0](3).IceVolume), - ]; Index: sm/trunk-jpl/jenkins/javascript/karma/scripts/test112.js =================================================================== --- /issm/trunk-jpl/jenkins/javascript/karma/scripts/test112.js (revision 21203) +++ (revision ) @@ -1,16 +1,0 @@ -//Test Name: SquareShelfConstrainedSurfSlop2d -var md = new model(); -triangle(md,square[0],150000.); -setmask(md,'all',''); -parameterize(md); -setflowequation(md,'SSA','all'); -//md.cluster=generic('name',oshostname(),'np',3); -md=solve(md,SurfaceSlopeSolutionEnum()); - -//Fields and tolerances to track changes -field_names =['SurfaceSlopeX','SurfaceSlopeY']; -field_tolerances=[1e-13,1e-13]; -field_values=[ - (md.results.SurfaceSlopeSolution[0].SurfaceSlopeX), - (md.results.SurfaceSlopeSolution[0].SurfaceSlopeY), - ]; Index: sm/trunk-jpl/jenkins/javascript/karma/scripts/test113.js =================================================================== --- /issm/trunk-jpl/jenkins/javascript/karma/scripts/test113.js (revision 21203) +++ (revision ) @@ -1,17 +1,0 @@ -//Test Name: SquareShelfConstrainedSurfSlope3d -var md = new model(); -triangle(md,square[0],150000.); -setmask(md,'all',''); -parameterize(md); -md.extrude(md,5,1.); -setflowequation(md,'SSA','all'); -//md.cluster=generic('name',oshostname(),'np',3); -md=solve(md,SurfaceSlopeSolutionEnum()); - -//Fields and tolerances to track changes -field_names =['SurfaceSlopeX','SurfaceSlopeY']; -field_tolerances=[1e-13,1e-13]; -field_values=[ - (md.results.SurfaceSlopeSolution[0].SurfaceSlopeX), - (md.results.SurfaceSlopeSolution[0].SurfaceSlopeY), - ]; Index: sm/trunk-jpl/jenkins/javascript/karma/scripts/test114.js =================================================================== --- /issm/trunk-jpl/jenkins/javascript/karma/scripts/test114.js (revision 21203) +++ (revision ) @@ -1,16 +1,0 @@ -//Test Name: SquareShelfConstrainedBedSlop2d -var md = new model(); -triangle(md,square[0],150000.); -setmask(md,'all',''); -parameterize(md); -setflowequation(md,'SSA','all'); -//md.cluster=generic('name',oshostname(),'np',3); -md=solve(md,BedSlopeSolutionEnum()); - -//Fields and tolerances to track changes -field_names =['BedSlopeX','BedSlopeY']; -field_tolerances=[1e-13,1e-13]; -field_values=[ - (md.results.BedSlopeSolution[0].BedSlopeX), - (md.results.BedSlopeSolution[0].BedSlopeY), - ]; Index: sm/trunk-jpl/jenkins/javascript/karma/scripts/test115.js =================================================================== --- /issm/trunk-jpl/jenkins/javascript/karma/scripts/test115.js (revision 21203) +++ (revision ) @@ -1,17 +1,0 @@ -//Test Name: SquareShelfConstrainedBedSlop3d -var md = new model(); -triangle(md,square[0],150000.); -setmask(md,'all',''); -parameterize(md); -md.extrude(md,5,1.); -setflowequation(md,'SSA','all'); -//md.cluster=generic('name',oshostname(),'np',3); -md=solve(md,BedSlopeSolutionEnum()); - -//Fields and tolerances to track changes -field_names =['BedSlopeX','BedSlopeY']; -field_tolerances=[1e-13,1e-13]; -field_values=[ - (md.results.BedSlopeSolution[0].BedSlopeX), - (md.results.BedSlopeSolution[0].BedSlopeY), - ]; Index: sm/trunk-jpl/jenkins/javascript/karma/scripts/test116.js =================================================================== --- /issm/trunk-jpl/jenkins/javascript/karma/scripts/test116.js (revision 21203) +++ (revision ) @@ -1,26 +1,0 @@ -//Test Name: SquareShelfConstrainedBalThic2d -var md = new model(); -triangle(md,square[0],150000.); -setmask(md,'all',''); -parameterize(md); -//Add boundary conditions on thickness on the border -pos=[]; -for (var i = 0; i < md.mesh.vertexonboundary.length; ++i) { - if ((md.mesh.vertexonboundary[i] !== 0)) { - pos.push(i); - }; -} - -for (var i = 0; i < pos.length; ++i) { - md.balancethickness.spcthickness[pos[i]] = md.geometry.thickness[pos[i]]; -} -setflowequation(md,'SSA','all'); -//md.cluster=generic('name',oshostname(),'np',3); -md=solve(md,BalancethicknessSolutionEnum()); - -//Fields and tolerances to track changes -field_names =['Thickness']; -field_tolerances=[1e-13]; -field_values=[ - (md.results.BalancethicknessSolution[0].Thickness), - ]; Index: sm/trunk-jpl/jenkins/javascript/karma/scripts/test120.js =================================================================== --- /issm/trunk-jpl/jenkins/javascript/karma/scripts/test120.js (revision 21203) +++ (revision ) @@ -1,31 +1,0 @@ -//Test Name: SquareShelfConstrainedEnthalpyStea -var md = new model(); -function zeros(...args) { - var array = []; - for (var i = 0; i < args[0]; ++i) { - array.push(args.length == 1 ? 0 : zeros(args.slice(1))); - } - return array; -} -triangle(md,square[0],180000.); -setmask(md,'all',''); -parameterize(md); -md.extrude(md,3,1.); -setflowequation(md,'SSA','all'); -md.timestepping.time_step=0; -md.initialization.waterfraction=zeros(md.mesh.numberofvertices,1); -md.initialization.watercolumn=zeros(md.mesh.numberofvertices,1); -md.thermal.isenthalpy = 1; -md.thermal.isdynamicbasalspc = 1; - -//md.cluster=generic('name',oshostname(),'np',3); -md=solve(md,ThermalSolutionEnum()); - -//Fields and tolerances to track changes -field_names =['Enthalpy','Waterfraction','Temperature']; -field_tolerances=[1e-13,3e-10,1e-13]; -field_values=[ - (md.results.ThermalSolution[0].Enthalpy), - (md.results.ThermalSolution[0].Waterfraction), - (md.results.ThermalSolution[0].Temperature), - ]; Index: sm/trunk-jpl/jenkins/javascript/karma/scripts/test121.js =================================================================== --- /issm/trunk-jpl/jenkins/javascript/karma/scripts/test121.js (revision 21203) +++ (revision ) @@ -1,42 +1,0 @@ -//Test Name: SquareShelfConstrainedEnthalpyTran -function zeros(...args) { - var array = []; - for (var i = 0; i < args[0]; ++i) { - array.push(args.length == 1 ? 0 : zeros(args.slice(1))); - } - return array; -} -var md = new model(); -triangle(md,square[0],180000.); -setmask(md,'all',''); -parameterize(md); -md.extrude(md,3,1.); -setflowequation(md,'SSA','all'); -//md.cluster=generic('name',oshostname(),'np',3); -md.initialization.waterfraction=zeros(md.mesh.numberofvertices,1); -md.initialization.watercolumn=zeros(md.mesh.numberofvertices,1); -md.transient.isstressbalance=0; -md.transient.ismasstransport=0; -md.transient.issmb=1; -md.transient.isthermal=1; -md.transient.isgroundingline=0; -md.thermal.isenthalpy=1; -md.thermal.isdynamicbasalspc=1; -md=solve(md,TransientSolutionEnum()); - -//Fields and tolerances to track changes -field_names =['Enthalpy1','Waterfraction1','Temperature1', - 'Enthalpy2','Waterfraction2','Temperature2', - 'Enthalpy3','Waterfraction3','Temperature3']; -field_tolerances=[1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-9,1e-13]; -field_values=[ - (md.results.TransientSolution[0](1).Enthalpy), - (md.results.TransientSolution[0](1).Waterfraction), - (md.results.TransientSolution[0](1).Temperature), - (md.results.TransientSolution[0](2).Enthalpy), - (md.results.TransientSolution[0](2).Waterfraction), - (md.results.TransientSolution[0](2).Temperature), - (md.results.TransientSolution[0](3).Enthalpy), - (md.results.TransientSolution[0](3).Waterfraction), - (md.results.TransientSolution[0](3).Temperature), - ]; Index: sm/trunk-jpl/jenkins/javascript/karma/scripts/test122.js =================================================================== --- /issm/trunk-jpl/jenkins/javascript/karma/scripts/test122.js (revision 21203) +++ (revision ) @@ -1,65 +1,0 @@ -//Test Name: SquareShelfConstrainedTransHOEnth -function zeros(...args) { - var array = []; - for (var i = 0; i < args[0]; ++i) { - array.push(args.length == 1 ? 0 : zeros(args.slice(1))); - } - return array; -} -var md = new model(); -triangle(md,square[0],200000.); -setmask(md,'all',''); -parameterize(md); -md.extrude(md,3,1.); -setflowequation(md,'HO','all'); -md.initialization.waterfraction=zeros(md.mesh.numberofvertices,1); -md.initialization.watercolumn=zeros(md.mesh.numberofvertices,1); -md.thermal.isenthalpy=1; -md.thermal.isdynamicbasalspc=1; -md.thermal.stabilization=2; -//md.cluster=generic('name',oshostname(),'np',3); -md=solve(md,TransientSolutionEnum()); - -//Fields and tolerances to track changes -field_names =['Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1','Enthalpy1','Waterfraction1', - 'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2','Enthalpy2','Waterfraction2', - 'Vx3','Vy3','Vz3','Vel3','Pressure3','Bed3','Surface3','Thickness3','Temperature3','Enthalpy3','Waterfraction3']; -field_tolerances=[ - 1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09, - 1e-09,1e-09,1e-08,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09, - 1e-09,1e-09,1e-08,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09]; -field_values=[ - (md.results.TransientSolution[0](1).Vx), - (md.results.TransientSolution[0](1).Vy), - (md.results.TransientSolution[0](1).Vz), - (md.results.TransientSolution[0](1).Vel), - (md.results.TransientSolution[0](1).Pressure), - (md.results.TransientSolution[0](1).Base), - (md.results.TransientSolution[0](1).Surface), - (md.results.TransientSolution[0](1).Thickness), - (md.results.TransientSolution[0](1).Temperature), - (md.results.TransientSolution[0](1).Enthalpy), - (md.results.TransientSolution[0](1).Waterfraction), - (md.results.TransientSolution[0](2).Vx), - (md.results.TransientSolution[0](2).Vy), - (md.results.TransientSolution[0](2).Vz), - (md.results.TransientSolution[0](2).Vel), - (md.results.TransientSolution[0](2).Pressure), - (md.results.TransientSolution[0](2).Base), - (md.results.TransientSolution[0](2).Surface), - (md.results.TransientSolution[0](2).Thickness), - (md.results.TransientSolution[0](2).Temperature), - (md.results.TransientSolution[0](2).Enthalpy), - (md.results.TransientSolution[0](2).Waterfraction), - (md.results.TransientSolution[0](3).Vx), - (md.results.TransientSolution[0](3).Vy), - (md.results.TransientSolution[0](3).Vz), - (md.results.TransientSolution[0](3).Vel), - (md.results.TransientSolution[0](3).Pressure), - (md.results.TransientSolution[0](3).Base), - (md.results.TransientSolution[0](3).Surface), - (md.results.TransientSolution[0](3).Thickness), - (md.results.TransientSolution[0](3).Temperature), - (md.results.TransientSolution[0](3).Enthalpy), - (md.results.TransientSolution[0](3).Waterfraction), - ]; Index: sm/trunk-jpl/jenkins/javascript/karma/scripts/test123.js =================================================================== --- /issm/trunk-jpl/jenkins/javascript/karma/scripts/test123.js (revision 21203) +++ (revision ) @@ -1,23 +1,0 @@ -//Test Name: SquareShelfConstrainedTranMisfitSurface -var md = new model(); -triangle(md,square[0],180000.); -setmask(md,'all',''); -parameterize(md); -setflowequation(md,'SSA','all'); -//md.cluster=generic('name',oshostname(),'np',3); - -fake_surface=[[md.geometry.surface+100,1.1],[md.geometry.surface+200,2.1],[md.geometry.surface+300,2.5]]; - -md.transient.requested_outputs=['default','SurfaceMisfit']; -md.outputdefinition.definitions=[misfit('name','SurfaceMisfit', 'definitionenum',Outputdefinition1Enum,'model_enum',SurfaceEnum,'observation_enum',SurfaceObservationEnum,'observation',fake_surface,'timeinterpolation','nearestneighbor','weights',ones(md.mesh.numberofvertices,1),'weights_enum',WeightsSurfaceObservationEnum)]; - -md=solve(md,TransientSolutionEnum()); - -//Fields and tolerances to track changes -field_names =['SurfaceMisfitFirstStep','SurfaceMisfitSecondStep','SurfaceMisfitThirdStep']; -field_tolerances=[1e-13,1e-13,1e-13]; -field_values=[ - (md.results.TransientSolution[0](1).SurfaceMisfit), - (md.results.TransientSolution[0](2).SurfaceMisfit), - (md.results.TransientSolution[0](3).SurfaceMisfit) - ]; Index: sm/trunk-jpl/jenkins/javascript/karma/scripts/test124.js =================================================================== --- /issm/trunk-jpl/jenkins/javascript/karma/scripts/test124.js (revision 21203) +++ (revision ) @@ -1,49 +1,0 @@ -//Test Name: SquareShelfConstrainedTranFSFreeSurface -var md = new model(); -triangle(md,square[0],150000.); -setmask(md,'all',''); -parameterize(md); -md.extrude(md,3,1.); -setflowequation(md,'FS','all'); - -//Free surface -md.masstransport.isfreesurface=1; -md.timestepping.time_step=0.00001; -md.timestepping.final_time=0.00005; - -//Go solve -//md.cluster=generic('name',oshostname(),'np',3); -md=solve(md,TransientSolutionEnum()); - -//Fields and tolerances to track changes -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=[ - 2e-09,3e-9,3e-9,3e-9,1e-13,1e-12,1e-12, - 2e-09,3e-9,3e-9,3e-9,1e-10,1e-10,1e-10, - 2e-09,3e-9,3e-9,3e-9,1e-10,1e-10,1e-10]; -field_values=[ - (md.results.TransientSolution[0](1).Vx), - (md.results.TransientSolution[0](1).Vy), - (md.results.TransientSolution[0](1).Vel), - (md.results.TransientSolution[0](1).Pressure), - (md.results.TransientSolution[0](1).Base), - (md.results.TransientSolution[0](1).Surface), - (md.results.TransientSolution[0](1).Thickness), - (md.results.TransientSolution[0](2).Vx), - (md.results.TransientSolution[0](2).Vy), - (md.results.TransientSolution[0](2).Vel), - (md.results.TransientSolution[0](2).Pressure), - (md.results.TransientSolution[0](2).Base), - (md.results.TransientSolution[0](2).Surface), - (md.results.TransientSolution[0](2).Thickness), - (md.results.TransientSolution[0](3).Vx), - (md.results.TransientSolution[0](3).Vy), - (md.results.TransientSolution[0](3).Vel), - (md.results.TransientSolution[0](3).Pressure), - (md.results.TransientSolution[0](3).Base), - (md.results.TransientSolution[0](3).Surface), - (md.results.TransientSolution[0](3).Thickness), - ]; Index: sm/trunk-jpl/jenkins/javascript/karma/scripts/test125.js =================================================================== --- /issm/trunk-jpl/jenkins/javascript/karma/scripts/test125.js (revision 21203) +++ (revision ) @@ -1,49 +1,0 @@ -//Test Name: SquareShelfConstrainedRestartTranSSA2d -var md = new model(); -triangle(md,square[0],150000.); -setmask(md,'all',''); -parameterize(md); -setflowequation(md,'SSA','all'); -//md.cluster=generic('name',oshostname(),'np',1); -md.transient.requested_outputs=['IceVolume','TotalSmb']; - -//md.verbose=verbose('solution',true); -md.settings.recording_frequency=5; - -// time steps and resolution -md.timestepping.final_time=8; - -md=solve(md,TransientSolutionEnum()); -md2=solve(md,TransientSolutionEnum(),'restart',1); - -//Fields and tolerances to track changes -field_names =['Vx1','Vy1','Vel1','TotalSmb1','Bed1','Surface1','Thickness1','Volume1','Vx2','Vy2','Vel2','TotalSmb2','Bed2','Surface2','Thickness2','Volume2','Vx3','Vy3','Vel3','TotalSmb3','Bed3','Surface3','Thickness3','Volume3']; -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=[ - (md.results.TransientSolution(6).Vx)-(md2.results.TransientSolution[0](1).Vx), - (md.results.TransientSolution(6).Vy)-(md2.results.TransientSolution[0](1).Vy), - (md.results.TransientSolution(6).Vel)-(md2.results.TransientSolution[0](1).Vel), - (md.results.TransientSolution(6).TotalSmb)-(md2.results.TransientSolution[0](1).TotalSmb), - (md.results.TransientSolution(6).Base)-(md2.results.TransientSolution[0](1).Base), - (md.results.TransientSolution(6).Surface)-(md2.results.TransientSolution[0](1).Surface), - (md.results.TransientSolution(6).Thickness)-(md2.results.TransientSolution[0](1).Thickness), - (md.results.TransientSolution(6).IceVolume)-(md2.results.TransientSolution[0](1).IceVolume), - (md.results.TransientSolution(7).Vx)-(md2.results.TransientSolution[0](2).Vx), - (md.results.TransientSolution(7).Vy)-(md2.results.TransientSolution[0](2).Vy), - (md.results.TransientSolution(7).Vel)-(md2.results.TransientSolution[0](2).Vel), - (md.results.TransientSolution(7).TotalSmb)-(md2.results.TransientSolution[0](2).TotalSmb), - (md.results.TransientSolution(7).Base)-(md2.results.TransientSolution[0](2).Base), - (md.results.TransientSolution(7).Surface)-(md2.results.TransientSolution[0](2).Surface), - (md.results.TransientSolution(7).Thickness)-(md2.results.TransientSolution[0](2).Thickness), - (md.results.TransientSolution(7).IceVolume)-(md2.results.TransientSolution[0](2).IceVolume), - (md.results.TransientSolution(8).Vx)-(md2.results.TransientSolution[0](3).Vx), - (md.results.TransientSolution(8).Vy)-(md2.results.TransientSolution[0](3).Vy), - (md.results.TransientSolution(8).Vel)-(md2.results.TransientSolution[0](3).Vel), - (md.results.TransientSolution(8).TotalSmb)-(md2.results.TransientSolution[0](3).TotalSmb), - (md.results.TransientSolution(8).Base)-(md2.results.TransientSolution[0](3).Base), - (md.results.TransientSolution(8).Surface)-(md2.results.TransientSolution[0](3).Surface), - (md.results.TransientSolution(8).Thickness)-(md2.results.TransientSolution[0](3).Thickness), - (md.results.TransientSolution(8).IceVolume)-(md2.results.TransientSolution[0](3).IceVolume), - ]; Index: sm/trunk-jpl/jenkins/javascript/karma/scripts/test126.js =================================================================== --- /issm/trunk-jpl/jenkins/javascript/karma/scripts/test126.js (revision 21203) +++ (revision ) @@ -1,57 +1,0 @@ -//Test Name: SquareShelfConstrainedRestartTranHO3d -var md = new model(); -triangle(md,square[0],150000.); -setmask(md,'all',''); -parameterize(md); -//md.cluster=generic('name',oshostname(),'np',1); -md.extrude(md,3,1.); -setflowequation(md,'HO','all'); -md.transient.requested_outputs=['IceVolume','TotalSmb']; - -//md.verbose=verbose('solution',true); -md.settings.recording_frequency=5; - -// time steps and resolution -md.timestepping.final_time=8; - -md=solve(md,TransientSolutionEnum()); -md2=solve(md,TransientSolutionEnum(),'restart',1); - -//Fields and tolerances to track changes -field_names =['Vx1','Vy1','Vel1','TotalSmb1','Bed1','Surface1','Thickness1','Volume1','Temperature1','Pressure1','Vx2','Vy2','Vel2','TotalSmb2','Bed2','Surface2','Thickness2','Volume2','Temperature2','Pressure2','Vx3','Vy3','Vel3','TotalSmb3','Bed3','Surface3','Thickness3','Volume3','Temperature3','Pressure3']; -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]; -field_values=[ - (md.results.TransientSolution(6).Vx)-(md2.results.TransientSolution[0](1).Vx), - (md.results.TransientSolution(6).Vy)-(md2.results.TransientSolution[0](1).Vy), - (md.results.TransientSolution(6).Vel)-(md2.results.TransientSolution[0](1).Vel), - (md.results.TransientSolution(6).TotalSmb)-(md2.results.TransientSolution[0](1).TotalSmb), - (md.results.TransientSolution(6).Base)-(md2.results.TransientSolution[0](1).Base), - (md.results.TransientSolution(6).Surface)-(md2.results.TransientSolution[0](1).Surface), - (md.results.TransientSolution(6).Thickness)-(md2.results.TransientSolution[0](1).Thickness), - (md.results.TransientSolution(6).IceVolume)-(md2.results.TransientSolution[0](1).IceVolume), - (md.results.TransientSolution(6).Temperature)-(md2.results.TransientSolution[0](1).Temperature), - (md.results.TransientSolution(6).Pressure)-(md2.results.TransientSolution[0](1).Pressure), - (md.results.TransientSolution(7).Vx)-(md2.results.TransientSolution[0](2).Vx), - (md.results.TransientSolution(7).Vy)-(md2.results.TransientSolution[0](2).Vy), - (md.results.TransientSolution(7).Vel)-(md2.results.TransientSolution[0](2).Vel), - (md.results.TransientSolution(7).TotalSmb)-(md2.results.TransientSolution[0](2).TotalSmb), - (md.results.TransientSolution(7).Base)-(md2.results.TransientSolution[0](2).Base), - (md.results.TransientSolution(7).Surface)-(md2.results.TransientSolution[0](2).Surface), - (md.results.TransientSolution(7).Thickness)-(md2.results.TransientSolution[0](2).Thickness), - (md.results.TransientSolution(7).IceVolume)-(md2.results.TransientSolution[0](2).IceVolume), - (md.results.TransientSolution(7).Temperature)-(md2.results.TransientSolution[0](2).Temperature), - (md.results.TransientSolution(7).Pressure)-(md2.results.TransientSolution[0](2).Pressure), - (md.results.TransientSolution(8).Vx)-(md2.results.TransientSolution[0](3).Vx), - (md.results.TransientSolution(8).Vy)-(md2.results.TransientSolution[0](3).Vy), - (md.results.TransientSolution(8).Vel)-(md2.results.TransientSolution[0](3).Vel), - (md.results.TransientSolution(8).TotalSmb)-(md2.results.TransientSolution[0](3).TotalSmb), - (md.results.TransientSolution(8).Base)-(md2.results.TransientSolution[0](3).Base), - (md.results.TransientSolution(8).Surface)-(md2.results.TransientSolution[0](3).Surface), - (md.results.TransientSolution(8).Thickness)-(md2.results.TransientSolution[0](3).Thickness), - (md.results.TransientSolution(8).IceVolume)-(md2.results.TransientSolution[0](3).IceVolume), - (md.results.TransientSolution(8).Temperature)-(md2.results.TransientSolution[0](3).Temperature), - (md.results.TransientSolution(8).Pressure)-(md2.results.TransientSolution[0](3).Pressure), - ]; Index: sm/trunk-jpl/jenkins/javascript/karma/scripts/test201.js =================================================================== --- /issm/trunk-jpl/jenkins/javascript/karma/scripts/test201.js (revision 21203) +++ (revision ) @@ -1,18 +1,0 @@ -//Test Name: SquareShelfStressSSA2d -var md = new model(); -triangle(md,square[0],150000.); -setmask(md,'all',''); -parameterize(md); -setflowequation(md,'SSA','all'); -//md.cluster=generic('name',oshostname(),'np',3); -md=solve(md,StressbalanceSolutionEnum()); - -//Fields and tolerances to track changes -field_names =['Vx','Vy','Vel','Pressure']; -field_tolerances=[1e-13,1e-13,1e-13,1e-13]; -field_values=[ - (md.results.StressbalanceSolution[0].Vx), - (md.results.StressbalanceSolution[0].Vy), - (md.results.StressbalanceSolution[0].Vel), - (md.results.StressbalanceSolution[0].Pressure), - ]; Index: sm/trunk-jpl/jenkins/javascript/karma/scripts/test202.js =================================================================== --- /issm/trunk-jpl/jenkins/javascript/karma/scripts/test202.js (revision 21203) +++ (revision ) @@ -1,20 +1,0 @@ -//Test Name: SquareShelfStressSSA3d -var md = new model(); -triangle(md,square[0],180000.); -setmask(md,'all',''); -parameterize(md); -md.extrude(md,3,2.); -setflowequation(md,'SSA','all'); -//md.cluster=generic('name',oshostname(),'np',3); -md=solve(md,StressbalanceSolutionEnum()); - -//Fields and tolerances to track changes -field_names =['Vx','Vy','Vz','Vel','Pressure']; -field_tolerances=[1e-13,1e-13,1e-13,1e-13,1e-13]; -field_values=[ - (md.results.StressbalanceSolution[0].Vx), - (md.results.StressbalanceSolution[0].Vy), - (md.results.StressbalanceSolution[0].Vz), - (md.results.StressbalanceSolution[0].Vel), - (md.results.StressbalanceSolution[0].Pressure), - ]; Index: sm/trunk-jpl/jenkins/javascript/karma/scripts/test203.js =================================================================== --- /issm/trunk-jpl/jenkins/javascript/karma/scripts/test203.js (revision 21203) +++ (revision ) @@ -1,20 +1,0 @@ -//Test Name: SquareShelfStressHO -var md = new model(); -triangle(md,square[0],180000.); -setmask(md,'all',''); -parameterize(md); -md.extrude(md,3,2.); -setflowequation(md,'HO','all'); -//md.cluster=generic('name',oshostname(),'np',3); -md=solve(md,StressbalanceSolutionEnum()); - -//Fields and tolerances to track changes -field_names =['Vx','Vy','Vz','Vel','Pressure']; -field_tolerances=[2e-09,1e-09,1e-09,1e-09,1e-09]; -field_values=[ - (md.results.StressbalanceSolution[0].Vx), - (md.results.StressbalanceSolution[0].Vy), - (md.results.StressbalanceSolution[0].Vz), - (md.results.StressbalanceSolution[0].Vel), - (md.results.StressbalanceSolution[0].Pressure), - ]; Index: sm/trunk-jpl/jenkins/javascript/karma/scripts/test204.js =================================================================== --- /issm/trunk-jpl/jenkins/javascript/karma/scripts/test204.js (revision 21203) +++ (revision ) @@ -1,29 +1,0 @@ -//Test Name: SquareShelfStressFS -var md = new model(); -triangle(md,square[0],180000.); -setmask(md,'all',''); -parameterize(md); -md.extrude(md,3,2.); -setflowequation(md,'FS','all'); -//md.cluster=generic('name',oshostname(),'np',1); -md.stressbalance.shelf_dampening=1; -md.timestepping.time_step=0; -md1=solve(md,StressbalanceSolutionEnum()); -md.stressbalance.shelf_dampening=0; -md=solve(md,StressbalanceSolutionEnum()); - -//Fields and tolerances to track changes -field_names =['Vx','Vy','Vz','Vel','Pressure','Vx_damp','Vy_damp','Vz_damp','Vel_damp','Pressure_damp']; -field_tolerances=[1e-08,1e-08,2e-06,1e-08,1e-08,1e-08,1e-08,2e-07,1e-08,1e-08]; -field_values=[ - (md.results.StressbalanceSolution[0].Vx), - (md.results.StressbalanceSolution[0].Vy), - (md.results.StressbalanceSolution[0].Vz), - (md.results.StressbalanceSolution[0].Vel), - (md.results.StressbalanceSolution[0].Pressure), - (md1.results.StressbalanceSolution.Vx), - (md1.results.StressbalanceSolution.Vy), - (md1.results.StressbalanceSolution.Vz), - (md1.results.StressbalanceSolution.Vel), - (md1.results.StressbalanceSolution.Pressure), - ]; Index: sm/trunk-jpl/jenkins/javascript/karma/scripts/test205.js =================================================================== --- /issm/trunk-jpl/jenkins/javascript/karma/scripts/test205.js (revision 21203) +++ (revision ) @@ -1,20 +1,0 @@ -//Test Name: SquareShelfStressMHOPenalties -var md = new model(); -triangle(md,square[0],150000.); -setmask(md,'all',''); -parameterize(md); -md.extrude(md,3,2.); -setflowequation(md,'HO','../Exp/SquareHalfRight.exp','fill','SSA','coupling','penalties'); -//md.cluster=generic('name',oshostname(),'np',3); -md=solve(md,StressbalanceSolutionEnum()); - -//Fields and tolerances to track changes -field_names =['Vx','Vy','Vz','Vel','Pressure']; -field_tolerances=[2e-05,2e-05,1e-05,1e-05,1e-05]; -field_values=[ - (md.results.StressbalanceSolution[0].Vx), - (md.results.StressbalanceSolution[0].Vy), - (md.results.StressbalanceSolution[0].Vz), - (md.results.StressbalanceSolution[0].Vel), - (md.results.StressbalanceSolution[0].Pressure), - ]; Index: sm/trunk-jpl/jenkins/javascript/karma/scripts/test206.js =================================================================== --- /issm/trunk-jpl/jenkins/javascript/karma/scripts/test206.js (revision 21203) +++ (revision ) @@ -1,18 +1,0 @@ -//Test Name: SquareShelfTherStea -var md = new model(); -triangle(md,square[0],180000.); -setmask(md,'all',''); -parameterize(md); -md.extrude(md,3,1.); -setflowequation(md,'SSA','all'); -md.timestepping.time_step=0; -//md.cluster=generic('name',oshostname(),'np',3); -md=solve(md,ThermalSolutionEnum()); - -//Fields and tolerances to track changes -field_names =['Temperature','BasalforcingsGroundediceMeltingRate']; -field_tolerances=[1e-13,5e-6]; -field_values=[ - (md.results.ThermalSolution[0].Temperature), - (md.results.ThermalSolution[0].BasalforcingsGroundediceMeltingRate), - ]; Index: sm/trunk-jpl/jenkins/javascript/karma/scripts/test207.js =================================================================== --- /issm/trunk-jpl/jenkins/javascript/karma/scripts/test207.js (revision 21203) +++ (revision ) @@ -1,26 +1,0 @@ -//Test Name: SquareShelfTherTran -var md = new model(); -triangle(md,square[0],180000.); -setmask(md,'all',''); -parameterize(md); -md.extrude(md,3,1.); -setflowequation(md,'SSA','all'); -//md.cluster=generic('name',oshostname(),'np',3); -md.transient.isstressbalance=0; -md.transient.ismasstransport=0; -md.transient.issmb=1; -md.transient.isthermal=1; -md.transient.isgroundingline=0; -md=solve(md,TransientSolutionEnum()); - -//Fields and tolerances to track changes -field_names =['Temperature1','BasalforcingsGroundediceMeltingRate1','Temperature2','BasalforcingsGroundediceMeltingRate2','Temperature3','BasalforcingsGroundediceMeltingRate3']; -field_tolerances=[1e-13,1e-6,1e-13,1e-6,1e-13,1e-6]; -field_values=[ - (md.results.TransientSolution[0](1).Temperature), - (md.results.TransientSolution[0](1).BasalforcingsGroundediceMeltingRate), - (md.results.TransientSolution[0](2).Temperature), - (md.results.TransientSolution[0](2).BasalforcingsGroundediceMeltingRate), - (md.results.TransientSolution[0](3).Temperature), - (md.results.TransientSolution[0](3).BasalforcingsGroundediceMeltingRate), - ]; Index: sm/trunk-jpl/jenkins/javascript/karma/scripts/test208.js =================================================================== --- /issm/trunk-jpl/jenkins/javascript/karma/scripts/test208.js (revision 21203) +++ (revision ) @@ -1,49 +1,0 @@ -//Test Name: SquareShelfTranSSA2d -var md = new model(); -triangle(md,square[0],150000.); -setmask(md,'all',''); -parameterize(md); -setflowequation(md,'SSA','all'); -//md.cluster=generic('name',oshostname(),'np',3); -md.transient.requested_outputs=['default','FloatingArea','GroundedArea','TotalGroundedBmb','TotalFloatingBmb']; -for (var i = 0; i < md.basalforcings.floatingice_melting_rate.length; ++i) { - md.basalforcings.floatingice_melting_rate[i] = 1; -} -md=solve(md,TransientSolutionEnum()); - -//Fields and tolerances to track changes -field_names =['Vx1','Vy1','Vel1','Pressure1','Bed1','Surface1','Thickness1','TotalGroundedBmb1','TotalFloatingBmb1', - 'Vx2','Vy2','Vel2','Pressure2','Bed2','Surface2','Thickness2','TotalGroundedBmb2','TotalFloatingBmb2', - 'Vx3','Vy3','Vel3','Pressure3','Bed3','Surface3','Thickness3','TotalGroundedBmb3','TotalFloatingBmb3']; -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]; -field_values=[ - (md.results.TransientSolution[0](1).Vx), - (md.results.TransientSolution[0](1).Vy), - (md.results.TransientSolution[0](1).Vel), - (md.results.TransientSolution[0](1).Pressure), - (md.results.TransientSolution[0](1).Base), - (md.results.TransientSolution[0](1).Surface), - (md.results.TransientSolution[0](1).Thickness), - (md.results.TransientSolution[0](1).TotalGroundedBmb), - (md.results.TransientSolution[0](1).TotalFloatingBmb), - (md.results.TransientSolution[0](2).Vx), - (md.results.TransientSolution[0](2).Vy), - (md.results.TransientSolution[0](2).Vel), - (md.results.TransientSolution[0](2).Pressure), - (md.results.TransientSolution[0](2).Base), - (md.results.TransientSolution[0](2).Surface), - (md.results.TransientSolution[0](2).Thickness), - (md.results.TransientSolution[0](2).TotalGroundedBmb), - (md.results.TransientSolution[0](2).TotalFloatingBmb), - (md.results.TransientSolution[0](3).Vx), - (md.results.TransientSolution[0](3).Vy), - (md.results.TransientSolution[0](3).Vel), - (md.results.TransientSolution[0](3).Pressure), - (md.results.TransientSolution[0](3).Base), - (md.results.TransientSolution[0](3).Surface), - (md.results.TransientSolution[0](3).Thickness), - (md.results.TransientSolution[0](3).TotalGroundedBmb), - (md.results.TransientSolution[0](3).TotalFloatingBmb), - ]; Index: sm/trunk-jpl/jenkins/javascript/karma/scripts/test209.js =================================================================== --- /issm/trunk-jpl/jenkins/javascript/karma/scripts/test209.js (revision 21203) +++ (revision ) @@ -1,50 +1,0 @@ -//Test Name: SquareShelfTranSSA3d -var md = new model(); -triangle(md,square[0],180000.); -setmask(md,'all',''); -parameterize(md); -md.extrude(md,3,1.); -setflowequation(md,'SSA','all'); -//md.cluster=generic('name',oshostname(),'np',3); -md=solve(md,TransientSolutionEnum()); - -//Fields and tolerances to track changes -field_names =['Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1','BasalforcingsGroundediceMeltingRate1', - 'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2','BasalforcingsGroundediceMeltingRate2', - 'Vx3','Vy3','Vz3','Vel3','Pressure3','Bed3','Surface3','Thickness3','Temperature3','BasalforcingsGroundediceMeltingRate3']; -field_tolerances=[ - 1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-8, - 1e-13,1e-13,3e-9,1e-13,3.8e-11,3.85e-11,3.85e-11,3.85e-11,1e-13,2e-8, - 5e-11,6e-12,1e-08,6e-12,3.9e-11,1e-10,1e-10,1e-10,9e-12,5e-8]; -field_values=[ - (md.results.TransientSolution[0](1).Vx), - (md.results.TransientSolution[0](1).Vy), - (md.results.TransientSolution[0](1).Vz), - (md.results.TransientSolution[0](1).Vel), - (md.results.TransientSolution[0](1).Pressure), - (md.results.TransientSolution[0](1).Base), - (md.results.TransientSolution[0](1).Surface), - (md.results.TransientSolution[0](1).Thickness), - (md.results.TransientSolution[0](1).Temperature), - (md.results.TransientSolution[0](1).BasalforcingsGroundediceMeltingRate), - (md.results.TransientSolution[0](2).Vx), - (md.results.TransientSolution[0](2).Vy), - (md.results.TransientSolution[0](2).Vz), - (md.results.TransientSolution[0](2).Vel), - (md.results.TransientSolution[0](2).Pressure), - (md.results.TransientSolution[0](2).Base), - (md.results.TransientSolution[0](2).Surface), - (md.results.TransientSolution[0](2).Thickness), - (md.results.TransientSolution[0](2).Temperature), - (md.results.TransientSolution[0](2).BasalforcingsGroundediceMeltingRate), - (md.results.TransientSolution[0](3).Vx), - (md.results.TransientSolution[0](3).Vy), - (md.results.TransientSolution[0](3).Vz), - (md.results.TransientSolution[0](3).Vel), - (md.results.TransientSolution[0](3).Pressure), - (md.results.TransientSolution[0](3).Base), - (md.results.TransientSolution[0](3).Surface), - (md.results.TransientSolution[0](3).Thickness), - (md.results.TransientSolution[0](3).Temperature), - (md.results.TransientSolution[0](3).BasalforcingsGroundediceMeltingRate), - ]; Index: /issm/trunk-jpl/jenkins/javascript/karma/scripts/test231.js =================================================================== --- /issm/trunk-jpl/jenkins/javascript/karma/scripts/test231.js (revision 21204) +++ /issm/trunk-jpl/jenkins/javascript/karma/scripts/test231.js (revision 21204) @@ -0,0 +1,73 @@ +//Test Name: SquareShelfTranForcePos3d +var md = new model(); +triangle(md,square[0],350000.); +setmask(md,'all',''); +parameterize(md); +md.extrude(md,3,1.); +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); +for (var i = 0; i < smb.length; ++i) { + smb[i][0] *= 3.6; + smb[i].push(smb[i][0]*2); +} + +md.smb.mass_balance=smb.slice(); +md.smb.mass_balance[md.smb.mass_balance.length-1] = [1.5, 3.]; +md.trans.isthermal=0; + +md=solve(md,'Transient'); + +//Fields and tolerances to track changes +field_names =['Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','SmbMassBalance1', + 'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','SmbMassBalance2', + 'Vx3','Vy3','Vz3','Vel3','Pressure3','Bed3','Surface3','Thickness3','SmbMassBalance3', + 'Vx4','Vy4','Vz4','Vel4','Pressure4','Bed4','Surface4','Thickness4','SmbMassbalance4']; +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](1).Vx), + (md.results.TransientSolution[0](1).Vy), + (md.results.TransientSolution[0](1).Vz), + (md.results.TransientSolution[0](1).Vel), + (md.results.TransientSolution[0](1).Pressure), + (md.results.TransientSolution[0](1).Base), + (md.results.TransientSolution[0](1).Surface), + (md.results.TransientSolution[0](1).Thickness), + (md.results.TransientSolution[0](1).SmbMassBalance), + (md.results.TransientSolution[0](2).Vx), + (md.results.TransientSolution[0](2).Vy), + (md.results.TransientSolution[0](2).Vz), + (md.results.TransientSolution[0](2).Vel), + (md.results.TransientSolution[0](2).Pressure), + (md.results.TransientSolution[0](2).Base), + (md.results.TransientSolution[0](2).Surface), + (md.results.TransientSolution[0](2).Thickness), + (md.results.TransientSolution[0](2).SmbMassBalance), + (md.results.TransientSolution[0](3).Vx), + (md.results.TransientSolution[0](3).Vy), + (md.results.TransientSolution[0](3).Vz), + (md.results.TransientSolution[0](3).Vel), + (md.results.TransientSolution[0](3).Pressure), + (md.results.TransientSolution[0](3).Base), + (md.results.TransientSolution[0](3).Surface), + (md.results.TransientSolution[0](3).Thickness), + (md.results.TransientSolution[0](3).SmbMassBalance), + (md.results.TransientSolution[0](4).Vx), + (md.results.TransientSolution[0](4).Vy), + (md.results.TransientSolution[0](4).Vz), + (md.results.TransientSolution[0](4).Vel), + (md.results.TransientSolution[0](4).Pressure), + (md.results.TransientSolution[0](4).Base), + (md.results.TransientSolution[0](4).Surface), + (md.results.TransientSolution[0](4).Thickness), + (md.results.TransientSolution[0](4).SmbMassBalance), + ]; Index: /issm/trunk-jpl/jenkins/javascript/karma/scripts/test232.js =================================================================== --- /issm/trunk-jpl/jenkins/javascript/karma/scripts/test232.js (revision 21204) +++ /issm/trunk-jpl/jenkins/javascript/karma/scripts/test232.js (revision 21204) @@ -0,0 +1,31 @@ +//Test Name: SquareShelfTherTranForcTemp +var md = new model(); +triangle(md,square[0],180000.); +setmask(md,'all',''); +parameterize(md); +md.extrude(md,3,1.); +setflowequation(md,'SSA','all'); +//md.cluster=generic('name',oshostname(),'np',3); +md.thermal.spctemperature=[[md.thermal.spctemperature, md.thermal.spctemperature+5., md.thermal.spctemperature+10., md.thermal.spctemperature+15.], [1.5, 2.5, 3.5, 4.]]; +md.timestepping.time_step=1; +md.timestepping.final_time=4; +md.trans.isstressbalance=0; +md.trans.ismasstransport=0; +md.trans.issmb=1; +md.trans.isthermal=1; +md.trans.isgroundingline=0; +md=solve(md,'Transient'); + +//Fields and tolerances to track changes +field_names =['Temperature1','BasalforcingsGroundediceMeltingRate1','Temperature2','BasalforcingsGroundediceMeltingRate2','Temperature3','BasalforcingsGroundediceMeltingRate3','Temperature4','BasalforcingsGroundediceMeltingRate4']; +field_tolerances=[1e-13,1e-6,1e-13,1e-6,1e-13,1e-6,1e-13,1e-6]; +field_values=[ + (md.results.TransientSolution[0](1).Temperature), + (md.results.TransientSolution[0](1).BasalforcingsGroundediceMeltingRate), + (md.results.TransientSolution[0](2).Temperature), + (md.results.TransientSolution[0](2).BasalforcingsGroundediceMeltingRate), + (md.results.TransientSolution[0](3).Temperature), + (md.results.TransientSolution[0](3).BasalforcingsGroundediceMeltingRate), + (md.results.TransientSolution[0](4).Temperature), + (md.results.TransientSolution[0](4).BasalforcingsGroundediceMeltingRate), + ]; Index: /issm/trunk-jpl/jenkins/javascript/karma/scripts/test233.js =================================================================== --- /issm/trunk-jpl/jenkins/javascript/karma/scripts/test233.js (revision 21204) +++ /issm/trunk-jpl/jenkins/javascript/karma/scripts/test233.js (revision 21204) @@ -0,0 +1,64 @@ +//Test Name: SquareShelfTranHOForcTemp +var md = new model(); +triangle(md,square[0],200000.); +setmask(md,'all',''); +parameterize(md); +md.extrude(md,3,1.); +setflowequation(md,'HO','all'); +//md.cluster=generic('name',oshostname(),'np',3); +md.thermal.spctemperature=[[md.thermal.spctemperature, md.thermal.spctemperature+5.], [1., 2.]]; +md.timestepping.time_step=0.5; +md.timestepping.final_time=2.; +md=solve(md,'Transient'); + +//Fields and tolerances to track changes +field_names =['Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1','BasalforcingsGroundediceMeltingRate1', + 'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2','BasalforcingsGroundediceMeltingRate2', + 'Vx3','Vy3','Vz3','Vel3','Pressure3','Bed3','Surface3','Thickness3','Temperature3','BasalforcingsGroundediceMeltingRate3', + 'Vx4','Vy4','Vz4','Vel4','Pressure4','Bed4','Surface4','Thickness4','Temperature4','BasalforcingsGroundediceMeltingRate4']; +field_tolerances=[1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09, + 1e-09,1e-09,1e-08,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-06, + 1e-08,1e-09,1e-08,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-06, + 1e-09,1e-09,1e-08,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-06]; +field_values=[ + (md.results.TransientSolution[0](1).Vx), + (md.results.TransientSolution[0](1).Vy), + (md.results.TransientSolution[0](1).Vz), + (md.results.TransientSolution[0](1).Vel), + (md.results.TransientSolution[0](1).Pressure), + (md.results.TransientSolution[0](1).Base), + (md.results.TransientSolution[0](1).Surface), + (md.results.TransientSolution[0](1).Thickness), + (md.results.TransientSolution[0](1).Temperature), + (md.results.TransientSolution[0](1).BasalforcingsGroundediceMeltingRate), + (md.results.TransientSolution[0](2).Vx), + (md.results.TransientSolution[0](2).Vy), + (md.results.TransientSolution[0](2).Vz), + (md.results.TransientSolution[0](2).Vel), + (md.results.TransientSolution[0](2).Pressure), + (md.results.TransientSolution[0](2).Base), + (md.results.TransientSolution[0](2).Surface), + (md.results.TransientSolution[0](2).Thickness), + (md.results.TransientSolution[0](2).Temperature), + (md.results.TransientSolution[0](2).BasalforcingsGroundediceMeltingRate), + (md.results.TransientSolution[0](3).Vx), + (md.results.TransientSolution[0](3).Vy), + (md.results.TransientSolution[0](3).Vz), + (md.results.TransientSolution[0](3).Vel), + (md.results.TransientSolution[0](3).Pressure), + (md.results.TransientSolution[0](3).Base), + (md.results.TransientSolution[0](3).Surface), + (md.results.TransientSolution[0](3).Thickness), + (md.results.TransientSolution[0](3).Temperature), + (md.results.TransientSolution[0](3).BasalforcingsGroundediceMeltingRate), + (md.results.TransientSolution[0](4).Vx), + (md.results.TransientSolution[0](4).Vy), + (md.results.TransientSolution[0](4).Vz), + (md.results.TransientSolution[0](4).Vel), + (md.results.TransientSolution[0](4).Pressure), + (md.results.TransientSolution[0](4).Base), + (md.results.TransientSolution[0](4).Surface), + (md.results.TransientSolution[0](4).Thickness), + (md.results.TransientSolution[0](4).Temperature), + (md.results.TransientSolution[0](4).BasalforcingsGroundediceMeltingRate), + ]; Index: /issm/trunk-jpl/jenkins/javascript/karma/scripts/test234.js =================================================================== --- /issm/trunk-jpl/jenkins/javascript/karma/scripts/test234.js (revision 21204) +++ /issm/trunk-jpl/jenkins/javascript/karma/scripts/test234.js (revision 21204) @@ -0,0 +1,93 @@ +//Test Name: SquareShelfTranForceNeg2dDakotaSamp +var md = new model(); +triangle(md,square[0],180000.); +setmask(md,'all',''); +parameterize(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; + +smb=ones(md.mesh.numberofvertices,1); +for (var i = 0; i < smb.length; ++i) { + smb[i][0] *= 3.6; + smb[i].push(smb[i][0]*-1); +} + +md.smb.mass_balance= smb.slice(); +md.smb.mass_balance[md.smb.mass_balance.length-1] = [1.5, 3.]; +md.trans.isthermal=0; +//Dakota options + +//dakota version +version=IssmConfig('_DAKOTA_VERSION_'); version=version.toString().slice(0,2); + +//partitioning +md.qmu.numberofpartitions=20; +partitioner(md,'package','chaco','npart',md.qmu.numberofpartitions,'weighting','on'); +md.qmu.partition=md.qmu.partition-1; + +//variables +md.qmu.variables.surface_mass_balance=normal_uncertain('scaled_SmbMassBalance',1,0.1); + +//responses +md.qmu.responses.MaxVel=response_function('MaxVel',[],[0.0001, 0.001, 0.01, 0.25, 0.5, 0.75, 0.99, 0.999, 0.9999]); +md.qmu.responses.IceVolume=response_function('IceVolume',[],[0.0001, 0.001, 0.01, 0.25, 0.5, 0.75, 0.99, 0.999, 0.9999]); +md.qmu.responses.MassFlux1=response_function('indexed_MassFlux_1',[],[0.0001, 0.001, 0.01, 0.25, 0.5, 0.75, 0.99, 0.999, 0.9999]); +md.qmu.responses.MassFlux2=response_function('indexed_MassFlux_2',[],[0.0001, 0.001, 0.01, 0.25, 0.5, 0.75, 0.99, 0.999, 0.9999]); +md.qmu.responses.MassFlux3=response_function('indexed_MassFlux_3',[],[0.0001, 0.001, 0.01, 0.25, 0.5, 0.75, 0.99, 0.999, 0.9999]); +md.qmu.responses.MassFlux4=response_function('indexed_MassFlux_4',[],[0.0001, 0.001, 0.01, 0.25, 0.5, 0.75, 0.99, 0.999, 0.9999]); +md.qmu.responses.MassFlux5=response_function('indexed_MassFlux_5',[],[0.0001, 0.001, 0.01, 0.25, 0.5, 0.75, 0.99, 0.999, 0.9999]); +md.qmu.responses.massFlux6=response_function('indexed_MassFlux_6',[],[0.0001, 0.001, 0.01, 0.25, 0.5, 0.75, 0.99, 0.999, 0.9999]); + +//mass flux profiles +md.qmu.mass_flux_profiles=['../Exp/MassFlux1.exp','../Exp/MassFlux2.exp','../Exp/MassFlux3.exp','../Exp/MassFlux4.exp','../Exp/MassFlux5.exp','../Exp/MassFlux6.exp']; +md.qmu.mass_flux_profile_directory=pwd; + +//// nond_sampling study +md.qmu.method=dakota_method('nond_samp'); +md.qmu.method[md.qmu.method.length-1]=dmeth_params_set(md.qmu.method[md.qmu.method.length-1],'seed',1234,'samples',20,'sample_type','lhs'); +dver=textscan(IssmConfig('_DAKOTA_VERSION_'),'//[0123456789].//[0123456789].//[0123456789]'); +if (((str2num(dver[1][1])==4 && str2num(dver[2][1])>2) || str2num(dver[1][1])>4)) { + md.qmu.method[md.qmu.method.length-1]=dmeth_params_set(md.qmu.method(end),'rng','rnum2'); +} + +//parameters +md.qmu.params.direct=true; +md.qmu.params.analysis_components=''; +md.qmu.params.interval_type='forward'; +md.qmu.params.tabular_graphics_data=true; +md.qmu.isdakota=1; + +if (version>=6) { + md.qmu.params.analysis_driver='matlab'; + md.qmu.params.evaluation_scheduling='master'; + md.qmu.params.processors_per_evaluation=2; +} else { + md.qmu.params.analysis_driver='stressbalance'; + md.qmu.params.evaluation_concurrency=1; +} + + +md.stressbalance.reltol=Math.pow(10,-5); //tighten for qmu analyses +md.trans.requested_outputs=['IceVolume']; + +//solve +md=solve(md,'Transient','overwrite','y'); +md.qmu.results=md.results.dakota; + +//Fields and tolerances to track changes +md.results.dakota.moments=[]; +for (var i = 0; i < 8; ++i) { + md.results.dakota.moments.push(md.results.dakota.dresp_out[i].mean); +} +for (var i = 0; i < 8; ++i) { + md.results.dakota.moments.push(md.results.dakota.dresp_out[i].stddev); +} +field_names =['moments']; +field_tolerances=[1e-11]; +field_values=[ + md.results.dakota.moments, + ];