Reorganisation

Author: dadevlee

Megasplice construction/Basin_corrections_Cramer_carbon.csv

@@ -0,0 +1,85 @@
+setwd("CHANGE TO YOUR DIRECTORY")
+library(astrochron)
+library(quantmod)
+rm(list=ls())
+
+Basin_correction=read.csv("Basin_corrections_Cramer_carbon.csv")
+Pacific_d13C_trend=read.csv("Pacific_d13C_trend_carbon.csv",header = T)
+
+
+Site926=read.csv("Site_926.csv")
+colnames(Site926) <- c("Age (Ma)", "Benthic d13C")
+#Site982_Andersson=read.csv("Site_982_Andersson2003_ETP.csv")
+#colnames(Site982_Andersson) <- c("Age (Ma)", "Benthic d13C")
+Site982=read.csv("Site_982_Drury_wuellstorfi.csv",header = T)
+colnames(Site982) <- c("Age (Ma)", "Benthic d13C")
+Site1090=read.csv("Site_1090.csv")
+colnames(Site1090) <- c("Age (Ma)", "Benthic d13C")
+Site1146=read.csv("Site_1146_Holbourn2018_AJD.csv")
+colnames(Site1146) <- c("Age (Ma)", "Benthic d13C")
+# Site1146=read.csv("Site_1146_Holbourn2018.csv")
+# colnames(Site1146) <- c("Age (Ma)", "Benthic d13C")
+Site1218=read.csv("Site_1218.csv")
+colnames(Site1218) <- c("Age (Ma)", "Benthic d13C")
+idx=which(Site1218$`Age (Ma)`<35)
+Site1218=Site1218[idx,]
+Site1337=read.csv("Site_1337.csv")
+colnames(Site1337) <- c("Age (Ma)", "Benthic d13C")
+Site1338=read.csv("Site_1338.csv")
+colnames(Site1338) <- c("Age (Ma)", "Benthic d13C")
+Site1264=read.csv("Site_1264.csv")
+colnames(Site1264) <- c("Age (Ma)", "Benthic d13C")
+Site1267=read.csv("Site_1267.csv")
+colnames(Site1267)<- c("Age (Ma)", "Benthic d13C")
+Site1143=read.csv("Site_1143.csv")
+colnames(Site1143)<- c("Age (Ma)", "Benthic d13C")
+Site1337_Drury=read.csv("Site_1337_Drury.csv")
+colnames(Site1337_Drury)<- c("Age (Ma)", "Benthic d13C")
+
+###########################
+# Bring to Pacific trend
+###########################
+
+for (name in c("Site926", "Site982", "Site1090","Site1146","Site1218","Site1337","Site1338","Site1264","Site1267","Site1143","Site1337_Drury")){
+
+b=c()
+eval(parse(text=paste('output=mwStats(sortNave(',name,'),win=0.5,conv = 1,genplot=F)')))
+for (i in 1:nrow(output)){
+  a=which(abs(Pacific_d13C_trend$Age..GTS2004.-output[i,1])==min(abs(Pacific_d13C_trend$Age..GTS2004.-output[i,1])))  
+  b[i]=output[i,2]-Pacific_d13C_trend$Pacific.d13Ctrend[a[1]]
+}
+b=cbind(output$Center_win,b)
+colnames(b)<- c("Center_win","Offset_record_Pacific_trend")
+eval(parse(text=paste('N=nrow(',name,')')))
+for (j in 1:N) {
+  eval(parse(text=paste('a=which(abs(b[,1]-',name,'[j,1])==min(abs(b[,1]-',name,'[j,1])))')))  
+  eval(parse(text=paste(name,'[j,3]=',name,'[j,2]-b[a[1],2]')))
+  }
+eval(parse(text=paste('colnames(',name,') <- c("Age (Ma)", "Benthic d13C","Benthic d13C Pacific")')))
+eval(parse(text=paste0('write.csv(',name,',"Pacific_corrected_data/',name,'_Pacific_corr.csv",row.names = F)')))
+}
+
+####################################
+# Concatenation of records
+####################################
+
+BP=c(3.2844,5.174,6.677,7.891,12.834,15.8,20.025,20.215,22.56,23.9)
+
+pos <- which(Site1267[,1]>=0 & Site1267[,1]<=BP[1])
+pos1 <- which(Site1264[,1]>=BP[1] & Site1264[,1]<=BP[2])
+pos2 <- which(Site982[,1]>=BP[2] & Site982[,1]<=BP[3])
+pos3 <- which(Site1337_Drury[,1]>=BP[3] & Site1337_Drury[,1]<=BP[4])
+#pos4 <- which(Site1143[,1]>=BP[4] & Site1143[,1]<=BP[5])
+pos5 <- which(Site1146[,1]>=BP[4] & Site1146[,1]<=BP[5])
+pos6 <- which(Site1338[,1]>=BP[5] & Site1338[,1]<=BP[6])
+pos7 <- which(Site1337[,1]>=BP[6] & Site1337[,1]<=BP[7])
+pos8 <- which(Site1090[,1]>=BP[7] & Site1090[,1]<=BP[8])
+pos9 <- which(Site926[,1]>=BP[8] & Site926[,1]<=BP[9])
+pos10 <- which(Site1090[,1]>=BP[9] & Site1090[,1]<=BP[10])
+pos11 <- which(Site1218[,1]>=BP[10] & Site1218[,1]<=35)
+
+Splice=rbind(Site1267[pos,c(1,3)],Site1264[pos1,c(1,3)],Site982[pos2,c(1,3)], Site1337_Drury[pos3,c(1,3)],Site1146[pos5,c(1,3)], Site1338[pos6,c(1,3)], Site1337[pos7,c(1,3)],Site1090[pos8,c(1,3)],Site926[pos9,c(1,3)],Site1090[pos10,c(1,3)],Site1218[pos11,c(1,3)])
+Splice=Splice[is.finite(Splice[,2]), ]
+plot(Splice,type="l")
+
+write.csv(Splice,"d13C_megasplice.csv",row.names=F)

Megasplice construction/MegaSplice_carbon_construction.R

@@ -0,0 +1 @@
+Age (GTS2004),Pacific d13Ctrend0,-0.1840.1,-0.1910.2,-0.1920.3,-0.1940.4,-0.20.5,-0.2120.6,-0.230.7,-0.2510.8,-0.270.9,-0.2791,-0.2751.1,-0.2621.2,-0.2441.3,-0.2261.4,-0.2151.5,-0.2171.6,-0.231.7,-0.2481.8,-0.2681.9,-0.2842,-0.2962.1,-0.32.2,-0.2962.3,-0.2812.4,-0.2542.5,-0.2132.6,-0.1682.7,-0.1382.8,-0.122.9,-0.1073,-0.0963.1,-0.0863.2,-0.0763.3,-0.0643.4,-0.053.5,-0.0373.6,-0.0273.7,-0.023.8,-0.0243.9,-0.0444,-0.0744.1,-0.1074.2,-0.144.3,-0.174.4,-0.1954.5,-0.2184.6,-0.244.7,-0.2624.8,-0.2844.9,-0.2845,-0.2285.1,-0.1415.2,-0.0665.3,-0.0145.4,0.0215.5,0.0415.6,0.0585.7,0.075.8,0.075.9,0.0586,0.0426.1,0.0276.2,0.0186.3,0.0136.4,0.0096.5,0.0216.6,0.056.7,0.096.8,0.1376.9,0.1787,0.2067.1,0.2327.2,0.2567.3,0.2817.4,0.3087.5,0.3397.6,0.3737.7,0.4087.8,0.4397.9,0.4648,0.4828.1,0.4968.2,0.5048.3,0.5058.4,0.5098.5,0.5258.6,0.568.7,0.6128.8,0.6678.9,0.7129,0.7419.1,0.7549.2,0.7549.3,0.7449.4,0.7339.5,0.7219.6,0.7049.7,0.6849.8,0.6679.9,0.65310,0.64110.1,0.62910.2,0.61410.3,0.59610.4,0.57710.5,0.5610.6,0.54910.7,0.54810.8,0.55710.9,0.5711,0.58511.1,0.60211.2,0.62611.3,0.65811.4,0.69811.5,0.74211.6,0.7811.7,0.80811.8,0.82711.9,0.83612,0.8412.1,0.8512.2,0.86812.3,0.88812.4,0.90412.5,0.91412.6,0.91512.7,0.91512.8,0.92112.9,0.93513,0.95213.1,0.97613.2,1.00813.3,1.04713.4,1.09513.5,1.15613.6,1.21913.7,1.27413.8,1.31913.9,1.34914,1.36414.1,1.37314.2,1.38514.3,1.39914.4,1.41614.5,1.43814.6,1.45914.7,1.46814.8,1.46814.9,1.46515,1.4615.1,1.4615.2,1.46415.3,1.4715.4,1.47315.5,1.4715.6,1.46115.7,1.44915.8,1.43215.9,1.40916,1.38616.1,1.36416.2,1.34416.3,1.32516.4,1.30416.5,1.27216.6,1.22516.7,1.16316.8,1.09416.9,1.02417,0.9617.1,0.90517.2,0.85217.3,0.79417.4,0.74217.5,0.70617.6,0.68717.7,0.68517.8,0.69917.9,0.71618,0.72318.1,0.72618.2,0.73118.3,0.74118.4,0.75118.5,0.75818.6,0.75718.7,0.75118.8,0.74718.9,0.74919,0.7619.1,0.77519.2,0.78919.3,0.80219.4,0.81619.5,0.82719.6,0.83619.7,0.83619.8,0.82219.9,0.79920,0.77820.1,0.76520.2,0.76320.3,0.76920.4,0.77620.5,0.78720.6,0.80220.7,0.82420.8,0.85620.9,0.89421,0.9221.1,0.92821.2,0.91921.3,0.90421.4,0.90421.5,0.92921.6,0.97621.7,1.03321.8,1.07821.9,1.10822,1.12422.1,1.12722.2,1.12422.3,1.12722.4,1.13222.5,1.1422.6,1.16222.7,1.18922.8,1.20322.9,1.19423,1.14923.1,1.06723.2,0.97623.3,0.89523.4,0.82723.5,0.77823.6,0.7523.7,0.73323.8,0.71823.9,0.70824,0.70324.1,0.724.2,0.69224.3,0.6824.4,0.66424.5,0.6524.6,0.63724.7,0.63124.8,0.63324.9,0.63925,0.64625.1,0.65425.2,0.65925.3,0.66325.4,0.66725.5,0.67225.6,0.67925.7,0.68825.8,0.69525.9,0.70126,0.70426.1,0.70326.2,0.70226.3,0.70126.4,0.726.5,0.69726.6,0.68826.7,0.67326.8,0.65426.9,0.63427,0.61627.1,0.60327.2,0.59727.3,0.59827.4,0.59927.5,0.59827.6,0.58827.7,0.5627.8,0.51627.9,0.4728,0.42928.1,0.40228.2,0.39428.3,0.39828.4,0.40528.5,0.41228.6,0.4228.7,0.43128.8,0.44428.9,0.45229,0.45429.1,0.44729.2,0.43229.3,0.41429.4,0.39729.5,0.38429.6,0.37929.7,0.38829.8,0.4129.9,0.44330,0.47930.1,0.50830.2,0.52530.3,0.5330.4,0.5330.5,0.52930.6,0.53630.7,0.55330.8,0.58230.9,0.62131,0.66131.1,0.68831.2,0.69931.3,0.69931.4,0.69231.5,0.68831.6,0.69331.7,0.70331.8,0.71731.9,0.73332,0.75132.1,0.77432.2,0.832.3,0.82532.4,0.84432.5,0.85432.6,0.85332.7,0.84732.8,0.84332.9,0.8533,0.87533.1,0.92133.2,0.98433.3,1.05433.4,1.11533.5,1.15133.6,1.15533.7,1.13333.8,1.09533.9,1.04234,0.97934.1,0.91534.2,0.85834.3,0.81234.4,0.78734.5,0.77834.6,0.7734.7,0.74534.8,0.70934.9,0.69635,0.67835.1,0.6635.2,0.65135.3,0.65535.4,0.66935.5,0.68735.6,0.70435.7,0.72435.8,0.79935.9,0.8736,0.93136.1,0.97336.2,1.00536.3,0.99336.4,0.9736.5,0.93336.6,0.89336.7,0.8636.8,0.84136.9,0.83737,0.84337.1,0.8537.2,0.85537.3,0.84737.4,0.82537.5,0.79937.6,0.77937.7,0.76237.8,0.7437.9,0.70138,0.64438.1,0.57738.2,0.51938.3,0.4938.4,0.49338.5,0.51238.6,0.53838.7,0.56438.8,0.58738.9,0.61339,0.65339.1,0.71739.2,0.79239.3,0.82839.4,0.80639.5,0.75339.6,0.69939.7,0.66239.8,0.64539.9,0.64740,0.6640.1,0.67640.2,0.6940.3,0.69240.4,0.67440.5,0.64140.6,0.60440.7,0.5740.8,0.54840.9,0.54441,0.55141.1,0.5641.2,0.56741.3,0.57341.4,0.57741.5,0.58441.6,0.59741.7,0.61341.8,0.62941.9,0.64442,0.65642.1,0.66142.2,0.66342.3,0.66942.4,0.67842.5,0.69242.6,0.71542.7,0.73942.8,0.7642.9,0.77743,0.78743.1,0.78843.2,0.78443.3,0.76543.4,0.7343.5,0.68543.6,0.6343.7,0.57643.8,0.5443.9,0.51344,0.48944.1,0.47444.2,0.4744.3,0.47344.4,0.47644.5,0.47244.6,0.46244.7,0.45244.8,0.45144.9,0.46145,0.47145.1,0.47445.2,0.47145.3,0.45945.4,0.44445.5,0.43445.6,0.42945.7,0.43345.8,0.45245.9,0.50446,0.59346.1,0.66846.2,0.70746.3,0.72246.4,0.72546.5,0.72746.6,0.73746.7,0.75846.8,0.79146.9,0.83347,0.86847.1,0.86647.2,0.80647.3,0.70447.4,0.61347.5,0.54747.6,0.51347.7,0.51147.8,0.54747.9,0.59448,0.65948.1,0.73948.2,0.82948.3,0.91248.4,0.94648.5,0.94548.6,0.9448.7,0.93448.8,0.92648.9,0.92249,0.92949.1,0.94449.2,0.9649.3,0.97649.4,0.98549.5,0.98349.6,0.97549.7,0.96549.8,0.95249.9,0.93750,0.91850.1,0.89750.2,0.87950.3,0.87250.4,0.85750.5,0.80750.6,0.72850.7,0.63650.8,0.53850.9,0.46351,0.42151.1,0.39851.2,0.39151.3,0.3951.4,0.38251.5,0.35351.6,0.30651.7,0.25551.8,0.21951.9,0.20252,0.19652.1,0.19752.2,0.19652.3,0.1952.4,0.1852.5,0.17352.6,0.17652.7,0.19652.8,0.2352.9,0.26753,0.29153.1,0.28653.2,0.25553.3,0.21153.4,0.16653.5,0.13553.6,0.11953.7,0.10253.8,0.05953.9,-0.00554,-0.05954.1,-0.06654.2,0.01454.3,0.17254.4,0.35254.5,0.51954.6,0.65554.7,0.74654.8,0.80754.9,0.86955,0.95355.1,1.05355.2,1.14955.3,1.22955.4,1.28255.5,1.30355.6,1.30655.7,1.31355.8,1.33555.9,1.38656,1.46756.1,1.54256.2,1.58756.3,1.6156.4,1.6256.5,1.6356.6,1.67856.7,1.77456.8,1.88756.9,1.9957,2.07657.1,2.1457.2,2.1857.3,2.20557.4,2.22157.5,2.22357.6,2.21657.7,2.20257.8,2.17957.9,2.14658,2.10358.1,2.05258.2,2.00358.3,1.95958.4,1.92758.5,1.91258.6,1.91458.7,1.92358.8,1.93258.9,1.92259,1.8959.1,1.85159.2,1.81459.3,1.78659.4,1.77159.5,1.76459.6,1.74159.7,1.68159.8,1.58959.9,1.4860,1.37160.1,1.28560.2,1.23460.3,1.19560.4,1.16360.5,1.13760.6,1.11460.7,1.09160.8,1.0760.9,1.04361,1.01461.1,0.98761.2,0.96861.3,0.9661.4,0.96261.5,0.97461.6,0.99761.7,1.03361.8,1.08661.9,1.15162,1.21262.1,1.26162.2,1.29262.3,1.29462.4,1.26462.5,1.20262.6,1.11162.7,1.01962.8,0.96962.9,0.96563,0.95963.1,0.9463.2,0.91163.3,0.88463.4,0.86663.5,0.85963.6,0.85563.7,0.84163.8,0.81463.9,0.77964,0.75364.1,0.76964.2,0.84664.3,0.96164.4,1.07864.5,1.15764.6,1.18764.7,1.19164.8,1.1964.9,1.19465,1.26765.1,1.40565.2,1.5465.3,1.6465.4,1.60365.5,1.53565.6,1.4365.7,1.33965.8,1.26465.9,1.19566,1.14366.1,1.12566.2,1.15566.3,1.22966.4,1.31166.5,1.38866.6,1.44566.7,1.47166.8,1.46266.9,1.40967,1.34667.1,1.30467.2,1.28867.3,1.29967.4,1.33467.5,1.37667.6,1.41767.7,1.4567.8,1.46967.9,1.47968,1.4968.1,1.50568.2,1.52468.3,1.54568.4,1.54868.5,1.52568.6,1.48768.7,1.44668.8,1.40968.9,1.37869,1.34769.1,1.30769.2,1.25569.3,1.20169.4,1.15169.5,1.10469.6,1.0669.7,1.01869.8,0.97269.9,0.92770,0.88970.1,0.8670.2,0.83870.3,0.81970.4,0.80170.5,0.79770.6,0.79770.7,0.80370.8,0.81370.9,0.82471,0.83771.1,0.85171.2,0.86371.3,0.87371.4,0.87771.5,0.87571.6,0.86971.7,0.8671.8,0.85671.9,0.86172,0.87572.1,0.89272.2,0.90372.3,0.90672.4,0.90772.5,0.90572.6,0.972.7,0.90772.8,0.93372.9,0.96573,1.00173.1,1.03473.2,1.0673.3,1.07673.4,1.08973.5,1.10273.6,1.12373.7,1.16373.8,1.2373.9,1.30774,1.33974.1,1.31574.2,1.27374.3,1.23374.4,1.20574.5,1.20674.6,1.20674.7,1.19874.8,1.18474.9,1.16575,1.1475.1,1.11275.2,1.11775.3,1.12775.4,1.13975.5,1.1575.6,1.15175.7,1.13375.8,1.09475.9,1.04676,1.00376.1,0.97776.2,0.97476.3,0.98376.4,0.99776.5,1.01176.6,1.02276.7,1.03576.8,1.05976.9,1.07677,1.09377.1,1.11577.2,1.13877.3,1.15377.4,77.5,77.6,77.7,77.8,77.9,78,78.1,78.2,78.3,78.4,78.5,78.6,78.7,78.8,78.9,79,79.1,79.2,79.3,79.4,79.5,79.6,79.7,79.8,79.9,