#!/usr/bin/env python
# Wenchang Yang (wenchang@princeton.edu)
# Tue Sep 15 12:39:19 EDT 2020
if __name__ == '__main__':
    from misc.timer import Timer
    tt = Timer(f'start {__file__}')
import sys, os.path, os, glob
import xarray as xr, numpy as np, pandas as pd
#import matplotlib.pyplot as plt
#more imports
#
if __name__ == '__main__':
    tt.check('end import')
#
#start from here
ds_ctl = xr.open_dataset('../HiRAM_CTL1990/fig_scatter_ntc_vs_all.nc')
ds = xr.open_dataset('fig_scatter_ntc_vs_all.nc')
 
if __name__ == '__main__':
    from wyconfig import * #my plot settings
    figname = __file__.replace('.py', f'_{tt.today()}.png')
    ax = ds_ctl.ntc.to_dataframe(name='ntc_ctl').plot.bar(align='center')
    ds.ntc.to_dataframe(name='ntc_plus2k_2xco2').plot.bar(ax=ax, color='C1', align='edge')

    ax.set_ylabel('# per month')
    ax.set_title('HiRAM TC frequency annual cycle')

    plt.savefig(figname)
    print('[saved]:', figname)

    tt.check(f'**Done**')
    plt.show()