convert-CRISPCAv0-L1B_PCS.sh

#!/bin/bash -e
BASE="$( cd -P "$( dirname "${BASH_SOURCE[0]}" )" && pwd )"

l1a_nc_copy_vars() {
  if [[ -d "$CRIS_L1B_HOME" ]]; then
    "$CRIS_L1B_HOME/lib/ShellB3/bin/python3" -m sips_cris_l0l1a.netcdf "$@"
  elif [[ -x "$(which nc_copy_vars 2>/dev/null)" ]]; then
    nc_copy_vars "$@"
  else
    echo "${Oops:?Unable to find nc_copy_vars utility; is CRIS_L1B_HOME set?}"
  fi
}

py_ncrename() {
  if [[ -d "$CRIS_L1B_HOME" ]]; then
    "$CRIS_L1B_HOME/lib/ShellB3/bin/python3" "${BASE}/nc_rename.py" "$@"
  elif [[ -x "$(which python3 2>/dev/null)" ]]; then
    python3 "${BASE}/nc_rename.py" "$@"
  else
    echo "${Oops:?Unable to find python3 with netCDF4 library; is CRIS_L1B_HOME set?}"
  fi
}

usage() {
  cat <<EOF
Create a L1B_PCS file with PC Scores, given an L1B input and a CRISPCAv0 input.

Usage:
convert-CRISPCAv0-L1B_PCS.sh L1B-file CRISPCAv0-file {optional-output-file}

Example:
/path/to/convert-CRISPCAv0-L1B_PCS.sh \\
  input/SNDR.J1.CRIS.20210411T0448.m06.g049.L1B.std.v03_08.G.210411135350.nc \\
  input/CRISPCAv0.J1.CRIS.20210411T0448.m06.g049.L1B.std.v03_08.G.210411135350.nc \\
  SNDR.J1.CRIS.20210411T0448.m06.g049.L1B_PCS.std.v03_08.G.210411135350.nc
EOF
  exit 1
}

default_output_name() {
  local fn="$(basename "$1")"
  echo "${fn/.L1B./.L1B_PCS.}"
}


FN_L1B="${1}"
FN_PCS="${2}"
FN_L1BPCS="${3:-$(default_output_name "$FN_L1B")}"
FN_GUIDE=${BASE}/anc/cris_l1b_pcs_inherit.nc
FN_TEMPLATE=${BASE}/anc/template.nc

if [[ ! -f "$FN_PCS" ]]; then
  usage
  exit 1
fi

cp "$FN_TEMPLATE" "${FN_L1BPCS}-"

# reformat the pcs file
cp -f "$FN_PCS" "${FN_L1BPCS}-pcs-"
py_ncrename \
    v:globalCoef:rad_global_pcs \
    v:localCoef:rad_local_pcs \
    v:U:rad_local_pcs_U \
    v:M:rad_local_pcs_M \
    v:D:rad_local_pcs_D \
    d:nFORs:xtrack \
    d:nFOVs:fov \
    d:nScans:atrack \
    d:gpcs:globalpc \
    d:pcs:localpc \
    d:channels:wnum \
    "${FN_L1BPCS}-pcs-"

test ! -f "$FN_L1BPCS" || rm -f "$FN_L1BPCS"
l1a_nc_copy_vars -vv "${FN_L1BPCS}-" / "$FN_L1B" / "$FN_GUIDE"
l1a_nc_copy_vars -vv "${FN_L1BPCS}-" / "${FN_L1BPCS}-pcs-" / "${FN_L1BPCS}-pcs-"

rm -f "${FN_L1BPCS}-pcs-"
mv "${FN_L1BPCS}-" "$FN_L1BPCS"
echo "> $FN_L1BPCS"