Skip to content
Snippets Groups Projects
Commit 4d27a624 authored by Johannes Hörner's avatar Johannes Hörner
Browse files

Reworked scripts, removed unnecessary clutter code

parent 3004de6c
Branches
Tags
No related merge requests found
Expand all Hide whitespace changes
Inline Side-by-side
Showing
with 9327 additions and 10327 deletions
python_scripts /bifurcation_Semtner_unlim-lim5.ipynb
+ 6
− 139
View file @ 4d27a624
%% Cell type:markdown id: tags:
# Create CO2 - ice-edge latitude diagrams and calculate stability
Stability is assessed by reading simulation data and applying the stability criterion. CO2 - ice-edge latitude diagrams are created by manual input of stability and equilibrium ice-edge latitude. The actual bifurcation diagram is added afterwards as a best guess.
%% Cell type:code id: tags:
``` python
%matplotlib inline
import xarray as xr
import numpy as np
import matplotlib.pyplot as plt
from scipy import stats
def load_experiment(expname): #loads the dataset of a simulation
fname = expname +"_atm_2d_ml.mm.gm.nc" #filename of global yearly mean
dpath = "/work/bb1092/pp_JH/" +expname +"/" #simulation path
DS = xr.open_dataset(dpath +fname, decode_times=False) #loading of dataset
print(dpath +fname)
return expname, DS # returns the name of the experiment & the actual dataset
def get_var(dataset, varname, offsettime): #gets the dataarray with one specific variable
da=getattr(dataset,varname) #read dataarray
da.squeeze() #squeeze dataarray (time is the only dimension)
if offsettime:
da=da.assign_coords(time=((da.time-da.time[0])/360)) #change time units from days to years & move the origin to 0
else:
da=da.assign_coords(time=((da.time)/360)) #change time units from days to years
return da
def plot_simulation(axes,co2, startlat, endlat, col, stable): #plot a simulation into the bifurcation diagram
handle, =axes.plot([co2,co2],[np.sin(np.radians(startlat)),np.sin(np.radians(endlat))],color=col,linestyle='--') #plot the line
#plt.plot(co2,np.sin(np.radians(startlat)),'bo',fillstyle='none')#
if stable:
if stable==2:
axes.plot(co2,np.sin(np.radians(endlat)),markeredgecolor=col,marker='o',markerfacecolor="none")
else:
axes.plot(co2,np.sin(np.radians(endlat)),color=col,marker='o')
else:
if endlat<startlat:
axes.plot(co2,np.sin(np.radians(endlat)),color=col,marker='v')
else:
axes.plot(co2,np.sin(np.radians(endlat)),color=col,marker='^')
return handle
def plot_simulation_merged(axes,co2, startlat, endlat, col, stable, offset=0): #plot a simulation into the bifurcation diagram
#plt.plot(co2,np.sin(np.radians(startlat)),'bo',fillstyle='none')#
handle = []
if stable==2: # metastable
axes.plot(co2 *(1+offset),np.sin(np.radians(endlat)),markeredgecolor=col,marker='o',markerfacecolor="none",clip_on=False)
elif stable==1: # stable from warm
axes.plot(co2 *(1+offset),np.sin(np.radians(endlat)),color=col,marker='o',clip_on=False)
elif stable==3: # towards Snowball
if endlat==0:
axes.plot(co2 *(1+offset),np.sin(np.radians(endlat)),color=col,marker='v',clip_on=False)
else:
axes.plot(co2 *(1+offset),np.sin(np.radians(endlat)),color=col,markerfacecolor='none',marker='v',clip_on=False)
handle, =axes.plot([co2 *(1+offset),co2 *(1+offset)],[np.sin(np.radians(startlat)),np.sin(np.radians(endlat))],color=col,linestyle='-',alpha=0.2,clip_on=False) #plot the line
elif stable==4: # towards icefree
axes.plot(co2 *(1+offset),np.sin(np.radians(endlat)),color=col,markerfacecolor='none',marker='^',clip_on=False)
handle, =axes.plot([co2 *(1+offset),co2 *(1+offset)],[np.sin(np.radians(startlat)),np.sin(np.radians(endlat))],color=col,linestyle='-',alpha=0.2,clip_on=False) #plot the line
return handle
def plot_co2_equivalent(co2, startlat, endlat, col, stable, heatflux):
CO2_fac=np.exp(heatflux/5.35)
handle, =plt.plot([co2*CO2_fac,co2*CO2_fac],[np.sin(np.radians(startlat)),np.sin(np.radians(endlat))],color=col,linestyle='--') #plot the line
#plt.plot(co2,np.sin(np.radians(startlat)),'bo',fillstyle='none')#
if stable:
plt.plot(co2*CO2_fac,np.sin(np.radians(endlat)),color=col,marker='o')
else:
if endlat<startlat:
plt.plot(co2*CO2_fac,np.sin(np.radians(endlat)),color=col,marker='v')
else:
plt.plot(co2*CO2_fac,np.sin(np.radians(endlat)),color=col,marker='^')
return handle
def legend_color(ax, handle_array, pos):
legend = ax.legend(handle_array,handlelength=0, handletextpad=0, edgecolor='none', facecolor='none', markerscale=0, loc=pos)
for item in legend.legendHandles:
item.set_visible(False)
for text in legend.get_texts():
if text.get_text()=='Winton':
text.set_color('C1')
if text.get_text()=='3L-Winton':
text.set_color('C1')
elif text.get_text()=='Semtner':
text.set_color('C0')
elif text.get_text()=='0L-Semtner':
text.set_color('C0')
elif text.get_text()=='Semtner_5m':
text.set_color('C2')
elif text.get_text()=='0L-Semtner-lim5':
text.set_color('C2')
elif text.get_text()=='0L-Semtner-limited':
text.set_color('C2')
elif text.get_text()=='1438ppmv':
text.set_color('C0')
elif text.get_text()=='1500ppmv':
text.set_color('C1')
elif text.get_text()=='3000ppmv':
text.set_color('C2')
elif text.get_text()=='5000ppmv':
text.set_color('C3')
return legend
def stability_index(DS):
LW = get_var(DS,"rlut",True)
SW = get_var(DS,"rsdt",True)-get_var(DS,"rsut",True)
T = get_var(DS,'ts',True)
N = np.squeeze(SW[-50*12:]-LW[-50*12:])
Nparm = stats.linregress(LW.time[-50*12:]-LW.time[-50*12],N)
f=(Nparm.slope) #*np.sign(np.squeeze(SW-LW).values[-1])
return f
def is_stable(f):
if f>0: return False
else: return True
```
%% Cell type:code id: tags:
``` python
experiment1, DS1 = load_experiment("mlo_aqua_1500ppmv_hice_unlim")
experiment2, DS2 = load_experiment("mlo_aqua_1594ppmv_hice_unlim")
experiment3, DS3 = load_experiment("mlo_aqua_1688ppmv_hice_unlim")
experiment4, DS4 = load_experiment("mlo_aqua_1875ppmv_hice_unlim")
experiment5, DS5 = load_experiment("mlo_aqua_2250ppmv_hice_unlim")
experiment1_Jor, DS1_Jor = load_experiment("mlo_aqua_3000ppmv_77sic_hice_unlim")
experiment2_Jor, DS2_Jor = load_experiment("mlo_aqua_3000ppmv_61sic_hice_unlim")
experiment3_Jor, DS3_Jor = load_experiment("mlo_aqua_3750ppmv_69sic_hice_unlim")
experiment4_Jor, DS4_Jor = load_experiment("mlo_aqua_3750ppmv_77sic_hice_unlim")
experiment5_Jor, DS5_Jor = load_experiment("mlo_aqua_3907ppmv_77sic_hice_unlim")
experiment6_Jor, DS6_Jor = load_experiment("mlo_aqua_4063ppmv_77sic_hice_unlim")
experiment7_Jor, DS7_Jor = load_experiment("mlo_aqua_4219ppmv_77sic_hice_unlim")
experiment8_Jor, DS8_Jor = load_experiment("mlo_aqua_4375ppmv_77sic_hice_unlim")
experiment9_Jor, DS9_Jor = load_experiment("mlo_aqua_5000ppmv_77sic_hice_unlim")
experiment10_Jor, DS10_Jor = load_experiment("mlo_aqua_5000ppmv_83sic_hice_unlim")
experiment11_Jor, DS11_Jor = load_experiment("mlo_aqua_10000ppmv_77sic_hice_unlim")
experiment1_5m, DS1_5m = load_experiment("mlo_aqua_1500ppmv")
experiment2_5m, DS2_5m = load_experiment("mlo_aqua_3000ppmv")
experiment3_5m, DS3_5m = load_experiment("mlo_aqua_5000ppmv")
DS_array=[DS1, DS2, DS3, DS4, DS5] #, DS1_Jor, DS2_Jor, DS3_Jor, DS4_Jor, DS5_Jor, DS6_Jor
# , DS7_Jor, DS8_Jor, DS9_Jor, DS10_Jor, DS11_Jor
```
%% Output
/work/bb1092/pp_JH/mlo_aqua_1500ppmv_hice_unlim/mlo_aqua_1500ppmv_hice_unlim_atm_2d_ml.mm.gm.nc
/work/bb1092/pp_JH/mlo_aqua_1594ppmv_hice_unlim/mlo_aqua_1594ppmv_hice_unlim_atm_2d_ml.mm.gm.nc
/work/bb1092/pp_JH/mlo_aqua_1688ppmv_hice_unlim/mlo_aqua_1688ppmv_hice_unlim_atm_2d_ml.mm.gm.nc
/work/bb1092/pp_JH/mlo_aqua_1875ppmv_hice_unlim/mlo_aqua_1875ppmv_hice_unlim_atm_2d_ml.mm.gm.nc
/work/bb1092/pp_JH/mlo_aqua_2250ppmv_hice_unlim/mlo_aqua_2250ppmv_hice_unlim_atm_2d_ml.mm.gm.nc
/work/bb1092/pp_JH/mlo_aqua_3000ppmv_77sic_hice_unlim/mlo_aqua_3000ppmv_77sic_hice_unlim_atm_2d_ml.mm.gm.nc
/work/bb1092/pp_JH/mlo_aqua_3000ppmv_61sic_hice_unlim/mlo_aqua_3000ppmv_61sic_hice_unlim_atm_2d_ml.mm.gm.nc
/work/bb1092/pp_JH/mlo_aqua_3750ppmv_69sic_hice_unlim/mlo_aqua_3750ppmv_69sic_hice_unlim_atm_2d_ml.mm.gm.nc
/work/bb1092/pp_JH/mlo_aqua_3750ppmv_77sic_hice_unlim/mlo_aqua_3750ppmv_77sic_hice_unlim_atm_2d_ml.mm.gm.nc
/work/bb1092/pp_JH/mlo_aqua_3907ppmv_77sic_hice_unlim/mlo_aqua_3907ppmv_77sic_hice_unlim_atm_2d_ml.mm.gm.nc
/work/bb1092/pp_JH/mlo_aqua_4063ppmv_77sic_hice_unlim/mlo_aqua_4063ppmv_77sic_hice_unlim_atm_2d_ml.mm.gm.nc
/work/bb1092/pp_JH/mlo_aqua_4219ppmv_77sic_hice_unlim/mlo_aqua_4219ppmv_77sic_hice_unlim_atm_2d_ml.mm.gm.nc
/work/bb1092/pp_JH/mlo_aqua_4375ppmv_77sic_hice_unlim/mlo_aqua_4375ppmv_77sic_hice_unlim_atm_2d_ml.mm.gm.nc
/work/bb1092/pp_JH/mlo_aqua_5000ppmv_77sic_hice_unlim/mlo_aqua_5000ppmv_77sic_hice_unlim_atm_2d_ml.mm.gm.nc
/work/bb1092/pp_JH/mlo_aqua_5000ppmv_83sic_hice_unlim/mlo_aqua_5000ppmv_83sic_hice_unlim_atm_2d_ml.mm.gm.nc
/work/bb1092/pp_JH/mlo_aqua_10000ppmv_77sic_hice_unlim/mlo_aqua_10000ppmv_77sic_hice_unlim_atm_2d_ml.mm.gm.nc
/work/bb1092/pp_JH/mlo_aqua_1500ppmv/mlo_aqua_1500ppmv_atm_2d_ml.mm.gm.nc
/work/bb1092/pp_JH/mlo_aqua_3000ppmv/mlo_aqua_3000ppmv_atm_2d_ml.mm.gm.nc
/work/bb1092/pp_JH/mlo_aqua_5000ppmv/mlo_aqua_5000ppmv_atm_2d_ml.mm.gm.nc
%% Cell type:code id: tags:
``` python
### %matplotlib inline
import numpy as np
import matplotlib.pyplot as plt
save_plot = True
plot_unlim = True
plot_5m = True
offset=0.01
labelsize=10
ticksize=10
fig, ax = plt.subplots(figsize=(6,4))
ax.spines['left'].set_position(('outward',5))
ax.spines['bottom'].set_position(('outward',5))
# 5m simulations
if plot_5m:
color_5m='C2'
#b=plot_simulation_merged(ax,1250, 90, 0, color_5m, 1)
#plot_simulation_merged(ax,1375, 90, 0, color_5m, 1)
plot_simulation_merged(ax,1438, 90, 0, color_5m, 1)
plot_simulation_merged(ax,1500, 90, 46.6, color_5m, 1)
plot_simulation_merged(ax,1594, 90, 50.3, color_5m, 1)
#plot_simulation_merged(ax,1250, 90, 0, color_5m, 1)
plot_simulation_merged(ax,3000, 90, 65, color_5m, 1)
plot_simulation_merged(ax,5000, 90, 90, color_5m, 1)
# Jormungand
#mlo_aqua_2000ppmv_Jor2
#plot_simulation_merged(ax,2000, 32.68, 29.34, color_5m, 3)
#mlo_aqua_2000ppmv_Jor5
plot_simulation_merged(ax,2000, 13.48, 0, color_5m, 3)
#mlo_aqua_2000ppmv_Jor5
#plot_simulation_merged(ax,2000, 32.68, 29.34, color_5m, 3)
#mlo_aqua_2500ppmv_Jor2
plot_simulation_merged(ax,2500, 13.48, 0, color_5m, 3)
#mlo_aqua_2500ppmv_Jor5
#plot_simulation_merged(ax,2500, 32.68, 36.16, color_5m, 4)
#mlo_aqua_3000ppmv_Jor
plot_simulation_merged(ax,3000, 28.03, 36.87, color_5m, 4)
#mlo_aqua_3000ppmv_Jor2
plot_simulation_merged(ax,3000, 13.48, 0, color_5m, 3)
#mlo_aqua_3000ppmv_Jor3
#plot_simulation_merged(ax,3000, 16.26, 22.95, color_5m, 4)
#mlo_aqua_5000ppmv_Jor
#plot_simulation_merged(ax,5000, 28.03, 34.75, color_5m, 4)
#mlo_aqua_5000ppmv_Jor2
plot_simulation_merged(ax,5000, 13.29, 30.00, color_5m, 4)
#mlo_aqua_5000ppmv_Jor4
plot_simulation_merged(ax,5000, 7.47, 15.07, color_5m, 4, -offset)
#mlo_aqua_10000ppmv_Jor
#plot_simulation_merged(ax,10000, 28.03, 41.30, color_5m, 4)
#hice_unlim simulations
if plot_unlim:
color_unlim='C0'
offset=0.01
a=plot_simulation_merged(ax,1500, 90, 0, color_unlim, 1)
plot_simulation_merged(ax,1594, 90, 49.54, color_unlim, 1, offset)
plot_simulation_merged(ax,1688, 90, 51.07, color_unlim, 1)
plot_simulation_merged(ax,1875, 90, 53.20, color_unlim, 1)
plot_simulation_merged(ax,2250, 90, 56.85, color_unlim, 1)
plot_simulation_merged(ax,3000, 90, 65, color_unlim, 1, offset) #5m limited, but always under 5m
plot_simulation_merged(ax,3000, 14.46, 0, color_unlim, 3, offset)
#plot_simulation_merged(ax,3000, 22.05, 22.83, color_unlim, 4)
plot_simulation_merged(ax,3750, 14.46, 0, color_unlim, 3)
#plot_simulation_merged(ax,3750, 19.2, 18.7, color_unlim, 4)
#plot_simulation_merged(ax,3907, 14.46, 15.1, color_unlim, 2)
plot_simulation_merged(ax,4063, 14.46, 11.5, color_unlim, 3)
#plot_simulation_merged(ax,4219, 14.46, 16.23, color_unlim, 2)
plot_simulation_merged(ax,4375, 14.46, 18.2, color_unlim, 4)
plot_simulation_merged(ax,5000, 14.46, 24, color_unlim, 4, offset)
#plot_simulation_merged(ax,5000, 11.6, 14.3, color_unlim, 2)
plot_simulation_merged(ax,5000, 90, 90, color_unlim, 1, offset) # limited, but no ice at all
#plot_simulation_merged(ax,10000, 14.46, 18, color_unlim, 4)
ax.set_yticks([np.sin(np.radians(0)),np.sin(np.radians(10)),np.sin(np.radians(20)),np.sin(np.radians(30)),np.sin(np.radians(45)),np.sin(np.radians(60)),np.sin(np.radians(90))])
ax.set_ylim(-0.1,1.1)
ax.set_xlabel("Atmospheric CO$_2$ [ppmv]",size=labelsize)
ax.set_ylabel("Ice-Edge Latitude [°]",size=labelsize)
plt.tick_params(axis='both', which='major', labelsize=ticksize)
ax.spines['right'].set_color('none')
ax.spines['top'].set_color('none')
ax.yaxis.tick_left()
ax.xaxis.tick_bottom()
ax.set_xscale('log')
ax.set_xticks(np.linspace(1000,10000,10) )
ax.set_xticklabels([1000, 2000,3000,4000,5000,'','','','',10000] )
ax.set_xlim(1000,5000)
ax.set_ylim(0,1)
# extra ticks at bifurcation points
#semtner_bif=1641
#semtner_5m_bif=1547
semtner_bif=1547
semtner_5m_bif=1469
ax.vlines([semtner_5m_bif],-0.04,-0.01,color=color_5m,lw=1,clip_on=False)
ax.vlines([semtner_bif],-0.04,-0.01,color=color_unlim,lw=1,clip_on=False)
#ax.axvline(semtner_bif,-20,0.1)
ax.annotate(semtner_bif,(semtner_bif,0.02),color=color_unlim,ha='center',clip_on=False)
ax.annotate(semtner_5m_bif,(semtner_5m_bif,0.1),xytext=(semtner_5m_bif,-0.05),color=color_5m,ha='center',va='top',clip_on=False)
ax.set_yticklabels([0,10,20,30,45,60,90])
legend_color(ax,["0L-Semtner","0L-Semtner-limited"],2)
if save_plot:
plt.savefig("plots/bifurcation_5mlim_merged.pdf")
```
%% Output
%% Cell type:code id: tags:
``` python
```
%% Cell type:code id: tags:
``` python
```
%% Cell type:code id: tags:
``` python
```
......
python_scripts /bifurcation_Winton_Semtner_unlim.ipynb
+ 4
− 564
View file @ 4d27a624
This diff is collapsed.
python_scripts /energyflux_icelimit_comparison.ipynb
+ 57
− 332
View file @ 4d27a624
Source diff could not be displayed: it is too large. Options to address this: view the blob.
python_scripts /overview.ipynb
+ 187
− 227
View file @ 4d27a624
Source diff could not be displayed: it is too large. Options to address this: view the blob.
python_scripts /plots/bifurcation_5mlim_merged.pdf
+ 0
− 0
View file @ 4d27a624
No preview for this file type
python_scripts /plots/bifurcation_winton_merged.pdf
+ 0
− 0
View file @ 4d27a624
No preview for this file type
python_scripts /plots/heatflux_comparison.svg
+ 25
− 25
View file @ 4d27a624
......@@ -32,10 +32,10 @@ z
<defs>
<path d="M 0 0
L 0 3.5
" id="me7473f8e97" style="stroke:#000000;stroke-width:0.8;"/>
" id="mb8e67a1687" style="stroke:#000000;stroke-width:0.8;"/>
</defs>
<g>
<use style="stroke:#000000;stroke-width:0.8;" x="54" xlink:href="#me7473f8e97" y="252"/>
<use style="stroke:#000000;stroke-width:0.8;" x="54" xlink:href="#mb8e67a1687" y="252"/>
</g>
</g>
<g id="text_1">
......@@ -71,7 +71,7 @@ z
<g id="xtick_2">
<g id="line2d_2">
<g>
<use style="stroke:#000000;stroke-width:0.8;" x="93.857143" xlink:href="#me7473f8e97" y="252"/>
<use style="stroke:#000000;stroke-width:0.8;" x="93.857143" xlink:href="#mb8e67a1687" y="252"/>
</g>
</g>
<g id="text_2">
......@@ -135,7 +135,7 @@ z
<g id="xtick_3">
<g id="line2d_3">
<g>
<use style="stroke:#000000;stroke-width:0.8;" x="133.714286" xlink:href="#me7473f8e97" y="252"/>
<use style="stroke:#000000;stroke-width:0.8;" x="133.714286" xlink:href="#mb8e67a1687" y="252"/>
</g>
</g>
<g id="text_3">
......@@ -149,7 +149,7 @@ z
<g id="xtick_4">
<g id="line2d_4">
<g>
<use style="stroke:#000000;stroke-width:0.8;" x="173.571429" xlink:href="#me7473f8e97" y="252"/>
<use style="stroke:#000000;stroke-width:0.8;" x="173.571429" xlink:href="#mb8e67a1687" y="252"/>
</g>
</g>
<g id="text_4">
......@@ -174,7 +174,7 @@ z
<g id="xtick_5">
<g id="line2d_5">
<g>
<use style="stroke:#000000;stroke-width:0.8;" x="213.428571" xlink:href="#me7473f8e97" y="252"/>
<use style="stroke:#000000;stroke-width:0.8;" x="213.428571" xlink:href="#mb8e67a1687" y="252"/>
</g>
</g>
<g id="text_5">
......@@ -204,7 +204,7 @@ z
<g id="xtick_6">
<g id="line2d_6">
<g>
<use style="stroke:#000000;stroke-width:0.8;" x="253.285714" xlink:href="#me7473f8e97" y="252"/>
<use style="stroke:#000000;stroke-width:0.8;" x="253.285714" xlink:href="#mb8e67a1687" y="252"/>
</g>
</g>
<g id="text_6">
......@@ -219,7 +219,7 @@ z
<g id="xtick_7">
<g id="line2d_7">
<g>
<use style="stroke:#000000;stroke-width:0.8;" x="293.142857" xlink:href="#me7473f8e97" y="252"/>
<use style="stroke:#000000;stroke-width:0.8;" x="293.142857" xlink:href="#mb8e67a1687" y="252"/>
</g>
</g>
<g id="text_7">
......@@ -234,7 +234,7 @@ z
<g id="xtick_8">
<g id="line2d_8">
<g>
<use style="stroke:#000000;stroke-width:0.8;" x="333" xlink:href="#me7473f8e97" y="252"/>
<use style="stroke:#000000;stroke-width:0.8;" x="333" xlink:href="#mb8e67a1687" y="252"/>
</g>
</g>
<g id="text_8">
......@@ -249,7 +249,7 @@ z
<g id="xtick_9">
<g id="line2d_9">
<g>
<use style="stroke:#000000;stroke-width:0.8;" x="372.857143" xlink:href="#me7473f8e97" y="252"/>
<use style="stroke:#000000;stroke-width:0.8;" x="372.857143" xlink:href="#mb8e67a1687" y="252"/>
</g>
</g>
<g id="text_9">
......@@ -488,10 +488,10 @@ z
<defs>
<path d="M 0 0
L -3.5 0
" id="m313bc820b0" style="stroke:#000000;stroke-width:0.8;"/>
" id="ma0fcf8f027" style="stroke:#000000;stroke-width:0.8;"/>
</defs>
<g>
<use style="stroke:#000000;stroke-width:0.8;" x="54" xlink:href="#m313bc820b0" y="252"/>
<use style="stroke:#000000;stroke-width:0.8;" x="54" xlink:href="#ma0fcf8f027" y="252"/>
</g>
</g>
<g id="text_11">
......@@ -504,7 +504,7 @@ L -3.5 0
<g id="ytick_2">
<g id="line2d_11">
<g>
<use style="stroke:#000000;stroke-width:0.8;" x="54" xlink:href="#m313bc820b0" y="224.82"/>
<use style="stroke:#000000;stroke-width:0.8;" x="54" xlink:href="#ma0fcf8f027" y="224.82"/>
</g>
</g>
<g id="text_12">
......@@ -517,7 +517,7 @@ L -3.5 0
<g id="ytick_3">
<g id="line2d_12">
<g>
<use style="stroke:#000000;stroke-width:0.8;" x="54" xlink:href="#m313bc820b0" y="197.64"/>
<use style="stroke:#000000;stroke-width:0.8;" x="54" xlink:href="#ma0fcf8f027" y="197.64"/>
</g>
</g>
<g id="text_13">
......@@ -530,7 +530,7 @@ L -3.5 0
<g id="ytick_4">
<g id="line2d_13">
<g>
<use style="stroke:#000000;stroke-width:0.8;" x="54" xlink:href="#m313bc820b0" y="170.46"/>
<use style="stroke:#000000;stroke-width:0.8;" x="54" xlink:href="#ma0fcf8f027" y="170.46"/>
</g>
</g>
<g id="text_14">
......@@ -577,7 +577,7 @@ z
<g id="ytick_5">
<g id="line2d_14">
<g>
<use style="stroke:#000000;stroke-width:0.8;" x="54" xlink:href="#m313bc820b0" y="143.28"/>
<use style="stroke:#000000;stroke-width:0.8;" x="54" xlink:href="#ma0fcf8f027" y="143.28"/>
</g>
</g>
<g id="text_15">
......@@ -609,7 +609,7 @@ z
<g id="ytick_6">
<g id="line2d_15">
<g>
<use style="stroke:#000000;stroke-width:0.8;" x="54" xlink:href="#m313bc820b0" y="116.1"/>
<use style="stroke:#000000;stroke-width:0.8;" x="54" xlink:href="#ma0fcf8f027" y="116.1"/>
</g>
</g>
<g id="text_16">
......@@ -622,7 +622,7 @@ z
<g id="ytick_7">
<g id="line2d_16">
<g>
<use style="stroke:#000000;stroke-width:0.8;" x="54" xlink:href="#m313bc820b0" y="88.92"/>
<use style="stroke:#000000;stroke-width:0.8;" x="54" xlink:href="#ma0fcf8f027" y="88.92"/>
</g>
</g>
<g id="text_17">
......@@ -667,7 +667,7 @@ z
<g id="ytick_8">
<g id="line2d_17">
<g>
<use style="stroke:#000000;stroke-width:0.8;" x="54" xlink:href="#m313bc820b0" y="61.74"/>
<use style="stroke:#000000;stroke-width:0.8;" x="54" xlink:href="#ma0fcf8f027" y="61.74"/>
</g>
</g>
<g id="text_18">
......@@ -680,7 +680,7 @@ z
<g id="ytick_9">
<g id="line2d_18">
<g>
<use style="stroke:#000000;stroke-width:0.8;" x="54" xlink:href="#m313bc820b0" y="34.56"/>
<use style="stroke:#000000;stroke-width:0.8;" x="54" xlink:href="#ma0fcf8f027" y="34.56"/>
</g>
</g>
<g id="text_19">
......@@ -912,7 +912,7 @@ z
</g>
</g>
<g id="line2d_19">
<path clip-path="url(#pb047c6e5ce)" d="M 54 238.064221
<path clip-path="url(#p2b19aa5384)" d="M 54 238.064221
L 55.594286 236.749226
L 57.188571 236.521669
L 58.782857 236.072591
......@@ -1095,7 +1095,7 @@ L 372.857143 112.438198
" style="fill:none;stroke:#1f77b4;stroke-linecap:square;stroke-width:1.5;"/>
</g>
<g id="line2d_20">
<path clip-path="url(#pb047c6e5ce)" d="M 54 252
<path clip-path="url(#p2b19aa5384)" d="M 54 252
L 79.508571 252
L 81.102857 251.356772
L 82.697143 252
......@@ -1306,7 +1306,7 @@ L 433 238.859914
" style="fill:none;stroke:#ff7f0e;stroke-linecap:square;stroke-width:1.5;"/>
</g>
<g id="line2d_21">
<path clip-path="url(#pb047c6e5ce)" d="M 54 143.318592
<path clip-path="url(#p2b19aa5384)" d="M 54 143.318592
L 55.594286 144.163657
L 57.188571 144.232215
L 58.782857 143.89198
......@@ -1455,7 +1455,7 @@ L 312.274286 92.241839
" style="fill:none;stroke:#2ca02c;stroke-linecap:square;stroke-width:1.5;"/>
</g>
<g id="line2d_22">
<path clip-path="url(#pb047c6e5ce)" d="M 54 144.8917
<path clip-path="url(#p2b19aa5384)" d="M 54 144.8917
L 55.594286 147.211242
L 57.188571 148.655568
L 58.782857 148.855479
......@@ -1675,7 +1675,7 @@ L 83 91.970938
</g>
</g>
<defs>
<clipPath id="pb047c6e5ce">
<clipPath id="p2b19aa5384">
<rect height="217.44" width="334.8" x="54" y="34.56"/>
</clipPath>
</defs>
......
python_scripts /plots/heatflux_scatter_flux_icelat.pdf
+ 0
− 0
View file @ 4d27a624
No preview for this file type
python_scripts /plots/heatflux_scatter_flux_sit.svg
+ 9034
− 9034
View file @ 4d27a624
This diff is collapsed.
python_scripts /plots/icelat_comparison.pdf
+ 0
− 0
View file @ 4d27a624
No preview for this file type
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment