need to store cell nodes, not elem nodes

src/discretization.f90

@@ -51,7 +51,7 @@ subroutine discretization_init(homogenizationAt,microstructureAt,&
   integer, optional,           intent(in) :: &
     sharedNodesBeginn
 
-  write(6,'(/,a)')   ' <<<+-  discretization init  -+>>>'
+  write(6,'(/,a)') ' <<<+-  discretization init  -+>>>'; flush(6)
 
   discretization_nElem = size(microstructureAt,1)
   discretization_nIP   = size(IPcoords0,2)/discretization_nElem
@@ -81,15 +81,15 @@ subroutine discretization_results
 
   real(pReal), dimension(:,:), allocatable :: u
   
-  call results_closeGroup(results_addGroup(trim('current/geometry')))
+  call results_closeGroup(results_addGroup('current/geometry'))
   
   u = discretization_NodeCoords (1:3,:discretization_sharedNodesBeginn) &
     - discretization_NodeCoords0(1:3,:discretization_sharedNodesBeginn)
-  call results_writeDataset('current/geometry',u,'u_n','nodal displacements','m')
+  call results_writeDataset('current/geometry',u,'u_n','displacements of the nodes','m')
   
   u = discretization_IPcoords &
     - discretization_IPcoords0
-  call results_writeDataset('current/geometry',u,'u_c','cell center displacements','m')
+  call results_writeDataset('current/geometry',u,'u_p','displacements of the materialpoints','m')
 
 end subroutine discretization_results
 

src/mesh_marc.f90

@@ -102,7 +102,7 @@ subroutine mesh_init(ip,el)
 
   call discretization_init(microstructureAt,homogenizationAt,&
                            ip_reshaped,&
-                           node0_elem)
+                           node0_cell)
                            
   call writeGeometry(0,connectivity_elem,&
                      reshape(connectivity_cell,[elem%NcellNodesPerCell,elem%nIPs*nElems]),&
@@ -152,11 +152,11 @@ subroutine writeGeometry(elemType, &
                             
   coordinates_temp = coordinates_nodes
   call results_writeDataset('geometry',coordinates_temp,'x_n', &
-                            'coordinates of the nodes','m')
+                            'initial coordinates of the nodes','m')
                             
   coordinates_temp = coordinates_points
   call results_writeDataset('geometry',coordinates_temp,'x_p', &
-                            'coordinates of the material points','m')
+                            'initial coordinates of the materialpoints','m')
                       
   call results_closeJobFile