#!/usr/bin/env python
# -*- coding: utf-8 -*-

import numpy as np
import matplotlib.pyplot as plt
from scipy import stats
from matplotlib import ticker
from mpl_toolkits.axes_grid1 import AxesGrid
from matplotlib.backends.backend_pdf import PdfPages

nums = 2

dd1 = np.linspace(5,75,25)
dd2 = np.linspace(100,900,50)

conc1 = np.random.randn(nums*16,25) * 30
conc2 = np.random.randn(nums*16,50) * 10

# Creating PdfPages object
pdf = PdfPages('test_speed2.pdf') 

xx = 0 
for i in range(nums):
    fig = plt.figure(figsize=(11, 8.5))
    grid = AxesGrid(fig, 111, nrows_ncols=(4,4), axes_pad = 0.0, share_all=0, aspect=0)
    for j in range(xx*16, xx*16+16):
        grid[j%16].plot(dd1, conc1[j], color='r', lw=1.5)
        grid[j%16].plot(dd2, conc2[j], color='b', lw=1.5)
        grid[j%16].set_xlim(1,2000)
        grid[j%16].set_ylim(-50,50)
        grid[j%16].set_xscale('log')
        grid[j%16].set_xticks([10, 100, 1000])
        grid[j%16].set_yticks([-25, 0, 25])
        grid[j%16].yaxis.set_minor_locator(ticker.MultipleLocator(5))
    fig.text(0.5, 0.01, u"Xlabel", fontsize=16, transform=fig.transFigure)
    fig.text(0.01, 0.6, u"Ylabel", fontsize=16, rotation='90', transform=fig.transFigure)
    fig.subplots_adjust(left=0.06, right=0.98, wspace=0.2, bottom=0.07, top=0.98)
    pdf.savefig(fig)
    plt.close(fig)
    xx = xx + 1

pdf.close()

#plt.show()