Changeset 10658


Ignore:
Timestamp:
11/15/11 13:15:01 (13 years ago)
Author:
schlegel
Message:

newest dak test

File:
1 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk/test/NightlyRun/dak.m

    r10655 r10658  
    1212smb = ones(md.mesh.numberofvertices,1)*3.6;
    1313smb=[ smb smb*-1 ];
    14 %md.diagnostic.spcvx=ones(md.mesh.numberofvertices,1);
    15 %md.diagnostic.spcvy=ones(md.mesh.numberofvertices,1);
    1614
    1715md.surfaceforcings.mass_balance=smb;
     
    2624%Dakota inputs,outputs
    2725md.qmu.variables.surface_mass_balance=normal_uncertain('scaled_SurfaceforcingsMassBalance',1,0.1);
    28 %md.qmu.variables.drag_coefficient=normal_uncertain('scaled_FrictionCoefficient',1,0.5);
    29 %md.qmu.variables.rheology_B=normal_uncertain('scaled_MaterialsRheologyB',1,.5);
    3026md.qmu.responses.IceVolume=response_function('IceVolume',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
    3127md.qmu.responses.MaxVel=response_function('MaxVel',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
     
    4440
    4541%partitioning
    46 md.qmu.numberofpartitions=1;
     42md.qmu.numberofpartitions=20;
    4743md=partitioner(md,'package','chaco','npart',md.qmu.numberofpartitions,'weighting','on');
    4844md.qmu.partition=md.qmu.partition-1;
     
    5248
    5349md.diagnostic.reltol=10^-5; %tighten for qmu analyses
    54 md.debug.valgrind=1;
    55 %md.verbose=verbose('1111111');
     50%md.debug.valgrind=1;
    5651%md.verbose.solution=true;
    5752md=solve(md,TransientSolutionEnum,'overwrite','y');
Note: See TracChangeset for help on using the changeset viewer.