Changeset 23793 for issm/trunk-jpl/test/NightlyRun/test2001.py
- Timestamp:
- 03/13/19 03:17:46 (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/test/NightlyRun/test2001.py
r21585 r23793 1 1 #Test Name: SquareSheetConstrainedGia2d 2 2 import numpy as np 3 import copy4 3 from model import * 5 4 from socket import gethostname … … 11 10 12 11 13 #Define a model 14 md =model()15 md =triangle(md,'../Exp/Square.exp',100000.)16 md =setmask(md,'','')17 md =parameterize(md,'../Par/SquareSheetConstrained.py')12 #Define a model 13 md = model() 14 md = triangle(md, '../Exp/Square.exp', 100000.) 15 md = setmask(md, '', '') 16 md = parameterize(md, '../Par/SquareSheetConstrained.py') 18 17 19 #Indicate what you want to compute 20 md.gia.cross_section_shape =1 # for square-edged x-section18 #Indicate what you want to compute 19 md.gia.cross_section_shape = 1 # for square-edged x-section 21 20 22 21 #Define loading history (see test2001.m for the description) 23 md.timestepping.start_time =2400000 # 2,400 kyr24 md.timestepping.final_time =2500000 # 2,500 kyr25 md.geometry.thickness =np.vstack((np.hstack((md.geometry.thickness*0.0, 0.0)),26 np.hstack((md.geometry.thickness/2.0, 0.1)),27 28 29 22 md.timestepping.start_time = 2400000 # 2, 400 kyr 23 md.timestepping.final_time = 2500000 # 2, 500 kyr 24 md.geometry.thickness = np.vstack((np.hstack((md.geometry.thickness * 0.0, 0.0)), 25 np.hstack((md.geometry.thickness / 2.0, 0.1)), 26 np.hstack((md.geometry.thickness, 0.2)), 27 np.hstack((md.geometry.thickness, 1.0)), 28 np.hstack((md.geometry.thickness, md.timestepping.start_time)))).T 30 29 31 #Solve for GIA deflection 32 md.cluster =generic('name',gethostname(),'np',3)33 md =solve(md,'Gia')30 #Solve for GIA deflection 31 md.cluster = generic('name', gethostname(), 'np', 3) 32 md = solve(md, 'Gia') 34 33 35 34 #Fields and tolerances to track changes 36 field_names =['GiaW','GiadWdt'] 37 field_tolerances=[1e-13,1e-13] 38 field_values =[md.results.GiaSolution.GiaW, 39 md.results.GiaSolution.GiadWdt] 40 35 field_names = ['GiaW', 'GiadWdt'] 36 field_tolerances = [1e-13, 1e-13] 37 field_values = [md.results.GiaSolution.GiaW, 38 md.results.GiaSolution.GiadWdt]
Note:
See TracChangeset
for help on using the changeset viewer.