Exemple #1
0
def do_day(yyyymmdd='20200101', n=35,m=10):
    for i in range(n):
        date = ndates(i * 24, yyyymmdd, fmt='%Y-%m-%d')
        print(date)
        fin = f'{p_in}/*{date}*.nc'
        fout = f'{p_out}/ufs_land_output.{date}.nc'
        ncea(fin, fout)
        print('deleting...')
        run2(f'rm -f {fin}')
        time.sleep(m)
Exemple #2
0
def ncea(fin='/p/t*.nc', fout='/p/te.nc'):
    run2(f'/cpc/home/li.xu/anaconda3/bin/ncea {fin} -O {fout}')
def ufs(restart='20210821', n=35):
    date = ndates(0, restart, fmt='%Y-%m-%d')
    nml(date, n)
    run2('./ufsLand.exe')
    from mean2day import do_day2
    do_day2(restart, n, 1)
from xu import ndates,run2,last_day

#/cpc/drought/pdf/noahmp/conus_0p5/restart/ufs_land_restart.2021-01-31_23-00-00.nc
def keep(f):
    '''last day of month to keep'''
    dt=f[-22:-3] # totla 19char
    date=dt[:10]
    yyyy=date[:4]
    mm=date[5:7]
    dd=date[8:10]

    if dd == last_day(yyyy,mm):
        return True
    else:
        return False


p='/cpc/drought/pdf/noahmp/conus_0p5/restart/ufs_land_restart.*.nc'

from glob import glob

files=glob(p)

for f in files:
    if not keep(f):
        run2(f'mv {f} /cpc/drought/pdf/noahmp/conus_0p5/trash')
    else:
        print('keep:',f)