Changeset 8011 for issm/trunk/test/Miscellaneous/Bump/Bump.par
- Timestamp:
- 04/21/11 11:27:33 (14 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk/test/Miscellaneous/Bump/Bump.par
r3768 r8011 1 %Start defining model parameters here2 3 1 %Choose the radius and heigth of the bump: 4 2 height=100; 5 3 6 4 %dynamics 7 md.debug=2;8 md.dt=0.1*md.yts; %1 year9 md.ndt=md.dt*3;10 md.artificial_diffusivity=1;11 12 5 hmin=1000; 13 6 hmax=1000; … … 17 10 md.bed=-md.rho_ice/md.rho_water*md.thickness+10; 18 11 md.surface=md.bed+md.thickness; 19 md.bed=-md.rho_ice/md.rho_water*md.thickness+height*exp(-((md.x-50000).^2+(md.y-50000).^2)/( 2000)^2)+10;12 md.bed=-md.rho_ice/md.rho_water*md.thickness+height*exp(-((md.x-50000).^2+(md.y-50000).^2)/(4000)^2)+10; 20 13 md.thickness=md.surface-md.bed; 21 14 … … 39 32 40 33 %Parallel options 41 md.cluster='larsen';42 md.np=8;43 md.time=50;44 34 md.waitonlock=1; 45 35 md.connectivity=100; 46 47 %Control method options48 md.nsteps=15;49 md.fit=[2*ones(10,1);0*ones(5,1)];50 md.optscal=1000*ones(md.nsteps,1);51 md.maxiter=20*ones(md.nsteps,1);52 md.maxcontrolconstraint=500;53 md=solversettomumps(md);
Note:
See TracChangeset
for help on using the changeset viewer.