linux scripts

• Get vertical temperature at several nodes

set -e

nodes=(90 113 2243 2554 2867 6600)

for j in "${nodes[@]}"; do
echo "Processing node ${j}"
  for i in {1..366}; do
    printf -v nn "%04d" "$i"

    ncks -C -v h,art1,temp,siglay,Times \
         -F -d node,"$j" \
         "superior_${nn}.nc" \
         "tmp_${nn}.nc"
  done

  ncrcat tmp_????.nc "gls_ver_temp_2018_LowRes_${j}.nc"
  rm -f tmp_????.nc
done