#!/bin/bash
# run_get_ModisData.bsh
# get AQUA data for DUST input data
. ~/.bash_profile

dirpath=$(dirname $0)

#local_dir="/home/oper/Dust/OSPO_Package_Dust_GRIB/scripts"

. $dirpath/run_Dust_script.bsh

#pre_indir="/data/Dust/input"
#export pre_indir

echo 'run_get_ModisData.bsh started at ' $(date)

YR=$1
jd=$2

if [ ! $YR ]
then

YR=(`(date -d '1 days ago' '+%Y')`)
pwd
jd=(`(date -d '1 days ago' '+%j')`)
fi

echo 'YR jd ' $YR $jd

satid=(MYD)
skipftp=0
skipdo=0
skipplot=0
A=0
###
#while [ $A -le 0 ]

#do
  YJD=${YR[A]}${jd[A]}
  hour=16
  while [ $hour -le 22 ]
  do
    minute=0
    while [ $minute -le 55 ]
    do
      if test $minute -lt 10
      then
        granule=${hour}0$minute
      else
        granule=$hour$minute
      fi
#      echo $granule

      if test $skipftp == 0
      then
        cd $pre_indir
#        echo $pre_indir
#        rm -fr *.hdf
#       for MODNAME in 021KM 03 35_L2
#       for MODNAME in 021KM
#        for MODNAME in 03
#        do
          hdfname1=MYD021KM.P${YR[A]}${jd[A]}.$granule.hdf
          hdfname2=MYD03.P${YR[A]}${jd[A]}.$granule.hdf
          hdfname3=MOD35_L2.P${YR[A]}${jd[A]}.$granule.hdf
          hdfname4=MOD04_L2.P${YR[A]}${jd[A]}.$granule.hdf
          echo  $hdfname1  $hdfname2  $hdfname3
#scp oper@10.144.3.76:/data3p/modis/pub/modis/MOD021KM/$hdfname $pre_indir
sftp 10.144.3.76 <<ENDPOINT
cd /data3p/modis/pub/modis/MYD021KM
lcd $pre_indir
get /data3p/modis/pub/modis/MYD021KM/$hdfname1
cd /data3p/modis/pub/modis/MYD03
lcd $pre_indir
get  /data3p/modis/pub/modis/MYD03/$hdfname2
cd /data3p/modis/pub/modis/MOD35_L2
lcd $pre_indir
get /data3p/modis/pub/modis/MOD35_L2/$hdfname3
cd /data3p/modis/pub/modis/MOD04_L2
lcd $pre_indir/MYD04_L2
get /data3p/modis/pub/modis/MOD04_L2/$hdfname4
bye
ENDPOINT
cd $pre_indir
mv MYD021KM.P${YR[A]}${jd[A]}.$granule.hdf MOD021KM.P${YR[A]}${jd[A]}.$granule.hdf
mv MYD03.P${YR[A]}${jd[A]}.$granule.hdf MOD03.P${YR[A]}${jd[A]}.$granule.hdf
        
#        done
      fi
#      inc_minute
       minute=`expr $minute + 5`
     done
#      inc_hour
      hour=`expr $hour + 1`
    done
#      inc_A
#     A=`expr $A + 1`
# done