1 | #Test Name: SquareSheetHydrologyDCTwoLayers
|
---|
2 | import numpy
|
---|
3 | from model import *
|
---|
4 | from triangle import *
|
---|
5 | from setmask import *
|
---|
6 | from parameterize import *
|
---|
7 | from transient import *
|
---|
8 | from setflowequation import *
|
---|
9 | from EnumDefinitions import *
|
---|
10 | from solve import *
|
---|
11 | from MatlabFuncs import *
|
---|
12 |
|
---|
13 | from generic import generic
|
---|
14 |
|
---|
15 | md=triangle(model(),'../Exp/Square.exp',100000.)
|
---|
16 | md=setmask(md,'','')
|
---|
17 | md=parameterize(md,'../Par/SquareSheetConstrained.py')
|
---|
18 | md.transient=transient.setallnullparameters(md.transient)
|
---|
19 | md.transient.ishydrology=True
|
---|
20 | md.transient.issmb=True
|
---|
21 | md=setflowequation(md,'SSA','all')
|
---|
22 | md.cluster=generic('name',oshostname(),'np',1)
|
---|
23 | md.hydrology=hydrologydc()
|
---|
24 | md.hydrology=md.hydrology.initialize(md)
|
---|
25 |
|
---|
26 | md.hydrology.isefficientlayer=1
|
---|
27 | md.hydrology.sedimentlimit_flag=1
|
---|
28 | md.hydrology.sedimentlimit=800.0
|
---|
29 | md.initialization.sediment_head=numpy.zeros((md.mesh.numberofvertices,1))
|
---|
30 | md.hydrology.spcsediment_head=float('NaN')*numpy.ones((md.mesh.numberofvertices,1))
|
---|
31 |
|
---|
32 | md.basalforcings.groundedice_melting_rate = 2.0*numpy.ones((md.mesh.numberofvertices,1))
|
---|
33 | md.basalforcings.floatingice_melting_rate = 0.0*numpy.ones((md.mesh.numberofvertices,1))
|
---|
34 | md.hydrology.sediment_transmitivity= 3.0*numpy.ones((md.mesh.numberofvertices,1))
|
---|
35 |
|
---|
36 | md.initialization.epl_head=numpy.zeros((md.mesh.numberofvertices,1))
|
---|
37 | md.initialization.epl_thickness=numpy.ones((md.mesh.numberofvertices,1))
|
---|
38 | md.hydrology.spcepl_head=float('NaN')*numpy.ones((md.mesh.numberofvertices,1))
|
---|
39 | md.hydrology.mask_eplactive_node=numpy.zeros((md.mesh.numberofvertices,1))
|
---|
40 | md.hydrology.epl_conductivity=30
|
---|
41 | md.hydrology.epl_initial_thickness=1
|
---|
42 | md.hydrology.epl_colapse_thickness=1.0e-3
|
---|
43 | md.hydrology.epl_thick_comp=1
|
---|
44 | md.hydrology.epl_max_thickness=1
|
---|
45 | md.timestepping.time_step=0.2
|
---|
46 | md.timestepping.final_time=2.0
|
---|
47 |
|
---|
48 | md=solve(md,'TransientSolution')
|
---|
49 |
|
---|
50 | field_names =['SedimentWaterHead1','EplWaterHead1','SedimentHeadResidual1', \
|
---|
51 | 'SedimentWaterHead4','EplWaterHead4','SedimentHeadResidual4', \
|
---|
52 | 'SedimentWaterHead5','EplWaterHead5','SedimentHeadResidual5', \
|
---|
53 | 'SedimentWaterHead9','EplWaterHead9','SedimentHeadResidual9']
|
---|
54 | field_tolerances=[ \
|
---|
55 | 1e-13, 1e-13, 1e-13, \
|
---|
56 | 1e-13, 1e-13, 1e-13, \
|
---|
57 | 1e-13, 5e-12, 1e-11, \
|
---|
58 | 1e-13, 5e-12, 1e-11]
|
---|
59 | field_values=[md.results.TransientSolution[0].SedimentHead, \
|
---|
60 | md.results.TransientSolution[0].EplHead, \
|
---|
61 | md.results.TransientSolution[0].SedimentHeadResidual, \
|
---|
62 | md.results.TransientSolution[3].SedimentHead, \
|
---|
63 | md.results.TransientSolution[3].EplHead, \
|
---|
64 | md.results.TransientSolution[3].SedimentHeadResidual, \
|
---|
65 | md.results.TransientSolution[4].SedimentHead, \
|
---|
66 | md.results.TransientSolution[4].EplHead, \
|
---|
67 | md.results.TransientSolution[4].SedimentHeadResidual, \
|
---|
68 | md.results.TransientSolution[8].SedimentHead, \
|
---|
69 | md.results.TransientSolution[8].EplHead, \
|
---|
70 | md.results.TransientSolution[8].SedimentHeadResidual]
|
---|