UM/MesoscaleModelling: get_dump

#!/bin/bash
#
# W. McGinty
# 4 Apr 2018
#
# Unified script for getting start dumps
#
set -u

# Check if MASS is up
moo ls moose:/opfc > /dev/null 2>&1
RC=$?
if [[ $RC -ne 0 ]] ; then
 echo "MASS is down" 2>&1
 exit 1
fi

if [[ $# -ne 1 ]]; then
  echo "Usage: get_dump <date>" 2>&1
  echo "dates in the form yyyymmddhh" 2>&1
  exit 1
fi

startdate=$1
if [[ ${#startdate} -ne 10 ]]; then
 echo "Start date must be 10 chars e.g. 2011050112" 2>&1
 exit 1
fi

# Split the date
  YRMON=${startdate:0:6}
  Date=${startdate:0:8}
  hr=${startdate:8:2}

# Decide which extraction method to use 

if   [[ $startdate > 2015082606 ]] ; then
  moo get moose:/opfc/atm/global/rerun/${YRMON}.file/${Date}T${hr}00Z_glm_t+0 .
elif [[ $startdate > 2014020406 ]] ; then
  moo get moose:/opfc/atm/global/rerun/${YRMON}.file/${Date}${hr}_glm_t+0 .
elif [[ $startdate > 2009123118 ]] ; then
  moo get moose:/opfc/atm/global/rerun/${YRMON}.file/${Date}_qwqg${hr}.T+0 .
elif [[ $startdate > 2009062006 ]] ; then
  moo get moose:/opfc/atm/global/rerun/${YRMON}.file/${Date}_qwqu${hr}.T+3 .
elif [[ $startdate > 2006010100 ]] ; then
# It's really old, in packed archive file format ...

# Force hour to 06 as that's all that's available
  hr=06
  PAXFILE=coprr.udQU${hr}.${Date}.pax

  moo get moose:/opfc/atm/global/rerun/${YRMON}.file/${PAXFILE} .

  if [[ $? -ne 0 ]]; then
   echo "No such file for date $startdate" 2>&1
   exit 1
  fi

  pax -rf $PAXFILE                     # unpack it

  if [[ -e op ]]; then
    ROOTDIR=op/daily/datawgl
  else
    ROOTDIR=opdaily/datawgl
  fi

  for j in $ROOTDIR/*
  do
    fn=${j##*/}                        # remove the directory part
    mv $ROOTDIR/$fn ${Date}_$fn # Stamp it with the date
  done

  rm $PAXFILE
  rm -rf ${ROOTDIR%%/*}
else
  echo "This is prehistoric.  Earliest Met Office date is 2006-01-01"
  echo "You'll need to get a GRIB start dump from ECMWF"
fi