source: issm/trunk-jpl/test/NightlyRun/test333.py@ 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.8 KB
Line 
1#Test Name: SquareSheetHydrologyDCTwoLayers
2import numpy
3from model import *
4from triangle import *
5from setmask import *
6from parameterize import *
7from transient import *
8from setflowequation import *
9from EnumDefinitions import *
10from solve import *
11from MatlabFuncs import *
12
13from generic import generic
14
15md=triangle(model(),'../Exp/Square.exp',100000.)
16md=setmask(md,'','')
17md=parameterize(md,'../Par/SquareSheetConstrained.py')
18md.transient=transient.setallnullparameters(md.transient)
19md.transient.ishydrology=True
20md.transient.issmb=True
21md=setflowequation(md,'SSA','all')
22md.cluster=generic('name',oshostname(),'np',1)
23md.hydrology=hydrologydc()
24md.hydrology=md.hydrology.initialize(md)
25
26md.hydrology.isefficientlayer=1
27md.hydrology.sedimentlimit_flag=1
28md.hydrology.sedimentlimit=800.0
29md.initialization.sediment_head=numpy.zeros((md.mesh.numberofvertices,1))
30md.hydrology.spcsediment_head=float('NaN')*numpy.ones((md.mesh.numberofvertices,1))
31
32md.basalforcings.groundedice_melting_rate = 2.0*numpy.ones((md.mesh.numberofvertices,1))
33md.basalforcings.floatingice_melting_rate = 0.0*numpy.ones((md.mesh.numberofvertices,1))
34md.hydrology.sediment_transmitivity= 3.0*numpy.ones((md.mesh.numberofvertices,1))
35
36md.initialization.epl_head=numpy.zeros((md.mesh.numberofvertices,1))
37md.initialization.epl_thickness=numpy.ones((md.mesh.numberofvertices,1))
38md.hydrology.spcepl_head=float('NaN')*numpy.ones((md.mesh.numberofvertices,1))
39md.hydrology.mask_eplactive_node=numpy.zeros((md.mesh.numberofvertices,1))
40md.hydrology.epl_conductivity=30
41md.hydrology.epl_initial_thickness=1
42md.hydrology.epl_colapse_thickness=1.0e-3
43md.hydrology.epl_thick_comp=1
44md.hydrology.epl_max_thickness=1
45md.timestepping.time_step=0.2
46md.timestepping.final_time=2.0
47
48md=solve(md,'Transient')
49
50field_names =['SedimentWaterHead1','EplWaterHead1','SedimentHeadResidual1', \
51 'SedimentWaterHead4','EplWaterHead4','SedimentHeadResidual4', \
52 'SedimentWaterHead5','EplWaterHead5','SedimentHeadResidual5', \
53 'SedimentWaterHead9','EplWaterHead9','SedimentHeadResidual9']
54field_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]
59field_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]
Note: See TracBrowser for help on using the repository browser.