import os

for fieldIdx, field in enumerate(["SZ", "SAnomalyZ", "TZ", "TAnomalyZ", "OHCZ", "OHCAnomalyZ", "OHC", "OHCAnomaly"]):

    for baselineStartYear in [1, 11, 21]:

        for baselineLength in [1, 2, 3, 4, 5]:

            for movingAveragePoints in [1, 12, 24, 36, 48, 60, 72, 84, 96, 108, 120, 132, 144, 156, 168, 180, 192, 204, 216, 228, 240, 252, 264, 276, 288, 300, 312, 324, 336, 348, 360, 372, 384, 396, 408, 420, 432, 444, 456, 468, 480]:

                os.system('ln -sf %s_%d_global_20170502.beta1_05.A_WCYCL1850S.ne30_oECv3_ICG.edison.png image_files/%d/%d/%d_%d_global_20170502.beta1_05.A_WCYCL1850S.ne30_oECv3_ICG.edison.png' % (field, movingAveragePoints, baselineStartYear, baselineLength, fieldIdx, movingAveragePoints))