Changeset 16137 for issm/trunk/test/NightlyRun/test1206.py
- Timestamp:
- 09/16/13 09:43:55 (12 years ago)
- Location:
- issm/trunk
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk
- Property svn:ignore
-
old new 1 nightlylog 2 configure.sh 1 3 par 2 4 ad
-
- Property svn:mergeinfo changed
/issm/trunk-jpl merged: 15397-15401,15403-15487,15489-15701,15704-15735,15737-16076,16082-16133
- Property svn:ignore
-
issm/trunk/test
- Property svn:mergeinfo changed
-
issm/trunk/test/NightlyRun/test1206.py
r14310 r16137 10 10 11 11 """ 12 The aim of this program is to compare a model with an analytical solution given in MacAyealEISMINT : Lessons in Ice-Sheet Modeling.12 The aim of this program is to compare a model with an analytical solution given in SSA EISMINT : Lessons in Ice-Sheet Modeling. 13 13 """ 14 14 … … 32 32 #We extrude the model to have a 3d model 33 33 md.extrude(numlayers,1.) 34 md=setflowequation(md,' pattyn','all')34 md=setflowequation(md,'HO','all') 35 35 36 36 #Spc the nodes on the bed 37 37 pos=numpy.nonzero(md.mesh.vertexonbed) 38 md. diagnostic.spcvx[pos]=0.39 md. diagnostic.spcvy[pos]=0.40 md. diagnostic.spcvz[pos]=0.38 md.stressbalance.spcvx[pos]=0. 39 md.stressbalance.spcvy[pos]=0. 40 md.stressbalance.spcvz[pos]=0. 41 41 42 42 #Now we can solve the problem 43 43 md.cluster=generic('name',oshostname(),'np',8) 44 md=solve(md, DiagnosticSolutionEnum())44 md=solve(md,StressbalanceSolutionEnum()) 45 45 46 46 #Calculate the depth averaged velocity field (2d): 47 vx=md.results. DiagnosticSolution.Vx48 vy=md.results. DiagnosticSolution.Vy47 vx=md.results.StressbalanceSolution.Vx 48 vy=md.results.StressbalanceSolution.Vy 49 49 vel=numpy.zeros((md.mesh.numberofvertices2d,1)) 50 50 … … 94 94 pass 95 95 # set(gcf,'Color','w') 96 # printmodel(' pattynstatic','png','margin','on','marginsize',25,'frame','off','resolution',0.7,'hardcopy','off');97 # system(['mv pattynstatic.png ' ISSM_DIR '/website/doc_pdf/validation/Images/EISMINT/IceSheet']);96 # printmodel('HOstatic','png','margin','on','marginsize',25,'frame','off','resolution',0.7,'hardcopy','off'); 97 # system(['mv HOstatic.png ' ISSM_DIR '/website/doc_pdf/validation/Images/EISMINT/IceSheet']); 98 98 99 99 #Fields and tolerances to track changes
Note:
See TracChangeset
for help on using the changeset viewer.