-
Notifications
You must be signed in to change notification settings - Fork 2
/
a1101_00_00_Jahn2018.py
114 lines (88 loc) · 2.75 KB
/
a1101_00_00_Jahn2018.py
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
'''
Created on May 10, 2019
@author: keiin
'''
from pylab import *
from Savefig2 import *
from FigSetting2 import *
import matplotlib.ticker as ptick
#OOOOOOOOOOOOOOOOOOOOOOOOOOOO
# Data preparation
#OOOOOOOOOOOOOOOOOOOOOOOOOOOO
a = genfromtxt('C:\\Users\\keiin\\Google Drive\\0 Paper\\3 phytoplanktonmodel\\Data\\Jahn 2018\\Jahn2018Altogether.csv',delimiter=',')
MuLHC = a[:,1]
LHC = a[:,2]
MuRIB = a[:,4]
RIB = a[:,5]
MuCBM = a[:,7]
CBM = a[:,8]
#OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO
# Line preparation (Linear interpolation)
#OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO
MuFit = arange(0,0.12+0.01,0.01)
LHCfit = -1.7718*MuFit + 0.3152 #Jahn2018Altogether.xlsx
RIBfit = 0.9507*MuFit + 0.1012
CBMfit = 0.7206*MuFit + 0.0921
#OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO
# Preparing functrions
#OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO
def Dpi():
return 600
def SaveFolder():
return 'Phytoplankton model\\Jahn2018'
#OOOOOOOOOOOOOOO
# Figures
#OOOOOOOOOOOOOOO
ColorLHC = 'orange'
ColorRIB = '#00ff00'
ColorCBM = 'g'
# figure(1)
# plot(MuFit*24,LHCfit,color = ColorLHC)
# plot(MuLHC*24,LHC,'o',color = ColorLHC)
# ylim(ymin=0)
#
# figure(2)
# plot(MuFit*24,RIBfit,color = ColorRIB)
# plot(MuRIB*24,RIB,'o',color = ColorRIB)
# ylim(ymin=0)
#
# figure(3)
# plot(MuFit*24,CBMfit,color = ColorCBM)
# plot(MuCBM*24,CBM,'o',color = ColorCBM)
# ylim(ymin=0)
Dashes=(5,3)
figure(4)
plot(MuFit*24,LHCfit,dashes=Dashes,color = ColorLHC,label='LHC')
plot(MuLHC*24,LHC,'o',color = ColorLHC)
plot(MuFit*24,RIBfit,dashes=Dashes,color = ColorRIB,label='RIB')
plot(MuRIB*24,RIB,'o',color = ColorRIB)
plot(MuFit*24,CBMfit,dashes=Dashes,color = ColorCBM,label='CBM')
plot(MuCBM*24,CBM,'o',color = ColorCBM)
xlabel('$\mu_{max}$ (d$^{-1}$)')
ylabel('Protein mass fraction')
ylim(ymin=0)
title('Data: Synechocystis',y=1.02)
legend(loc=4,ncol=2,fontsize=20)
Savefig2(SaveFolder(),'DataPlot',Dpi())
# figure(1)
# xlabel('$\mu$ (d$^{-1}$)')
# ax1=gca(); ax2=ax1.twinx()
# ax1.set_zorder(ax2.get_zorder()+1) # put ax in front of ax2
# ax1.patch.set_visible(False) # hide the 'canvas'
#
# Ebar2(ax1,Carb,CarbMean,CarbSTDup,CarbSTDdown,'g','Carbohydrate')
# Ebar2(ax1,Protein,ProteinMean,ProteinSTDup,ProteinSTDdown,'r','Protein')
# Ebar2(ax2,D*nan,Ddata*nan,Ddata*nan,Ddata*nan,'r','Protein')
# Ebar2(ax2,D*nan,Ddata*nan,Ddata*nan,Ddata*nan,'g','Carb.')
# Ebar2(ax2,Lipid,LipidMean,LipidSTDup,LipidSTDdown,'b','Lipid')
#
# title('Prochlorococcus',y=1.02)
# ax1.set_ylim(18,80) #or (18,70)
# ax2.set_ylim(4,17.23) #or (4,15.1)
# ax1.set_ylabel('C allocation ($\%$)')
# ax2.set_ylabel('C allocation ($\%$): Lipid')
#
# ax2.legend(loc=2,fontsize=23)
#
# Savefig2(SaveFolder(),'DataPlot',Dpi())
show()