![(please configure the [header_logo] section in trac.ini)](/trac/issm/chrome/common/trac_banner.png)
Changeset 16171
- Timestamp:
- 09/18/13 13:49:15 (12 years ago)
- Location:
- issm/trunk-jpl/test/Par
- Files:
-
- 5 edited
-
79North.py (modified) (3 diffs)
-
Pig.py (modified) (3 diffs)
-
SquareSheetShelf.py (modified) (3 diffs)
-
SquareShelf.py (modified) (4 diffs)
-
SquareShelf2.py (modified) (4 diffs)
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/test/Par/79North.py
r16170 r16171 12 12 #Geometry and observation 13 13 f = netCDF4.Dataset('../Data/79North.nc','r') 14 x = reshape(f.variables['x'][:],(-1))15 y = reshape(f.variables['y'][:],(-1))14 x = numpy.reshape(f.variables['x'][:],(-1)) 15 y = numpy.reshape(f.variables['y'][:],(-1)) 16 16 vx = f.variables['vx'][:] 17 17 vy = f.variables['vy'][:] … … 28 28 29 29 #Materials 30 md.initialization.temperature=(273.-20.)* ones((md.mesh.numberofvertices,1))30 md.initialization.temperature=(273.-20.)*numpy.ones((md.mesh.numberofvertices,1)) 31 31 md.materials.rheology_B=paterson(md.initialization.temperature) 32 md.materials.rheology_n=3.* ones((md.mesh.numberofelements,1))32 md.materials.rheology_n=3.*numpy.ones((md.mesh.numberofelements,1)) 33 33 md.initialization.temperature=md.initialization.temperature 34 34 … … 37 37 38 38 #Friction 39 md.friction.coefficient=50.* ones((md.mesh.numberofvertices,1))39 md.friction.coefficient=50.*numpy.ones((md.mesh.numberofvertices,1)) 40 40 md.friction.coefficient[numpy.nonzero(md.mask.groundedice_levelset<0.)[0]]=0. 41 md.friction.p= ones((md.mesh.numberofelements,1))42 md.friction.q= ones((md.mesh.numberofelements,1))41 md.friction.p=numpy.ones((md.mesh.numberofelements,1)) 42 md.friction.q=numpy.ones((md.mesh.numberofelements,1)) 43 43 44 44 #Ice shelf melting and surface mass balance -
issm/trunk-jpl/test/Par/Pig.py
r16170 r16171 12 12 #Geometry and observation 13 13 f = netCDF4.Dataset('../Data/Pig.nc','r') 14 x = reshape(f.variables['x'][:],(-1))15 y = reshape(f.variables['y'][:],(-1))14 x = numpy.reshape(f.variables['x'][:],(-1)) 15 y = numpy.reshape(f.variables['y'][:],(-1)) 16 16 vx_obs = f.variables['vx_obs'][:] 17 17 vy_obs = f.variables['vy_obs'][:] … … 32 32 33 33 #Materials 34 md.initialization.temperature=(273.-20.)* ones((md.mesh.numberofvertices,1))34 md.initialization.temperature=(273.-20.)*numpy.ones((md.mesh.numberofvertices,1)) 35 35 md.materials.rheology_B=paterson(md.initialization.temperature) 36 md.materials.rheology_n=3.* ones((md.mesh.numberofelements,1))36 md.materials.rheology_n=3.*numpy.ones((md.mesh.numberofelements,1)) 37 37 md.initialization.temperature=md.initialization.temperature 38 38 … … 41 41 42 42 #Friction 43 md.friction.coefficient=50.* ones((md.mesh.numberofvertices,1))43 md.friction.coefficient=50.*numpy.ones((md.mesh.numberofvertices,1)) 44 44 md.friction.coefficient[numpy.nonzero(md.mask.groundedice_levelset<0.)[0]]=0. 45 md.friction.p= ones((md.mesh.numberofelements,1))46 md.friction.q= ones((md.mesh.numberofelements,1))45 md.friction.p=numpy.ones((md.mesh.numberofelements,1)) 46 md.friction.q=numpy.ones((md.mesh.numberofelements,1)) 47 47 48 48 #Numerical parameters -
issm/trunk-jpl/test/Par/SquareSheetShelf.py
r16170 r16171 26 26 #Initial velocity 27 27 f = netCDF4.Dataset('../Data/SquareSheetShelf.nc','r') 28 x = reshape(f.variables['x'][:],(-1))29 y = reshape(f.variables['y'][:],(-1))28 x = numpy.reshape(f.variables['x'][:],(-1)) 29 y = numpy.reshape(f.variables['y'][:],(-1)) 30 30 vx = f.variables['vx'][:] 31 31 vy = f.variables['vy'][:] … … 39 39 40 40 #Materials 41 md.initialization.temperature=(273.-20.)* ones((md.mesh.numberofvertices,1))41 md.initialization.temperature=(273.-20.)*numpy.ones((md.mesh.numberofvertices,1)) 42 42 md.materials.rheology_B=paterson(md.initialization.temperature) 43 md.materials.rheology_n=3.* ones((md.mesh.numberofelements,1))43 md.materials.rheology_n=3.*numpy.ones((md.mesh.numberofelements,1)) 44 44 45 45 #Damage … … 51 51 52 52 #Friction 53 md.friction.coefficient=20.* ones((md.mesh.numberofvertices,1))53 md.friction.coefficient=20.*numpy.ones((md.mesh.numberofvertices,1)) 54 54 md.friction.coefficient[numpy.nonzero(md.mask.groundedice_levelset<0.)[0]]=0. 55 md.friction.p= ones((md.mesh.numberofelements,1))56 md.friction.q= ones((md.mesh.numberofelements,1))55 md.friction.p=numpy.ones((md.mesh.numberofelements,1)) 56 md.friction.q=numpy.ones((md.mesh.numberofelements,1)) 57 57 58 58 #Numerical parameters -
issm/trunk-jpl/test/Par/SquareShelf.py
r16170 r16171 23 23 iVelF = netCDF4.Dataset('../Data/SquareShelf.nc','r') 24 24 25 x= reshape(iVelF.variables['x'][:],(-1))26 y= reshape(iVelF.variables['y'][:],(-1))25 x=numpy.reshape(iVelF.variables['x'][:],(-1)) 26 y=numpy.reshape(iVelF.variables['y'][:],(-1)) 27 27 vx=iVelF.variables['vx'][:] 28 28 vy=iVelF.variables['vy'][:] … … 55 55 56 56 #Materials 57 md.initialization.temperature = (273.-20.)* ones((md.mesh.numberofvertices,1))57 md.initialization.temperature = (273.-20.)*numpy.ones((md.mesh.numberofvertices,1)) 58 58 md.materials.rheology_B = paterson(md.initialization.temperature) 59 md.materials.rheology_n = 3.* ones((md.mesh.numberofelements,1))59 md.materials.rheology_n = 3.*numpy.ones((md.mesh.numberofelements,1)) 60 60 61 61 #Damage … … 63 63 64 64 #Friction 65 md.friction.coefficient = 20.* ones((md.mesh.numberofvertices,1))65 md.friction.coefficient = 20.*numpy.ones((md.mesh.numberofvertices,1)) 66 66 md.friction.coefficient[numpy.nonzero(md.mask.groundedice_levelset<0.)[0]]=0. 67 md.friction.p = ones((md.mesh.numberofelements,1))68 md.friction.q = ones((md.mesh.numberofelements,1))67 md.friction.p = numpy.ones((md.mesh.numberofelements,1)) 68 md.friction.q = numpy.ones((md.mesh.numberofelements,1)) 69 69 70 70 #Numerical parameters … … 77 77 md.steadystate.reltol = 0.02 78 78 md.stressbalance.reltol = 0.02 79 md.stressbalance.abstol = nan79 md.stressbalance.abstol = float('nan') 80 80 md.timestepping.time_step = 1. 81 81 md.timestepping.final_time = 3. -
issm/trunk-jpl/test/Par/SquareShelf2.py
r16170 r16171 23 23 iVelF = netCDF4.Dataset('../Data/SquareShelf.nc','r') 24 24 25 x= reshape(iVelF.variables['x'][:],(-1))26 y= reshape(iVelF.variables['y'][:],(-1))25 x=numpy.reshape(iVelF.variables['x'][:],(-1)) 26 y=numpy.reshape(iVelF.variables['y'][:],(-1)) 27 27 vx=iVelF.variables['vx'][:] 28 28 vy=iVelF.variables['vy'][:] … … 55 55 56 56 #Materials 57 md.initialization.temperature = (273.-20.)* ones((md.mesh.numberofvertices,1))57 md.initialization.temperature = (273.-20.)*numpy.ones((md.mesh.numberofvertices,1)) 58 58 md.materials.rheology_B = paterson(md.initialization.temperature) 59 md.materials.rheology_n = 3.* ones((md.mesh.numberofelements,1))59 md.materials.rheology_n = 3.*numpy.ones((md.mesh.numberofelements,1)) 60 60 61 61 #Damage … … 63 63 64 64 #Friction 65 md.friction.coefficient = 20.* ones((md.mesh.numberofvertices,1))65 md.friction.coefficient = 20.*numpy.ones((md.mesh.numberofvertices,1)) 66 66 md.friction.coefficient[numpy.nonzero(md.mask.groundedice_levelset<0.)[0]]=0. 67 md.friction.p = ones((md.mesh.numberofelements,1))68 md.friction.q = ones((md.mesh.numberofelements,1))67 md.friction.p = numpy.ones((md.mesh.numberofelements,1)) 68 md.friction.q = numpy.ones((md.mesh.numberofelements,1)) 69 69 70 70 #Numerical parameters … … 77 77 md.steadystate.reltol = 0.02 78 78 md.stressbalance.reltol = 0.02 79 md.stressbalance.abstol = nan79 md.stressbalance.abstol = float('nan') 80 80 md.timestepping.time_step = 1. 81 81 md.timestepping.final_time = 3.
Note:
See TracChangeset
for help on using the changeset viewer.
