pro fa_append_1d,array1,array2,index1,index2,array,index,max_index,map
ndimensions1=size(array1,/dimensions)
ndimensions2=size(array2,/dimensions)
if ndimensions1[0] NE ndimensions2[0] then begin
print,'Error: Dimensions Do Not Agree'
return
endif
if n_elements(ndimensions1) EQ 1 then ndimensions1=[ndimensions1[0],1]
if n_elements(ndimensions2) EQ 1 then ndimensions2=[ndimensions2[0],1]
max_index=ndimensions1[1]-1
mod_index2=index2
map=intarr(ndimensions2[1])
array=replicate(array1[0],ndimensions1[0],ndimensions1[1]+ndimensions2[1])
array[*,0:max_index]=array1
for iii=0,ndimensions2[1]-1 do begin
flag=1
for jjj=0,ndimensions1[1]-1 do begin
test=equal_arrays(array2[*,iii],array1[*,jjj],tolerance=.001,/silence)
if test EQ 1 then begin
map[iii]=jjj
wherearray=where(index2 EQ iii)
if wherearray[0] NE -1 then mod_index2[wherearray]=jjj
flag=0
break
endif
endfor
if flag then begin
++max_index
map[iii]=max_index
wherearray=where(index2 EQ iii)
if wherearray[0] NE -1 then mod_index2[wherearray]=max_index
array[*,max_index]=array2[*,iii]
endif
endfor
index=[index1,mod_index2]
array=array[*,0:max_index]
return
end