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

Last change on this file since 21049 was 21049, checked in by agscott1, 9 years ago

CHG: Replaced Enums with Strings in matlab and python. Updated corresponding cpp code.

File size: 2.2 KB
Line 
1md=triangle(model(),'../Exp/Square.exp',100000.);
2md=setmask(md,'','');
3md=parameterize(md,'../Par/SquareSheetConstrained.par');
4md.transient=deactivateall(md.transient);
5md.transient.ishydrology=1;
6md=setflowequation(md,'SSA','all');
7md.cluster=generic('name',oshostname(),'np',1);
8md.hydrology=(hydrologydc);
9md.hydrology=initialize(md.hydrology,md);
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);
16md.basalforcings.groundedice_melting_rate = 2.0*ones(md.mesh.numberofvertices,1);
17md.hydrology.sediment_transmitivity=1.5e-4*ones(md.mesh.numberofvertices,1);
18
19md.initialization.epl_head=0.0*ones(md.mesh.numberofvertices,1);
20md.initialization.epl_thickness=1.0*ones(md.mesh.numberofvertices,1);
21md.hydrology.spcepl_head=NaN*ones(md.mesh.numberofvertices,1);
22md.hydrology.mask_eplactive_node=0*ones(md.mesh.numberofvertices,1);
23
24md.hydrology.epl_conductivity=1.5e-2;
25md.hydrology.epl_initial_thickness=1.0;
26md.hydrology.epl_colapse_thickness=1.0e-6;
27md.hydrology.epl_thick_comp=1;
28md.hydrology.epl_max_thickness=5.0;
29
30md.hydrology.transfer_flag=1.0;
31md.hydrology.leakage_factor=3.9e-12;
32
33times=0:0.002:8.0;
34md.basalforcings.groundedice_melting_rate=ones(md.mesh.numberofvertices+1,length(times));
35
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
46md.basalforcings.groundedice_melting_rate(end,:)=times;
47
48md.timestepping.time_step=0.002;
49md.timestepping.final_time=8.0;
50
51md=solve(md,'TransientSolution');
52
53field_names ={'SedimentWaterHead5','EplWaterHead5','SedimentWaterHead40','EplWaterHead40'};
54field_tolerances={1e-13, 1e-13, 1e-13, 1e-13, 1e-13};
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.