Changeset 10658
- Timestamp:
- 11/15/11 13:15:01 (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk/test/NightlyRun/dak.m
r10655 r10658 12 12 smb = ones(md.mesh.numberofvertices,1)*3.6; 13 13 smb=[ smb smb*-1 ]; 14 %md.diagnostic.spcvx=ones(md.mesh.numberofvertices,1);15 %md.diagnostic.spcvy=ones(md.mesh.numberofvertices,1);16 14 17 15 md.surfaceforcings.mass_balance=smb; … … 26 24 %Dakota inputs,outputs 27 25 md.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);30 26 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]); 31 27 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]); … … 44 40 45 41 %partitioning 46 md.qmu.numberofpartitions= 1;42 md.qmu.numberofpartitions=20; 47 43 md=partitioner(md,'package','chaco','npart',md.qmu.numberofpartitions,'weighting','on'); 48 44 md.qmu.partition=md.qmu.partition-1; … … 52 48 53 49 md.diagnostic.reltol=10^-5; %tighten for qmu analyses 54 md.debug.valgrind=1; 55 %md.verbose=verbose('1111111'); 50 %md.debug.valgrind=1; 56 51 %md.verbose.solution=true; 57 52 md=solve(md,TransientSolutionEnum,'overwrite','y');
Note:
See TracChangeset
for help on using the changeset viewer.