#!/bin/sh
#################################################################
#  SCRIPT TO RETRIEVE FIELDs FROM GFS
#  SAVE IN ASCII FOR FORECAST
#################################################################

cd $DATA

set +x
echo "***************"
echo "* READPGB_6"
echo "***************"
set -x

f3=$1
f6=$2

#cat ${COMIN_GFS}/gfs.t${cycl}z.pgrb2f${f3} ${COMIN_GFS}/gfs.t${cycl}z.pgrb2bf${f3} > grib2f3
#cat ${COMIN_GFS}/gfs.t${cycl}z.pgrb2f${f6} ${COMIN_GFS}/gfs.t${cycl}z.pgrb2bf${f6} > grib2f6
#$utilexec/cnvgrib -g21 grib2f3 gribf3 &
#$utilexec/cnvgrib -g21 grib2f6 gribf6

#wait
   grib_in_3=/global/hires/glopara/com/gfs/para/gfs.${PDY}/gfs.t${cycl}z.master.grbf${f3}
   gribi_in_3=/global/hires/glopara/com/gfs/para/gfs.${PDY}/gfs.t${cycl}z.master.grbif${f3}
   gribf3=/ptmp/wx53cl/gfs.t${cycl}z.uv.grbf${f3}
   grib_in_6=/global/hires/glopara/com/gfs/para/gfs.${PDY}/gfs.t${cycl}z.master.grbf${f6}
   gribi_in_6=/global/hires/glopara/com/gfs/para/gfs.${PDY}/gfs.t${cycl}z.master.grbif${f6}
   gribf6=/ptmp/wx53cl/gfs.t${cycl}z.uv.grbf${f6}

   $utilexec/copygb -i0 -k'4*-1 10 200 0' -g4 $grib_in_3 $gribi_in_3 $gribf3
   $utilexec/copygb -i0 -k'4*-1 11 1 0' -a -g4 $grib_in_3 $gribi_in_3 $gribf3
   $utilexec/copygb -i0 -k'4*-1 84 1 0' -a -g4 $grib_in_3 $gribi_in_3 $gribf3
   $utilexec/copygb -i0 -k'4*-1 200 1 0' -a -g4 $grib_in_3 $gribi_in_3 $gribf3
   $utilexec/copygb -i0 -k'4*-1 201 1 0' -a -g4 $grib_in_3 $gribi_in_3 $gribf3
   $utilexec/copygb -i0 -k'4*-1 204 1 0' -a -g4 $grib_in_3 $gribi_in_3 $gribf3
   $utilexec/copygb -i0 -k'4*-1 211 1 0' -a -g4 $grib_in_3 $gribi_in_3 $gribf3


   $utilexec/copygb -i0 -k'4*-1 10 200 0' -g4 $grib_in_6 $gribi_in_6 $gribf6
   $utilexec/copygb -i0 -k'4*-1 11 1 0' -a -g4 $grib_in_6 $gribi_in_6 $gribf6
   $utilexec/copygb -i0 -k'4*-1 84 1 0' -a -g4 $grib_in_6 $gribi_in_6 $gribf6
   $utilexec/copygb -i0 -k'4*-1 200 1 0' -a -g4 $grib_in_6 $gribi_in_6 $gribf6
   $utilexec/copygb -i0 -k'4*-1 201 1 0' -a -g4 $grib_in_6 $gribi_in_6 $gribf6
   $utilexec/copygb -i0 -k'4*-1 204 1 0' -a -g4 $grib_in_6 $gribi_in_6 $gribf6
   $utilexec/copygb -i0 -k'4*-1 211 1 0' -a -g4 $grib_in_6 $gribi_in_6 $gribf6

echo 'READ_PGB_MASTER: gribf3 fhr=', $f3
echo 'READ_PGB_MASTER: gribf6 fhr=', $f6
#
# Create input files for uv_index_3.sh step
#
$utilexec/wgrib $gribf3  -s | grep "TOZNE" > $DATA/ozone.inv
$utilexec/wgrib $gribf3  -s | grep "ALBDO" > $DATA/albedo.inv
$utilexec/wgrib $gribf3  -s | grep ":DUVB" > $DATA/uvbcld.inv
$utilexec/wgrib $gribf3  -s | grep "CDUVB" > $DATA/uvbclr.inv
$utilexec/wgrib $gribf3  -s | grep "TMP:sfc" > $DATA/tsfc.inv
 
$utilexec/wgrib $gribf3  -i < $DATA/ozone.inv -text -o $DATA/uv_ozone.txt
$utilexec/wgrib $gribf3  -i < $DATA/albedo.inv -text -o $DATA/uv_albedo.txt
$utilexec/wgrib $gribf3  -i < $DATA/uvbcld.inv -text -o $DATA/uv_cloud.txt
$utilexec/wgrib $gribf3  -i < $DATA/uvbclr.inv -text -o $DATA/uv_clear.txt
$utilexec/wgrib $gribf3  -i < $DATA/tsfc.inv -text -o $DATA/uv_tsfc.txt
#
# Create input files for uv_index_6.sh step
#
$utilexec/wgrib $gribf3 -s | grep "DSWRF:sfc" > $DATA/swdown3.inv
$utilexec/wgrib $gribf3 -s | grep "USWRF:sfc" > $DATA/swup3.inv
$utilexec/wgrib $gribf3 -s | grep ":DUVB" > $DATA/uvbcld3.inv
$utilexec/wgrib $gribf3 -s | grep "CDUVB" > $DATA/uvbclr3.inv

$utilexec/wgrib $gribf6 -s | grep "TOZNE" > $DATA/ozone6.inv
$utilexec/wgrib $gribf6 -s | grep "DSWRF:sfc" > $DATA/swdown6.inv
$utilexec/wgrib $gribf6 -s | grep "USWRF:sfc" > $DATA/swup6.inv
$utilexec/wgrib $gribf6 -s | grep ":DUVB" > $DATA/uvbcld6.inv
$utilexec/wgrib $gribf6 -s | grep "CDUVB" > $DATA/uvbclr6.inv
$utilexec/wgrib $gribf6 -s | grep "TMP:sfc" > $DATA/tsfc6.inv

$utilexec/wgrib $gribf3 -i < $DATA/swdown3.inv -text -o $DATA/uv_swdn3.txt
$utilexec/wgrib $gribf3 -i < $DATA/swup3.inv -text -o $DATA/uv_swup3.txt
$utilexec/wgrib $gribf3 -i < $DATA/uvbcld3.inv -text -o $DATA/uv_cloud3.txt
$utilexec/wgrib $gribf3 -i < $DATA/uvbclr3.inv -text -o $DATA/uv_clear3.txt

$utilexec/wgrib $gribf6 -i < $DATA/ozone6.inv -text -o $DATA/uv_ozone6.txt
$utilexec/wgrib $gribf6 -i < $DATA/swdown6.inv -text -o $DATA/uv_swdn6.txt
$utilexec/wgrib $gribf6 -i < $DATA/swup6.inv -text -o $DATA/uv_swup6.txt
$utilexec/wgrib $gribf6 -i < $DATA/uvbcld6.inv -text -o $DATA/uv_cloud6.txt
$utilexec/wgrib $gribf6 -i < $DATA/uvbclr6.inv -text -o $DATA/uv_clear6.txt
$utilexec/wgrib $gribf6 -i < $DATA/tsfc6.inv -text -o $DATA/uv_tsfc6.txt

exit