anchor.diff:      write(unit=vtk(rf)%u,iostat=E_IO)repeat(' ',vtk(rf)%indent)//'<Cells>'//end_rec ; vtk(rf)%indent = vtk(rf)%indent + 2
       write(unit=vtk(rf)%u,iostat=E_IO)repeat(' ',vtk(rf)%indent)//&
                                        '<DataArray type="Int32" Name="connectivity" format="binary">'//end_rec
  +    allocate(cocp(1:size(transfer([int(offset(NC)*BYI4P,I4P),connect],cocp))))
       cocp = transfer([int(offset(NC)*BYI4P,I4P),connect],cocp)
       call b64_encode(nB=int(BYI1P,I4P),n=cocp,code=coc64)
  +    deallocate(cocp)
       write(unit=vtk(rf)%u,iostat=E_IO)repeat(' ',vtk(rf)%indent+2)//tochar(coc64)//end_rec
       write(unit=vtk(rf)%u,iostat=E_IO)repeat(' ',vtk(rf)%indent)//'</DataArray>'//end_rec
       write(unit=vtk(rf)%u,iostat=E_IO)repeat(' ',vtk(rf)%indent)//'<DataArray type="Int32" Name="offsets" format="binary">'//end_rec
  +    allocate(cocp(1:size(transfer([int(NC*BYI4P,I4P),offset],cocp))))
       cocp = transfer([int(NC*BYI4P,I4P),offset],cocp)
       call b64_encode(nB=int(BYI1P,I4P),n=cocp,code=coc64)
  +    deallocate(cocp)
       write(unit=vtk(rf)%u,iostat=E_IO)repeat(' ',vtk(rf)%indent+2)//tochar(coc64)//end_rec
       write(unit=vtk(rf)%u,iostat=E_IO)repeat(' ',vtk(rf)%indent)//'</DataArray>'//end_rec
       write(unit=vtk(rf)%u,iostat=E_IO)repeat(' ',vtk(rf)%indent)//'<DataArray type="Int8" Name="types" format="binary">'//end_rec
  
positive.code:     s_buffer=repeat(' ',vtk(rf)%indent)//'<DataArray type="Int32" Name="'//trim(varname)// &
               '" NumberOfComponents="1" format="binary">'
      write(vtk(rf)%u,iostat=E_IO)trim(s_buffer)//end_rec
      varp = transfer([int(NC_NN*BYI4P,I4P),var],varp)
      call b64_encode(nB=int(BYI1P,I4P),n=varp,code=var64) ; deallocate(varp)
      write(vtk(rf)%u,iostat=E_IO)repeat(' ',vtk(rf)%indent+2)//tochar(var64)//end_rec ; deallocate(var64)
  
negative.code:     do n1=1,NC_NN
        var(1+(n1-1)*3:1+(n1-1)*3+2)=[varX(n1),varY(n1),varZ(n1)]
      enddo
      allocate(varp(1:size(transfer([int(3*NC_NN*BYI4P,I4P),var],varp))))
      varp = transfer([int(3*NC_NN*BYI4P,I4P),var],varp) ; deallocate(var)
      call b64_encode(nB=int(BYI1P,I4P),n=varp,code=var64) ; deallocate(varp)
      write(vtk(rf)%u,iostat=E_IO)repeat(' ',vtk(rf)%indent+2)//tochar(var64)//end_rec ; deallocate(var64)
  