#!/bin/bash

echo -n "What format is the data in? (netcdf/binary): "
read nb

if [ $nb = 'binary' ]; then

  echo -n "Would you like to use default settings? (y/n): "
  read yn

  if [ $yn = 'n' ]; then

    rm -f .darwinview/binfilenames
    rm -f .darwinview/*.datbin

    echo -n "Please enter screen resolution: "
    read win
    echo $win > .darwinview/binconfig
    echo -n "Please enter dimensions of data: "
    read dim
    echo $dim >> .darwinview/binconfig
    echo -n "Please enter dimensions of data sets: "
    read dimset
    echo $dimset >> .darwinview/binconfig
    echo .darwinview/binfilenames >> .darwinview/binconfig
    echo -n "Please enter directory containing data: "
    read dir
  
    echo -n "Please enter name of file containing desired species: "
    read fn 
 
    vn=( `cat $fn` ) 

    for name in ${vn[@]}; do
      touch .darwinview/$name.datbin
      echo .darwinview/$name.datbin >> .darwinview/binfilenames 
      for i in `ls $dir`; do
        echo $dir/$i/`ls $dir/$i | grep $name` >> .darwinview/$name.datbin
      done
    done

  fi

elif [ $nb = 'netcdf' ]; then

  echo -n "Would you like to use default settings? (y/n): "
  read yn
 
  if [ $yn = 'n' ]; then

    rm -f .darwinview/ncfilenames
    rm -f .darwinview/*.datnc

    echo -n "Please enter screen resolution: "
    read win
    echo $win > .darwinview/ncconfig
    echo -n "Please enter dimensions of data: "
    read dim
    echo $dim >> .darwinview/ncconfig
    echo -n "Please enter dimensions of data sets: "
    read dimset
    echo $dimset >> .darwinview/ncconfig
    echo .darwinview/ncfilenames >> .darwinview/ncconfig

    echo -n "Please enter directory containing data: "
    read tmp 
    echo $tmp > .darwinview/netcdf
    echo -n "Please enter directory to write to: "
    read tmp 
    echo $tmp >> .darwinview/netcdf

    echo -n "Please enter name of file containing desired species: "
    read tmp 
    echo $tmp >> .darwinview/netcdf
    echo -n "Please enter name of file containing desired tiles: "
    read tmp 
    echo $tmp >> .darwinview/netcdf
    echo -n "Please enter name of file containing desired time steps: "
    read tmp
    echo $tmp >> .darwinview/netcdf

  fi

  indir=( `cat .darwinview/netcdf | head -1` )
  outdir=( `cat .darwinview/netcdf | head -2 | tail -1` )
  fn=( `cat .darwinview/netcdf | head -3 | tail -1` )
  vn=( `cat $fn` )
  fn=( `cat .darwinview/netcdf | head -4 | tail -1` )
  tlist=( `cat $fn` )
  fn=( `cat .darwinview/netcdf | tail -1` )
  itlist=( `cat $fn` )

  initdir=`pwd`
  gcc ~/netcdf.c -I/usr/include/netcdf-3 -L/usr/lib/netcdf-3 -lnetcdf
  for tstep in ${itlist[@]}; do 
    cd $initdir
    cd $outdir
    mkdir -p $tstep
    cd $tstep
    for name in ${vn[@]}; do
      for tile in ${tlist[@]}; do
        /home/marissa/a.out ~/${indir}/ptr_tave.${tstep}.${tile}.nc $name
      done
    done
  done

  cd $initdir

  for name in ${vn[@]}; do
    touch .darwinview/$name.datnc
    echo .darwinview/$name.datnc >> .darwinview/ncfilenames 
    for i in `ls $outdir`; do
      echo $outdir/$i/`ls $outdir/$i | grep $name` >> .darwinview/$name.datnc
    done
  done

fi

gcc darwin.c -lglut
./a.out $nb