def do_file(file, src_grd, dst_grd):
dst_dir = './SODA/'
zeta = remap_clm(file, 'ssh', src_grd, dst_grd, dst_dir=dst_dir)
dst_grd = pyroms.grid.get_ROMS_grid('ARCTIC2', zeta=zeta)
remap_clm(file, 'temp', src_grd, dst_grd, dst_dir=dst_dir)
remap_clm(file, 'salt', src_grd, dst_grd, dst_dir=dst_dir)
remap_clm_uv(file, src_grd, dst_grd, dst_dir=dst_dir)
# merge file
clim_file = dst_dir + file.rsplit(
'/')[-1][:-3] + '_clim_' + dst_grd.name + '.nc'
out_file = dst_dir + file.rsplit(
'/')[-1][:-3] + '_ssh_clim_' + dst_grd.name + '.nc'
command = ('ncks', '-a', '-O', out_file, clim_file)
# subprocess.check_call(command)
subprocess.call(command)
os.remove(out_file)
out_file = dst_dir + file.rsplit(
'/')[-1][:-3] + '_temp_clim_' + dst_grd.name + '.nc'
command = ('ncks', '-a', '-A', out_file, clim_file)
# subprocess.check_call(command)
subprocess.call(command)
os.remove(out_file)
out_file = dst_dir + file.rsplit(
'/')[-1][:-3] + '_salt_clim_' + dst_grd.name + '.nc'
command = ('ncks', '-a', '-A', out_file, clim_file)
# subprocess.check_call(command)
subprocess.call(command)
os.remove(out_file)
out_file = dst_dir + file.rsplit(
'/')[-1][:-3] + '_u_clim_' + dst_grd.name + '.nc'
command = ('ncks', '-a', '-A', out_file, clim_file)
# subprocess.check_call(command)
subprocess.call(command)
os.remove(out_file)
out_file = dst_dir + file.rsplit(
'/')[-1][:-3] + '_v_clim_' + dst_grd.name + '.nc'
command = ('ncks', '-a', '-A', out_file, clim_file)
# subprocess.check_call(command)
subprocess.call(command)
os.remove(out_file)
def do_file(file, src_grd, dst_grd):
dst_dir='./SODA/'
zeta = remap_clm(file, 'ssh', src_grd, dst_grd, dst_dir=dst_dir)
dst_grd = pyroms.grid.get_ROMS_grid('ARCTIC2', zeta=zeta)
remap_clm(file, 'temp', src_grd, dst_grd, dst_dir=dst_dir)
remap_clm(file, 'salt', src_grd, dst_grd, dst_dir=dst_dir)
remap_clm_uv(file, src_grd, dst_grd, dst_dir=dst_dir)
# merge file
clim_file = dst_dir + file.rsplit('/')[-1][:-3] + '_clim_' + dst_grd.name + '.nc'
out_file = dst_dir + file.rsplit('/')[-1][:-3] + '_ssh_clim_' + dst_grd.name + '.nc'
command = ('ncks', '-a', '-O', out_file, clim_file)
# subprocess.check_call(command)
subprocess.call(command)
os.remove(out_file)
out_file = dst_dir + file.rsplit('/')[-1][:-3] + '_temp_clim_' + dst_grd.name + '.nc'
command = ('ncks', '-a', '-A', out_file, clim_file)
# subprocess.check_call(command)
subprocess.call(command)
os.remove(out_file)
out_file = dst_dir + file.rsplit('/')[-1][:-3] + '_salt_clim_' + dst_grd.name + '.nc'
command = ('ncks', '-a', '-A', out_file, clim_file)
# subprocess.check_call(command)
subprocess.call(command)
os.remove(out_file)
out_file = dst_dir + file.rsplit('/')[-1][:-3] + '_u_clim_' + dst_grd.name + '.nc'
command = ('ncks', '-a', '-A', out_file, clim_file)
# subprocess.check_call(command)
subprocess.call(command)
os.remove(out_file)
out_file = dst_dir + file.rsplit('/')[-1][:-3] + '_v_clim_' + dst_grd.name + '.nc'
command = ('ncks', '-a', '-A', out_file, clim_file)
# subprocess.check_call(command)
subprocess.call(command)
os.remove(out_file)
lst_file = lst_file + lst
print('Build CLM file from the following file list:')
print(lst_file)
print(' ')
src_grd = pyroms_toolbox.Grid_HYCOM.get_nc_Grid_HYCOM('/archive/u1/uaf/kate/HYCOM/Svalbard/HYCOM_GLBa0.08_North_grid2.nc')
dst_grd = pyroms.grid.get_ROMS_grid('ARCTIC2')
for file in lst_file:
# remap
zeta = remap_clm(file, 'ssh', src_grd, dst_grd, dst_dir=dst_dir)
dst_grd = pyroms.grid.get_ROMS_grid('ARCTIC2', zeta=zeta)
remap_clm(file, 'temp', src_grd, dst_grd, dst_dir=dst_dir)
remap_clm(file, 'salt', src_grd, dst_grd, dst_dir=dst_dir)
remap_clm_uv(file, src_grd, dst_grd, dst_dir=dst_dir)
# merge file
clim_file = dst_dir + file.rsplit('/')[-1][:-3] + '_clim_' + dst_grd.name + '.nc'
out_file = dst_dir + file.rsplit('/')[-1][:-3] + '_ssh_clim_' + dst_grd.name + '.nc'
command = ('ncks', '-a', '-O', out_file, clim_file)
print(command)
subprocess.check_call(command)
os.remove(out_file)
out_file = dst_dir + file.rsplit('/')[-1][:-3] + '_temp_clim_' + dst_grd.name + '.nc'
command = ('ncks', '-a', '-A', out_file, clim_file)
print(command)
subprocess.check_call(command)
os.remove(out_file)
out_file = dst_dir + file.rsplit('/')[-1][:-3] + '_salt_clim_' + dst_grd.name + '.nc'
lst_file = lst_file + lst
print 'Build CLM file from the following file list:'
print lst_file
print ' '
src_grd = pyroms_toolbox.Grid_HYCOM.get_nc_Grid_HYCOM('/archive/u1/uaf/kate/HYCOM/Svalbard/HYCOM_GLBa0.08_North_grid2.nc')
dst_grd = pyroms.grid.get_ROMS_grid('ARCTIC2')
for file in lst_file:
# remap
zeta = remap_clm(file, 'ssh', src_grd, dst_grd, dst_dir=dst_dir)
dst_grd = pyroms.grid.get_ROMS_grid('ARCTIC2', zeta=zeta)
remap_clm(file, 'temp', src_grd, dst_grd, dst_dir=dst_dir)
remap_clm(file, 'salt', src_grd, dst_grd, dst_dir=dst_dir)
remap_clm_uv(file, src_grd, dst_grd, dst_dir=dst_dir)
# merge file
clim_file = dst_dir + file.rsplit('/')[-1][:-3] + '_clim_' + dst_grd.name + '.nc'
out_file = dst_dir + file.rsplit('/')[-1][:-3] + '_ssh_clim_' + dst_grd.name + '.nc'
command = ('ncks', '-a', '-O', out_file, clim_file)
print command
subprocess.check_call(command)
os.remove(out_file)
out_file = dst_dir + file.rsplit('/')[-1][:-3] + '_temp_clim_' + dst_grd.name + '.nc'
command = ('ncks', '-a', '-A', out_file, clim_file)
print command
subprocess.check_call(command)
os.remove(out_file)
out_file = dst_dir + file.rsplit('/')[-1][:-3] + '_salt_clim_' + dst_grd.name + '.nc'