#!/bin/bash
set -xv

gribfile=$1
gblllnc=G2N_${gribfile%%.g*}.nc

rm -f $gblllnc

rh=RH

cat << EOF > g2n.nml
&splice
outputunit  = "pres", "hPa"
levinvert   = .TRUE.
badecmwf    = .FALSE.
ensnum      = 0
splicename  = "HGT", "TMP", "UGRD", "VGRD", "$rh", "TMP", "PRMSL", "HGT", "LAND", "PRES", "ICEC"
spliceunit  = "ISBL", "ISBL", "ISBL", "ISBL", "ISBL", "SFC", "MSL", "SFC", "SFC", "SFC", "SFC"
newname     = "hgt", "temp", "u", "v", "rh", "tss", "mslp", "zs", "land", "ps", "fracice"
/
EOF

g2n -f nc -i $gribfile -o $gblllnc < g2n.nml > log.g2n
[ ! -e $gblllnc ] && exit 0

##########################################################################

unittime=`ncdump -c $gblllnc | grep time:unit | cut -d\" -f2`
units=`echo $unittime | cut -d" " -f1`
unitdate=`echo $unittime | cut -d" " -f 3`
kdate_s=`echo $unitdate | cut -c1-10|sed 's/-//g'`
ktime_s=`echo $unittime | cut -d" " -f 4 | sed 's/://g' | cut -c1-4`
unit=$unitdate":"$ktime_s

mm=`echo $unitdate | cut -c6-7`

echo "================================================> mm=$mm"
echo "================================================> kdate_s=$kdate_s"
echo "================================================> ktime_s=$ktime_s"

mkdir -p output
mv -f ${gblllnc} output/
rm -f g2n.nml log.g2n