search

lanl / LaGriT

Los Alamos Grid Toolbox (LaGriT) is a library of user callable tools that provide mesh generation, mesh optimization and dynamic mesh maintenance in two and three dimensions.
https://lanl.github.io/LaGriT/
Other
116 stars 48 forks source link

mesh.fehm nodes running together dump_fehm('mesh') #191

Closed shane-mckinney closed 4 years ago

shane-mckinney commented 4 years ago

This is the result of grepping for the node number that FEHM is trying to read when producing its first RTE.

$ grep "1856605" mesh.fehmn
   1856605 -3.891200000000E+03 -6.826666667000E+01  1.000000000000E+02
  11856605  2.184533333330E+03  3.754666666670E+03  6.500000000000E+02
   8053911766005117660041185660511856305
  888475 1947205 1856605 1856604 1856303
 1856605 4335801 4245501 4336101 4245500
 434504211856907118566051185660611766005
 434654011856906118566051194750711947207
 437717111946905118566051194720511856604
 441002611856303118566051185630411946905

As you can see, it looks like when all integers are 8 digits or bigger, they string together without a space in between!

This is the pylagrit script I am running to generate this mesh:

from pylagrit import PyLaGriT
import numpy

lg = PyLaGriT(lagrit_exe='/project/eesdev/tam/LaGriT/src/mylagrit')

x = y = numpy.linspace(-5120,5120,301)
z = numpy.linspace(0,1000,201)
mtet = lg.gridder(x,y,z,elem_type='tet',connect=True)

mtet.setatt('imt',1)
pcap = mtet.pset_geom([-5000,5000,50],[-5000,5000,1000])
pcap.setatt('imt',2)

mtet.resetpts_itp()

mtet.dump_fehm('mesh')
mtet.dump('mesh.inp')

I believe this is happening in mtet.dump_fehm('mesh') as the problem is seen in mesh.fehm

shane-mckinney commented 4 years ago

You can check to see if the file has these ran-together instance in them by doing: awk -F: '{if($1>10000000000)print$1}' <mesh.stor where 'mesh.stor' is any file you want to check for a number greater than 10000000000