#!/bin/bash

#-----------------------------------------------------
# Run the sfcupdate program on the WCOSS Dell
#
# Invoke this script as follows:
# "cat $script | bsub"
#
# Do not run with more than one mpi task.
#-----------------------------------------------------

#BSUB -oo log
#BSUB -eo log
#BSUB -M 3000
#BSUB -q debug
#BSUB -P NAM-T2O
#BSUB -J sfcupdate
#BSUB -W 0:05
#BSUB -x                  # not shared
#BSUB -n 1                # total tasks
#BSUB -R "span[ptile=1]"  # tasks per node
#BSUB -R "affinity[core(1):distribute=balance]"

set -x

module purge
module load EnvVars/1.0.2
module load ips/18.0.1.163
module load impi/18.0.1
module load lsf/10.1
module list

# path/name of the sfcupdate executable.
RUN_DIR=${LS_SUBCWD}
EXE_DIR=${RUN_DIR}/../../exec
EXE=$EXE_DIR/emcsfc_sfcupdate

WORK_DIR="/gpfs/dell1/stmp/$LOGNAME/sfcupdate/"
rm -fr $WORK_DIR
mkdir -p $WORK_DIR
cd $WORK_DIR

# get nems file from nco's nmmb parallel.
DATA_DIR=/gpfs/dell2/emc/modeling/noscrub/George.Gayno/landutil.git/sfcupdate
ln -fs ${DATA_DIR}/ndas.t06z.input_nemsio.d01.tm12_presfcupdate   ./
cp ndas.t06z.input_nemsio.d01.tm12_presfcupdate ndas.t06z.input_nemsio.d01.tm12_afterupdate
chmod 644  ndas.t06z.input_nemsio.d01.tm12_afterupdate
cp ${DATA_DIR}/seaice.2014103018.nam.grb     ./seaice.2014103018.nam.grb
cp ${DATA_DIR}/snow.2014103018.nam.grb       ./snow.2014103018.nam.grb
cp ${DATA_DIR}/sst.2014103018.nam.grb        ./sst.2014103018.nam.grb

# directory where fixed fields for nmmb grid are located.  includes landmask, albedo, greenness.
FIXED_DIR=$DATA_DIR

cat > ${WORK_DIR}/fort.41 << !
 &grid_info
  domain_name="nam"
 /
 &cycle
  cycle_year=2014
  cycle_month=10
  cycle_day=30
  cycle_hour=18
  fcst_hour=0.0
 /
 &model_flds
  first_guess_file="${WORK_DIR}/ndas.t06z.input_nemsio.d01.tm12_presfcupdate"
  first_guess_file_type="nems"
 /
 &fixed_flds
  lsmask_file="${FIXED_DIR}/nam_slmask.grb"
 /
 &time_varying_analysis_flds
  snowfree_albedo_climo_file="${FIXED_DIR}/nam_snowfree_albedo.grb"
  greenfrc_climo_file="${FIXED_DIR}/nam_vegfrac.grb"
  z0_climo_file=""
  seaice_anal_dir="${WORK_DIR}"
  seaice_climo_file=""
  sst_anal_dir="${WORK_DIR}"
  sst_climo_file=""
  merge_coeff_sst=0.
  sst_anlcyc_update=.false.
  soilm_climo_file=""
  merge_coeff_soilm=99999.
  snow_depth_anal_dir="${WORK_DIR}"
  snow_depth_climo_file=""
  merge_coeff_snow_depth=0.
  fire_anal_dir="${WORK_DIR}"
 /
 &output
  output_file="${WORK_DIR}/ndas.t06z.input_nemsio.d01.tm12_afterupdate"
  output_file_type="nems"
 /
 &settings
  thinned=.false.
  nsoil=4
 /

 &soil_parameters
  soil_type_src="statsgo"
  smclow = 0.5
  smchigh = 3.0
  smcmax = 0.395, 0.421, 0.434, 0.476, 0.476, 0.439,
           0.404, 0.464, 0.465, 0.406, 0.468, 0.457,
           0.464, -9.99,  0.20, 0.421
  beta =  4.05,  4.26,  4.74,  5.33,  5.33,  5.25,
          6.77,  8.72,  8.17, 10.73, 10.39, 11.55,
          5.25, -9.99,  4.05,  4.26
  psis =  0.0350, 0.0363, 0.1413, 0.7586, 0.7586, 0.3548,
          0.1349, 0.6166, 0.2630, 0.0977, 0.3236, 0.4677,
          0.3548,  -9.99, 0.0350, 0.0363
  satdk = 1.7600e-4, 1.4078e-5, 5.2304e-6, 2.8089e-6, 2.8089e-6,
          3.3770e-6, 4.4518e-6, 2.0348e-6, 2.4464e-6, 7.2199e-6,
          1.3444e-6, 9.7394e-7, 3.3770e-6,     -9.99, 1.4078e-5,
          1.4078e-5
 /
 &veg_parameters
  veg_type_src="igbp"
  salp = 4.0
  snup = 0.080, 0.080, 0.080, 0.080, 0.080, 0.020,
         0.020, 0.060, 0.040, 0.020, 0.010, 0.020,
         0.020, 0.020, 0.013, 0.013, 0.010, 0.020,
         0.020, 0.020
 /
!

mpirun $EXE 
rc=$?

exit $rc