drc_pwd=${PWD}

unset IFS

spt_src="${BASH_SOURCE[0]}"

[[ -z "${spt_src}" ]] && spt_src="${0}" # Use ${0} when BASH_SOURCE is unavailable (e.g., dash)

while [ -h "${spt_src}" ]; do # Recursively resolve ${spt_src} until file is no longer a symlink

drc_spt="$( cd -P "$( dirname "${spt_src}" )" && pwd )"

spt_src="$(readlink "${spt_src}")"

[[ ${spt_src} != /* ]] && spt_src="${drc_spt}/${spt_src}" # If ${spt_src} was relative symlink, resolve it relative to path where symlink file was located

cmd_ln="${spt_src} ${@}"

drc_spt="$( cd -P "$( dirname "${spt_src}" )" && pwd )"

spt_nm=$(basename ${spt_src}) # [sng] Script name (unlike $0, ${BASH_SOURCE[0]} works well with 'source <script>')

spt_pid=$$ # [nbr] Script PID (process ID)

if [ -z "${HDF5_USE_FILE_LOCKING}" ]; then

export HDF5_USE_FILE_LOCKING='FALSE'

fi # HDF5_USE_FILE_LOCKING

if [ -z "${HOSTNAME}" ]; then

if [ -f /bin/hostname ] && [ -x /bin/hostname ]; then

export HOSTNAME=`/bin/hostname`

elif [ -f /usr/bin/hostname ] && [ -x /usr/bin/hostname ]; then

export HOSTNAME=`/usr/bin/hostname`

fi # !hostname

fi # HOSTNAME

if [ -z "${DATA}" ]; then

case "${HOSTNAME}" in

acme1* ) DATA="/home/${USER}" ; ;; # LLNL acme1

andes* | titan* ) DATA="/gpfs/alpine/world-shared/cli115/${USER}" ; ;; # OLCF andes compute nodes named andesNNN, 256 GB/node

blues* | blogin* | b[0123456789][0123456789][0123456789] ) DATA="/lcrc/project/ACME/${USER}" ; ;; # ANL/LCRC blues compute nodes named bNNN, 36|64 cores|GB/node

chrysalis* | chrlogin* | chr-[0123456789][0123456789][0123456789][0123456789] ) DATA="/lcrc/project/ACME/${USER}" ; ;; # ANL/LCRC chrysalis compute nodes named chr-NNNN, 64|256 cores|GB/node

*cheyenne* ) DATA="/glade/p/work/${USER}" ; ;; # NCAR cheyenne compute nodes named, e.g., r8i0n8, r5i3n16, r12i5n29 ... 18|(64/128) cores|GB/node (cheyenne login nodes 256 GB)

compy* ) DATA="/qfs/people/${USER}/data" ; ;; # PNNL compy compute nodes all named nNNNN, 40|192 cores|GB/node (compy login nodes also 192 GB)

constance* | node* ) DATA='/scratch' ; ;; # PNNL

cooley* | cc[0123456789][0123456789][0123456789] | mira* ) DATA="/projects/OceanClimate_2/${USER}" ; ;; # ALCF cooley compute nodes named ccNNN, 384 GB/node

cori* ) DATA="${SCRATCH}" ; ;; # NERSC cori compute nodes named nidNNNNN (Haswell/KNL) or cmemNN (AMD) with (32/68/64)|(96/128/1024) cores|GB/node (knl/haswell/amd) (login nodes 512 GB)

perlmutter* | login[0123456789][0123456789] ) DATA="${SCRATCH}" ; ;; # NERSC perlmutter compute nodes named nidNNNNNN (CPU) with (64)|(512) cores|GB/node (cpu) (login nodes 512 GB)

theta* ) DATA="/projects/ClimateEnergy_3/${USER}" ; ;; # ALCF theta compute nodes named fxm, 64|192 cores|GB/node

* ) DATA='/tmp' ; ;; # Other

esac # !HOSTNAME

fi # DATA

OMP_PROC_BIND=false

hrd_pth='Yes' # [sng] Hard-code machine-dependent paths/modules if HOSTNAME in database

if [ "${hrd_pth}" = 'Yes' ] && [ "${NCO_PATH_OVERRIDE}" = 'Yes' ]; then

# If HOSTNAME is not in database, change hrd_pth_fnd to 'No' in case-statement default fall-through

hrd_pth_fnd='Yes' # [sng] Machine-dependent paths/modules for HOSTNAME found in database

case "${HOSTNAME}" in

acme1* )

if [ ${spt_nm} = 'ncremap' ]; then

E3SMU_ROOT='/usr/local/e3sm_unified/envs/base/envs/e3sm_unified_latest'

fi # !ncremap

export PATH='/home/zender1/bin'\:${PATH}

export LD_LIBRARY_PATH='/home/zender1/lib'\:${LD_LIBRARY_PATH} ; ;;

andes* )

# 20190827: Must guarantee finding mpirun

source ${MODULESHOME}/init/sh # 20150607: PMC Ensures find module commands will be found

if [ ${spt_nm} = 'ncremap' ]; then

E3SMU_ROOT='/ccs/proj/cli900/sw/rhea/e3sm-unified/base/envs/e3sm_unified_latest'

fi # !ncremap

export PATH='/ccs/home/zender/bin_andes'\:${PATH}

export LD_LIBRARY_PATH='/ccs/home/zender/lib_andes'\:${LD_LIBRARY_PATH} ; ;;

blues* | blogin* | b[0123456789][0123456789][0123456789] )

if [ ${spt_nm} = 'ncremap' ]; then

E3SMU_ROOT='/lcrc/soft/climate/e3sm-unified/base/envs/e3sm_unified_1.6.0_anvil'

fi # !ncremap

export PATH='/home/zender/bin_blues'\:${PATH}

export LD_LIBRARY_PATH='/home/zender/lib_blues'\:${LD_LIBRARY_PATH} ; ;;

chrysalis* | chrlogin* | chr-[0123456789][0123456789][0123456789][0123456789] )

if [ ${spt_nm} = 'ncremap' ]; then

E3SMU_ROOT='/lcrc/soft/climate/e3sm-unified/base/envs/e3sm_unified_1.6.0_chrysalis'

fi # !ncremap

export PATH='/gpfs/fs1/soft/chrysalis/spack/opt/spack/linux-centos8-x86_64/gcc-9.2.0/netcdf-c-4.7.4-a4uk6zy/bin:/home/zender/bin_chrysalis'\:${PATH}

export LD_LIBRARY_PATH='/gpfs/fs1/soft/chrysalis/spack/opt/spack/linux-centos8-x86_64/gcc-9.2.0/netcdf-c-4.7.4-a4uk6zy/lib:/home/zender/lib_chrysalis'\:${LD_LIBRARY_PATH} ; ;;

*cheyenne* )

# 20180112: Cheyenne support not yet tested in batch mode

if [ ${spt_nm} = 'ncremap' ]; then

# On cheyenne, module load ncl installs ERWG in /glade/u/apps/ch/opt/ncl/6.4.0/intel/17.0.1/bin (i.e., ${NCARG_ROOT}/bin)

module load ncl

fi # !ncremap

if [ -n "${NCARG_ROOT}" ]; then

export PATH="${PATH}:/glade/u/apps/ch/opt/ncl/6.6.2/gnu/8.3.0/bin"

fi # !NCARG_ROOT

export PATH='/glade/u/home/zender/bin'\:${PATH}

export LD_LIBRARY_PATH='/glade/u/apps/ch/opt/netcdf/4.6.3/gnu/9.1.0/lib:/glade/u/apps/ch/opt/udunits/2.2.26/gnu/9.1.0/lib:/glade/u/apps/ch/opt/gsl/2.4/gnu/6.3.0/lib:/glade/u/home/zender/lib'\:${LD_LIBRARY_PATH} ; ;;

compy* )

if [ ${spt_nm} = 'ncremap' ]; then

# 20210519: This script takes significant time (5-10 seconds) to load

source /compyfs/software/mbtempest.envs.sh

E3SMU_ROOT='/share/apps/E3SM/conda_envs/base/envs/e3sm_unified_1.6.0_compy'

fi # !ncremap

export PATH='/qfs/people/zender/bin'\:${PATH}

export LD_LIBRARY_PATH='/qfs/people/zender/lib'\:${LD_LIBRARY_PATH} ; ;;

cooley* | cc[0123456789][0123456789][0123456789] )

# 20160421: Split cooley from mira binary locations to allow for different system libraries

# http://www.mcs.anl.gov/hs/software/systems/softenv/softenv-intro.html

soft add +mvapich2

export PBS_NUM_PPN=12 # Spoof PBS on Soft (which knows nothing about node capabilities)

export PATH='/home/zender/bin_cooley'\:${PATH}

export LD_LIBRARY_PATH='/home/zender/lib_cooley'\:${LD_LIBRARY_PATH} ; ;;

cori* )

if [ ${spt_nm} = 'ncremap' ]; then

MOAB_DIR=/project/projectdirs/e3sm/software/moab

TEMPESTREMAP_DIR=/project/projectdirs/e3sm/software/tempestremap

E3SMU_ROOT='/global/common/software/e3sm/anaconda_envs/e3sm_unified_latest'

fi # !ncremap

if [ -n "${NCARG_ROOT}" ]; then

export PATH="${PATH}:${NCARG_ROOT}/bin"

export LD_LIBRARY_PATH="${LD_LIBRARY_PATH}:${NCARG_ROOT}/lib"

fi # !NCARG_ROOT

export PATH='/global/homes/z/zender/bin_cori'\:${PATH}

export LD_LIBRARY_PATH='/global/homes/z/zender/lib_cori'\:${LD_LIBRARY_PATH} ; ;;

mira* )

export PATH='/home/zender/bin_mira'\:${PATH}

export LD_LIBRARY_PATH='/soft/libraries/netcdf/current/library:/home/zender/lib_mira'\:${LD_LIBRARY_PATH} ; ;;

theta* )

export PATH='/opt/cray/pe/netcdf/4.6.1.2/gnu/7.1/bin'\:${PATH}

export LD_LIBRARY_PATH='/opt/cray/pe/netcdf/4.6.1.2/gnu/7.1/lib'\:${LD_LIBRARY_PATH} ; ;;

titan* )

source ${MODULESHOME}/init/sh # 20150607: PMC Ensures find module commands will be found

module load gcc

if [ ${spt_nm} = 'ncremap' ]; then

# 20170831: Use module load ncl (6.3.0 lacks ERWG)

module load ncl # 20170916 OK

fi # !ncremap

if [ -n "${NCARG_ROOT}" ]; then

export PATH="${PATH}:${NCARG_ROOT}/bin"

fi # !NCARG_ROOT

export PATH='/ccs/home/zender/bin_titan'\:${PATH}

export LD_LIBRARY_PATH='/opt/cray/netcdf/4.4.1.1/GNU/49/lib:/sw/xk6/udunits/2.1.24/sl_gcc4.5.3/lib:/ccs/home/zender/lib_titan'\:${LD_LIBRARY_PATH} ; ;;

* ) # Default fall-through

hrd_pth_fnd='No' ; ;;

esac # !HOSTNAME

if [ -n "${MOAB_DIR}" ]; then

export PATH="${PATH}:${MOAB_DIR}/bin"

export LD_LIBRARY_PATH="${LD_LIBRARY_PATH}:${MOAB_DIR}/lib"

fi # !MOAB_DIR

if [ -n "${TEMPESTREMAP_DIR}" ]; then

export PATH="${PATH}:${TEMPESTREMAP_DIR}/bin"

export LD_LIBRARY_PATH="${LD_LIBRARY_PATH}:${TEMPESTREMAP_DIR}/lib"

fi # !TEMPESTREMAP_DIR

# 20210519: E3SM-U supplies many uncommon executables, e.g., ESMF_RegridWeightGen

# Append E3SM-U to end of PATH so NCO binaries not redirected, e.g., from CSZ's development directory to E3SM-U

if [ -n "${E3SMU_ROOT}" ]; then

export PATH="${PATH}:${E3SMU_ROOT}/bin"

fi # !MOAB_DIR

fi # !hrd_pth && !NCO_PATH_OVERRIDE

if [ "${E3SMU_MPI}" = 'NOMPI' ] && [ -n "${E3SMU_SCRIPT}" ] && [ -n "${CONDA_PREFIX}" ]; then

export PATH="${CONDA_PREFIX}/bin"\:${PATH}

fi # !E3SMU_MPI

nco_exe=`which ncks`

if [ -z "${nco_exe}" ]; then

echo "ERROR: Unable to find NCO, nco_exe = ${nco_exe}"

echo "${spt_nm}: HINT Carefully examine your environment setup (e.g., .bashrc) to avoid inadvertently overriding (with, e.g., conda-initialization) paths intended to be provided by an analysis-package environment (e.g., E3SM-Unified)"

exit 1

fi # !nco_exe

while [ -h "${nco_exe}" ]; do

drc_nco="$( cd -P "$( dirname "${nco_exe}" )" && pwd )"

nco_exe="$(readlink "${nco_exe}")"

[[ ${nco_exe} != /* ]] && nco_exe="${drc_nco}/${nco_exe}"

drc_nco="$( cd -P "$( dirname "${nco_exe}" )" && pwd )"

nco_vrs=$(ncks --version 2>&1 > /dev/null | grep NCO | awk '{print $5}')

nco_sng=$(ncks --version 2>&1 > /dev/null | grep NCO | awk -F '"' '{print $2}')

if [ -z "${nco_vrs}" ]; then

echo "${spt_nm}: ERROR ${nco_exe} dies with error message on next line:"

$(ncks --version)

if [ "${NCO_PATH_OVERRIDE}" != 'Yes' ]; then

printf "HINT: Run-time errors due to link issues (e.g., libraries missing or not found) might be solved at supported national labs (ALCF, NCAR, NERSC, OLCF, PNNL) by employing NCO machine-dependent hardcoded paths/modules. To try this, re-run command after setting \"export NCO_PATH_OVERRIDE=Yes\".\n"

fi # !NCO_PATH_OVERRIDE

exit 1

fi # !nco_vrs

lbr_vrs=$(ncks --library 2>&1 > /dev/null | awk '{print $6}')

if [ "${drc_spt}" != "${drc_nco}" ]; then

echo "INFO: Mixture of NCO scripts and binaries from different locations. Script ${spt_nm} is from directory ${drc_spt} while NCO binaries are from directory ${drc_nco}. Normally the script and binaries are from the same executables directory. This INFO may be safely ignored for customized scripts and/or binaries that the user has intentionally split into different directories."

echo "HINT: Conflicting script and binary directories may result from 1) Hardcoding an NCO script and/or binary pathnames, 2) Incomplete NCO installations in one or more directories in the \$PATH environment variable, 3) (Re-)Installing or (re-)building NCO without issuing a \"hash -r\" command afterward to update the executable pathnames that the shell remembers from before."

fi # drc_spt

if [ -n "${TERM}" ]; then

# Set fonts for legibility

if [ -x /usr/bin/tput ] && tput setaf 1 &> /dev/null; then

fnt_bld=`tput bold` # Bold

fnt_nrm=`tput sgr0` # Normal

fnt_rvr=`tput smso` # Reverse

fnt_tlc=`tput sitm` # Italic

else

fnt_bld="\e[1m" # Bold

fnt_nrm="\e[0m" # Normal

fnt_rvr="\e[07m" # Reverse

fnt_tlc="\e[3m" # Italic

fi # !tput

fi # !TERM

par_bck='background' # [sng] Parallelism: background

par_mpi='mpi' # [sng] Parallelism: MPI

par_srl='serial' # [sng] Parallelism: serial

vrb_0=0 # [enm] Verbosity level: Quiet

vrb_1=1 # [enm] Verbosity level: Standard, minimal file I/O

vrb_2=2 # [enm] Verbosity level: All file I/O

vrb_3=3 # [enm] Verbosity level: English

vrb_4=4 # [enm] Verbosity level: Pedantic

fmt_sng_nc3='CLASSIC (netCDF3 classic CDF1)'

fmt_sng_nc4='NETCDF4 (netCDF4 extended HDF5)'

fmt_sng_nc5='64BIT_DATA (netCDF3/pnetCDF CDF5)'

fmt_sng_nc6='64BIT_OFFSET (netCDF3 64bit CDF2)'

fmt_sng_nc7='NETCDF4_CLASSIC (netCDF4 classic HDF5)'

a2o_flg='No' # [flg] Atmosphere-to-ocean (only used by Tempest mesh generator)

add_fll='No' # [flg] Add _FillValue to fields with empty destination cells

alg_typ='nco_con' # [sng] Algorithm for remapping (bilinear|conserve|conserve2nd|nearestdtos|neareststod|patch|tempest|se2fv_flx|se2fv_stt|se2fv_alt|se2se|fv2se_flx|fv2se_stt|fv2se_alt|fv2fv_flx|fv2fv_stt|nco_con|nco_dwe)

att_flg='No' # [flg] Add/alter attributes (e.g., _FillValue) before regridding

bch_pbs='No' # [sng] PBS batch (non-interactive) job

bch_slr='No' # [sng] SLURM batch (non-interactive) job

chk_map='Yes' # [flg] Check map-file quality

cln_flg='Yes' # [flg] Clean-up (remove) intermediate files before exiting

clm_flg='No' # [flg] Invoked by ncclimo script

cmp_sng='' # [sng] Compression string

cnv_exe_mbt='mbconvert' # [sng] MOAB-Tempest SCRIP/Exodus conversion to H5M executable

cnv_exe_tps='ConvertSCRIPToExodus' # [sng] TempestRemap SCRIP-to-Exodus conversion executable

cnv_opt_e2m='-o PARALLEL=WRITE_PART -O PARALLEL=BCAST_DELETE -O PARTITION=TRIVIAL -O PARALLEL_RESOLVE_SHARED_ENTS' # [sng] Recommended options to mbconvert for (high-resolution) Exodus to MOAB grids https://acme-climate.atlassian.net/wiki/spaces/ED/pages/932380820/Offline+remapping+workflow+with+mbtempest

d2f_flg='No' # [flg] Convert double-precision fields to single-precision

d2f_opt='-M dbl_flt' # [sng] Option string to convert double-precision fields to single-precision

dbg_lvl=0 # [nbr] Debugging level

dfl_lvl='' # [enm] Deflate level

dpt_exe_mpas='add_depth.py' # [sng] Depth coordinate addition command for MPAS

dpt_flg='No' # [flg] Add depth coordinate to MPAS files

dpt_fl='' # [sng] Depth file with refBottomDepth for MPAS ocean

drc_in='' # [sng] Input file directory

drc_in_xmp='~/drc_in' # [sng] Input file directory for examples

drc_out="${drc_pwd}" # [sng] Output file directory

drc_out_xmp="~/rgr" # [sng] Output file directory for examples

dst_fl='' # [sng] Destination file

dst_xmp='dst.nc' # [sng] Destination file for examples

dt_sng=`date +%Y%m%d` # [sng] Date string for MWF map names

dvn_flg='Yes' # [flg] Send internal NCO regridder output to /dev/null

esmf_typ='' # [sng] ESMF extrapolation type (nearestidavg|neareststod|none)

fl_fmt='' # [enm] Output file format

fl_nbr=0 # [nbr] Number of files to remap

flg_grd_only='No' # [flg] Create grid then exit before regridding

flg_hrz='No' # [flg] Perform horizontal regridding

gaa_sng="--gaa remap_script=${spt_nm} --gaa remap_hostname=${HOSTNAME} --gaa remap_version=${nco_vrs}" # [sng] Global attributes to add

grd_dst='' # [sng] Destination grid-file

grd_dst_xmp='grd_dst.nc' # [sng] Destination grid-file for examples

grd_sng='' # [sng] Grid string

grd_src='' # [sng] Source grid-file

grd_src_xmp='grd_src.nc' # [sng] Source grid-file for examples

grd_vrt_xmp='grd_vrt.nc' # [sng] Vertical grid-file for examples

hdr_pad='10000' # [B] Pad at end of header section

hnt_dst='' # [sng] ERWG hint for destination grid

hnt_src='' # [sng] ERWG hint for source grid

in_fl='' # [sng] Input file

in_xmp='in.nc' # [sng] Input file for examples

inp_aut='No' # [sng] Input file list automatically generated (in ncclimo, or specified with -i in ncremap)

inp_glb='No' # [sng] Input file list from globbing directory

inp_psn='No' # [sng] Input file list from positional arguments

inp_std='No' # [sng] Input file list from stdin

job_nbr=2 # [nbr] Job simultaneity for parallelism

map_fl='' # [sng] Map-file

map_nfo_sng='' # [sng] Informational string about map-file

map_rsl_fl='' # [sng] File containing results of weight-generation command (i.e., map_fl or map_trn_fl for monotr)

map_trn_fl='' # [sng] Map-file transpose (for Tempest)

map_mk='No' # [flg] Generate map-file (i.e., map does not yet exist)

map_usr_flg='No' # [flg] User supplied argument to --map option

map_xmp='map.nc' # [sng] Map-file for examples

mbt_typ_grd=0 # [flg] mbtempest mesh type (CS=0, RLL=1, ICO=2, OVERLAP_FILES=3, OVERLAP_MEMORY=4, OVERLAP_MOAB=5])

mbt_typ_map=5 # [flg] mbtempest mesh type (CS=0, RLL=1, ICO=2, OVERLAP_FILES=3, OVERLAP_MEMORY=4, OVERLAP_MOAB=5])

mem_mb=0 # [MB] Megabytes of RAM per srun job in Cori SLURM in MPI mode (0 indicates unlimited RAM)

mlt_map_flg='Yes' # [sng] Multi-map flag

mpi_flg='No' # [sng] Parallelize over nodes

mpi_nbr=1 # [sng] Number of tasks-per-node that MPI-enabled weight generators will request (should be factor of prt_nbr)

msh_exe_tps='GenerateOverlapMesh' # [sng] TempestRemap mesh-generation executable

msh_fl='' # [sng] Mesh-file (for Tempest)

msk_apl='No' # [flg] Apply msk_out to variables after regridding

msk_dst='' # [sng] Mask-template variable in destination file

msk_out='No' # [flg] Put mask variable in regridded file

msk_src='' # [sng] Mask-template variable in source file

mss_val='-9.99999979021476795361e+33' # [frc] Missing value for MPAS (ocean+seaice)

nco_opt='--no_tmp_fl' # [sng] NCO options (e.g., '-6 -t 1')

nd_nbr=1 # [nbr] Number of nodes

no_stg_grd='Yes' # [flg] Omit staggered grid variables

nsx_alg='' # [enm] Intersection algorithm (mbtempest only)

out_fl='' # [sng] Output file

out_xmp='out.nc' # [sng] Output file for examples

par_typ="${par_srl}" # [sng] Parallelism type

pdq_flg='Yes' # [sng] Permute dimensions before regridding

pdq_opt='' # [sng] ncpdq dimension permutation option

ppc_prc='' # [nbr] Precision-preserving compression number of significant digits

prc_sgs='No' # [sng] SGS requested

prc_typ='' # [sng] Procedure type

prt_alg='zoltan' # [sng] MOAB-Tempest H5M partition algorithm

prt_exe_mbt='mbpart' # [sng] MOAB-Tempest H5M partition executable

prt_nbr=8 # [nbr] Partition number for mbpart (should be multiple of mpi_nbr)

prs_stt='' # [sng] Preserved statistic ('integral' or 'mean')

qea_flg='No' # [sng] Quasi-equal area mode for TempestRemap

rgr_opt='--rgr lat_nm_out=lat --rgr lon_nm_out=lon' # [sng] Regridding options

rnr_thr='' # [frc] Renormalization option

se_np_nbr=4 # [nbr] Spectral Element polynomial order of discretization (for TempestRemap SE grids)

sgs_frc='' # [sng] Sub-grid fraction variable

sgs_frc_usr='' # [sng] Sub-grid fraction variable

sgs_msk='' # [sng] Sub-grid mask variable

sgs_nrm='1.0' # [frc] Sub-grid normalization

skl_fl='' # [sng] Skeleton file

std_chk='Yes' # [sng] Check stdin for input file list

std_flg='No' # [sng] Input available from pipe to stdin

thr_nbr=2 # [nbr] Thread number for regridder

trn_map='No' # [flg] Tempest transpose map (i.e., fv2se_flx == monotr)

ugrid_fl='' # [sng] UGRID file

unq_sfx=".pid${spt_pid}" # [sng] Unique suffix

var_lst='' # [sng] Variables to process (empty means all)

var_rgr='' # [sng] CF template variable

var_xmp='FSNT' # [sng] Variable list for examples

vrb_lvl=${vrb_2} # [enm] Verbosity level

vrs_prn='No' # [sng] Print version information

vrt_fl='' # [sng] Vertical coordinate file

vrt_nm='' # [sng] Vertical coordinate name

vrt_opt='' # [sng] Vertical options (e.g., '--rgr xtr_mth=mss_val')

vrt_ntp='' # [sng] Vertical interpolation type (lin|log)

vrt_xtr='' # [sng] Vertical extrapolation type (mss_val|nrs_ngh)

wgt_exe_cdo='cdo' # [sng] ESMF executable (weight generation only)

wgt_exe_esmf='ESMF_RegridWeightGen' # [sng] ESMF executable (weight generation only)

wgt_exe_mbt='mbtempest' # [sng] MOAB-Tempest executable (both mesh- and weight-generation)

wgt_exe_nco='ncks' # [sng] NCO executable (both mesh- and weight-generation)

wgt_exe_tps='GenerateOfflineMap' # [sng] TempestRemap weight-generation executable

wgt_typ='' # [sng] Weight-generator program ('esmf', 'mbtempest', 'nco', or 'tempest')

wgt_opt='' # [sng] Weight-generator options

wgt_opt_cdo='' # [sng] CDO weight-generation options

wgt_opt_esmf='--no_log --ignore_unmapped' # [sng] ESMF_RegridWeightGen weight-generation options (ESMF < 7.1.0r)

wgt_opt_nco='-O --dmm_in_mk' # [sng] NCO weight-generation options

wgt_opt_tps='' # [sng] TempestRemap weight-generation options

msh_opt_tps='--allow_no_overlap' # [sng] TempestRemap mesh-generation options (20210421: add allow_no_overlap required by regional-to-regional grids like GrIS)

xcl_flg='No' # [sng] Exclude rather than extract variable list

xtn_var='' # [sng] Extensive variables (e.g., 'TSurfStd_ct')

xtr_nsp=8 # [nbr] Extrapolation number of source points (for ESMF & NCO only)

xtr_xpn=2.0 # [frc] Extrapolation exponent (absolute value) (for ESMF & NCO only)

if [ -d "${TMPDIR}" ]; then

# Fancy %/ syntax removes trailing slash (e.g., from $TMPDIR)

drc_tmp="${TMPDIR%/}"

elif [ -d '/tmp' ]; then

drc_tmp='/tmp'

drc_tmp=${PWD}

fi # !gpfs

function fnc_usg_prn { # NB: dash supports fnc_nm (){} syntax, not function fnc_nm{} syntax

} # end fnc_usg_prn()

function ncvarlst { ncks --trd -m ${1} | grep -E ': type' | cut -f 1 -d ' ' | sed 's/://' | sort ; }

function ncdmnlst { ncks --cdl -m ${1} | cut -d ':' -f 1 | cut -d '=' -s -f 1 ; }

function dst_is_grd {

} # end dst_is_grd()

arg_nbr=$#

if [ ${arg_nbr} -eq 0 ]; then

fnc_usg_prn

fi # !arg_nbr

while getopts :34567a:CD:d:f:g:G:h:I:i:j:L:Mm:n:O:o:P:p:R:r:s:T:t:Uu:V:v:W:w:x:-: OPT; do

case ${OPT} in

3) fl_fmt='3' ;; # File format

4) fl_fmt='4' ;; # File format

5) fl_fmt='5' ;; # File format

6) fl_fmt='6' ;; # File format

7) fl_fmt='7' ;; # File format

a) alg_typ="${OPTARG}" ;; # Algorithm

C) clm_flg='Yes' ;; # Climo flag (undocumented)

D) dbg_lvl="${OPTARG}" ;; # Debugging level

d) dst_fl="${OPTARG}" ;; # Destination file

g) grd_dst="${OPTARG}" ;; # Destination grid-file

G) grd_sng="${OPTARG}" ;; # Grid generation string

I) drc_in="${OPTARG}" ;; # Input directory

i) in_fl="${OPTARG}" ;; # Input file

j) job_usr="${OPTARG}" ;; # Job simultaneity

L) dfl_lvl="${OPTARG}" ;; # Deflate level

M) mlt_map_flg='No' ;; # Multi-map flag

m) map_fl="${OPTARG}" ;; # Map-file

n) nco_opt="${OPTARG} ${nco_opt}" ;; # NCO options

O) drc_usr="${OPTARG}" ;; # Output directory

o) out_fl="${OPTARG}" ;; # Output file

P) prc_typ="${OPTARG}" ;; # Procedure type

p) par_typ="${OPTARG}" ;; # Parallelism type

r) rnr_thr="${OPTARG}" ;; # Renormalization threshold

R) rgr_opt="${OPTARG}" ;; # Regridding options

s) grd_src="${OPTARG}" ;; # Source grid-file

T) tmp_usr="${OPTARG}" ;; # Temporary directory

t) thr_usr="${OPTARG}" ;; # Thread number

U) pdq_opt='-U' ;; # Unpack input

u) unq_usr="${OPTARG}" ;; # Unique suffix

V) var_rgr="${OPTARG}" ;; # CF template variable

v) var_lst="${OPTARG}" ;; # Variables

W) wgt_opt_usr="${OPTARG}" ;; # Weight-generator options

w) wgt_usr="${OPTARG}" ;; # Weight-generator command

x) xtn_var="${OPTARG}" ;; # Extensive variables

-) LONG_OPTARG="${OPTARG#*=}"

case ${OPTARG} in

# Hereafter ${OPTARG} is long argument key, and ${LONG_OPTARG}, if any, is long argument value

# Long options with no argument, no short option counterpart

# Long options with argument, no short option counterpart

# Long options with short counterparts, ordered by short option key

a2o | atm2ocn | b2l | big2ltl | l2s | lrg2sml ) a2o_flg='Yes' ;; # # Atmosphere-to-ocean

a2o=?* | atm2ocn=?* | b2l=?* | big2ltl=?* | l2s=?* | lrg2sml=?* ) echo "No argument allowed for --${OPTARG switch}" >&2; exit 1 ;; # # Atmosphere-to-ocean

add_fll | add_fill_value | fll_mpt | fill_empty ) add_fll='Yes' ;; # # Add _FillValue to fields with empty destination cells

add_fll=?* | add_fill_value=?* | fll_mpt=?* | fill_empty=?* ) echo "No argument allowed for --${OPTARG switch}" >&2; exit 1 ;; # # Add _FillValue to fields with empty destination cells

alg_typ=?* | algorithm=?* | regrid_algorithm=?* ) alg_typ="${LONG_OPTARG}" ;; # -a # Algorithm

area_dgn | dgn_area | area_diagnose | diagnose_area ) area_dgn='Yes' ;; # # Diagnose (rather than copy) grid_area to inferred grid-file

area_dgn=?* | dgn_area=?* | area_diagnose=?* | diagnose_area=?* ) echo "No argument allowed for --${OPTARG switch}" >&2; exit 1 ;; # # Diagnose (rather than copy) grid_area in inferred grids

clm_flg | climatology_flag ) clm_flg='Yes' ;; # -C # Climo flag (undocumented)

clm_flg=?* | climatology_flag=?* ) echo "No argument allowed for --${OPTARG switch}" >&2; exit 1 ;; # # clm_flg

cmp=?* | cmp_sng=?* | compression=?* | cdc=?* | codec=?* ) cmp_sng="${LONG_OPTARG}" ;; # # Compression string

d2f | d2s | dbl_flt | dbl_sgl | double_float ) d2f_flg='Yes' ;; # # Convert double-precision fields to single-precision

d2f=?* | d2s=?* | dbl_flt=?* | dbl_sgl=?* | double_float=?* ) echo "No argument allowed for --${OPTARG switch}" >&2; exit 1 ;; # # D2F

dbg_lvl=?* | dbg=?* | debug=?* | debug_level=?* ) dbg_lvl="${LONG_OPTARG}" ;; # -d # Debugging level

devnull=?* | dev_nll=?* | dev_null=?* | dvn_flg=?* ) dvn_flg="${LONG_OPTARG}" ;; # # Send NCO weight-generator messages to /dev/null

dfl_lvl=?* | deflate=?* | dfl=?* ) dfl_lvl="${LONG_OPTARG}" ;; # -L # Deflate level

dpt | depth | add_dpt | add_depth ) dpt_flg='Yes' ;; # # Add depth coordinate to MPAS files

dpt=?* | depth=?* | add_dpt=?* | add_depth=?* ) echo "No argument allowed for --${OPTARG switch}" >&2; exit 1 ;; # # DPT

dpt_fl=?* | mpas_fl=?* | mpas_file=?* | depth_file=?* ) dpt_fl="${LONG_OPTARG}" ;; # # Depth file with refBottomDepth for MPAS ocean

dst_fl=?* | destination_file=?* | tpl=?* | tpl_fl=?* | template=?* | template_file=?* ) dst_fl="${LONG_OPTARG}" ;; # -d # Destination file

grd_dst=?* | grid_dest=?* | dst_grd=?* | dest_grid=?* | destination_grid=?* ) grd_dst="${LONG_OPTARG}" ;; # -g # Destination grid-file

grd_sng=?* | grid_generation=?* | grid_gen=?* | grid_string=?* ) grd_sng="${LONG_OPTARG}" ;; # -G # Grid generation string

drc_in=?* | in_drc=?* | dir_in=?* | in_dir=?* | input=?* ) drc_in="${LONG_OPTARG}" ;; # -i # Input directory

dt_sng=?* | date_string=?* ) dt_sng="${LONG_OPTARG}" ;; # # Date string for MWF map names

esmf_typ=?* | esmf_mth=?* | esmf_extrap_type=?* | esmf_extrap_method=?* ) esmf_typ="${LONG_OPTARG}" ;; # # ESMF extrapolation type (nearestidavg|neareststod|none)

fl_fmt=?* | fmt_out=?* | file_format=?* | format_out=?* ) fl_fmt="${LONG_OPTARG}" ;; # # Output file format

hrd_pth=?* | hard_path=?* | csz_exe=?* | csz_bin_lib=?* ) echo "No argument allowed for --${OPTARG switch}" >&2; exit 1 ;; # # Use hard-coded paths on known machines

in_fl=?* | in_file=?* | input_file=?* ) in_fl="${LONG_OPTARG}" ;; # -i # Input file

job_nbr=?* | job_number=?* | jobs=?* ) job_usr="${LONG_OPTARG}" ;; # -j # Job simultaneity

map_fl=?* | map=?* | map_file=?* | rgr_map=?* | regrid_map=?* ) map_fl="${LONG_OPTARG}" ;; # -m # Map-file

mem_mb=?* ) mem_mb="${LONG_OPTARG}" ;; # # Megabytes of RAM per srun job in Cori SLURM in MPI mode

mlt_map | multimap | no_multimap | nomultimap ) mlt_map_flg='No' ;; # -M # Multi-map flag

mlt_map=?* | multimap=?* | no_multimap=?* | nomultimap=?* ) echo "No argument allowed for --${OPTARG switch}" >&2; exit 1 ;; # -M # Multi-map flag

mpi_nbr=?* | mpi_number=?* | tsk_nbr=?* | task_number=?* ) mpi_nbr_usr="${LONG_OPTARG}" ;; # # Number of tasks-per-node that MPI-enabled weight generators will request

mpi_pfx=?* | mpi_prefix=?* | srun_cmd=?* | srun_command=?* ) mpi_pfx_usr="${LONG_OPTARG}" ;; # # Prefix for MPI-enabled weight generators

msh_fl=?* | msh=?* | mesh=?* | mesh_file=?* ) msh_fl="${LONG_OPTARG}" ;; # # Mesh file

msk_apl | mask_apply | msk_app ) msk_apl='Yes' ;; # # Apply msk_out to variables after regridding

msk_apl=?* | mask_apply=?* | msk_app=?* ) echo "No argument allowed for --${OPTARG switch}" >&2; exit 1 ;; # # Apply msk_out to variables after regridding

msk_dst=?* | dst_msk=?* | mask_destination=?* | mask_dst=?* ) msk_dst="${LONG_OPTARG}" ;; # # Mask-template variable in destination file

msk_out=?* | out_msk=?* | mask_output=?* | mask_out=?* ) msk_out="${LONG_OPTARG}" ;; # # Put mask variable in regridded file, argument is name

msk_src=?* | src_msk=?* | mask_source=?* | mask_src=?* ) msk_src="${LONG_OPTARG}" ;; # # Mask-template variable in source file

mss_val=?* | fll_val=?* | missing_value=?* | fill_value=?* ) mss_val_usr="${LONG_OPTARG}" ;; # # Missing value for MPAS

nco_opt=?* | nco=?* | nco_options=?* ) nco_opt="${LONG_OPTARG} ${nco_opt}" ;; # -n # NCO options

nm_dst=?* | name_dst=?* | nm_sht_dst=?* | short_name_destination=?* ) nm_dst="${LONG_OPTARG}" ;; # # Short name of destination grid

nm_src=?* | name_src=?* | nm_sht_src=?* | short_name_source=?* ) nm_src="${LONG_OPTARG}" ;; # # Short name of source grid

no_area | no_cll_msr | no_cell_measures ) no_cll_msr='Yes' ;; # # Omit cell_measures variables

no_area=?* | no_cell_msr=?* | no_cell_measures=?* ) echo "No argument allowed for --${OPTARG switch}" >&2; exit 1 ;; # # Omit cell_measures variables

no_frm_trm | no_frm | no_formula_terms ) no_frm_trm='Yes' ;; # # Omit formula_terms variables

no_frm_trm=?* | no_frm=?* | no_formula_terms=?* ) echo "No argument allowed for --${OPTARG switch}" >&2; exit 1 ;; # # Omit formula_terms variables

no_permute | no_prm | no_pdq | no_ncpdq ) pdq_flg='No' ;; # # Do not permute dimensions before regridding

no_permute=?* | no_prm=?* | no_pdq=?* | no_ncpdq=?* ) echo "No argument allowed for --${OPTARG switch}" >&2; exit 1 ;; # # Do not permute dimensions before regridding

nsx_alg=?* | intersection_algorithm=?* ) nsx_alg="${LONG_OPTARG}" ;; # # Intersection algorithm used by mbtempest

drc_out=?* | out_drc=?* | dir_out=?* | out_dir=?* | output=?* ) drc_usr="${LONG_OPTARG}" ;; # -O # Output directory

out_fl=?* | output_file=?* | out_file=?* ) out_fl="${LONG_OPTARG}" ;; # -o # Output file

prc_typ=?* | prm_typ=?* | procedure=?* ) prc_typ="${LONG_OPTARG}" ;; # -P # Procedure type

par_typ=?* | par_md=?* | parallel_type=?* | parallel_mode=?* | parallel=?* ) par_typ="${LONG_OPTARG}" ;; # -p # Parallelism type

pdq=?* | prm=?* | prm_opt=?* | permute=?* ) pdq_opt="-a ${LONG_OPTARG}" ;; # # ncpdq dimension permutation option

ppc=?* | ppc_prc=?* | precision=?* | quantize=?* ) ppc_prc="${LONG_OPTARG}" ;; # # Precision-preserving compression

prs_stt=?* | preserve=?* | preserved_statistic=?* ) prs_stt_usr="${LONG_OPTARG}" ;; # # Preserved statistic

prt=?* | prt_nbr=?* | partition_number=?* ) prt_nbr="${LONG_OPTARG}" ;; # # Partition number for MOAB parallelism

qea | qea_flg | np2 | quasi-equal-area ) qea_flg='Yes' ;; # # Quasi-equal area mode for TempestRemap

qea=?* | qea_flg=?* | np2=?* | quasi-equal-area=?* ) echo "No argument allowed for --${OPTARG switch}" >&2; exit 1 ;; # # Quasi-equal area mode for TempestRemap

rgr_opt=?* | regrid_options=?* ) rgr_opt="${LONG_OPTARG}" ;; # -R # Regridding options

rnr_thr=?* | thr_rnr=?* | rnr=?* | renormalize=?* | renormalization_threshold=?* ) rnr_thr="${LONG_OPTARG}" ;; # -r # Renormalization threshold

rgn_dst=?* | dst_rgn=?* | regional_destination=?* ) hnt_dst='--dst_regional' ;; # # Regional destination grid

rgn_dst=?* | dst_rgn=?* | regional_destination=?* ) echo "No argument allowed for --${OPTARG switch}" >&2; exit 1 ;; # # Regional destination grid

rgn_src=?* | src_rgn=?* | regional_source=?* ) hnt_src='--src_regional' ;; # # Regional source grid

rgn_src=?* | src_rgn=?* | regional_source=?* ) echo "No argument allowed for --${OPTARG switch}" >&2; exit 1 ;; # # Regional source grid

rrg_bb_wesn=?* | bb_wesn=?* | bb=?* | bounding_box=?* ) bb_wesn="${LONG_OPTARG}" ; prc_typ='rrg' ; ;; # # Regional regridding bounding-box WESN order

rrg_dat_glb=?* | dat_glb=?* | data_global=?* | global_data=?* ) dat_glb="${LONG_OPTARG}" ; prc_typ='rrg' ; ;; # # Regional regridding global data file

rrg_grd_glb=?* | grd_glb=?* | grid_global=?* | global_grid=?* ) grd_glb="${LONG_OPTARG}" ; prc_typ='rrg' ; ;; # # Regional regridding global grid file

rrg_grd_rgn=?* | grd_rgn=?* | grid_regional=?* | regional_grid=?* ) grd_rgn="${LONG_OPTARG}" ; prc_typ='rrg' ; ;; # # Regional regridding regional grid file

rrg_rnm_sng=?* | rnm_sng=?* | rename_string=?* ) rnm_sng="${LONG_OPTARG}" ; prc_typ='rrg' ; ;; # # Regional regridding rename string