#!/usr/bin/env tcsh
#
# usage: $ sh ./make-figure.sh ?????.nc
#
# history: 2005/12/09 Hiroshi Taniguchi
#

FILESUF=`expr "$1" : "\(.*\).nc"` # 拡張子 .nc だけを除いた名前
FILESUF=`expr "$FILESUF" : "\(exp.....\).*"` # _beta???... だけを除いた名前

gpview $1@ek --wsn 2
mv -f dcl.ps $FILESUF@ek.ps

gpview $1@ektot --wsn 2
mv -f dcl.ps $FILESUF@ektot.ps

gpview $1@ens --wsn 2
mv -f dcl.ps $FILESUF@ens.ps

gpview $1@enstot --wsn 2
mv -f dcl.ps $FILESUF@enstot.ps

gpview $1@hvisc --aspect=1 --exch --wsn 2 
mv -f dcl.ps $FILESUF@hvisc.ps

    for i in velx
    do
	for t in 0.00 0.20 0.40 0.60 0.80 1.00 1.20 1.40 1.60 1.80 2.00 2.20 2.40 2.60 2.80 3.00 3.20 3.40 3.60 3.80 4.00 4.20 4.40 4.60 4.80 5.00 5.20 5.40 5.60 5.80 6.00 6.20 6.40 6.60 6.80 7.00 7.20 7.40 7.60 7.80 8.00 8.20 8.40 8.60 8.80 9.00 9.20 9.40 9.60 9.80 10.0
#	for t in 0.00 0.02 0.04 0.06 0.08 0.10 0.12 0.14 0.16 0.18 0.20 0.22 0.24 0.26 0.28 0.30 0.32 0.34 0.36 0.38 0.40 0.42 0.44 0.46 0.48 0.50
	do
	    gpview $1@$i,t=$t --exch --mean x --aspect=1 --range -1:1 --wsn 2 
	    mv -f dcl.ps $FILESUF@$i\_t=$t.ps
	done
    done

    for i in eksp enssp
    do
	for t in 0.00 0.20 0.40 0.60 0.80 1.00 1.20 1.40 1.60 1.80 2.00 2.20 2.40 2.60 2.80 3.00 3.20 3.40 3.60 3.80 4.00 4.20 4.40 4.60 4.80 5.00 5.20 5.40 5.60 5.80 6.00 6.20 6.40 6.60 6.80 7.00 7.20 7.40 7.60 7.80 8.00 8.20 8.40 8.60 8.80 9.00 9.20 9.40 9.60 9.80 10.0
#	for t in 0.00 0.02 0.04 0.06 0.08 0.10 0.12 0.14 0.16 0.18 0.20 0.22 0.24 0.26 0.28 0.30 0.32 0.34 0.36 0.38 0.40 0.42 0.44 0.46 0.48 0.50
    do
	gpview $1@$i,t=$t --exch --aspect=1 --wsn 2 --nocont --slevels=0,1e-7,1e-6,1e-5,1e-4,1e-3,1e-2,1e-1,999
	    mv -f dcl.ps $FILESUF@$i\_t=$t.ps
#	    mv -f dcl.ps $1@$i\_t=$t.ps
	done
    done

    for i in vor
    do
	for t in 0.00 0.20 0.40 0.60 0.80 1.00 1.20 1.40 1.60 1.80 2.00 2.20 2.40 2.60 2.80 3.00 3.20 3.40 3.60 3.80 4.00 4.20 4.40 4.60 4.80 5.00 5.20 5.40 5.60 5.80 6.00 6.20 6.40 6.60 6.80 7.00 7.20 7.40 7.60 7.80 8.00 8.20 8.40 8.60 8.80 9.00 9.20 9.40 9.60 9.80 10.0
#	for t in 0.00 0.02 0.04 0.06 0.08 0.10 0.12 0.14 0.16 0.18 0.20 0.22 0.24 0.26 0.28 0.30 0.32 0.34 0.36 0.38 0.40 0.42 0.44 0.46 0.48 0.50
    do
	gpview $1@$i,t=$t --aspect=1 --wsn 2 --nocont --slevels -60,-16,-14,-12,-10,-8,-6,-4,-2,0,2,4,6,8,10,12,14,16,60
	    mv -f dcl.ps $FILESUF@$i\_t=$t.ps
#	    mv -f dcl.ps $1@$i\_t=$t.ps
	done
    done

    for i in strfunc
    do
	for t in 0.00 0.20 0.40 0.60 0.80 1.00 1.20 1.40 1.60 1.80 2.00 2.20 2.40 2.60 2.80 3.00 3.20 3.40 3.60 3.80 4.00 4.20 4.40 4.60 4.80 5.00 5.20 5.40 5.60 5.80 6.00 6.20 6.40 6.60 6.80 7.00 7.20 7.40 7.60 7.80 8.00 8.20 8.40 8.60 8.80 9.00 9.20 9.40 9.60 9.80 10.0
#	for t in 0.00 0.02 0.04 0.06 0.08 0.10 0.12 0.14 0.16 0.18 0.20 0.22 0.24 0.26 0.28 0.30 0.32 0.34 0.36 0.38 0.40 0.42 0.44 0.46 0.48 0.50
	do
	gpview $1@$i,t=$t --aspect=1 --wsn 2 --nocont --slevels -2.5,-0.5,-0.45,-0.4,-0.35,-0.3,-0.25,-0.2,-0.15,-0.1,-0.08,-0.06,-0.04,-0.02,0,0.02,0.04,0.06,0.08,0.1,0.15,0.2,0.25,0.3,0.35,0.4,0.45,0.5,2.5
	    mv -f dcl.ps $FILESUF@$i\_t=$t.ps
#	    mv -f dcl.ps $1@$i\_t=$t.ps
	done
    done

for i in `ls *.ps`
do
    ps2epsi $i
    convert -density 180 -geometry 44% -rotate 90 `basename $i .ps`.epsi `basename $i .ps`.gif
    rm -f `basename $i .ps`.epsi
#    convert -density 180 -geometry 1000x750 -rotate 90 `basename $i .ps`.epsi `basename $i .ps`.gif
#    convert -density 180 -geometry 800x600 -rotate 90 `basename $i .ps`.epsi `basename $i .ps`.gif
done

DIR=`expr "$1" : "\(.*\).nc"` # 拡張子 .nc だけを除いたディレクトリ名
if [ -d $DIR ]; then 
    echo "directory $DIR exists."
    else
    echo "mkdir $DIR"
    mkdir $DIR
fi
mv $FILESUF*.gif $DIR
mv $FILESUF*.ps $DIR
rm -f $FILESUF*.epsi
echo "rm -f $FILESUF*.epsi"