import xarray as xr
from glob import glob
import os

gcms = [
    "access1.0_RCP45",
    "access1.0_RCP85",
    "access1.3_RCP85",
    "bcc-csm1.1_RCP85",
    "canesm2_RCP85",
    "ccsm4_RCP85",
    "csiro-mk3.6.0_RCP85",
    "fgoals-g2_RCP85",
    "gfdl-cm3_RCP85",
    "giss-e2-h_RCP85",
    "miroc5_RCP85",
    "mri-cgcm3_RCP85",
    "noresm1-m_RCP85"
]
gcms = ["pnnl"]

locs = ['eugene', 'omak', 'salem', 'vancouver']
idir = "/home/disk/rocinante/DATA/temp/crystal_fire/slp/data_daily/"
odir = "/home/disk/rocinante/DATA/temp/crystal_fire/slp/data_merge/"

for g in gcms:
    gidir = "{}/{}".format(idir, g)
    godir = "{}/{}".format(odir, g)
    os.makedirs(godir, exist_ok=True)
    for l in locs:
        print(g, l)
        #print('{}/*{}*.nc'.format(gidir, l))
        data = glob('{}/*{}*.nc'.format(gidir, l))
        ds = xr.open_mfdataset(data)
        df = ds.to_dataframe()
        df.to_csv('{}/{}_{}.csv'.format(godir, g, l))