#!/bin/tcsh

set model = "mm5p"

# note: winter and years at beginning and end of run more complicated;
# special-case code has been omitted for clarity

foreach y (1981 1986 1991 1996)

  set z = `expr $y + 4`
  
  foreach year (`seq $y $z`)
    ncra -d time,"$year-03-01 00:00:0.0","$year-06-01 00:00:0.0" \
	pr_"$model"_$y*nc spring-$year.nc
    ncra -d time,"$year-06-01 00:00:0.0","$year-09-01 00:00:0.0" \
	pr_"$model"_$y*nc summer-$year.nc
    ncra -d time,"$year-09-01 00:00:0.0","$year-12-01 00:00:0.0" \
	pr_"$model"_$y*nc autumn-$year.nc
    echo $year; date
  end

end

echo "Years done!"; date



foreach season (spring summer autumn)

  ncea $season-*.nc avg-$season.nc

end

echo "All done!"; date