#!/bin/ksh
set -x

#module load libjpeg/9c
#module load grib_util/1.2.2

module purge
module load envvar/1.0
module load intel/19.1.3.304
module load PrgEnv-gnu/8.1.0
#module load ips/18.0.1.163
module load prod_util/2.0.9
module load python/3.8.6
module load wgrib2/2.0.8
module load libpng/1.6.37
module load g2/3.4.1
module load netcdf/3.6.3

module use -a /u/Jun.Du/modulefiles
module load python3/test

if [ $# -eq 1 ]; then
  date0=$1
else                      
  date0=`date -u +%Y%m%d%H`
  minute=`date +%M`
  if [ $minute -lt 33 ]; then
    date0=`$NDATE -1 $date0`
  fi
fi
#. ~/dot.for.python
#source  ~/dot.for.bashrc

hr0=`echo $date0 | cut -c 9-10`
export HOMEurma=/lfs/h2/emc/vpppg/noscrub/Jun.Du/pcpanl/urma.v2.9.0 
export COMINprod=/lfs/h2/emc/ptmp/Jun.Du/pcpanl
export COMINpara=/lfs/h2/emc/ptmp/Jun.Du/pcpanl
export COMINpcpanl=/lfs/h2/emc/ptmp/Jun.Du/pcpanl
#export COMOUT=/gpfs/dell2/ptmp/Jun.Du/com/urma/prod
export COMOUT=/lfs/h2/emc/vpppg/noscrub/Jun.Du/pcpurma_plots
export DATA=/lfs/h2/emc/ptmp/Jun.Du/pcpanl
export RUN=pcpurma

cd $DATA
mkdir pyplot_prod
cd pyplot_prod

day0=`echo $date0 | cut -c 1-8`
hr0=`echo $date0 | cut -c 9-10`

mkdir -p $COMOUT/$RUN.$day0

grep 06h $COMINpcpanl/pcpanl.$day0/todo_urma.$date0 | grep conus > $COMOUT/$RUN.$day0/todo_urma.$date0

todourma=$COMOUT/$RUN.$day0/todo_urma.$date0

######### processing pcpurma files #################

for item in `cat $todourma`
do
  date=`echo $item | cut -c 1-10`
  day=`echo $item | cut -c 1-8`
  acc=`echo $item | cut -c 12-14`
  region=`echo $item | awk -F"." '{print $3}'`
  if [ $region = conus ]; then
    regnam=wexp
  else
    regnam=$region
  fi
 
  mkdir -p $COMOUT/$RUN.$day
  
  if [ $region = conus ]; then
    urmafile=pcpurma_${regnam}.${date}.$acc.grb2
    cp $COMINprod/$RUN.$day/$urmafile .
    urmapng=`echo $urmafile | sed -s 's/grb2/png/'`
    python $HOMEurma/util.dev/plt_pcpurma.py $urmafile $region
    cp $urmapng $COMOUT/$RUN.$day/.

    urmafile_para=pcpurma2p8_${regnam}.${date}.$acc.grb2
    cp $COMINpara/$RUN.$day/$urmafile $urmafile_para
    urmapng_para=`echo $urmafile_para | sed -s 's/grb2/png/'`
    python $HOMEurma/util.dev/plt_pcpurma.py $urmafile_para $region
    cp $urmapng_para $COMOUT/$RUN.$day/.

    urmamask=pcpurma_mask.${date}.$acc.grb2
    if [ -s $COMINprod/$RUN.$day/$urmamask ]; then
      cp $COMINprod/$RUN.$day/$urmamask .
      python $HOMEurma/util.dev/plt_urmamask.py $urmamask
      maskpng=`echo $urmamask | sed -s 's/grb2/png/'`
      cp $maskpng $COMOUT/$RUN.$day/.
    fi 
  fi

  RZDMDIR=/home/people/emc/www/htdocs/mmb/SREF/pcpanl/$day/URMA/$RUN.$day
  ssh wd20jd@emcrzdm "mkdir -p $RZDMDIR/plot"
  
  montage -geometry 480x425 $COMOUT/${RUN}.$day/$urmapng $COMOUT/pcpurma.$day/$urmapng_para $COMOUT/${RUN}.$day/pcpurma_comp.${date}.$acc.png

  scp $COMOUT/${RUN}.$day/pcpurma_comp.${date}.$acc.png wd20jd@emcrzdm:$RZDMDIR/plot/.
  scp $COMOUT/${RUN}.$day/$maskpng wd20jd@emcrzdm:$RZDMDIR/plot/.
  scp $COMOUT/${RUN}.$day/$urmapng wd20jd@emcrzdm:$RZDMDIR/plot/.
  scp $COMOUT/${RUN}.$day/$urmapng_para wd20jd@emcrzdm:$RZDMDIR/plot/.
  
done

exit