dataset = "icesat4_basins"
# resolution for the output file, set from config.
# resolution = cf.resolution # in km
# conservative regridding for bedmap2 and albmap data. does
# not yet work for the other datasets.
use_conservative_regridding = False
data_path = os.path.join(cf.output_data_path, dataset)
inputfile = os.path.join(data_path, 'icesat4_basins_1km_input.nc')
# pi.prepare_ncfile_for_cdo(inputfile)
cdo_targetgrid_file, regridded_file = pi.get_filenames_for_cdo(
cf.cdo_remapgridpath, data_path, dataset, cf.grid_id)
# Create a bash script that handles the regridding.
# Regridding can be a CPU-heavy task. Choose cluster_regridding=True in config.py
# If you want to submit to the cluster using SLURM.
# use 'sbatch cdo_remap.sh' to submit your job.
# Conservative regridding does not work for all datasets yet, use it for bedmap2 or albmap.
# We use cdo, see https://code.zmaw.de/projects/cdo/embedded/index.html
# FIXME: use remapnn instead of remapbil in script.
pi.write_regrid_command_file(cf, data_path, dataset, inputfile, cf.grid_id,
cdo_targetgrid_file, regridded_file,
cf.regridding_method)
# ## Integer regridding for basin values, can be run interactively.
# subprocess.check_call("cdo remapnn,"+cdo_targetgrid_file+" "+inputfile+" "+
# regridded_file, shell=True)
reload(cf)
import pism_input.pism_input as pi
reload(pi)
dataset = "accum"
# conservative regridding for bedmap2 and albmap data. does
# not yet work for the other datasets.
use_conservative_regridding = True
data_path = os.path.join(cf.output_data_path, dataset)
# prepare the input file for cdo remapping
# this step takes a while for high resolution data (i.e. 1km)
inputfile = os.path.join(data_path, 'accum_1km_input.nc')
pi.prepare_ncfile_for_cdo(inputfile)
cdo_targetgrid_file, regridded_file = pi.get_filenames_for_cdo(
cf.cdo_remapgridpath, data_path, dataset, cf.grid_id)
# Create a bash script that handles the regridding.
# Regridding can be a CPU-heavy task. Choose cluster_regridding=True in config.py
# If you want to submit to the cluster using SLURM.
# use 'sbatch cdo_remap.sh' to submit your job.
# Conservative regridding does not work for all datasets yet, use it for bedmap2 or albmap.
# We use cdo, see https://code.zmaw.de/projects/cdo/embedded/index.html
pi.write_regrid_command_file(cf, data_path, dataset, inputfile, cf.grid_id,
cdo_targetgrid_file, regridded_file,
use_conservative_regridding)
