source: issm/trunk-jpl/test/NightlyRun/test3300.m@ 21056

Last change on this file since 21056 was 21056, checked in by Mathieu Morlighem, 9 years ago

CHG: do not request solution in solution string

File size: 2.1 KB
RevLine 
[17376]1md=triangle(model(),'../Exp/Square.exp',100000.);
2md=setmask(md,'','');
3md=parameterize(md,'../Par/SquareSheetConstrained.par');
[17734]4md.transient=deactivateall(md.transient);
5md.transient.ishydrology=1;
[17376]6md=setflowequation(md,'SSA','all');
7md.cluster=generic('name',oshostname(),'np',1);
8md.hydrology=(hydrologydc);
[17422]9md.hydrology=initialize(md.hydrology,md);
[17376]10md.hydrology.isefficientlayer=1;
11md.hydrology.sedimentlimit_flag=1;
12md.hydrology.sedimentlimit=400.0;
13md.hydrology.sediment_thickness=20.0;
14md.initialization.sediment_head=0.0*ones(md.mesh.numberofvertices,1);
15md.hydrology.spcsediment_head=NaN*ones(md.mesh.numberofvertices,1);
[18582]16md.basalforcings.groundedice_melting_rate = 2.0*ones(md.mesh.numberofvertices,1);
[17376]17md.hydrology.sediment_transmitivity=1.5e-4*ones(md.mesh.numberofvertices,1);
18
19md.initialization.epl_head=0.0*ones(md.mesh.numberofvertices,1);
[17524]20md.initialization.epl_thickness=1.0*ones(md.mesh.numberofvertices,1);
[17376]21md.hydrology.spcepl_head=NaN*ones(md.mesh.numberofvertices,1);
22md.hydrology.mask_eplactive_node=0*ones(md.mesh.numberofvertices,1);
23
[18974]24md.hydrology.epl_conductivity=1.5e-2;
[17376]25md.hydrology.epl_initial_thickness=1.0;
[18974]26md.hydrology.epl_colapse_thickness=1.0e-6;
[18983]27md.hydrology.epl_thick_comp=1;
[17738]28md.hydrology.epl_max_thickness=5.0;
[17376]29
[18646]30md.hydrology.transfer_flag=1.0;
[18974]31md.hydrology.leakage_factor=3.9e-12;
[17376]32
[19031]33times=0:0.002:8.0;
[18582]34md.basalforcings.groundedice_melting_rate=ones(md.mesh.numberofvertices+1,length(times));
[17376]35
[18974]36for i=1:length(times)
37 if(times(i)<=1.0)
38 md.basalforcings.groundedice_melting_rate(:,i)=1.0;
39 elseif(times(i)<=6.0)
40 md.basalforcings.groundedice_melting_rate(:,i)=-0.2;
41 else
42 md.basalforcings.groundedice_melting_rate(:,i)=0.0;
43 end
44end
45
[18582]46md.basalforcings.groundedice_melting_rate(end,:)=times;
[17376]47
[19031]48md.timestepping.time_step=0.002;
[18974]49md.timestepping.final_time=8.0;
[17376]50
[21056]51md=solve(md,'Transient');
[17376]52
53field_names ={'SedimentWaterHead5','EplWaterHead5','SedimentWaterHead40','EplWaterHead40'};
54field_tolerances={1e-13, 1e-13, 1e-13, 1e-13, 1e-13};
[17734]55field_values={md.results.TransientSolution(5).SedimentHead, ...
56 md.results.TransientSolution(5).EplHead,...
57 md.results.TransientSolution(40).SedimentHead,...
58 md.results.TransientSolution(40).EplHead};
Note: See TracBrowser for help on using the repository browser.