#!/usr/bin/env python
# Wenchang Yang (wenchang@princeton.edu)
# Tue Mar 12 23:22:32 EDT 2019
import os, os.path, sys
from climdata import run_shell

data_name = 'PRECT'
#ifile_template = f'/tigress/wenchang/data/cesm/LME/{data_name}/b.e11.BLMTRC5CN.f19_g16.VOLC_GRA.iii.cam.h0.{data_name}.085001-184912.nc'
#ens = range(1,6)
ifile_template = f'/tigress/wenchang/data/cesm/LME/{data_name}/b.e11.BLMTRC5CN.f19_g16.iii.cam.h0.{data_name}.085001-184912.nc'
ens = range(1,14)

for i in ens:
    s = str(i).zfill(3)
    ifile1 = ifile_template.replace('iii', s)
    ifile2 = ifile1.replace('085001-184912', '185001-200512')
    ofile = ifile1.replace('085001-184912', '085001-200512').split('/')[-1]
    cmd = f'ncrcat -v {data_name} {ifile1} {ifile2} {ofile}'
    run_shell(cmd)