source: issm/trunk-jpl/test/NightlyRun/test902.m

Last change on this file was 24570, checked in by bdef, 5 years ago

BUG:NR fix
File size: 3.2 KB
Line 
1%Test Name: SquareSheetHydrologyDCTwoLayers
2md=triangle(model(),'../Exp/Square.exp',100000.);
3md=setmask(md,'','');
4md=parameterize(md,'../Par/SquareNoDyn.par');
5md.cluster=generic('name',oshostname(),'np',1);
6
7md.transient.ishydrology=1;
8md.hydrology=(hydrologydc);
9md.hydrology=initialize(md.hydrology,md);
10
11md.hydrology.isefficientlayer=1;
12md.hydrology.sedimentlimit_flag=1;
13md.hydrology.sedimentlimit=800.0;
14md.hydrology.transfer_flag = 0;
15md.hydrology.mask_thawed_node=ones(md.mesh.numberofvertices,1);
16md.initialization.sediment_head=0.0*ones(md.mesh.numberofvertices,1);
17md.hydrology.spcsediment_head=NaN*ones(md.mesh.numberofvertices,1);
18
19md.basalforcings.groundedice_melting_rate = 2.0*ones(md.mesh.numberofvertices,1);
20md.basalforcings.floatingice_melting_rate = 0.0*ones(md.mesh.numberofvertices,1);
21md.hydrology.sediment_transmitivity=3*ones(md.mesh.numberofvertices,1);
22
23md.initialization.epl_head=0.0*ones(md.mesh.numberofvertices,1);
24md.initialization.epl_thickness=1.0*ones(md.mesh.numberofvertices,1);
25md.hydrology.spcepl_head=NaN*ones(md.mesh.numberofvertices,1);
26md.hydrology.mask_eplactive_node=0*ones(md.mesh.numberofvertices,1);
27md.hydrology.epl_conductivity=30;
28md.hydrology.epl_initial_thickness=1;
29md.hydrology.epl_colapse_thickness=1.0e-3;
30md.hydrology.epl_thick_comp=1;
31md.hydrology.epl_max_thickness=1;
32md.hydrology.steps_per_step=10;
33md.timestepping.time_step=2.0;
34md.timestepping.final_time=2.0;
35
36md=solve(md,'Transient');
37
38%re-run with no substeps
39mdfine=md;
40mdfine.results=struct();
41mdfine.hydrology.steps_per_step=1;
42mdfine.timestepping.time_step=0.2;
43mdfine=solve(mdfine,'Transient');
44
45%Fields and tolerances to track changes
46field_names ={'SedimentWaterHead1','EplWaterHead1','SedimentHeadResidual1',...
47 'SedimentWaterHead4','EplWaterHead4','SedimentHeadResidual4',...
48 'SedimentWaterHead5','EplWaterHead5','SedimentHeadResidual5',...
49 'SedimentWaterHead9','EplWaterHead9','SedimentHeadResidual9',...
50 'EplWaterHead10', 'EplWaterHeadSubstep10', 'SedimentWaterHead10',...
51 'SedimentWaterHeadSubstep10'};
52field_tolerances={...
53 1e-13, 1e-13, 1e-13,...
54 1e-13, 1e-13, 1e-13,...
55 1e-13, 5e-12, 1e-11,...
56 1e-13, 5e-12, 1e-11,...
57 1e-13, 1e-13, 1e-13,...
58 1e-13};
59field_values={mdfine.results.TransientSolution(1).SedimentHead, ...
60 mdfine.results.TransientSolution(1).EplHead,...
61 mdfine.results.TransientSolution(1).SedimentHeadResidual,...
62 mdfine.results.TransientSolution(4).SedimentHead,...
63 mdfine.results.TransientSolution(4).EplHead,...
64 mdfine.results.TransientSolution(4).SedimentHeadResidual, ...
65 mdfine.results.TransientSolution(5).SedimentHead,...
66 mdfine.results.TransientSolution(5).EplHead,...
67 mdfine.results.TransientSolution(5).SedimentHeadResidual, ...
68 mdfine.results.TransientSolution(9).SedimentHead,...
69 mdfine.results.TransientSolution(9).EplHead,...
70 mdfine.results.TransientSolution(9).SedimentHeadResidual,...
71 md.results.TransientSolution(1).EplHead,...
72 md.results.TransientSolution(1).EplHeadSubstep,...
73 md.results.TransientSolution(1).SedimentHead,...
74 md.results.TransientSolution(1).SedimentHeadSubstep
75 };
Note: See TracBrowser for help on using the repository browser.