Create a geometry in Python

jesantos · October 2, 2021, 8:52pm

Dear Palabeacens and friends,

Do you guys know how to print a geometry using python?
I have not been successful. Here is the code to reproduce the issue:

github.com

je-santos/MPLBM-UT/blob/406865fc5513e4ac7dbf1f3bb567063b2d3fed2a/examples/pythonIO/create_geom.py#L49


    if args.num_slices:
        erock = np.pad(erock,[(args.num_slices,args.num_slices),(0,0),(0,0)])
        
    
    if args.print_size:
        size = erock.shape
        geom_name = f'{args.name}_{size[0]}_{size[1]}_{size[2]}'
    else:
        geom_name = args.name
        
    #np.savetxt(f'input/{geom_name}.dat', erock.flatten())
    #with open(f'input/{geom_name}.dat', 'w') as file:
    #    for i in range(erock.shape[0]):
    #        file.write(erock[i,:,:])
    #erock.astype('int8').tofile(f'input/{geom_name}.dat')
    np.savetxt(f'input/{geom_name}.dat', erock.reshape([252,-1]), fmt='%i', newline='/n', delimiter=' ')
    return erock






# read-in file

Many thanks