![(please configure the [header_logo] section in trac.ini)](/trac/issm/chrome/common/trac_banner.png)
Changeset 22287
- Timestamp:
- 12/01/17 08:14:15 (7 years ago)
- Location:
- issm/trunk-jpl/test
- Files:
-
- 3 edited
-
Archives/Archive333.arch (modified) ( previous)
-
NightlyRun/test333.m (modified) (2 diffs)
-
NightlyRun/test333.py (modified) (2 diffs)
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/test/NightlyRun/test333.m
r21480 r22287 34 34 35 35 md=solve(md,'Transient'); 36 %re-run with no substeps 37 mdfine=md; 38 mdfine.result=struct(); 39 mdfine.hydrology.steps_per_step=1; 40 mdfine.timestepping.time_step=0.2; 41 mdfine=solve(md,'Transient'); 36 42 37 43 %store=md.constants.g*md.hydrology.sediment_porosity*md.materials.rho_freshwater*((md.hydrology.sediment_compressibility/md.hydrology.sediment_porosity)+md.hydrology.water_compressibility); … … 43 49 %Fields and tolerances to track changes 44 50 field_names ={'SedimentWaterHead1','EplWaterHead1','SedimentHeadResidual1',... 45 'SedimentWaterHead4','EplWaterHead4','SedimentHeadResidual4',... 46 'SedimentWaterHead5','EplWaterHead5','SedimentHeadResidual5',... 47 'SedimentWaterHead9','EplWaterHead9','SedimentHeadResidual9'}; 51 'SedimentWaterHead4','EplWaterHead4','SedimentHeadResidual4',... 52 'SedimentWaterHead5','EplWaterHead5','SedimentHeadResidual5',... 53 'SedimentWaterHead9','EplWaterHead9','SedimentHeadResidual9',... 54 'EplWaterHead9'}; 48 55 field_tolerances={... 49 56 1e-13, 1e-13, 1e-13,... 50 57 1e-13, 1e-13, 1e-13,... 51 58 1e-13, 5e-12, 1e-11,... 52 1e-13, 5e-12, 1e-11}; 53 field_values={md.results.TransientSolution(1).SedimentHead, ... 54 md.results.TransientSolution(1).EplHead,... 55 md.results.TransientSolution(1).SedimentHeadResidual,... 56 md.results.TransientSolution(4).SedimentHead,... 57 md.results.TransientSolution(4).EplHead,... 58 md.results.TransientSolution(4).SedimentHeadResidual, ... 59 md.results.TransientSolution(5).SedimentHead,... 60 md.results.TransientSolution(5).EplHead,... 61 md.results.TransientSolution(5).SedimentHeadResidual, ... 62 md.results.TransientSolution(9).SedimentHead,... 63 md.results.TransientSolution(9).EplHead,... 64 md.results.TransientSolution(9).SedimentHeadResidual}; 59 1e-13, 5e-12, 1e-11,... 60 1e-13}; 61 field_values={mdfine.results.TransientSolution(1).SedimentHead, ... 62 mdfine.results.TransientSolution(1).EplHead,... 63 mdfine.results.TransientSolution(1).SedimentHeadResidual,... 64 mdfine.results.TransientSolution(4).SedimentHead,... 65 mdfine.results.TransientSolution(4).EplHead,... 66 mdfine.results.TransientSolution(4).SedimentHeadResidual, ... 67 mdfine.results.TransientSolution(5).SedimentHead,... 68 mdfine.results.TransientSolution(5).EplHead,... 69 mdfine.results.TransientSolution(5).SedimentHeadResidual, ... 70 mdfine.results.TransientSolution(9).SedimentHead,... 71 mdfine.results.TransientSolution(9).EplHead,... 72 mdfine.results.TransientSolution(9).SedimentHeadResidual,... 73 md.results.TransientSolution(0).EplHead}; -
issm/trunk-jpl/test/NightlyRun/test333.py
r21638 r22287 42 42 md.hydrology.epl_thick_comp=1 43 43 md.hydrology.epl_max_thickness=1 44 md.timestepping.time_step=0.2 44 md.hydrology.steps_per_step=10 45 md.timestepping.time_step=2.0 45 46 md.timestepping.final_time=2.0 46 47 47 48 md=solve(md,'Transient') 49 50 #re-run with no substeps 51 mdfine=md 52 mdfine.result=[] 53 mdfine.hydrology.steps_per_step=1 54 mdfine.timestepping.time_step=0.2 55 mdfine=solve(md,'Transient') 56 48 57 49 58 # sedvol=np.zeros(10) … … 59 68 'SedimentWaterHead4','EplWaterHead4','SedimentHeadResidual4', 60 69 'SedimentWaterHead5','EplWaterHead5','SedimentHeadResidual5', 61 'SedimentWaterHead9','EplWaterHead9','SedimentHeadResidual9'] 70 'SedimentWaterHead9','EplWaterHead9','SedimentHeadResidual9', 71 'EplWaterHead9'] 62 72 field_tolerances=[1e-13, 1e-13, 1e-13, 63 73 1e-13, 1e-13, 1e-13, 64 74 1e-13, 5e-12, 1e-11, 65 1e-13, 5e-12, 1e-11] 66 field_values=[md.results.TransientSolution[0].SedimentHead, 67 md.results.TransientSolution[0].EplHead, 68 md.results.TransientSolution[0].SedimentHeadResidual, 69 md.results.TransientSolution[3].SedimentHead, 70 md.results.TransientSolution[3].EplHead, 71 md.results.TransientSolution[3].SedimentHeadResidual, 72 md.results.TransientSolution[4].SedimentHead, 73 md.results.TransientSolution[4].EplHead, 74 md.results.TransientSolution[4].SedimentHeadResidual, 75 md.results.TransientSolution[8].SedimentHead, 76 md.results.TransientSolution[8].EplHead, 77 md.results.TransientSolution[8].SedimentHeadResidual] 75 1e-13, 5e-12, 1e-11, 76 1e-13] 77 field_values=[mdfine.results.TransientSolution[0].SedimentHead, 78 mdfine.results.TransientSolution[0].EplHead, 79 mdfine.results.TransientSolution[0].SedimentHeadResidual, 80 mdfine.results.TransientSolution[3].SedimentHead, 81 mdfine.results.TransientSolution[3].EplHead, 82 mdfine.results.TransientSolution[3].SedimentHeadResidual, 83 mdfine.results.TransientSolution[4].SedimentHead, 84 mdfine.results.TransientSolution[4].EplHead, 85 mdfine.results.TransientSolution[4].SedimentHeadResidual, 86 mdfine.results.TransientSolution[8].SedimentHead, 87 mdfine.results.TransientSolution[8].EplHead, 88 mdfine.results.TransientSolution[8].SedimentHeadResidual, 89 md.results.TransientSolution[0].EplHead] 90
Note:
See TracChangeset
for help on using the changeset viewer.
