anchor.diff:        if(E_IO == 0) then      
           read(aux,*, iostat=E_IO) nx1,nx2,ny1,ny2,nz1,nz2
           if(E_IO == 0) then      
  -          NN = (nx2-nx1+1)*(ny2-ny1+1)*(nz2-nz1+1)
             E_IO = search(from=pos, inside='Coordinates', to_find='DataArray', with_attribute='Name', of_value='X', &
                           buffer=s_buffer,content=data)
             call get_char(buffer=s_buffer, attrib='format', val=fmt,  E_IO=E_IO)
  
positive.code:               trim(adjustlt(Upper_Case(typeZ)))/='FLOAT64') then
              E_IO = -1_I4P 
            else
              allocate(X(NN), Y(NN), Z(NN), stat=E_IO)        
              read(unit=vtk(rf)%u, iostat=E_IO, pos = vtk(rf)%ioffset+offsX) N_Byte, (X(i), i=1,NN) !get appended array
              print*, N_Byte
              read(unit=vtk(rf)%u, iostat=E_IO, pos = vtk(rf)%ioffset+offsY) N_Byte, (Y(i), i=1,NN) !get appended array
              print*, N_Byte
              read(unit=vtk(rf)%u, iostat=E_IO, pos = vtk(rf)%ioffset+offsZ) N_Byte, (Z(i), i=1,NN) !get appended array
              print*, N_Byte
            endif
          endif
        endif
  
negative.code:               trim(adjustlt(Upper_Case(type)))/='FLOAT64') then
              E_IO = -1_I4P 
            else
              allocate(Y(ny1:ny2), stat=E_IO)
              read(data, fmt=*, iostat=E_IO) (Y(i), i=ny1,ny2) !get Y ascii array
            endif
            E_IO = search(from=pos, inside='Coordinates', to_find='DataArray', with_attribute='Name', of_value='Z', &
                          buffer=s_buffer,content=data)
  