#!/bin/bash

set -ax 

export PDY=$1
export cyc=$2

machine=`hostname | cut -c1`

source ../versions/run.ver

#module use ${PACKAGEROOT}/${model}.${aqm_ver}/modulefiles

 module use $(realpath ../modulefiles)


module load aqm/6.1.1

#TODAY=`date +%Y%m%d`

export SENDDBN=NO

export Yr=`echo ${PDY} | cut -c1-4`
export Mn=`echo ${PDY} | cut -c5-6`
export Yrmn=`echo ${PDY}  | cut -c1-6`
export Yr1=`echo ${PDYm1} | cut -c1-4`
export Yrmn1=`echo ${PDYm1} | cut -c1-6`
export Yr2=`echo ${PDYm2} | cut -c1-4`
export Yrmn2=`echo ${PDYm2} | cut -c1-6`

export stdm1=`$NDATE -6 ${PDY}${cyc}`
export PDY1=$(echo $stdm1 | cut -c 1-8)
export cyc1=$(echo $stdm1 | cut -c 9-10)
export stdm2=`$NDATE -12 ${PDY}${cyc}`
export PDY2=$(echo $stdm2 | cut -c 1-8)
export cyc2=$(echo $stdm2 | cut -c 9-10)
export stdm3=`$NDATE -18 ${PDY}${cyc}`
export PDY3=$(echo $stdm3 | cut -c 1-8)
export cyc3=$(echo $stdm3 | cut -c 9-10)
export stdm4=`$NDATE -24 ${PDY}${cyc}`
export PDY4=$(echo $stdm4 | cut -c 1-8)
export cyc4=$(echo $stdm4 | cut -c 9-10)

BASE=`pwd`
export HOMEaqm="$(dirname ${BASE})"
#export case=aqm.v7.0.c18
export usr_tmp=/lfs/h2/emc/ptmp/kai.wang/emc.para/com/aqm/v7.0/$case
export DATAROOT=${usr_tmp}/tmpnwprd_p2
export job=aqm_pm25_bias_cor
export DATA=${DATAROOT}/${job}.$$
export COMINbicordat=/lfs/h2/emc/physics/noscrub/kai.wang/Bias_correction/aqmv7.0
export COMINbicor=${COMINbicordat}/bcdata.${Yrmn}
export COMOUTbicordat=/lfs/h2/emc/physics/noscrub/kai.wang/Bias_correction/aqmv7.0
export COMOUTbicor=${COMINbicordat}/bcdata.${Yrmn}

export COMIN=/lfs/h2/emc/ptmp/kai.wang/emc.para/com/aqm/v7.0/$case.${PDY}/${cyc}
export COMINm1=/lfs/h2/emc/ptmp/kai.wang/emc.para/com/aqm/v7.0/$case.${PDY1}/${cyc1}
export COMINm2=/lfs/h2/emc/ptmp/kai.wang/emc.para/com/aqm/v7.0/$case.${PDY2}/${cyc2}
export COMINm3=/lfs/h2/emc/ptmp/kai.wang/emc.para/com/aqm/v7.0/$case.${PDY3}/${cyc3}
export COMINm4=/lfs/h2/emc/ptmp/kai.wang/emc.para/com/aqm/v7.0/$case.${PDY4}/${cyc4}
export COMOUT=/lfs/h2/emc/ptmp/kai.wang/emc.para/com/aqm/v7.0/$case.${PDY}/${cyc}

#export COMINbicor=/lfs/h2/emc/physics/noscrub/Jianping.Huang/prod/com/aqm/v6.1/bcdata.${Yrmn}
#export COMINbicorm1=/lfs/h2/emc/physics/noscrub/Jianping.Huang/prod/com/aqm/v6.1/bcdata.${Yrmn1}
#export COMINbicorm2=/lfs/h2/emc/physics/noscrub/Jianping.Huang/prod/com/aqm/v6.1/bcdata.${Yrmn2}

qsub  ${HOMEaqm}/working/jaqm_post_pm25_bicor_793.ecf

err=$?

if [ ${err} -ne 0 ]
then
  echo "FATAL - cmaq post bias correction cs failed with error code ${err}"
  exit 101
fi
echo "$0"
pwd
#