-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathFig.2_new.R
345 lines (313 loc) · 14 KB
/
Fig.2_new.R
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
setwd("C:/Users/Ginevra/Dropbox/BlackSea2/nuovi_script")
Clean_data<-read.table('Dataset_Finale_new_no_bosforo.txt', header=T)
str(Clean_data)
setwd("C:/Users/Ginevra/Dropbox/BlackSea2/implementazione/BlackSea_IIDraft/submission_figures")
tiff('f2.tiff', height = 25, width = 23, units = 'cm',
compression = "lzw", res = 300)
# mai make margins smaller mai=c(0.1,0.1,0.2,0.1),
par(mfrow=c(2,3),mar=c(2.3,4.5,3.3,1.4))
##par(mar=c( bottom , left , top, right )
#oma=c(1,1,1,4)
summary(Clean_data$Temp)
a<-seq(6,27, by=3)
par(bty='n')
#fig=c(x1, x2, y1, y2)
plot(Clean_data$Temp, Clean_data$sigma,xlim=c(6,27),
xaxt="n",cex=2.2,yaxt='n',
ylim=rev(range(Clean_data$sigma)), col="#fee08b00",
ylab='', cex.lab=2,
xlab=" ", pch=17)
rect(6, 10.5 , 27,15.64 , density = NULL, angle = 45,
col = '#abd9e922', border = '#abd9e9', lty = NULL, lwd = par("lwd"),
xpd = NULL)
rect(6, 14.25 , 27,15.285 , density = NULL, angle = 45,
col = '#ffffd922', border = '#e8e888', lty = NULL, lwd = par("lwd"),
xpd = NULL)
rect(6, 16.2 , 27,17.04 , density = NULL, angle = 45,
col = '#eadede22', border ='#d2c7c7', lty = NULL, lwd = par("lwd"),
xpd = NULL)
rect(6, 15.64 , 27,16.2 , density = NULL, angle = 45,
col = '#1d932722', border = '#35978f', lty = NULL, lwd = par("lwd"),
xpd = NULL)
rect(6, 17.04 , 27,17.2 , density = NULL, angle = 45,
col = '#a39b9b22', border = '#a39b9b', lty = NULL, lwd = par("lwd"),
xpd = NULL)
rect(6, 17.23 , 27,17.25 , density = NULL, angle = 45,
col = '#25232322', border = '#252323', lty = NULL, lwd = par("lwd"),
xpd = NULL)
par(new=T)
plot(Clean_data$Temp, Clean_data$sigma,xlim=c(6,27),
xaxt="n",cex=1.3,yaxt='n',
ylim=rev(range(Clean_data$sigma)), col="#fee08b",
ylab=' ', cex.lab=2,
xlab=" ", pch=24)
axis(3, at=a,line=0.3, col="darkgoldenrod1", cex.axis=1.6, cex.axis=1.6)
mtext("T \n (°C)", 3, line=0.3,at=1.5,col="darkgoldenrod1", cex=1.4)
text(16,12,'OL', cex=.9, font=2)
text(22,14.7,'CIL', cex=.9, font=2)
text(14,15.9,'SOL', cex=.9, font=2)
text(14,16.7,'AOL', cex=.9, font=2)
#text(13,22,"1", col="red", cex=1.4)
par(new=TRUE)
summary(Clean_data$Sal)
s<-seq(17,23,by=1)
plot(Clean_data$Sal, cex=2.2, Clean_data$sigma, ylim=rev(range(Clean_data$sigma)),
pch=4, col="#80cdc1",frame=FALSE, xaxt='n',yaxt='n',xlab="",ylab="",
xlim=c(17,23))
axis(1, at=s,line=.3, col="#35978f", cex.axis=1.6)
mtext("S",1, line=.5, col="#35978f", cex=1.4, at=15.5)
y2<-c(10.5,11.5,12.5,13.5,14.5,15.5,
16.5,17, 17.5)
y<-c(11,12,13,14,15,16,17)
axis(2, at=y,line=0, col="black", cex.axis=1.6 )
axis(2, at=y2,line=0,labels=F, col="black" , tck =-.02, cex.axis=1.6)
mtext(expression(paste(sigma[theta]*' (kg m'^-3*')')), 2, line=2.3,at=13.5,
col="black", cex=1.4)
#mtext("Sal",1, line=3, col="#35978f", cex=1.4, at=24)
#secondo---II-O2 HS mar=c(5.5,1,.5,1)
summary(Clean_data$Ox_uM)
Ox<-Clean_data[(Clean_data$Ox_uM>1),]
d <- seq(0,350,by=50)
plot(Ox$Ox_uM, Ox$sigma, ylim=rev(range(Clean_data$sigma)),
col="#abd9e900", frame=FALSE, xlim=c(0,350),cex=2.2,
xaxt='n', yaxt='n',xlab="", ylab='',
cex.lab=2, pch=8)
axis(3,line=0.3, col="#72b3c9",at=d, cex.axis=1.6)
#Lines <- list(bquote(), bquote())
mtext(expression(paste('O'[2]))
,side=3,
line=c(1),at=-90,col="#72b3c9", cex=1.4)
mtext("(µM)",side=3,
line=c(-0.3),at=-90,col="#72b3c9", cex=1.4)
rect(0, 10.5 , 350,15.64 , density = NULL, angle = 45,
col = '#abd9e922', border = '#abd9e9', lty = NULL, lwd = par("lwd"),
xpd = NULL)
rect(0, 14.25 , 350,15.285 , density = NULL,angle = 45,
col = '#ffffd933', border = '#e8e888', lty = NULL, lwd = par("lwd"),
xpd = NULL)
rect(0, 16.2 , 350,17.04 , density = NULL, angle = 45,
col = '#eadede22', border = '#d2c7c7', lty = NULL, lwd = par("lwd"),
xpd = NULL)
rect(0, 15.64 , 350,16.2 , density = NULL, angle = 45,
col = '#1d932722', border = '#35978f', lty = NULL, lwd = par("lwd"),
xpd = NULL)
rect(0, 17.04 , 350,17.2 , density = NULL, angle = 45,
col = '#a39b9b22', border = '#a39b9b', lty = NULL, lwd = par("lwd"),
xpd = NULL)
rect(0, 17.23 , 350,17.25 , density = NULL, angle = 45,
col = '#25232322', border = '#252323', lty = NULL, lwd = par("lwd"),
xpd = NULL)
par(new=T)
plot(Ox$Ox_uM, Ox$sigma, ylim=rev(range(Clean_data$sigma)),
col="#abd9e9",frame=FALSE, xlim=c(0,350),cex=2.2,
xaxt='n',yaxt='n', xlab="", ylab=' ', cex.lab=2,pch=8)
text(50,12,'OL', cex=.9, font=2)
text(50,14.6,'CIL', cex=.9, font=2)
text(220,15.9,'SOL', cex=.9, font=2)
text(250,16.6,'AOL', cex=.9, font=2)
#title(sub="µM", line=2.5)
par(new=TRUE)
HS<-Clean_data[(Clean_data$HS>0.1),]
summary(HS$HS)5
g <- seq(0,500,by=100)
plot(HS$HS, HS$sigma, ylim=rev(range(Clean_data$sigma)),xlim=c(0,500),cex=1.3,
col="#fb6a4a",frame=FALSE,xaxt='n', yaxt='n', xlab="", pch=25,ylab="")
axis(1,at=g,line=.3, col="#fb6a4a", cex.axis=1.6)
Lines <- list(bquote(), bquote("(µM)"))
mtext(expression(paste('HS'^-)),side=1,line=1,at=-90,col="#fb6a4a", cex=1.4)
mtext(expression(paste('µM')),side=1,line=-0.3,at=-90,col="#fb6a4a", cex=1.4))
# line=c(-0.3,1),at=-90,col="#fb6a4a", cex=1.4)
axis(2, at=y,line=0, col="black" , cex.axis=1.6)
axis(2, at=y2,line=0,labels=F, col="black" , tck =-.02, cex.axis=1.6)#mtext(expression(paste(sigma[theta]*' (kg m'^-3*')')),1, line=3, col="#35978f", cex=1.4, at=24)
mtext(expression(paste(sigma[theta]*' (kg m'^-3*')')), 2,
line=2.3,at=13.5,col="black", cex=1.4)
#--III --BAC + Mn,mar=c(5.5,1,.5,5)
e <- seq(0,5000,by=1000)
Mn<-Clean_data[(Clean_data$Mn>0.1),]
plot(Mn$Mn, Mn$sigma, ylim=rev(range(Clean_data$sigma)), xlim=c(0,5000),
col="#88419d00", bg="#d4b9da00", cex.lab=2, cex=2.2,
frame=FALSE,xaxt='n',yaxt='n', pch=24,xlab='',
ylab='')
axis(1,at=e,line=.3, col="#e1e22f", cex.axis=1.6)
Lines <- list(bquote('Mn'[D]), bquote("(nM)"))
mtext(do.call(expression, Lines),side=1,
line=c(-0.3,1),at=-1400,col="#e1e22f", cex=1.4)
rect(0, 10.5 , 5000,15.64 , density = NULL, angle = 45,
col = '#abd9e922', border = '#abd9e9', lty = NULL, lwd = par("lwd"),
xpd = NULL)
rect(0, 14.25 , 5000,15.285 , density = NULL, angle = 45,
col = '#ffffd933', border = '#e8e888', lty = NULL, lwd = par("lwd"),
xpd = NULL)
rect(0, 16.2 , 5000,17.04 , density = NULL, angle = 45,
col = '#eadede22', border = '#d2c7c7', lty = NULL, lwd = par("lwd"),
xpd = NULL)
rect(0, 15.64 , 5000,16.2 , density = NULL, angle = 45,
col = '#1d932722', border = '#35978f', lty = NULL, lwd = par("lwd"),
xpd = NULL)
rect(0, 17.04 , 5000,17.2 , density = NULL, angle = 45,
col = '#a39b9b22', border = '#a39b9b', lty = NULL, lwd = par("lwd"),
xpd = NULL)
rect(0, 17.23 , 5000,17.25 , density = NULL, angle = 45,
col = '#25232322', border = '#252323', lty = NULL, lwd = par("lwd"),
xpd = NULL)
par(new=TRUE)
ba<-seq(0.,0.5,by=0.1)
plot(Clean_data$BAC, Clean_data$sigma, ylim=rev(range(Clean_data$sigma)),
ylab=" ", cex=1.5, xlim=c(0,0.5), yaxt='n',lwd=1.5,
pch=21,col="#006837", frame=FALSE, xaxt='n',xlab=" ")
axis(3, at=ba,line=0.3, col="#006837", cex.axis=1.6)
mtext(expression(paste(' BAC \n (m'^-1*')')),3,line=0,at=-.13, col="#006837", cex=1.4)
text(0.45,12,'OL', cex=1.4, font=2)
text(0.45,14.7,'CIL', cex=.9, font=2)
text(0.45,15.9,'SOL', cex=.9, font=2)
text(.15,16.7,'AOL', cex=.9, font=2)
par(new=T)
plot(Mn$Mn, Mn$sigma, ylim=rev(range(Clean_data$sigma)), xlim=c(0,5000),
col="#e1e22f", bg="#d5f41900", cex.lab=2, cex=2.2,
frame=FALSE,xaxt='n',yaxt='n', pch=24, xlab='',lwd=1.5,
ylab='')
#mtext("dMn (nM)",1,line=3, at=c(7500),col="darkorange", cex=1)
axis(2, at=y,line=0, col="black", cex.axis=1.6 )
axis(2, at=y2,line=0,labels=F, col="black" , tck =-.02, cex.axis=1.6)
mtext(expression(paste(sigma[theta]*' (kg m'^-3*')')), 2,
line=2.3,at=13.5,col="black", cex=1.4)
#-IV ---NO3 + Fe
summary(Clean_data$NO3)
no3<-Clean_data[(Clean_data$NO3>0.1),]
plot(no3$NO3, no3$sigma, cex=2.2,xlim=c(0,5),
col="#66c2a400", lwd=1.3,frame=F, ylim=rev(range(Clean_data$sigma)),
ylab='', cex.lab=2,xaxt='n',yaxt='n', xlab="", pch=24)
rect(0, 10.5, 5, 15.64 , density = NULL, angle = 45,
col = '#abd9e922', border = '#abd9e9', lty = NULL, lwd = par("lwd"),
xpd = NULL)
rect(0, 14.25 , 5,15.28 , density = NULL, angle = 45,
col = '#ffffd922', border = '#e8e888', lty = NULL, lwd = par("lwd"),
xpd = NULL)
rect(0, 16.2 , 5,17.04 , density = NULL, angle = 45,
col = '#eadede22', border = '#d2c7c7', lty = NULL, lwd = par("lwd"),
xpd = NULL)
rect(0, 15.64 , 5,16.2 , density = NULL, angle = 45,
col = '#1d932722', border = '#35978f', lty = NULL, lwd = par("lwd"),
xpd = NULL)
rect(0, 17.04 ,5, 17.2 , density = NULL, angle = 45,
col = '#a39b9b22', border = '#a39b9b', lty = NULL, lwd = par("lwd"),
xpd = NULL)
rect(0, 17.235 , 5,17.25 , density = NULL, angle = 45,
col = '#25232322', border = '#252323', lty = NULL, lwd = par("lwd"),
xpd = NULL)
par(new=T)
plot(no3$NO3, no3$sigma, cex=2.1,xlim=c(0,5),
col="#66bd63", lwd=1.3,frame=F, ylim=rev(range(Clean_data$sigma)),
ylab='', cex.lab=2, xaxt='n',yaxt='n', xlab="", pch=4)
no3a<-seq(0,5,by=1)
axis(3,at=no3a,line=.3, col="#66bd63", cex.axis=1.6)
Lines <- list(bquote('NO'[3]), bquote("(µM)"))
mtext(do.call(expression, Lines),side=3,
line=c(1,-.4),at=c(-1.3), col="#66bd63", cex=1.4)
par(new=TRUE)
summary(Clean_data$Fe)
fe <- seq(0,500,by=100)
Fe<-Clean_data[(Clean_data$Fe>0.1),]
plot(Fe$Fe, Fe$sigma, ylim=rev(range(Clean_data$sigma)), cex=2.1,
lwd=1.4, xlim=c(0,500),
col="#f49b42", frame=FALSE, xaxt='n',yaxt='n',xlab="", pch=23,ylab="")
axis(1,at=fe,line=.3, col="#f49b42", cex.axis=1.6)
Lines <- list(bquote('Fe'[D]), bquote("(nM)"))
mtext(do.call(expression, Lines),side=1,
line=c(-.7,.8),at=c(-110), col="#f49b42", cex=1.4)
axis(2, at=y,line=0, col="black", cex.axis=1.6)
axis(2, at=y2,line=0,labels=F, col="black" , tck =-.02, cex.axis=1.6)
mtext(expression(paste(sigma[theta]*' (kg m'^-3*')')), 2,
line=2.3,at=14.5,col="black", cex=1.4)
text(100,12,'OL', cex=.9, font=2)
text(440,14.7,'CIL', cex=.9, font=2)
text(440,15.9,'SOL', cex=.9, font=2)
text(440,16.8,'AOL', cex=.9, font=2)
#-- V --- Hg
summary(Clean_data$Hg)
par(bty="n")
plot(Clean_data$Hg, Clean_data$sigma, xlim=c(0,5),
ylim=rev(range(Clean_data$sigma)), pch=15, xaxt='n', yaxt='n',
type="p",lty=4, col="#3690c000", cex=2.2,
ylab='', main=" ",
cex.lab=2, xlab=" " )
rect(0, 10.5 , 5,15.64 , density = NULL, angle = 45,
col = '#abd9e922', border = '#abd9e9', lty = NULL, lwd = par("lwd"),
xpd = NULL)
rect(0, 14.25 , 5,15.28 , density = NULL, angle = 45,
col = '#ffffd922', border = '#e8e888', lty = NULL, lwd = par("lwd"),
xpd = NULL)
rect(0, 16.2 , 5,17.04 , density = NULL, angle = 45,
col = '#eadede22', border = '#d2c7c7', lty = NULL, lwd = par("lwd"),
xpd = NULL)
rect(0, 15.64 , 5,16.2 , density = NULL, angle = 45,
col = '#1d932722', border = '#35978f', lty = NULL, lwd = par("lwd"),
xpd = NULL)
rect(0, 17.04 , 5,17.2 , density = NULL, angle = 45,
col = '#a39b9b22', border = '#a39b9b', lty = NULL, lwd = par("lwd"),
xpd = NULL)
rect(0, 17.23 , 5,17.25 , density = NULL, angle = 45,
col = '#25232322', border = '#252323', lty = NULL, lwd = par("lwd"),
xpd = NULL)
par(new=T)
plot(Clean_data$Hg, Clean_data$sigma, xlim=c(0,5),
ylim=rev(range(Clean_data$sigma)), pch=21, xaxt='n', yaxt='n',
type="p",lty=4, bg="#3690c088",col='#256486', cex=2.2,
ylab='', main=" ",
cex.lab=2, xlab=" " )
m<-seq(0,5,by=1)
axis(3,at=m,line=0.3, col="#3690c0", cex.axis=1.6)
Lines <- list(bquote('Hg'[D]), bquote("(pM)"))
mtext(do.call(expression, Lines),side=3,
line=c(1,-0.3),at=c(-1.2),col="#3690c0", cex=1.4)
axis(2, at=y,line=0, col="black", cex.axis=1.6)
axis(2, at=y2,line=0,labels=F, col="black" , tck =-.02, cex.axis=1.6)
mtext(expression(paste(sigma[theta]*' (kg m'^-3*')')), 2, line=2.3,
at=13.5,col="black", cex=1.4)
text(4,12,'OL', cex=.9, font=2)
text(4,14.7,'CIL', cex=.9, font=2)
text(4,15.9,'SOL', cex=.9, font=2)
text(1,16.6,'AOL', cex=.9, font=2)
# vI ---- Mehg----- mar=c(5.5,1,.5,5)
par(bty='n')
mm<-seq(0,1.2, by=0.2)
plot(Clean_data$MeHg, Clean_data$sigma, xlim=c(0,1.2),xaxt='n', yaxt='n',
ylim=rev(range(Clean_data$sigma)), pch=23, cex=2.2,
type="p",lty=4, col="#df65b000", bg="#df65b000", main=" ",
ylab='',cex.lab=2,xlab=" ")
rect(0, 10.5 , 1.2,15.64 , density = NULL, angle = 45,
col = '#abd9e922', border = '#abd9e9', lty = NULL, lwd = par("lwd"),
xpd = NULL)
rect(0, 14.25 , 1.2,15.28 , density = NULL, angle = 45,
col = '#ffffd922', border = '#e8e888', lty = NULL, lwd = par("lwd"),
xpd = NULL)
rect(0, 16.2 , 1.2,17.04 , density = NULL, angle = 45,
col = '#eadede22', border = '#d2c7c7', lty = NULL, lwd = par("lwd"),
xpd = NULL)
rect(0, 15.64 , 1.2,16.2 , density = NULL, angle = 45,
col = '#1d932722', border = '#35978f', lty = NULL, lwd = par("lwd"),
xpd = NULL)
rect(0, 17.04 , 1.2,17.2 , density = NULL, angle = 45,
col = '#a39b9b22', border = '#a39b9b', lty = NULL, lwd = par("lwd"),
xpd = NULL)
rect(0, 17.23 , 1.2,17.25 , density = NULL, angle = 45,
col = '#25232322', border = '#252323', lty = NULL, lwd = par("lwd"),
xpd = NULL)
par(new=T)
plot(Clean_data$MeHg, Clean_data$sigma, xlim=c(0,1.2),xaxt='n', yaxt='n',
ylim=rev(range(Clean_data$sigma)), pch=22, cex=2.2,
type="p",lty=4, col="#e7298a", bg="#df65b088", main=" ",
ylab='', cex.lab=2, xlab=" ")
axis(3,at=mm,line=0.3, col="#df65b0", cex.axis=1.6)
Lines <- list(bquote('MeHg'[D]), bquote(" \n (pM)"))
mtext(do.call(expression, Lines),side=3,
line=c(1,-0.3),at=c(-.35),col="#df65b0", cex=1.4)
axis(2, at=y,line=0, col="black", cex.axis=1.6)
axis(2, at=y2,line=0,labels=F, col="black" , tck =-.02, cex.axis=1.6)
mtext(expression(paste(sigma[theta]*' (kg m'^-3*')')), 2, line=2.3,
at=13.5,col="black", cex=1.4)
text(.8,12,'OL', cex=.9, font=2)
text(.8,14.7,'CIL', cex=.9, font=2)
text(.95,15.9,'SOL', cex=.9, font=2)
text(.15,16.6,'AOL', cex=.9, font=2)
dev.off()