#!/bin/sh
#####################################################################
# This script obtains the AFWA, CLAVR, GFS and HIRESW cloud data
# and interpolate them onto 212 grid as for the NAM grid
# History:   Binbin Zhou  -- May 2010, original version
#            Julia Zhu  -- June 30th, second version
# Usage: verf_g2g_get_cloud.sh DATA_TYPE(MODEL NAME)
#####################################################################
set -x

model_name=$1
vday=$2

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

if [ $model_name = afwa ]; then    # AFWA data
  AFWADIR=${AFWADIR:-/dcom/us007003}
  OUTDIR=$COMOUT
  for HH in 00 06 12 18
  do
    $copygb -B-1 -K-1 -x -g212 $AFWADIR/$vday/wgrbbul/cloud/grbclouds.${vday}${HH}00 $COMOUT/afwa.t${HH}z.grd212.f00
    echo "afwa.t${HH}z.grd212.f00 done"
  done

elif [ $model_name = clavr ]; then   # CLAVR data
  CLAVRDIR=${CLAVRDIR:-/dcom/us007003}
  for HH in 00 06 12 18
  do
    $copygb -B-1 -K-1  -x -g212 -i0 $CLAVRDIR//$vday/wgrbbul/clavrx/clvrxgrb.${vday}${HH} $COMOUT/clavr.t${HH}z.grd212.f00
     echo "clavr.t${HH}z.grd212.f00 done"
  done

elif [ $model_name = gfs ]; then   # GFS data
  COMGFS=${COMGFS:-/com/gfs/prod}
  for HH in 00 06 12 18
  do
    fhr=06
    while [ $fhr -le 84 ]
    do
       $cnvgrib -g21 $COMGFS/gfs.$vday/gfs.t${HH}z.pgrb2f${fhr} gfs.t${HH}z.pgrb1f${fhr}
       $copygb -B-1 -K-1  -x -g212 gfs.t${HH}z.pgrb1f${fhr} $COMOUT/gfs.t${HH}z.grd212.f${fhr} 
    
       rm -f gfs.t${HH}z.pgrb1f${fhr}
 
       fhr=`expr $fhr + 6`
       if [ $fhr -lt 10 ]; then fhr=0$fhr; fi
    done
  done
  
  echo "copying of GFS done"

elif [ $model_name = hiresw ]; then   # HIRESW data
  COMHIRESW=${COMHIRESW:-/com/hiresw/prod}
  # EAST region only runs iat 00/12Z
  for HH in 00 12
  do
    fhr=03
    while [ $fhr -le 48 ]
    do
      file1=$COMHIRESW/hiresw.$vday/eastnmm.t${HH}z.awpreg${fhr}.tm00
      $wgrib $file1 |grep "kpds5=71:kpds6=200:kpds7=0"|$wgrib  -i -grib $file1 -o retrived1
      $copygb -i3 -x -g212 retrived1 $COMOUT/eastnmm.t${HH}z.grd212.f${fhr}
      rm -rf retrived1

      file2=$COMHIRESW/hiresw.$vday/eastarw.t${HH}z.awpreg${fhr}.tm00
      $wgrib $file2 |grep "kpds5=71:kpds6=200:kpds7=0"|$wgrib  -i -grib $file2 -o retrived2
      $copygb -i3 -x -g212 retrived2 $COMOUT/eastarw.t${HH}z.grd212.f${fhr}
      rm -rf retrived2

      fhr=`expr $fhr + 3`
      if [ $fhr -lt 10 ]; then fhr=0$fhr; fi
    done
  done

  # WEST region only runs at 06Z
  for HH in 06
  do
    fhr=03
    while [ $fhr -le 48 ]
    do
      file3=$COMHIRESW/hiresw.$vday/westnmm.t${HH}z.awpreg${fhr}.tm00
      $wgrib $file3 |grep "kpds5=71:kpds6=200:kpds7=0"|$wgrib  -i -grib $file3 -o retrived3
      $copygb -i3 -x -g212 retrived3 $COMOUT/westnmm.t${HH}z.grd212.f${fhr}
      rm -rf retrivedwest3

      file4=$COMHIRESW/hiresw.$vday/westarw.t${HH}z.awpreg${fhr}.tm00
      $wgrib $file4 |grep "kpds5=71:kpds6=200:kpds7=0"|$wgrib  -i -grib $file4 -o retrived4
      $copygb -i3 -x -g212 retrived4 $COMOUT/westarw.t${HH}z.grd212.f${fhr}
      rm -rf retrived4

      fhr=`expr $fhr + 3`
      if [ $fhr -lt 10 ]; then fhr=0$fhr; fi
    
    done
  done
fi