Changeset 21477 for issm/trunk-jpl/test/NightlyRun/test333.py
- Timestamp:
- 01/11/17 00:26:07 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/test/NightlyRun/test333.py
r21411 r21477 48 48 md=solve(md,'Transient') 49 49 50 sedvol=np.zeros(10)51 eplvol=np.zeros(10)52 time=np.arange(0.2,2.1,0.2)53 store=md.constants.g*md.hydrology.sediment_porosity*md.materials.rho_freshwater*((md.hydrology.sediment_compressibility/md.hydrology.sediment_porosity)+md.hydrology.water_compressibility)54 sedstore=20.0*store55 for i in range(0,10):56 sedvol[i]=np.mean(md.results.TransientSolution[i].SedimentHead)*sedstore57 eplvol[i]=np.mean(md.results.TransientSolution[i].EplHead)*store*np.mean(md.results.TransientSolution[i].HydrologydcEplThickness)50 # sedvol=np.zeros(10) 51 # eplvol=np.zeros(10) 52 # time=np.arange(0.2,2.1,0.2) 53 # store=md.constants.g*md.hydrology.sediment_porosity*md.materials.rho_freshwater*((md.hydrology.sediment_compressibility/md.hydrology.sediment_porosity)+md.hydrology.water_compressibility) 54 # sedstore=20.0*store 55 # for i in range(0,10): 56 # sedvol[i]=np.mean(md.results.TransientSolution[i].SedimentHead)*sedstore 57 # eplvol[i]=np.mean(md.results.TransientSolution[i].EplHead)*store*np.mean(md.results.TransientSolution[i].HydrologydcEplThickness) 58 58 59 59 field_names=['SedimentWaterHead1','EplWaterHead1','SedimentHeadResidual1',
Note:
See TracChangeset
for help on using the changeset viewer.