#!/bin/bash

# Field                     Min            Max
# Dark_Matter_Density       0              3.374807e+15
# Density                   5.089015e+11   3.189541e+14              
# X_Ray_Luminosity          2.454139e-02   1.179837e+11
# Temp_X_Ray_Weighted       1.000000e+00   2.360228e+08
# SZ_Y_Effect               7.556865e-15   3.576121e-0
# SZ_Kinetic                -2.837694e-05  2.379747e-05

plotarray=/home/ux455215/Projects/rpwagner/projects/cadac/rcidata/scripts/plotarray

datadir=/oasis/scratch/ux455215/temp_project/lightcone

dset_numbers=`seq -w 9 36`

for dset_number in $dset_numbers
do
    cd $datadir/RD00$dset_number
    for axis in x y z
    do
	for half in front back
	do
	  echo "Plotting fields for RD00$dset_number.$axis.l2.$half.project"
            # Dark matter
	    $plotarray -c hot -m 0 -M 3.374807e+15 Dark_Matter_Density RD00$dset_number.$axis.l2.$half.project
	    
            # Baryon density
	    $plotarray -l -c jet -m 5.089015e+11 -M 3.189541e+14 Density RD00$dset_number.$axis.l2.$half.project

            # X-ray luminosity
	    $plotarray -l -c gist_heat -m 2.454139e-02 -M 1.179837e+11 X_Ray_Luminosity RD00$dset_number.$axis.l2.$half.project
	    
            # X-ray-weighted temperature
            # Lower limit raised for aesthetic appeal
	    $plotarray -l -c spectral -m 1e2 -M 2.360228e+08 Temp_X_Ray_Weighted RD00$dset_number.$axis.l2.$half.project
	    
            # SZ y
	    $plotarray -l -c Paired -m 7.556865e-15 -M 3.576121e-0 SZ_Y_Effect RD00$dset_number.$axis.l2.$half.project
	    
            # SZ kinetic
	    $plotarray -c redshift -m -2.837694e-06 -M 2.379747e-06 SZ_Kinetic RD00$dset_number.$axis.l2.$half.project
	    
	done
    done
done