#! /bin/bash

#============================================================
# <isrc> : base source directory of input SSMI raw data files
#============================================================
isrc=$CMORPH_1/GPROF

#=====================================
# <data> : base directory on $DATA_OUT
#=====================================
dout=$DATA_OUT/observations/land_air/long_range/global/CMAP_PEN_RT
#========================= CF2
#if [ ! -d ${dout} ]; then
#  mkdir -p ${dout}
#fi
#== CF2

#--------------------------------------------------
# <yyyymmdd>    : target day (yesterday)
# <yyyy>        : year of target day
# <mm>          : month of target day
# <yyjjj>       : year and Julian day of target day
# <yyyymmdd_10> : 10 days before target day
#--------------------------------------------------
cdate=$(date -u)
yyyymmdd=$(date -u "+%Y%m%d" --date="-1 days")
yyyy=$(echo ${yyyymmdd} | cut -c1-4)
mm=$(echo ${yyyymmdd} | cut -c5-6)
yyjjj=$(date -u "+%y%j" --date="${yyyymmdd}")
yyyymmdd_10=$(date -u "+%Y%m%d" --date="${yyyymmdd} -10 days")

#--------------------------------------------
# <base> : base directory of CMAP_PEN project
#--------------------------------------------
base=$CMAP_PEN_RT

#----------------------------------------------
# <bin>    : sub-directory of executables
# <input>  : sub-directory of input data files
# <output> : sub-directory of output data files
# <dlog>   : sub-directory of log files
#----------------------------------------------
  prog=${base}/prog
 input=${base}/input
output=${base}/output
  dlog=${base}/log
if [ ! -d ${input} ]; then
  mkdir -p ${input}
fi
if [ ! -d ${output} ]; then
  mkdir -p ${output}
fi
if [ ! -d ${dlog} ]; then
  mkdir -p ${dlog}
fi

#------------------------
# <log>    : log file
# <log_10> : old log file
#------------------------
   ymd=$(date -u "+%Y%m%d" --date="${cdate}")
ymd_rm=$(date -u "+%Y%m%d" --date="${cdate} -30 days")
   log=${dlog}/SSMI.log_${ymd}
log_rm=${dlog}/SSMI.log_${ymd_rm}
rm -f ${log_rm}

#-------------------------------------------------------------------------
# set up output data sub-directories
#-------------------------------------------------------------------------
# <d_ssmi> : sub-directory of SSMI rain & num data files on $DATA_OUT
# <o_ssmi> : sub-directory of SSMI rain & num data files
# <o_ppt>  : sub-directory of SSMI 2.5deg daily & pentad data yearly files
#-------------------------------------------------------------------------
d_ssmi=${dout}/SSMI/${yyyy}/${yyyy}${mm}
o_ssmi=${output}/SSMI/${yyyy}/${yyyy}${mm}
 o_ppt=${output}/SSMI/${yyyy}/ppt
#=========================== CF2
#if [ ! -d ${d_ssmi} ]; then
#  mkdir -p ${d_ssmi}
#fi
 #== CF2
 
if [ ! -d ${o_ssmi} ]; then
  mkdir -p ${o_ssmi}
fi
if [ ! -d ${o_ppt} ]; then
  mkdir -p ${o_ppt}
fi
###############################################################################
#set oo_ssmi=/cpc/prod_tst/cpcsat/CMAP_PEN_RT/PRODUCTS/SSMI/${yyyy}/${yyyy}${mm}
#if( ! -d ${oo_ssmi} ) then
#  mkdir -p ${oo_ssmi}
#endif

#---------------------------------------
# <f_list> : list of SSMI raw data files
# <f_rain> : SSMI rain data file
# <f_num>  : SSMI num data file
#---------------------------------------
f_list=dmsp_raw_file_list_${yyyymmdd}
f_rain=rain_daily_F13F14F15_${yyjjj}.dat
 f_num=num_F13F14F15_${yyjjj}.dat

#---------
# log time
#---------
echo "############################ [Start]" >> ${log}
date -u >> ${log}
echo "----------------------------" >> ${log}
echo >> ${log}

#-------------------------------
# clean up <input> sub-directory
#-------------------------------
cd ${input}
rm -f ${f_list} ${f_rain} ${f_num} DM??-raw1-${yyyymmdd}????

#--------------------------------------------------------------
# copy SSMI raw data files from source to <input> sub-directory
#--------------------------------------------------------------
cd ${input}
for id in 16 17 18; do
  cp -p ${isrc}/DM${id}/work/DM${id}-raw1-${yyyymmdd}???? ${input}
done
ls -1 DM??-raw1-${yyyymmdd}???? > ${f_list}

#------------------------------------
# generate SSMI rain & num data files
#------------------------------------
cd ${prog}
./SSMI_cal_rain_num_from_raw.x ${yyyymmdd} >> ${log}
echo >> ${log}

#------------------------------------------
# update SSMI 2.5deg daily data yearly file
#------------------------------------------
cd ${prog}
./SSMI_def_dly_2.5deg.x ${yyyymmdd} >> ${log}
echo >> ${log}

#-------------------------------------------
# update SSMI 2.5deg pentad data yearly file
#-------------------------------------------
cd ${prog}
./SSMI_def_pen_2.5deg.x ${yyyymmdd} >> ${log}
echo >> ${log}

#-------------------------------
# clean up <input> sub-directory
#-------------------------------
cd ${input}
rm -f ${f_list} DM??-raw1-${yyyymmdd}????

#-----------------------------------------------------------
# copy/move SSMI rain & num data files to $DATA_OUT/<output>
#-----------------------------------------------------------
bzip2 ${input}/${f_rain}
#====================================== CF2
#cp -p ${input}/${f_rain}.bz2 ${d_ssmi}
#== CF2
#######################################
#cp -p ${input}/${f_rain}.bz2 ${oo_ssmi}
mv    ${input}/${f_rain}.bz2 ${o_ssmi}

bzip2 ${input}/${f_num}
#====================================== CF2
#cp -p ${input}/${f_num}.bz2  ${d_ssmi}
#== CF2
#######################################
#cp -p ${input}/${f_num}.bz2  ${oo_ssmi}
mv    ${input}/${f_num}.bz2  ${o_ssmi}

echo "----------------------------" >> ${log}
date -u >> ${log}
echo "############################ [Done]" >> ${log}
echo >> ${log}

exit