#! /bin/bash

# Set environment...
# ml GCC ParaStationMPI      # for MPI runs
# export LD_LIBRARY_PATH=../libs/build/lib:$LD_LIBRARY_PATH # removed ../ at the start
# export OMP_NUM_THREADS=4

# Setup...
# working from mptrac/trial directory
trac=../src

# Create directories...
rm -rf data plots && mkdir -p data plots

# Set timestep and timerange of simulation...
t0=$($trac/time2jsec 2005 6 20 0 0 0 0)
t1=$($trac/time2jsec 2005 6 20 21 0 0 0)

# Create control parameter file...
cat > data/trac.ctl <<EOF
NQ = 9
QNT_NAME[0] = t
QNT_NAME[1] = u
QNT_NAME[2] = v
QNT_NAME[3] = w
QNT_NAME[4] = z
QNT_NAME[5] = pv
QNT_NAME[6] = ps
QNT_NAME[7] = pt
QNT_NAME[8] = m
METBASE = data_in/ea
MET_DT_OUT = 10800.0 # this used to be 86400
#CLIM_OH_FILENAME = ../../data/clams_radical_species.nc
#SPECIES = SO2
BOUND_MASS = 0.0
BOUND_DPS = 100.0
CONV_CAPE = 0.0
TDEC_TROP = 259200.0
TDEC_STRAT = 259200.0
DT_MET = 10800.0 # this gives timestep information (10800 = 3hrs which is the spacing bw timesteps for this trial example). used to be 86400
T_STOP = $t1
GRID_LON0 = -30
GRID_LON1 = 98.9
GRID_LAT0 = 5
GRID_LAT1 = 69.8
GRID_NX = 300
GRID_NY = 90
EOF

# Set initial air parcel positions...
$trac/atm_init data/trac.ctl data/atm_init.tab \
	       INIT_T0 "$t0" INIT_T1 "$t0" \
	       INIT_Z0 10.0 INIT_Z1 10.0 \
	       INIT_LON0 75 INIT_LON1 75 \
	       INIT_LAT0 25 INIT_LAT1 25

# Split air parcels...
$trac/atm_split data/trac.ctl data/atm_init.tab data/atm_split.tab \
		SPLIT_N 10000 SPLIT_M 1e9 SPLIT_DX 30.0 SPLIT_DZ 1.0

# Calculate trajectories...
echo "data" > data/dirlist
$trac/trac data/dirlist trac.ctl atm_split.tab \
	   ATM_BASENAME atm GRID_BASENAME grid