Ignore:
Timestamp:
07/26/17 15:39:33 (8 years ago)
Author:
Mathieu Morlighem
Message:

CHG: no more geometrical metric

File:
1 edited

Legend:

Unmodified
Added
Removed
  • TabularUnified issm/trunk-jpl/test/NightlyRun/test514.py

    r21411 r21877  
    1414hVertices=10000.*np.ones((27))
    1515hVertices[0:5]=1000.
    16 md=bamg(model(),'domain','../Exp/Pig.exp','hmax',20000.,'hVertices',hVertices,'gradation',3.,'geometricalmetric',1)
     16md=bamg(model(),'domain','../Exp/Pig.exp','hmax',20000.,'hVertices',hVertices,'gradation',3.)
    1717x1=md.mesh.x
    1818y1=md.mesh.y
     
    2929metric=ComputeMetric(hessian,2./9.,1.,1000.,25.*10.**3,[])
    3030md.miscellaneous.dummy=metric
    31 md2=bamg(copy.deepcopy(md),'metric',md.miscellaneous.dummy,'hmin',1000.,'hmax',20000.,'gradation',3.,'geometricalmetric',1)
     31md2=bamg(copy.deepcopy(md),'metric',md.miscellaneous.dummy,'hmin',1000.,'hmax',20000.,'gradation',3.)
    3232x3=md2.mesh.x
    3333y3=md2.mesh.y
    3434
    3535#refine existing mesh 2
    36 md2=bamg(copy.deepcopy(md),'metric',md.miscellaneous.dummy,'hmin',1000.,'hmax',20000.,'gradation',3.,'geometricalmetric',1,'anisomax',1.)
     36md2=bamg(copy.deepcopy(md),'metric',md.miscellaneous.dummy,'hmin',1000.,'hmax',20000.,'gradation',3.,'anisomax',1.)
    3737x4=md2.mesh.x
    3838y4=md2.mesh.y
     
    4141hVertices=np.nan*np.ones((md.mesh.numberofvertices))
    4242hVertices[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)
     43md2=bamg(copy.deepcopy(md),'metric',md.miscellaneous.dummy,'hmin',1000.,'hmax',20000.,'gradation',3.,'anisomax',1.,'hVertices',hVertices)
    4444x5=md2.mesh.x
    4545y5=md2.mesh.y
    4646
    4747#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.]]))
     48md2=bamg(copy.deepcopy(md),'field',md.inversion.vy_obs,'hmin',1000.,'hmax',20000.,'gradation',3.,'Hessiantype',0,'err',np.array([[1.]]))
    4949x6=md2.mesh.x
    5050y6=md2.mesh.y
    5151
    5252#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.]]))
     53md2=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.]]))
    5454x7=md2.mesh.x
    5555y7=md2.mesh.y
Note: See TracChangeset for help on using the changeset viewer.