--- MITgcm_contrib/darwinview/src/darwin.bash	2007/07/19 17:03:18	1.2
+++ MITgcm_contrib/darwinview/src/darwin.bash	2007/07/24 16:54:05	1.3
@@ -1,45 +1,120 @@
 #!/bin/bash
 
-echo -n "Would you like to use default settings? (y/n): "
-read yn
+echo -n "What format is the data in? (netcdf/binary): "
+read nb
 
-if echo $yn | grep -q n 
-then
+if [ $nb = 'binary' ]; then
 
-  rm .darwinview/filenames
-  rm .darwinview/*.dat
+  echo -n "Would you like to use default settings? (y/n): "
+  read yn
 
-  echo -n "Please enter screen resolution: "
-  read win
-  echo $win > .darwinview/config
-  echo -n "Please enter dimensions of data: "
-  read dim
-  echo $dim >> .darwinview/config
-  echo -n "Please enter dimensions of data sets: "
-  read dimset
-  echo $dimset >> .darwinview/config
-  echo .darwinview/filenames >> .darwinview/config
-  echo -n "Please enter directory containing data: "
-  read dir
-
-  name="start"
-
-  while [ $name != 'q' ] 
-  do
-   echo -n "Please enter name of species: "
-   read name
-   if [ $name != 'q' ]
-   then
-    echo .darwinview/$name.dat >> .darwinview/filenames 
-    for i in `ls $dir`
-    do
-     touch .darwinview/$name.dat
-     echo $dir/$i/`ls $dir/$i | grep $name` >> .darwinview/$name.dat
+  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
-   fi
   done
 
 fi
 
-make
-./a.out
+gcc darwin.c -lglut
+./a.out $nb