import os
import glob
import numpy as np

for iyear in [2021,2022]:
    for imonth in range(1,13):
#for iyear in [2022]:
#    for imonth in range(1,2):
        for iday in range(1,31):
            in_path = '/gpfs/fs1/Data/ccmc/RT_DATA_TREE_ROOT/ccmc_rt_data_tree/model/ionosphere/WAM-IPE/postprocessing/post_files/{:04d}{:02d}/'.format(iyear,imonth)
            out_path = '/gpfs/fs1/Data/ccmc/RT_DATA_TREE_ROOT/ccmc_rt_data_tree/model/ionosphere/WAM-IPE/WAM-IPE-July2021-22/{:04d}/{:02d}/{:02d}/'.format(iyear,imonth,iday)
            fin_paths = np.sort(glob.glob(in_path + 'wam_fixed_height.gsm10.{:04d}{:02d}{:02d}*.nc'.format(iyear,imonth,iday)))
            fnames = [os.path.basename(x) for x in fin_paths]
        
            if len(fnames) !=0:
                #print(out_path + 'postprocessed/')
                #os.system('mkdir {}'.format(out_path + 'postprocessed/'))
                for fin_path, fname in zip(fin_paths,fnames):
                    os.system('\cp {} {}'.format(fin_path,out_path + 'postprocessed/' + fname)) 

        print('{:04d}{:02d} done.'.format(iyear,imonth))