Changeset 16609
- Timestamp:
- 11/04/13 09:39:44 (11 years ago)
- Location:
- issm/trunk-jpl/test
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/test/NightlyRun/test1502.m
r16523 r16609 1 1 printingflag = false; 2 2 3 md=triangle(model(),'../Exp/Square.exp', 450000.);3 md=triangle(model(),'../Exp/Square.exp',350000.); 4 4 md=setmask(md,'all',''); 5 5 md=parameterize(md,'../Par/SquareShelf.par'); 6 6 md=setflowequation(md,'SSA','all'); 7 %md=extrude(md,3,1.);7 md=extrude(md,3,1.); 8 8 md.cluster=generic('name',oshostname(),'np',1); 9 9 md.transient.isthermal=0; 10 10 11 %md.timestepping.time_step=1.;12 %md.settings.output_frequency=1;13 %md.timestepping.final_time=2000.;11 md.timestepping.time_step=1.; 12 md.settings.output_frequency=1; 13 md.timestepping.final_time=2000.; 14 14 15 15 %Solve for thinning rate -> -1 * surface mass balance 16 16 smb= 2.*ones(md.mesh.numberofvertices,1); 17 %md.surfaceforcings.mass_balance= smb;18 %md.basalforcings.melting_rate= smb;17 md.surfaceforcings.mass_balance= smb; 18 md.basalforcings.melting_rate= smb; 19 19 md.verbose=verbose('all'); 20 20 21 21 md=solve(md,MasstransportSolutionEnum()); 22 22 23 error24 23 for i=1:10 25 24 md=solve(md,MasstransportSolutionEnum());
Note:
See TracChangeset
for help on using the changeset viewer.