/[escript]/trunk/downunder/py_src/domainbuilder.py
ViewVC logotype

Diff of /trunk/downunder/py_src/domainbuilder.py

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 4433 by gross, Fri May 31 12:09:58 2013 UTC revision 4434 by caltinay, Mon Jun 3 00:43:29 2013 UTC
# Line 442  class DomainBuilder(object): Line 442  class DomainBuilder(object):
442          else:          else:
443              # this should give us about meter-accuracy with lat/lon coords              # this should give us about meter-accuracy with lat/lon coords
444              self._dom_origin = [1e-5*np.floor(oi*1e5) for oi in origin]              self._dom_origin = [1e-5*np.floor(oi*1e5) for oi in origin]
445              f=[1.,1.,  1./self.getReferenceSystem().getHeightUnit() ]              if self.__dim==3:
446                    f=[1.,1.,  1./self.getReferenceSystem().getHeightUnit()]
447                else:
448                    f=[1., 1./self.getReferenceSystem().getHeightUnit()]
449          # cell size / point spacing          # cell size / point spacing
450          spacing = DX + [np.floor((self._v_depth+self._v_air_layer)/self._v_num_cells)]          spacing = DX + [np.floor((self._v_depth+self._v_air_layer)/self._v_num_cells)]
451          #self._spacing = [float(np.floor(si)) for si in spacing]          #self._spacing = [float(np.floor(si)) for si in spacing]
# Line 454  class DomainBuilder(object): Line 457  class DomainBuilder(object):
457          if self.__dim==3:          if self.__dim==3:
458              dom=Brick(*NE, l0=lo[0], l1=lo[1], l2=lo[2])              dom=Brick(*NE, l0=lo[0], l1=lo[1], l2=lo[2])
459          else:          else:
460              dom=Rectangle(*NE, l0=lo[0], l1=lo[2])              dom=Rectangle(*NE, l0=lo[0], l1=lo[1])
461    
462          # ripley may internally adjust NE and length, so recompute          # ripley may internally adjust NE and length, so recompute
463          self._dom_len=[sup(dom.getX()[i])-inf(dom.getX()[i]) for i in range(self.__dim)]          self._dom_len=[sup(dom.getX()[i])-inf(dom.getX()[i]) for i in range(self.__dim)]

Legend:
Removed from v.4433  
changed lines
  Added in v.4434

  ViewVC Help
Powered by ViewVC 1.1.26