#!/usr/bin/env python

# lock processing
import os
from sadcat.tools import lockproc
lockFileName = lockproc()

# headless matplotlib (for johansen)
import matplotlib
matplotlib.use('Agg')

from sadcat import *

# add downloaded MODIS images to catalog
cic = CommonImagesCatalog()
cic.fill_database('modis', '/Data/sat/downloads/MODIS/NorthSea/', 'A*.L2_LAC.NorthNorwegianSeas.hdf', ModisL2Image)
cic = None

# L2-processing of MODIS
#modis_std: standard processing - only image generation
habCatalogStd = L2Catalog('modis_std')
habCatalogStd.process(ModisL2Image)
habCatalogStd = None
# L3 processing: temporal and spatial binning, visualisation
l3pModisStd = L3PROC('modis_std_NorthSea', 'modis_std', ModisL2Image, 20)
l3pModisStd.l3('daily')
l3pModisStd.l3('weekly')

#modis_boreali: boreali processing - retrieval of chl, tsm, doc with BOREALI
#habCatalogBor = L2Catalog('modis_boreali')
#habCatalogBor.process(ModisL2Image)

# L3 processing: temporal and spatial binning, visualisation
#l3pModisBor = L3PROC('modis_boreali_NorthSea', 20)
#l3pModisBor.l3(habCatalogBor, 'daily')
#l3pModisBor.l3(habCatalogBor, 'weekly')

# unlock processing
os.remove(lockFileName)