#!/bin/ksh
#############################################################################
# script: verf_g2g_mdl2wrf_grid.sh $model
#         to convert NAM, RAP and SREF-NMM/ARW grids to Hi-res WRF grid
# Author: Binbin Zhou /SAIC
#         March 25, 2010
#############################################################################
set -x

modnam=$1

export vday=${vday:-$PDYm1}
export vdate=${vdate:-$vday$cyc}


export wgrb=${wgrb:-/nwprod/util/exec/wgrib}
export cpgb=${cpgb:-/nwprod/util/exec/copygb}
export cnvgrib=${cnvgrib:-/nwprod/util/exec/cnvgrib}

west='255 3 884 614 24500 -129200 8 -108000 5000 5000 0 64 0 40500 40500 0 0'   #west
east='255 3 884 614 22100 -109800 8  -89000 5000 5000 0 64 0 38000 38000 0 0'   #east

#############################################
#1:Get convert mosaic to Hires-WRF grid 
#############################################
if [ $modnam = mosaic ]; then
  COMMOSAIC=${COMMOSAIC:-/com/hourly/prod/radar}
  cycles="00 03 06 09 12 15 18 21"
  for cyc in $cycles
  do
    mosaic=$COMMOSAIC.$vday/refd3d.t${cyc}z.grbf00
    $cpgb -g"$west" -i2,1 -x $mosaic $COMOUT/refd3d.t${cyc}z.westgrid.f00
    $cpgb -g"$east" -i2,1 -x $mosaic $COMOUT/refd3d.t${cyc}z.eastgrid.f00
    echo 'copygb mosaic ' $cyc ' done!'
  done
fi

########################################
#2:Get convert NAM to Hires-WRF grid
########################################
if [ $modnam = nam ]; then
  COMNAM=${COMNAM:-/com/nam/prod/nam}
  cycles="00 06 12 18"
  for cyc in  $cycles
  do
    for fhr in 06 12 18 24 30 36 42 48 54 60 66 72 78 84
    do
      nam=$COMNAM.$vday/nam.t${cyc}z.awphys${fhr}.tm00
      $wgrb $nam |grep "kpds5=212:kpds6=200" |$wgrb -i -grib $nam -o $DATA/temp.${cyc}.${fhr}
      $cpgb -g"$west" -i2,1 -x  $DATA/temp.${cyc}.${fhr} $COMOUT/nam.t${cyc}z.westgrid.f${fhr}
      $cpgb -g"$east" -i2,1 -x  $DATA/temp.${cyc}.${fhr} $COMOUT/nam.t${cyc}z.eastgrid.f${fhr}
      echo 'copygb nam ' $cyc $fhr ' done!'
    done
  done
  rm -f $DATA/temp*
fi

#######################################
#3:Get convert RAP to Hires-WRF grid
#######################################
if [ $modnam = rap ]; then
  COMRAP=${COMRAP:-/com/rap/prod/rap}
  cycles="00 06 12 18"
  for cyc in $cycles
  do
    for fhr in 06 12 18
    do 
      rap2=$COMRAP.$vday/rap.t${cyc}z.awp130pgrbf${fhr}.grib2
      rap1=$DATA/rap.t${cyc}z.pgrb13.f${fhr}
      $cnvgrib -g21 $rap2 $rap1 
      $wgrb $rap1 |grep "kpds5=212:kpds6=200" |$wgrb -i -grib $rap1 -o $DATA/temp.${cyc}.${fhr}
      $cpgb -g"$west" -i2,1 -x  $DATA/temp.${cyc}.${fhr} $COMOUT/rap.t${cyc}z.westgrid.f${fhr}
      $cpgb -g"$east" -i2,1 -x  $DATA/temp.${cyc}.${fhr} $COMOUT//rap.t${cyc}z.eastgrid.f${fhr}
      echo 'copygb rap ' $cyc $fhr ' done!'
    done
  done
  rm -f $DATA/temp*
  rm -f $DATA/*pgrb13*
fi

#####################################################
#4: Get convert SREF nmm control grid
#####################################################
if [ $modnam = srefnmm ]; then 
  COMSREF=${COMSREF:-/com/sref/prod/sref}
  cycles="03 09 15 21"
  mdl=nmm
  for cyc in $cycles
  do
    for fhr in 06 12 18 24 30 36 42 48 54 60 66 72 78 84
    do
      srefmdl=$COMSREF.$vday/$cyc/pgrb/sref_${mdl}.t${cyc}z.pgrb212.ctl.f${fhr}
      $wgrb $srefmdl |grep "kpds5=212:kpds6=200" |$wgrb -i -grib $srefmdl -o $DATA/temp.${cyc}.${fhr}
      $cpgb -g"$west" -i2,1 -x $DATA/temp.${cyc}.${fhr} $COMOUT/${modnam}.t${cyc}z.westgrid.f${fhr}
      $cpgb -g"$east" -i2,1 -x $DATA/temp.${cyc}.${fhr} $COMOUT/${modnam}.t${cyc}z.eastgrid.f${fhr}
      echo 'copygb srefnmm ' $cyc $fhr ' done!'
    done
  done
  rm -f $DATA/temp*
  rm -f $DATA/temp*
fi

#####################################################
#4: Get convert SREF arw control grid
#####################################################
if [ $modnam = srefarw ]; then 
  COMSREF=${COMSREF:-/com/sref/prod/sref}
  cycles="03 09 15 21"
  mdl=em
  for cyc in $cycles
  do
    for fhr in 06 12 18 24 30 36 42 48 54 60 66 72 78 84
    do
      srefmdl=$COMSREF.$vday/$cyc/pgrb/sref_${mdl}.t${cyc}z.pgrb212.ctl.f${fhr}
      $wgrb $srefmdl |grep "kpds5=212:kpds6=200" |$wgrb -i -grib $srefmdl -o $DATA/temp.${cyc}.${fhr}
      $cpgb -g"$west" -i2,1 -x $DATA/temp.${cyc}.${fhr} $COMOUT/${modnam}.t${cyc}z.westgrid.f${fhr}
      $cpgb -g"$east" -i2,1 -x $DATA/temp.${cyc}.${fhr} $COMOUT/${modnam}.t${cyc}z.eastgrid.f${fhr}
      echo 'copygb srefarw ' $cyc $fhr ' done!'
    done
  done
  rm -f $DATA/temp*
  rm -f $DATA/temp*
fi

exit