Index: /issm/trunk-jpl/test/NightlyRun/test230.js =================================================================== --- /issm/trunk-jpl/test/NightlyRun/test230.js (revision 21215) +++ /issm/trunk-jpl/test/NightlyRun/test230.js (revision 21215) @@ -0,0 +1,73 @@ +//Test Name: SquareShelfTranForceNeg3d +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]*-1); +} + +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/test/NightlyRun/test231.js =================================================================== --- /issm/trunk-jpl/test/NightlyRun/test231.js (revision 21215) +++ /issm/trunk-jpl/test/NightlyRun/test231.js (revision 21215) @@ -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/test/NightlyRun/test232.js =================================================================== --- /issm/trunk-jpl/test/NightlyRun/test232.js (revision 21215) +++ /issm/trunk-jpl/test/NightlyRun/test232.js (revision 21215) @@ -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/test/NightlyRun/test233.js =================================================================== --- /issm/trunk-jpl/test/NightlyRun/test233.js (revision 21215) +++ /issm/trunk-jpl/test/NightlyRun/test233.js (revision 21215) @@ -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/test/NightlyRun/test234.js =================================================================== --- /issm/trunk-jpl/test/NightlyRun/test234.js (revision 21215) +++ /issm/trunk-jpl/test/NightlyRun/test234.js (revision 21215) @@ -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, + ]; Index: /issm/trunk-jpl/test/NightlyRun/test235.js =================================================================== --- /issm/trunk-jpl/test/NightlyRun/test235.js (revision 21215) +++ /issm/trunk-jpl/test/NightlyRun/test235.js (revision 21215) @@ -0,0 +1,88 @@ +//Test Name: SquareShelfTranForceNeg2dDakotaLocal +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,100); + +//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; + +//method +md.qmu.method =dakota_method('nond_l'); + +//parameters +md.qmu.params.direct=true; +md.qmu.params.analysis_components=''; +md.qmu.params.interval_type='forward'; +md.qmu.params.fd_gradient_step_size='0.1'; +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, + ];