/[MITgcm]/MITgcm_contrib/darwinview/src/darwin.bash
ViewVC logotype

Contents of /MITgcm_contrib/darwinview/src/darwin.bash

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph


Revision 1.6 - (show annotations) (download)
Fri Aug 3 19:14:33 2007 UTC (16 years, 8 months ago) by marissa
Branch: MAIN
CVS Tags: HEAD
Changes since 1.5: +21 -2 lines
final bash and c versions

1 #!/bin/bash
2
3 mkdir -p .darwinview
4
5 echo -n "What format is the data in? (binary/netcdf): "
6 read nb
7
8 if [ $nb = 'binary' ]; then
9
10 echo -n "Would you like to use default settings? (y/n): "
11 read yn
12
13 if [ $yn = 'n' ]; then
14
15 rm -f .darwinview/binfilenames
16 rm -f .darwinview/*.datbin
17
18 xrandr | grep '*' | awk -F' ' '{print $1}' > .darwinview/binconfig
19 echo -n "Please enter dimensions of data: "
20 read dim
21 echo $dim >> .darwinview/binconfig
22 echo -n "Please enter dimensions of data sets: "
23 read dimset
24 echo $dimset >> .darwinview/binconfig
25 echo .darwinview/binfilenames >> .darwinview/binconfig
26
27 echo -n "Please enter directory containing data: "
28 read dir
29
30 echo -n "Please enter name of file containing desired species: "
31 read fn
32
33 vn=( `cat $fn` )
34
35 for name in ${vn[@]}; do
36 touch .darwinview/$name.datbin
37 echo .darwinview/$name.datbin >> .darwinview/binfilenames
38 for i in `ls $dir`; do
39 echo $dir/$i/`ls $dir/$i | grep $name` >> .darwinview/$name.datbin
40 done
41 done
42
43 elif [ `echo $yn | grep -v y` ]; then
44 echo "Error: please enter 'y' or 'n.'"
45 exit
46 fi
47
48 elif [ $nb = 'netcdf' ]; then
49
50 echo -n "Would you like to use default settings? (y/n): "
51 read yn
52
53 if [ $yn = 'n' ]; then
54
55 rm -f .darwinview/ncfilenames
56 rm -f .darwinview/*.datnc
57
58 xrandr | grep '*' | awk -F' ' '{print $1}' > .darwinview/ncconfig
59 echo -n "Please enter dimensions of data: "
60 read dim
61 echo $dim >> .darwinview/ncconfig
62 echo -n "Please enter dimensions of data sets: "
63 read dimset
64 echo $dimset >> .darwinview/ncconfig
65 echo .darwinview/ncfilenames >> .darwinview/ncconfig
66
67 echo -n "Please enter directory containing data: "
68 read tmp
69 echo $tmp > .darwinview/netcdf
70
71
72 echo -n "Please enter directory to write to: "
73 read tmp
74 echo $tmp >> .darwinview/netcdf
75
76 echo -n "Please enter name of file containing desired species: "
77 read tmp
78 echo $tmp >> .darwinview/netcdf
79
80 echo -n "Please enter name of file containing desired tiles: "
81 read tmp
82 echo $tmp >> .darwinview/netcdf
83
84 echo -n "Please enter name of file containing desired time steps: "
85 read tmp
86 echo $tmp >> .darwinview/netcdf
87
88 elif [ `echo $yn | grep -v y` ]; then
89 echo "Error: please enter 'y' or 'n.'"
90 exit
91 fi
92
93 indir=( `cat .darwinview/netcdf | head -1` )
94 outdir=( `cat .darwinview/netcdf | head -2 | tail -1` )
95 fn=( `cat .darwinview/netcdf | head -3 | tail -1` )
96 vn=( `cat $fn` )
97 fn=( `cat .darwinview/netcdf | head -4 | tail -1` )
98 tlist=( `cat $fn` )
99 fn=( `cat .darwinview/netcdf | tail -1` )
100 itlist=( `cat $fn` )
101
102 initdir=`pwd`
103 gcc netcdf.c -I/usr/include/netcdf-3 -L/usr/lib/netcdf-3 -lnetcdf
104 for tstep in ${itlist[@]}; do
105 cd $initdir
106 cd $outdir
107 mkdir -p $tstep
108 cd $tstep
109 for name in ${vn[@]}; do
110 for tile in ${tlist[@]}; do
111 $initdir/a.out ~/${indir}/ptr_tave.${tstep}.${tile}.nc $name
112 done
113 done
114 done
115
116 cd $initdir
117
118 for name in ${vn[@]}; do
119 echo $name
120 touch .darwinview/$name.datnc
121 echo .darwinview/$name.datnc >> .darwinview/ncfilenames
122 for i in `ls $outdir`; do
123 echo $outdir/$i/`ls $outdir/$i | grep $name` >> .darwinview/$name.datnc
124 done
125 done
126
127 else
128 echo "Error: please enter 'binary' or 'netcdf.'"
129 exit
130 fi
131
132 if [ ! `ls | grep jet.dat` ]; then
133 echo "Error: jet.dat not found."
134 exit
135 fi
136
137 gcc darwin.c -lglut
138 ./a.out $nb

  ViewVC Help
Powered by ViewVC 1.1.22