Commit 0ad7bcce authored by Giovanni La Mura's avatar Giovanni La Mura
Browse files

Enable inertia.py script to parse FORTRAN and C++ floats

parent 17aa00c1

src/scripts/inertia.py

+6 −6
Original line number Diff line number Diff line
@@ -25,7 +25,7 @@ 


import math

import numpy as np

#import pdb

import pdb



from sys import argv
@@ -137,11 +137,11 @@ def get_centers(config): 
    for i in range(nsph):

        fline = geom_file.readline()

        elems = ingest_line(fline)

        value = float(elems[0])

        value = float(elems[0].replace('D', 'E').replace('d', 'E'))

        centers['x'].append(value)

        value = float(elems[1])

        value = float(elems[1].replace('D', 'E').replace('d', 'E'))

        centers['y'].append(value)

        value = float(elems[2])

        value = float(elems[2].replace('D', 'E').replace('d', 'E'))

        centers['z'].append(value)

    geom_file.close()

    return centers
@@ -403,8 +403,8 @@ def get_types(config): 
        while (read_layers < num_layers):

            fline = sc_file.readline()

            split_line = fline.split('(')[1].split(',')

            rval = float(split_line[0])

            ival = float(split_line[1][:-2])

            rval = float(split_line[0].replace('D', 'E').replace('d', 'E'))

            ival = float(split_line[1][:-2].replace('D', 'E').replace('d', 'E'))

            dc0 = (rval, ival)

            if not dc0 in found_dc0s:

                found_dc0s.append(dc0)