Script to crash mpiMcmc, presumably based on the seed

mpiMcmc.csh

#!/bin/csh
#
# run mpiMcmc with Girardi models on different photometric data subsets

set filters = (BV BVI BVIK JHK UB UBJHK UBVRI UBVRIJHK VI)  # first ran with thin=1, seed=73
set filters = (BVI BVIK JHK UBJHK)				# try these with thin = 3
set filters = (UBV)						# missed before, try thin=1 (ended up bimodal)

set MSRGB    = 0	# Girardi
set THIN     = 10
set SEED     = 703						# try UBV again, with a different seed (still bimodal)
set SEED     = 73						# try UBV again, with yet a different seed, and thin=10

set dir      = "photfiles"
set code_dir = "base/BUILD/bin"

foreach filter ($filters)
   set photfile = "ngc188.$filter.fromUBVRIJHK.noBS.aveprobs.phot"
   set outfile  = "mcmc.$photfile:r.t$THIN.s$SEED"
   $code_dir/mpiMcmc --photFile $dir/$photfile --outputFileBase $dir/$outfile --msRgbModel $MSRGB --thin $THIN --seed $SEED
end