diff --git a/VERSION b/VERSION
index 428961dd5..fa691fc66 100644
--- a/VERSION
+++ b/VERSION
@@ -1 +1 @@
-v3.0.0-alpha-203-g3e990489
+v3.0.0-alpha-216-gb51a1b7b
diff --git a/src/HDF5_utilities.f90 b/src/HDF5_utilities.f90
index faeff03e3..a13ac1a8a 100644
--- a/src/HDF5_utilities.f90
+++ b/src/HDF5_utilities.f90
@@ -88,7 +88,7 @@ subroutine HDF5_utilities_init
   integer         :: hdferr
   integer(SIZE_T) :: typeSize
 
-  write(6,'(/,a)') ' <<<+-  HDF5_Utilities init  -+>>>'
+  print'(/,a)', ' <<<+-  HDF5_Utilities init  -+>>>'
 
 !--------------------------------------------------------------------------------------------------
 !initialize HDF5 library and check if integer and float type size match
diff --git a/src/crystallite.f90 b/src/crystallite.f90
index bed36d79b..2ee85024b 100644
--- a/src/crystallite.f90
+++ b/src/crystallite.f90
@@ -1561,9 +1561,9 @@ subroutine crystallite_restartWrite
   integer(HID_T) :: fileHandle, groupHandle
   character(len=pStringLen) :: fileName, datasetName
 
-  print'(a)', ' writing field and constitutive data required for restart to file';flush(6)
+  print*, ' writing field and constitutive data required for restart to file';flush(6)
 
-  write(fileName,'(a,i0,a)', trim(getSolverJobName())//'_',worldrank,'.hdf5'
+  write(fileName,'(a,i0,a)') trim(getSolverJobName())//'_',worldrank,'.hdf5'
   fileHandle = HDF5_openFile(fileName,'a')
 
   call HDF5_write(fileHandle,crystallite_partionedF,'F')
@@ -1575,14 +1575,14 @@ subroutine crystallite_restartWrite
 
   groupHandle = HDF5_addGroup(fileHandle,'constituent')
   do i = 1,size(material_name_phase)
-    write(datasetName,'(i0,a)', i,'_omega_plastic'
+    write(datasetName,'(i0,a)') i,'_omega_plastic'
     call HDF5_write(groupHandle,plasticState(i)%state,datasetName)
   enddo
   call HDF5_closeGroup(groupHandle)
 
   groupHandle = HDF5_addGroup(fileHandle,'materialpoint')
   do i = 1, material_Nhomogenization
-    write(datasetName,'(i0,a)', i,'_omega_homogenization'
+    write(datasetName,'(i0,a)') i,'_omega_homogenization'
     call HDF5_write(groupHandle,homogState(i)%state,datasetName)
   enddo
   call HDF5_closeGroup(groupHandle)
@@ -1604,7 +1604,7 @@ subroutine crystallite_restartRead
 
   print'(/,a,i0,a)', ' reading restart information of increment from file'
 
-  write(fileName,'(a,i0,a)', trim(getSolverJobName())//'_',worldrank,'.hdf5'
+  write(fileName,'(a,i0,a)') trim(getSolverJobName())//'_',worldrank,'.hdf5'
   fileHandle = HDF5_openFile(fileName)
 
   call HDF5_read(fileHandle,crystallite_F0, 'F')
diff --git a/src/discretization.f90 b/src/discretization.f90
index 4520bf2ca..3bb905cf3 100644
--- a/src/discretization.f90
+++ b/src/discretization.f90
@@ -49,7 +49,7 @@ subroutine discretization_init(microstructureAt,&
   integer, optional,           intent(in) :: &
     sharedNodesBegin                                                                                !< index of first node shared among different processes (MPI)
 
-  write(6,'(/,a)') ' <<<+-  discretization init  -+>>>'; flush(6)
+  print'(/,a)', ' <<<+-  discretization init  -+>>>'; flush(6)
 
   discretization_nElem = size(microstructureAt,1)
   discretization_nIP   = size(IPcoords0,2)/discretization_nElem
diff --git a/src/element.f90 b/src/element.f90
index 836c8b2b1..3ccc2fb78 100644
--- a/src/element.f90
+++ b/src/element.f90
@@ -922,7 +922,7 @@ subroutine tElement_init(self,elemType)
 
   self%nIPneighbors = size(self%IPneighbor,1)
 
-  write(6,'(/,a)') ' <<<+-  element_init  -+>>>'; flush(6)
+  print'(/,a)', ' <<<+-  element_init  -+>>>'; flush(6)
 
   print*, 'element type:      ',self%elemType
   print*, '  geom type:       ',self%geomType