diff --git a/processing/post/spectral_post.py b/processing/post/spectral_post.py new file mode 100644 index 000000000..d1db46192 --- /dev/null +++ b/processing/post/spectral_post.py @@ -0,0 +1,18 @@ +print('post processing for mpie_spectral') + +results = open('results.out', 'rb') + +header = results.read(4) +end = results.read(1024).find(header) +results.seek(4) +loadcase = results.read(end) +print('loadcase:', loadcase) +# +begin = end +results.seek(end+4) +header = results.read(4) +print(header) +end = results.read(1024).find(header) +#results.seek(4) +#workingdir = results.read(end) +#print('workingdir:', workingdir) diff --git a/processing/pre/voronoi fast.f90 b/processing/pre/voronoi fast.f90 index 699dccbe7..0c676bd2e 100644 --- a/processing/pre/voronoi fast.f90 +++ b/processing/pre/voronoi fast.f90 @@ -171,8 +171,27 @@ program voronoi !write visualization files (in case wanted) if (choice == 'y' .or. choice == 'Y') then print*, 'for more information on gmsh: http://geuz.org/gmsh/' - + ! write full mesh out + open(20, file = ((trim(name))//'_3Dfull.msh')) + write(20, '(A, /, A, /, A, /, A, /, I10)'), '$MeshFormat', '2.1 0 8', '$EndMeshFormat', '$Nodes', abc + do i = 1, abc + write(20, '(I10, I10, I10, I10)'), i, mod((i-1), a) +1, mod(((i-1)/a), b) +1, mod(((i-1)/(ab)), c) +1 + end do + write(20, '(A, /, A, /, I10)'), '$EndNodes', '$Elements', abc + do i = 1, abc + write(20, '(I10, A, I10, A, I10)'), i, ' 15 2', grainMap(i), ' 2', i + end do + write(20, '(A)'), '$EndElements' + write(20, '(A, /, A, /, A, /, A, /, A, /, A, /, A, /, A, /, I10)'), '$NodeData', '1', '"Grain No."', '1', & + &'0.0', '3', '0', '1', abc + do i = 1, abc + write(20, '(I10, tr2, I10)'), i, grainMap(i) + end do + write(20, *), '$EndNodeData' + close(20) + +! write 3d skin out open(20, file = ((trim(name))//'_3D.msh')) write(20, '(A, /, A, /, A, /, A, /, I10)'), '$MeshFormat', '2.1 0 8', '$EndMeshFormat', '$Nodes', abc_Red do j = 1, abc_Red