#! /bin/csh
#
# Create 1-D files of all mag/plasma parameters from the 1-year time-dependent 
# solutions
#

cd /Volumes/data1-hd2/sumseq/ADAPT/adapt_old/adapt_one_year_scaled/

# find out the total number of time dumps in the directory
set n=0
foreach file (br*.hdf)
  @ n = $n + 1
end
echo n=$n

set i=1
while ($i < $n)
  @ seqn = 1000000 + $i
  set seqn=`echo $seqn | sed 's/^1//'`
  extract -x0 215 -x1 215 -y0 1.571 -y1 1.571 -z0 0.0 -z1 360. br$seqn.hdf /Users/pete/matt-1-year-eq-slices/br_1D_$seqn.hdf
  extract -x0 215 -x1 215 -y0 1.571 -y1 1.571 -z0 0.0 -z1 360. bt$seqn.hdf /Users/pete/matt-1-year-eq-slices/bt_1D_$seqn.hdf
  extract -x0 215 -x1 215 -y0 1.571 -y1 1.571 -z0 0.0 -z1 360. bp$seqn.hdf /Users/pete/matt-1-year-eq-slices/bp_1D_$seqn.hdf
  extract -x0 215 -x1 215 -y0 1.571 -y1 1.571 -z0 0.0 -z1 360. vr$seqn.hdf /Users/pete/matt-1-year-eq-slices/vr_1D_$seqn.hdf
  extract -x0 215 -x1 215 -y0 1.571 -y1 1.571 -z0 0.0 -z1 360. vt$seqn.hdf /Users/pete/matt-1-year-eq-slices/vt_1D_$seqn.hdf
  extract -x0 215 -x1 215 -y0 1.571 -y1 1.571 -z0 0.0 -z1 360. vp$seqn.hdf /Users/pete/matt-1-year-eq-slices/vp_1D_$seqn.hdf
  extract -x0 215 -x1 215 -y0 1.571 -y1 1.571 -z0 0.0 -z1 360. p$seqn.hdf /Users/pete/matt-1-year-eq-slices/p_1D_$seqn.hdf
  extract -x0 215 -x1 215 -y0 1.571 -y1 1.571 -z0 0.0 -z1 360. t$seqn.hdf /Users/pete/matt-1-year-eq-slices/t_1D_$seqn.hdf
  extract -x0 215 -x1 215 -y0 1.571 -y1 1.571 -z0 0.0 -z1 360. rho$seqn.hdf /Users/pete/matt-1-year-eq-slices/rho_1D_$seqn.hdf
  hdf2txt /Users/pete/matt-1-year-eq-slices/br_1D_$seqn
  hdf2txt /Users/pete/matt-1-year-eq-slices/bt_1D_$seqn
  hdf2txt /Users/pete/matt-1-year-eq-slices/bp_1D_$seqn
  hdf2txt /Users/pete/matt-1-year-eq-slices/vr_1D_$seqn
  hdf2txt /Users/pete/matt-1-year-eq-slices/vt_1D_$seqn
  hdf2txt /Users/pete/matt-1-year-eq-slices/vp_1D_$seqn
  hdf2txt /Users/pete/matt-1-year-eq-slices/t_1D_$seqn
  hdf2txt /Users/pete/matt-1-year-eq-slices/p_1D_$seqn
  hdf2txt /Users/pete/matt-1-year-eq-slices/rho_1D_$seqn
  @ i = $i + 1
end