4cubit.cmd('brick x 1 y 1 z 1')
5cubit.cmd('set duplicate block elements on')
6cubit.cmd('block 1 volume 1 ')
7cubit.cmd('block 1 attribute count 3')
9cubit.cmd('block 1 name "MAT_ELASTIC1"')
11cubit.cmd('block 1 attribute index 1 1.')
13cubit.cmd('block 1 attribute index 2 0.25')
15cubit.cmd('block 2 attribute index 3 1.0')
17cubit.cmd('block 2 volume 1 ')
18cubit.cmd('block 2 attribute count 2')
20cubit.cmd('block 2 name "MAT_THERMAL"')
21cubit.cmd('block 2 attribute count 2')
23cubit.cmd('block 2 attribute index 1 1.')
25cubit.cmd('block 2 attribute index 2 1.')
28cubit.cmd('create Displacement on vertex 7 dof 1 dof 2 dof 3 fix 0')
29cubit.cmd('create Displacement on vertex 4 dof 1 dof 2 fix 0')
30cubit.cmd('create Displacement on vertex 8 dof 1 fix 0')
33cubit.cmd('volume 1 size 1')
34cubit.cmd('volume 1 scheme Tetmesh')
35cubit.cmd('mesh volume 1')
37cubit.cmd('save as "/Users/likask/MyBuild/mofem-bitbucket/mofem/meshes/elasticity_atom_test_thermal_expansion.cub" overwrite')