#!/bin/bash
set -aeux

TMP_DIR=/scratchin/grupos/grpeta/projetos/tempo/oper/GEF_v1.0.0_oper/gef_trunk/PRP/gef_tmp
PRP_DIR=/scratchin/grupos/grpeta/projetos/tempo/oper/GEF_v1.0.0_oper/gef_trunk/PRP
ETA_RUN=/scratchout/grupos/grpeta/projetos/tempo/oper/GEF_v1.0.0_oper/2020082112
FInitBC=gfs0.25
Run_Date=2020082112
LM=${1}

ulimit -c unlimited

gribfile=${PRP_DIR}/INIT_DATA/${Run_Date}.bin
outfile=${ETA_RUN}/gfsdata01.dat
initout=${ETA_RUN}/init_vrbls01.dat
mpiout=${ETA_RUN}/vrbls01.

cd ${TMP_DIR}


${TMP_DIR}/get_gfs0.25_data.ksh ${Run_Date} 
${TMP_DIR}/${FInitBC}_deco.sh ${Run_Date} ${FInitBC}


echo ${Run_Date} | cut -c3-10 > ${PRP_DIR}/src/dprep/dgavn/adate.txt 
arq_grib=${ETA_RUN}/../INIT_DATA/${FInitBC}_${Run_Date}E.0000
echo $arq_grib > ${PRP_DIR}/src/dprep/dgavn/pas


rm -f ${PRP_DIR}/src/dprep/dgavn/InputModelInf.txt
cp ${PRP_DIR}/src/dprep/dgavn/InputModelInf_${FInitBC}.txt ${PRP_DIR}/src/dprep/dgavn/InputModelInf.txt
chmod 755 ${PRP_DIR}/src/dprep/dgavn/InputModelInf.txt

ln -sf ${PRP_DIR}/src/dprep/dgavn/InputModelInf.txt InputModelInf.txt
ln -sf ${PRP_DIR}/src/dprep/dgavn/adate.txt adate.txt


${PRP_DIR}/src/dprep/dgavn/read_${FInitBC}.x < ${PRP_DIR}/src/dprep/dgavn/pas > ${ETA_RUN}/../err/err.read_gfs 2>&1
mv ${TMP_DIR}/fort.15  ${ETA_RUN}/../INIT_DATA/${Run_Date}.bin


ln -sf ${ETA_RUN}/../INIT_DATA/${Run_Date}.bin  ${PRP_DIR}/INIT_DATA/${Run_Date}.bin


cd ${ETA_RUN}

ln -sf ${PRP_DIR}/data_in/static2/maxsnoalb.bin         fort.20
ln -sf ${PRP_DIR}/data_in/static2/alb1_ieee             fort.21
ln -sf ${PRP_DIR}/data_in/static2/alb2_ieee             fort.22
ln -sf ${PRP_DIR}/data_in/static2/alb3_ieee             fort.23
ln -sf ${PRP_DIR}/data_in/static2/alb4_ieee             fort.24
ln -sf ${PRP_DIR}/data_in/static2/ivgtyp_1d_ieee        fort.30
ln -sf ${PRP_DIR}/data_in/static2/isltyp_1d_ieee        fort.31
ln -sf ${PRP_DIR}/data_in/static2/islope_1d_ieee        fort.32
ln -sf ${PRP_DIR}/data_in/static2/rfusaflw_ieee         fort.42
ln -sf ${PRP_DIR}/data_in/static2/imsmask.ascii         fort.43

ln -sf ${PRP_DIR}/data_in/static2/eta_micro_lookup.dat  eta_micro_lookup.dat
ln -sf ${PRP_DIR}/INIT_DATA/feb1998/ICEC_WEASD.bin      ICEC_WEASD.bin
## sem mudanca nivel ln -sf ${PRP_DIR}/data_in/static2/co2.${LM}_25mb        co2.dat
ln -sf ${PRP_DIR}/data_in/static2/co2.${LM}_1mb         co2.dat
ln -sf ${PRP_DIR}/data_in/static2/snowdepth.grb         snowdepth.grb
ln -sf ${PRP_DIR}/data_in/static2/imssnow.grb           imssnow.grb
ln -sf ${PRP_DIR}/data_in/static2/veg.eta.grb           veg.eta.grb
ln -sf ${PRP_DIR}/data_in/grid/gridinit.dat             gridinit.dat
## sem mudanca nivel ln -sf ${PRP_DIR}/data_in/static2/deta_$LM              deta_$LM
ln -sf ${PRP_DIR}/data_in/static2/Deta_${LM}_1hpa       deta_$LM
ln -sf ${PRP_DIR}/data_in/grid/qintc_2.dat              qintc_2.dat
ln -sf ${PRP_DIR}/data_in/grid/grph.dat                 grph.dat
ln -sf ${PRP_DIR}/data_in/grid/topo.dat                 topo.dat
ln -sf ${PRP_DIR}/data_in/grid/zeff.dat                 zeff.dat
ln -sf ${PRP_DIR}/../cnstdata                           cnstdata


${PRP_DIR}/exe/const.exe > ${ETA_RUN}/../err/err.const 2>&1 

echo $gribfile > pasdg
echo $outfile >> pasdg
${PRP_DIR}/exe/dgeta2model_gbl.exe < pasdg > ${ETA_RUN}/../err/err.dgeta2model_gbl 2>&1

echo $outfile > pas
echo $initout >> pas
${PRP_DIR}/exe/init_vrb.exe < pas > ${ETA_RUN}/../err/err.init_vrb 2>&1

echo $initout > pas
echo $mpiout >> pas
${PRP_DIR}/exe/out4mpi.exe < pas > ${ETA_RUN}/../err/err.out4mpi 2>&1

rm -f fort.*

export PBS_SERVER=eslogin13
/opt/pbs/10.4.0.101257/bin/qsub ${TMP_DIR}/run_gef.sh

export PBS_SERVER=aux20-eth4
/opt/pbs/10.4.0.101257/bin/qsub ${TMP_DIR}/submit_post_aux.ksh

exit