if ( $model == 'prem' || $model == 'iasp91' ) then
set tt = taup_time -mod $model -h $depth -ph pP -deg $rec | grep " pP " | awk '{print $4}' |head -n 1 |grep -v "=="
if ( ${tt} == 1) then ; echo $tt $rec >> taup_pP_traveltime2.dat ;
mkdir TAUP
rm -f taup_2.dat
mv taup_.dat TAUP
these code cannot runing? no new dir(/TAUP)