<jupyter_start><jupyter_text>Plot Figure 3, panel a)<jupyter_code>pickle.HIGHEST_PROTOCOL
with open('../data/data_Fig3a.pickle', 'rb') as f:
    # The protocol version used is detected automatically, so we do not
    # have to specify it.
    data_Figure3a = pickle.load(f)
#-- plot everything together
fig=plt.figure(figsize=(13,7))
ax = fig.add_subplot(111)

#---- HR ---------------------------------
plt.plot(data_Figure3a['HR']['spinup-time'], data_Figure3a['HR']['spinup-data'], color='0.8',linewidth=3.) # label='HR spinup', '0.6'

# plot HR control
#
HR_contro_time=data_Figure3a['HR']['control-time']
HR_contro_data=data_Figure3a['HR']['control-data']
plt.plot(HR_contro_time, HR_contro_data, color='0.8',linewidth=3.,label='HR 1950-control') # original data

# plot HR scenario 
#
plt.plot(data_Figure3a['HR']['scenario-time'], data_Figure3a['HR']['scenario-data'], color='white',linewidth=6.,zorder=39) # original data
plt.plot(data_Figure3a['HR']['scenario-time'], data_Figure3a['HR']['scenario-data'], color='black',linewidth=3.,label='HR hist+scenario',zorder=40) # original data

plt.plot(np.zeros(1), np.zeros([1,1]), color='w', alpha=0, label=' '.ljust(28)) # dummy entry

#----- LR ---------------------------------
plt.plot(data_Figure3a['LR']['spinup-time'], data_Figure3a['LR']['spinup-data'], color='paleturquoise',linewidth=3.) # label='LR spinup', 'darkturquoise'

# plot LR control
plt.plot(data_Figure3a['LR']['control-time'], data_Figure3a['LR']['control-data'], color='paleturquoise',linewidth=3.,label='LR 1950-control') # original data
# plot LR scenario
plt.plot(data_Figure3a['LR']['scenario-time'], data_Figure3a['LR']['scenario-data'], color='darkcyan',linewidth=3.,label='LR hist+scenario')

#---- NSIDC ('50' is 1950, '79' is 1979) --
plt.plot(data_Figure3a['OBS']['time'], data_Figure3a['OBS']['data'], color='white', linewidth=4.,zorder=59)
plt.plot(data_Figure3a['OBS']['time'], data_Figure3a['OBS']['data'], color='red', linewidth=2., label='OBS (1979-2018)', zorder=60)


plt.ylabel('Extent (million square kilometres)',fontsize=18)
plt.xlabel('Year',fontsize=18)

plt.xticks(np.arange(50,200,10),np.arange(1950,2098,10),fontsize=18,rotation=45)
plt.axis([0.0,199,10.,25.]) # [xmin, xmax, ymin, ymax]
plt.yticks([10.0,12.0,14.0,16.0,18.0,20.,22.,24.], fontsize = 18)

plt.title('Antarctic Sea Ice Extent (September 1950-2099)',fontsize=22)
plt.legend(loc='lower left', ncol=2, fontsize=18)
plt.show()
saveFIG=True
if saveFIG==True:
   pp = PdfPages('./Figure3a.pdf')
   pp.savefig(fig,bbox_inches = 'tight')
   pp.close()<jupyter_output><empty_output>
{'fraction_non_alphanumeric': 0.2337962962962963, 'fraction_numerical': 0.059799382716049385, 'mean_word_length': 6.008108108108108, 'pattern_counts': {'":': 0, '<': 5, '<?xml version=': 0, '>': 5, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
{'chain_length': '1', 'dir': 'jupyter-structured-clean-dedup', 'id': '125147', 'license': 'MIT', 'n_tokens_mistral': 1073, 'n_tokens_neox': 961, 'n_words': 147, 'path': 'notebooks/Figure3a.ipynb', 'repo_name': 'trackow/AntarcticSeaIce_NatCommun', 'size': '2592'}