Ldir = Lfun.Lstart()

pth = os.path.abspath(Ldir['LO'] + 'plotting')
if pth not in sys.path:
    sys.path.append(pth)
import pfun

import numpy as np
import matplotlib.pyplot as plt
import pickle
from datetime import datetime, timedelta
import pandas as pd

import tef_fun
# get the DataFrame of all sections
sect_df = tef_fun.get_sect_df()

from warnings import filterwarnings
filterwarnings('ignore') # skip some warning messages

# choose input and organize output
Ldir = Lfun.Lstart()
indir0 = Ldir['LOo'] + 'tef/'
# choose the tef extraction to process
item = Lfun.choose_item(indir0)
indir0 = indir0 + item + '/'
indir = indir0 + 'bulk/'

in_fn = indir + 'sog5.p'

bulk = pickle.load(open(in_fn, 'rb'))
sys.path.append(os.path.abspath('../alpha'))
import Lfun
import zfun
Ldir = Lfun.Lstart()

sys.path.append(os.path.abspath(Ldir['LO'] + 'plotting'))
import pfun

import numpy as np
import matplotlib.pyplot as plt
import pickle
from datetime import datetime, timedelta

import tef_fun
# get the DataFrame of all sections
sect_df = tef_fun.get_sect_df()

from warnings import filterwarnings
filterwarnings('ignore')  # skip some warning messages

year = 2017
year_str = str(year)
item = 'cas6_v3_lo8b_2017.01.01_2017.12.31'

# choose input and organize output
Ldir = Lfun.Lstart()
indir00 = Ldir['LOo'] + 'tef/'
indir0 = indir00 + item + '/'
indir = indir0 + 'bulk/'

outdir = indir00 + 'bulk_plots_simple/'