# read sediment data
insed <- scan("../1.0/OC3_Coretop1.1.csv", sep = ",", skip=1, 
               list(lat=0,lon=0,dep=0,c13=0,dataset=0,coreflag=0,specflag=0))

rad <- pi/180
latseddeg <- insed$lat		# latitude in degrees
lonseddeg <- insed$lon		# longitude in degrees
latsedrad <- latseddeg*rad	# latitude in radians
lonsedrad <- lonseddeg*rad	# longitude in radians
depsedm <- insed$dep		# depth in m
depsedkm <- depsedm/1000.	# depth in km
c13sed <- insed$c13		# d13C in permil
datasetsed <- insed$dataset	# dataset number
coreflagsed <- insed$coreflag	# coreflag
specflagsed <- insed$specflag	# speciesflag

# use only the following datasets
isd <- which(datasetsed==1 | datasetsed==2 | datasetsed==3 | datasetsed==5 | datasetsed==6 | datasetsed==7 | datasetsed==8 | datasetsed==9 | datasetsed==10 | datasetsed==11 | datasetsed==12 | datasetsed==13 | datasetsed==14 | datasetsed==15 | datasetsed==16 | datasetsed==17 | datasetsed==18)
i123 <- which(datasetsed==1 | datasetsed==2 | datasetsed==3)
iwuel <- which(specflagsed==2)

# indicees for individual datasets
i1 <- which(datasetsed==1)	# Peterson et al. (2014; P14)
i2 <- which(datasetsed==2)	# Mix unpublished (Mix)
i3 <- which(datasetsed==3)	# Duplessy et al. (1984) not included in P14 (D84)
i4 <- which(datasetsed==4)	# Matsumoto (Indo-Pacific early Holocene dataset; not used)
i5 <- which(datasetsed==5)	# Mackensen (2012; M12) Southern Ocean
i6 <- which(datasetsed==6)	# Mackensen (2013; M13) Arctic
i7 <- which(datasetsed==7)	# Cartapanis 6 ka
i8 <- which(datasetsed==8)	# Cartapanis 10 cm
i9 <- which(datasetsed==9)	# Bostock
i10 <- which(datasetsed==10)	# Oppo
i11 <- which(datasetsed==11)	# Herguera
i12 <- which(datasetsed==12)	# Lund
i13 <- which(datasetsed==13)	# Martinez-Mendez
i14 <- which(datasetsed==14)	# Lynch-Stieglitz
i15 <- which(datasetsed==15)	# Waelbroeck & Michel
i16 <- which(datasetsed==16)	# Sikes
i17 <- which(datasetsed==17)	# Hoogakker
i18 <- which(datasetsed==18)	# Repschlaeger
la1 <- latseddeg[i1]
lo1 <- lonseddeg[i1]
la2 <- latseddeg[i2]
lo2 <- lonseddeg[i2]
la3 <- latseddeg[i3]
lo3 <- lonseddeg[i3]
la4 <- latseddeg[i4]
lo4 <- lonseddeg[i4]
la5 <- latseddeg[i5]
lo5 <- lonseddeg[i5]
la6 <- latseddeg[i6]
lo6 <- lonseddeg[i6]
la7 <- latseddeg[i7]
lo7 <- lonseddeg[i7]
la8 <- latseddeg[i8]
lo8 <- lonseddeg[i8]
la9 <- latseddeg[i9]
lo9 <- lonseddeg[i9]
la10 <- latseddeg[i10]
lo10 <- lonseddeg[i10]
la11 <- latseddeg[i11]
lo11 <- lonseddeg[i11]
la12 <- latseddeg[i12]
lo12 <- lonseddeg[i12]
la13 <- latseddeg[i13]
lo13 <- lonseddeg[i13]
la14 <- latseddeg[i14]
lo14 <- lonseddeg[i14]
la15 <- latseddeg[i15]
lo15 <- lonseddeg[i15]
la16 <- latseddeg[i16]
lo16 <- lonseddeg[i16]
la17 <- latseddeg[i17]
lo17 <- lonseddeg[i17]
la18 <- latseddeg[i18]
lo18 <- lonseddeg[i18]

iwue <- which(specflagsed == 2)
ikul <- which(specflagsed == 3)
ilob <- which(specflagsed == 4)
ipac <- which(specflagsed == 5)
iflo <- which(specflagsed == 6)