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):
        root_path = '/gpfs/fs1/Data/ccmc/RT_DATA_TREE_ROOT/ccmc_rt_data_tree/model/ionosphere/WAM-IPE/WAM-IPE-July2021-22/{:04d}/{:02d}/'.format(iyear,imonth)
        fpaths = glob.glob(root_path + '*/gsm10.*.nc')
        fnames = [os.path.basename(x) for x in fpaths]
        print(np.sort(fnames),len(fnames))
        if len(fnames) !=0:
            os.system('mkdir {:04d}{:02d}'.format(iyear,imonth))
            for fpath, fname in zip(fpaths,fnames):
                os.system('cp {} {:04d}{:02d}/{}'.format(fpath,iyear,imonth,fname))