Changeset 21877 for issm/trunk-jpl/test/NightlyRun/test514.py
- Timestamp:
- 07/26/17 15:39:33 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
TabularUnified issm/trunk-jpl/test/NightlyRun/test514.py ¶
r21411 r21877 14 14 hVertices=10000.*np.ones((27)) 15 15 hVertices[0:5]=1000. 16 md=bamg(model(),'domain','../Exp/Pig.exp','hmax',20000.,'hVertices',hVertices,'gradation',3. ,'geometricalmetric',1)16 md=bamg(model(),'domain','../Exp/Pig.exp','hmax',20000.,'hVertices',hVertices,'gradation',3.) 17 17 x1=md.mesh.x 18 18 y1=md.mesh.y … … 29 29 metric=ComputeMetric(hessian,2./9.,1.,1000.,25.*10.**3,[]) 30 30 md.miscellaneous.dummy=metric 31 md2=bamg(copy.deepcopy(md),'metric',md.miscellaneous.dummy,'hmin',1000.,'hmax',20000.,'gradation',3. ,'geometricalmetric',1)31 md2=bamg(copy.deepcopy(md),'metric',md.miscellaneous.dummy,'hmin',1000.,'hmax',20000.,'gradation',3.) 32 32 x3=md2.mesh.x 33 33 y3=md2.mesh.y 34 34 35 35 #refine existing mesh 2 36 md2=bamg(copy.deepcopy(md),'metric',md.miscellaneous.dummy,'hmin',1000.,'hmax',20000.,'gradation',3.,' geometricalmetric',1,'anisomax',1.)36 md2=bamg(copy.deepcopy(md),'metric',md.miscellaneous.dummy,'hmin',1000.,'hmax',20000.,'gradation',3.,'anisomax',1.) 37 37 x4=md2.mesh.x 38 38 y4=md2.mesh.y … … 41 41 hVertices=np.nan*np.ones((md.mesh.numberofvertices)) 42 42 hVertices[np.nonzero(md.mesh.vertexonboundary)]=500. 43 md2=bamg(copy.deepcopy(md),'metric',md.miscellaneous.dummy,'hmin',1000.,'hmax',20000.,'gradation',3.,' geometricalmetric',1,'anisomax',1.,'hVertices',hVertices)43 md2=bamg(copy.deepcopy(md),'metric',md.miscellaneous.dummy,'hmin',1000.,'hmax',20000.,'gradation',3.,'anisomax',1.,'hVertices',hVertices) 44 44 x5=md2.mesh.x 45 45 y5=md2.mesh.y 46 46 47 47 #refine existing mesh 4 48 md2=bamg(copy.deepcopy(md),'field',md.inversion.vy_obs,'hmin',1000.,'hmax',20000.,'gradation',3.,' geometricalmetric',1,'Hessiantype',0,'err',np.array([[1.]]))48 md2=bamg(copy.deepcopy(md),'field',md.inversion.vy_obs,'hmin',1000.,'hmax',20000.,'gradation',3.,'Hessiantype',0,'err',np.array([[1.]])) 49 49 x6=md2.mesh.x 50 50 y6=md2.mesh.y 51 51 52 52 #refine existing mesh 5 53 md2=bamg(copy.deepcopy(md),'field',np.hstack((md.inversion.vy_obs,md.geometry.thickness)),'hmin',1000.,'hmax',20000.,'gradation',3.,' geometricalmetric',1,'Hessiantype',1,'err',np.array([[10.,100.]]))53 md2=bamg(copy.deepcopy(md),'field',np.hstack((md.inversion.vy_obs,md.geometry.thickness)),'hmin',1000.,'hmax',20000.,'gradation',3.,'Hessiantype',1,'err',np.array([[10.,100.]])) 54 54 x7=md2.mesh.x 55 55 y7=md2.mesh.y
Note:
See TracChangeset
for help on using the changeset viewer.