Makefile.test

#  Copyright (C) 2002 Regents of the University of Michigan, 
#  portions used with permission 
#  For more information, see http://csem.engin.umich.edu/tools/swmf
#^CMP FILE TESTING

TESTDIR = run_test

TESTMACHINE = `hostname | sed -e 's/\..*//; s/[0-9]*$$//'`

test_help:
	@echo "    test                       (run all tests with ${MPIRUN})"
	@echo "    test MPIRUN=               (run tests serially)"
	@echo "    test NPFLAG=-n NP=3        (run tests with ${PARALLEL} -n 3)"
	@echo "    test PARALLEL='mpiexec'    (run tests with mpiexec)"
	@echo "    test8_compile              (configure and compile test8)"
	@echo "    test8_rundir               (create run directory for test8)"
	@echo "    test8_run NP=4             (run test8 with 4 MPI )"
	@echo "    test8_check                (check test8 results)"
	@echo "    test8_check BLESS=YES      (bless test8 results = update references)"
	@echo "    test8 TESTDIR=run_test08   (run test8 in run_test08 directory)"
	@echo "    test_gpu                   (run GPU related tests only)"
	@echo "    test_share_library         (run share/Library tests)"
	@echo "    test_util                  (run util/ tests)"
	@echo "    test_batl                  (run BATL tests if BATL/ is present)"
	@echo "    test_comp                  (run component tests)"
	@echo "    test_swmf                  (run SWMF tests below)"
	@(if [ -d PW/PWOM    ]; then echo "    test_pw                    (run PW/PWOM test)"; fi)
	@(if [ -d PS/DGCPM   ]; then echo "    test_ps                    (run PS/DGCPM test)"; fi)
	@(if [ -d SP/MFLAMPA ]; then echo "    test_sp                    (run SP/MFLAMPA test)"; fi)
	@(if [ -d IM/RCM2    ]; then echo "    test_swpc                  (run GM-IE-RCM2 test as used by SWPC)"; fi)
	@(if [ -d IM/RCM2    ]; then echo "    test_swpc_build            (create SWPC_build and run test_swpc)"; fi)
	@(if [ -d IM/RCM2    ]; then echo "    test_swpc_gpu              (run GM-IE-RCM2 test on GPU)"; fi)
	@(if [ -d IM/RCM2    ]; then echo "    test_swpc_large_gpu        (run large GM-IE-RCM2 test on GPU)"; fi)
	@(if [ -d IM/RCM2    ]; then echo "    test_swpc_multiion         (run multiion GM-IE-RCM2 SWPC test)"; fi)
	@(if [ -d IM/RCM2    ]; then echo "    test_swpc_multispecies     (run multispecies GM-IE-RCM2 SWPC test)"; fi)
	@(if [ -d IM/CIMI2   ]; then echo "    test_swpc_cimi             (run anisotropic GM-IE-CIMI2 SWPC test)"; fi)
	@(if [ -d PW/PWOM    ]; then echo "    test_swpc_pwom             (run GM-IE-RCM2-PWOM SWPC test)"; fi)
	@(if [ -d PW/PWOM    ]; then echo "    test_swpc_pwom_species     (run multispecies GM-IE-RCM2-PWOM SWPC test)"; fi)
	@(if [ -d PW/PWOM    ]; then echo "    test1                      (run GM-IE-PWOM test)"; fi)
	@(if [ -d IM/HEIDI   ]; then echo "    test4                      (run GM-IE-HEIDI test)"; fi)
	@(if [ -d GM/BATSRUS/srcUserExtra ]; then echo "    test5                      (run EE-SC test)"; fi)
	@(if [ -d GM/BATSRUS/srcUserExtra ]; then echo "    test7                      (run IH-OH test)"; fi)
	@(if [ -d GM/BATSRUS ]; then echo "    test8                      (run SC-IH coupling test with AWSoM-R)"; fi)
	@(if [ -d GM/BATSRUS ]; then echo "    test9                      (run SC-IH coupling test with AWSoM)"; fi)
	@(if [ -d PT/AMPS    ]; then echo "    test11                     (run GM-PT/AMPS coupler test)"; fi)
	@(if [ -d PC/AMPS    ]; then echo "    test12                     (run GM-PC/AMPS tests)"; fi)
	@(if [ -d PC/IPIC3D2 ]; then echo "    test12_ipic3d              (run GM-PC/IPIC3D2 tests)"; fi)
	@(if [ -d PT/AMPS    ]; then echo "    test13                     (run GM-PT/AMPS test: Europa)"; fi)
	@(if [ -d PT/AMPS    ]; then echo "    test14                     (run OH-PT/AMPS test)"; fi)
	@(if [ -d PT/AMPS    ]; then echo "    test14_1d                  (run OH-PT/AMPS test in 1D)"; fi)
	@(if [ -d PT/AMPS    ]; then echo "    test14pui                  (run OH-PT/AMPS test with PUI)"; fi)
	@(if [ -d PT/AMPS    ]; then echo "    test14pui1d                (run OH-PT/AMPS test with PUI in 1D)"; fi)
	@(if [ -d SP/MFLAMPA ]; then echo "    test15                     (run SC-IH-SP/MFLAMPA test)"; fi)
	@(if [ -d SP/MFLAMPA ]; then echo "    test15_build               (configure SC-IH-SP/MFLAMPA, run test15)"; fi)
	@(if [ -d PT/AMPS    ]; then echo "    test15_pt                  (run SC-IH-PT/AMPS test)"; fi)
	@(if [ -d PC/FLEKS   ]; then echo "    test16                     (run GM-PC/FLEKS test: periodic FLEKS)"; fi)
	@(if [ -d PC/FLEKS   ]; then echo "    test17                     (run GM-PC/FLEKS test: 2D Alfven wave)"; fi)
	@(if [ -d PC/FLEKS   ]; then echo "    test18                     (run GM-PC/FLEKS test: 3D reconnection)"; fi)
	@echo "    test_check                 (collect test results into test_swmf.res)"

.NOTPARALLEL:

test:
	@rm -f *.diff
	-@(${MAKE} test_share_library)
	-@(${MAKE} test_util)
	-@(${MAKE} test_comp)
	-@(${MAKE} test_swmf)
	-@(if [ -d BATL ]; then ${MAKE} test_batl; fi)
	-@(${MAKE} test_check)

test_gpu:
	@rm -f *.diff
	-@(${MAKE} test_swpc_gpu TESTDIR=run_test_swpc_gpu)
	-@(${MAKE} test_swpc_large_gpu TESTDIR=run_test_swpc_large_gpu)
	-@(cd GM/BATSRUS; ${MAKE} test_gpu)
	-@(${MAKE} test_check)

test_swmf:
	-@(if [ -d PS/DGCPM   ]; then ${MAKE} test_ps TESTDIR=run_test_ps; fi)
	-@(if [ -d SP/MFLAMPA ]; then ${MAKE} test_sp TESTDIR=run_test_sp; fi)
	-@(if [ -d IM/RCM2    ]; then ${MAKE} test_swpc_gpu TESTDIR=run_test_swpc_gpu; fi)
	-@(if [ -d IM/RCM2    ]; then ./Config.pl -noacc; fi)
	-@(if [ -d IM/RCM2    ]; then ${MAKE} test_swpc TESTDIR=run_test_swpc; fi)
	-@(if [ -d IM/RCM2    ]; then ${MAKE} test_swpc_multiion TESTDIR=run_test_swpc_multiion; fi)
	-@(if [ -d IM/RCM2    ]; then ${MAKE} test_swpc_multispecies TESTDIR=run_test_swpc_multispecies; fi)
	-@(if [ -d IM/CIMI2   ]; then ${MAKE} test_swpc_cimi TESTDIR=run_test_swpc_cimi; fi)
	-@(if [ -d PW/PWOM    ]; then ${MAKE} test_swpc_pwom TESTDIR=run_test_swpc_pwom; fi)
	-@(if [ -d PW/PWOM    ]; then ${MAKE} test_swpc_pwom_species TESTDIR=run_test_swpc_pwom_species; fi)
	-@(if [ -d PW/PWOM    ]; then ${MAKE} test1 TESTDIR=run_test01; fi)
	-@(if [ -d IM/HEIDI   ]; then ${MAKE} test4 TESTDIR=run_test04; fi)
	-@(if [ -d GM/BATSRUS/srcUserExtra ]; then ${MAKE} test5 TESTDIR=run_test05; fi)
	-@(if [ -d GM/BATSRUS/srcUserExtra ]; then ${MAKE} test7 TESTDIR=run_test07; fi)
	-@(if [ -d GM/BATSRUS ]; then ${MAKE} test8 TESTDIR=run_test08; fi)
	-@(if [ -d GM/BATSRUS ]; then ${MAKE} test9 TESTDIR=run_test09; fi)
	-@(${MAKE} test12)
	-@(if [ -d PT/AMPS    ]; then ${MAKE} test13 TESTDIR=run_test13; fi)
	-@(if [ -d PT/AMPS    ]; then ${MAKE} test14_1d TESTDIR=run_test14_1d; fi)
	-@(if [ -d PT/AMPS    ]; then ${MAKE} test14 TESTDIR=run_test14; fi)
	-@(if [ -d PT/AMPS    ]; then ${MAKE} test14pui TESTDIR=run_test14pui; fi)
	-@(if [ -d PT/AMPS    ]; then ${MAKE} test14pui1d TESTDIR=run_test14pui1d; fi)
	-@(if [ -d SP/MFLAMPA ]; then ${MAKE} test15 TESTDIR=run_test15; fi)
	-@(if [ -d PC/FLEKS   ]; then ${MAKE} test16 TESTDIR=run_test16; fi)
	-@(if [ -d PC/FLEKS   ]; then ${MAKE} test17 TESTDIR=run_test17; fi)
	-@(if [ -d PC/FLEKS   ]; then ${MAKE} test18 TESTDIR=run_test18; fi)
	-@(if [ -d SP/MFLAMPA ]; then ${MAKE} test19 TESTDIR=run_test19; fi)
	./Config.pl -noamrex

# Setting BLESS=YES or Y will copy the solution into the reference solution
BLESS=NO

DIFFNUM = ${SCRIPTDIR}/DiffNum.pl -BLESS=${BLESS}

DIFF_FILES = share/Library/test/*.diff util/*/src*/*.diff ??/*/*.diff *.diff

# This expression provides the number of AMPS processes used in a test
# Unfortunately AMPS results depend on the number of processors
NAMPS = `ls ${TESTDIR}/PT/thread* | wc -l | tr -d ' '`

# Note: the "cat ... > /dev/null" makes sure that file sizes are correct on nyx

test_check:
	-@cat ${DIFF_FILES} > /dev/null
	-@ls -ltr  ${DIFF_FILES} > test_swmf.res
	-@(if [ -f BATL/test11.diff ]; then ls -ltr BATL/*.diff >> test_swmf.res; fi)
	@cat test_swmf.res
	@echo '=============================================================='\
		>> test_swmf.res
	-@head -100 ${DIFF_FILES} >> test_swmf.res
	-@(if [ -f BATL/test11.diff ]; then head -100 BATL/*.diff >> test_swmf.res; fi)
	@echo ' '
	@echo 'See test_swmf.res and test_swmf.log for more detail...'
	@echo ' '
	@echo '=============================================================='\
		>> test_swmf.res
	share/Scripts/test_speed.pl . >> test_swmf.res

### Unit tests of the share/Library ###

test_share_library:
	rm -f share/Library/test/*.diff
	cd share/Library/src; ${MAKE} LIB
	cd share/Library/test; make test
	ls -l share/Library/test/*.diff

### Unit tests of the utilities ###

test_util:
	cd util; make test

### Unit tests of the BATL library ###

test_batl:
	cd BATL; ${MAKE} test

### Component tests ###

test_comp:
	for i in `ls -d [A-Z][A-Z]/*/ | grep -v /Empty/`; \
		do ( cd $$i; ${MAKE} test; ); done
	ls -ltr [A-Z][A-Z]/*/*.diff

### SWMF tests ###

### test_cz: convection zone model FSAM test ###

test_cz:
	@echo "starting..."        >  test_cz.diff
	@echo "test_cz_compile..." >> test_cz.diff
	${MAKE} test_cz_compile
	@echo "test_cz_rundir..."  >> test_cz.diff
	${MAKE} test_cz_rundir
	@echo "test_cz_run..."     >> test_cz.diff
	${MAKE} test_cz_run
	@echo "test_cz_check..."   >> test_cz.diff
	${MAKE} test_cz_check

test_cz_compile:
	./Config.pl -v=Empty,CZ/FSAM -fishpak
	cp CZ/FSAM/srcProblem/ModUserSetup.SWMF.f90 CZ/FSAM/src/ModUserSetup.f90
	cp CZ/FSAM/srcProblem/ModPar.SWMF.f90 CZ/FSAM/src/ModPar.f90
	${MAKE} SWMF

test_cz_rundir:
	rm -rf ${TESTDIR}
	${MAKE} rundir RUNDIR=${DIR}/${TESTDIR}
	cp Param/PARAM.in.test.CZ ${TESTDIR}/PARAM.in
	-Scripts/TestParam.pl -F ${TESTDIR}/PARAM.in
	cp ${CZDIR}/input/gong.l4b.14 ${TESTDIR}/

test_cz_run:
	cd ${TESTDIR}; ${PARALLEL} ${NPFLAG} 6 ./SWMF.exe > runlog

test_cz_check:
	-${DIFFNUM} -a=1e-5 -t -p="grep 'ref error='" \
		${TESTDIR}/runlog \
		${CZDIR}/output/runlog.ref \
		> test_cz.diff
	ls -ltr test_cz.diff

### test_pw: Polarwind test ###

test_pw:
	@echo "starting..." > test_pw.diff
	@echo "test_pw_compile..." >> test_pw.diff
	${MAKE} test_pw_compile
	@echo "test_pw_rundir..." >> test_pw.diff
	${MAKE} test_pw_rundir
	@echo "test_pw_run..." >> test_pw.diff
	${MAKE} test_pw_run
	@echo "test_pw_check..." >> test_pw.diff
	${MAKE} test_pw_check

test_pw_compile:
	./Config.pl -v=Empty,PW/PWOM -o=PW:Earth
	${MAKE} SWMF
	make PIDL

test_pw_rundir:
	rm -rf ${TESTDIR}
	${MAKE} rundir RUNDIR=${DIR}/${TESTDIR}
	cp Param/PARAM.in.test.PW ${TESTDIR}/PARAM.in
	-Scripts/TestParam.pl -F ${TESTDIR}/PARAM.in

test_pw_run:
	cd ${TESTDIR}/PW; rm -f log.out restartOUT/* plots/*
	cd ${TESTDIR}; ${MPIRUN} ./SWMF.exe > runlog

test_pw_check:
	-@(${DIFFNUM} -b -r=1e-8 \
		${TESTDIR}/PW/plots/north_plots_iline0001.out \
		${PWDIR}/data/output/Earth/north_plots_iline0001.out \
		> test_pw.diff)
	ls -ltr test_pw*.diff


### test_ps: Plasmasphere model DGCPM test ###

test_ps:
	@echo "starting..." > test_ps.diff
	@echo "test_ps_compile..." >> test_ps.diff
	${MAKE} test_ps_compile
	@echo "test_ps_rundir..." >> test_ps.diff
	${MAKE} test_ps_rundir
	@echo "test_ps_run..." >> test_ps.diff
	${MAKE} test_ps_run
	@echo "test_ps_check..." >> test_ps.diff
	${MAKE} test_ps_check	

test_ps_compile:
	./Config.pl -v=Empty,PS/DGCPM
	${MAKE} SWMF
	
test_ps_rundir:
	rm -rf ${TESTDIR}
	${MAKE} rundir RUNDIR=${DIR}/${TESTDIR}
	cp Param/PARAM.in.test.PS ${TESTDIR}/PARAM.in
	-Scripts/TestParam.pl -F ${TESTDIR}/PARAM.in
	cp ${PSDIR}/Input/kp_test.dat ${TESTDIR}/PS/Input/

test_ps_run:
	cd ${TESTDIR}; ${PARALLEL} ${NPFLAG} 1 ./SWMF.exe > runlog
	cd ${TESTDIR}; ./PostProc.pl -M RESULTS

test_ps_check:
	-@(${DIFFNUM} -t -r=1e-8 \
		${TESTDIR}/RESULTS/PS/dgcpm_t20011021_010000_000.dat \
		${PSDIR}/Output/dgcpm_t20011021_010000_000.ref \
		> test_ps.diff)
	-@(${DIFFNUM} -t -r=1e-8 \
		${TESTDIR}/RESULTS/PS/PS_t20011021_000000_000.log \
		${PSDIR}/Output/PS_t20011021_000000_000.ref \
		>> test_ps.diff)
	-@(${DIFFNUM} -t -r=1e-8 \
		${TESTDIR}/RESULTS/PS/slice_L6p60_t20011021_000000_000.dat \
		${PSDIR}/Output/slice_L6p60_t20011021_000000_000.ref \
		>> test_ps.diff)
	ls -ltr test_ps*.diff


### test_psie: Plasmasphere-Ionosphere Coupling: PS-IE/DGCPM-RIM ###

test_psie:
	@echo "starting..." > test_psie.diff
	@echo "test_psie_compile..." >> test_psie.diff
	${MAKE} test_psie_compile
	@echo "test_psie_rundir..." >> test_psie.diff
	${MAKE} test_psie_rundir
	@echo "test_psie_run..." >> test_psie.diff
	${MAKE} test_psie_run
	@echo "test_psie_check..." >> test_psie.diff
	${MAKE} test_psie_check	

test_psie_compile:
	./Config.pl -v=Empty,PS/DGCPM,IE/RIM
	${MAKE} SWMF

test_psie_rundir:
	rm -rf ${TESTDIR}
	${MAKE} rundir RUNDIR=${DIR}/${TESTDIR}
	cp Param/PARAM.in.test.PSIE  ${TESTDIR}/PARAM.in
	-Scripts/TestParam.pl -F ${TESTDIR}/PARAM.in
	cp ${PSDIR}/Input/IMF_NSturning_10nT.dat ${TESTDIR}/

test_psie_run:
	cd ${TESTDIR}; ${PARALLEL} ${NPFLAG} 1 ./SWMF.exe > runlog
	cd ${TESTDIR}/IE; ./pIE
	cd ${TESTDIR}; ./PostProc.pl -M RESULTS

test_psie_check:
	-@(${DIFFNUM} -t -r=1e-8                        \
	        ${TESTDIR}/RESULTS/IE/IE.log                         \
		${PSDIR}/Output/IE_test_psie.ref                     \
		> test_psie.diff)
	-@(${DIFFNUM} -t -r=1e-8                        \
	        ${TESTDIR}/RESULTS/PS/PS_t20000321_104500_000.log    \
		${PSDIR}/Output/PS_test_psie.ref                     \
		>> test_psie.diff)
	-@(${DIFFNUM} -t -r=1e-8                        \
	        ${TESTDIR}/RESULTS/PS/dgcpm_t20000321_105500_000.dat \
		${PSDIR}/Output/dgcpm_testpsie.ref                   \
		>> test_psie.diff)
	ls -ltr test_psie*.diff
	
### test1: PW+IE+GM ###

test1:
	@echo "test1_compile..." > test01_gm.diff
	@echo "test1_compile..." > test01_ie.diff
	@echo "test1_compile..." > test01_pw.diff
	${MAKE} test1_compile
	@echo "test1_rundir..." >> test01_gm.diff
	@echo "test1_rundir..." >> test01_ie.diff
	@echo "test1_rundir..." >> test01_pw.diff
	${MAKE} test1_rundir
	@echo "test1_run..." >> test01_gm.diff
	@echo "test1_run..." >> test01_ie.diff
	@echo "test1_run..." >> test01_pw.diff
	${MAKE} test1_run
	@echo "test1_check..." >> test01_gm.diff
	@echo "test1_check..." >> test01_ie.diff
	@echo "test1_check..." >> test01_pw.diff
	${MAKE} test1_check

test1_compile:
	./Config.pl -v=Empty,GM/BATSRUS,IE/Ridley_serial,PW/PWOM
	./Config.pl -o=GM:u=Default,e=Mhd,ng=2,g=8,8,8,PW:Earth,IE:g=91,181
	${MAKE} SWMF
	make PIDL

test1_rundir:
	rm -rf ${TESTDIR}
	${MAKE} rundir RUNDIR=${DIR}/${TESTDIR}
	cp Param/PARAM.in.test.GMIEPW ${TESTDIR}/PARAM.in
	-Scripts/TestParam.pl -F ${TESTDIR}/PARAM.in

test1_run:
	cd ${TESTDIR}/PW; rm -f log.out restartOUT/* plots/*
	cd ${TESTDIR}; ${MPIRUN} ./SWMF.exe > runlog
	cd ${TESTDIR}; ./PostProc.pl -M RESULTS

test1_check:
	-@(${DIFFNUM} -b -r=1e-5 \
		${TESTDIR}/RESULTS/PW/north_plots_iline0001.out \
		PW/PWOM/data/output/Earth/PW_north_plots_iline0001.out \
		> test01_pw.diff)
	-@(${DIFFNUM} -b -a=1e-20 -r=1e-7 \
		${TESTDIR}/RESULTS/PW/north_plots_iline0002.out \
		PW/PWOM/data/output/Earth/PW_north_plots_iline0002.out \
		>> test01_pw.diff)
	-@(${DIFFNUM} -b -r=2e-7 \
		${TESTDIR}/RESULTS/PW/north_plots_iline0003.out \
		PW/PWOM/data/output/Earth/PW_north_plots_iline0003.out \
		>> test01_pw.diff)
	-@(${DIFFNUM} -b -r=1e-7 \
		${TESTDIR}/RESULTS/PW/north_plots_iline0004.out \
		PW/PWOM/data/output/Earth/PW_north_plots_iline0004.out \
		>> test01_pw.diff)
	-@(${DIFFNUM} -b -r=1e-3 -a=1e-8 \
		${TESTDIR}/RESULTS/GM/log_n000001.log \
		output/test1/GM_log_n000001.log \
		> test01_gm.diff)
	-@(${DIFFNUM} -b -t -r=1e-5 \
		${TESTDIR}/RESULTS/IE/IE*.log \
		output/test1/IE.log \
		> test01_ie.diff)
	-@(${DIFFNUM} -b -r=1e-4 -a=1e-4 \
		${TESTDIR}/RESULTS/IE/it000321_104510_000.idl \
		output/test1/IE*.idl.gz \
		>> test01_ie.diff)
	ls -l test01_*.diff

### test2: GM+IE+RCM with multifluid GM-IM coupling ###

test2:  
	@echo "test2_compile..." > test02_gm_multi.diff
	@echo "test2_compile..." > test02_ie.diff	
	@echo "test2_compile..." > test02_im_rcm_multi.diff
	${MAKE} test2_compile
	@echo "test2_rundir..." >> test02_gm_multi.diff
	@echo "test2_rundir..." >> test02_ie.diff
	@echo "test2_rundir..." >> test02_im_rcm_multi.diff
	${MAKE} test2_rundir
	@echo "test2_run..." >> test02_gm_multi.diff
	@echo "test2_run..." >> test02_ie.diff
	@echo "test2_run..." >> test02_im_rcm_multi.diff
	${MAKE} test2_run
	@echo "test2_restart..." >> test02_gm_multi.diff
	@echo "test2_restart..." >> test02_ie.diff
	@echo "test2_restart..." >> test02_im_rcm_multi.diff
	${MAKE} test2_restart
	@echo "test2_check..." >> test02_gm_multi.diff
	@echo "test2_check..." >> test02_ie.diff
	@echo "test2_check..." >> test02_im_rcm_multi.diff
	${MAKE} test2_check

test2_compile:
	./Config.pl -v=Empty,GM/BATSRUS,IE/Ridley_serial,IM/RCM2
	./Config.pl -o=GM:u=Default,e=MhdIons,ng=2,g=8,8,8,IE:g=91,181
	${MAKE} SWMF
	make PIDL

test2_rundir:
	rm -rf ${TESTDIR}
	${MAKE} rundir RUNDIR=${DIR}/${TESTDIR}
	cp Param/PARAM.in.test.GMIERCM_Multi ${TESTDIR}/PARAM.in
	-Scripts/TestParam.pl -F ${TESTDIR}/PARAM.in

test2_run:
	cd ${TESTDIR}; ${MPIRUN} ./SWMF.exe > runlog
	cd ${TESTDIR}; ./Restart.pl

test2_restart:
	cd ${TESTDIR}; mv PARAM.in PARAM.in_start; rm -f PARAM.in_orig_
	cd ${TESTDIR}; cp Param/PARAM.in.test.restart.GMIERCM_Multi PARAM.in
	-Scripts/TestParam.pl -F ${TESTDIR}/PARAM.in
	cd ${TESTDIR}; ${MPIRUN} ./SWMF.exe > runlog_restart
	cd ${TESTDIR}/IM/plots/; cp 2d__max_t00000020.dat t20.txt
	cd ${TESTDIR}; ./PostProc.pl -M -cat RESULTS; \
		mv runlog_restart RESULTS

test2_check:
	-@(${DIFFNUM} -b -r=1e-5 -a=1e-8 \
		${TESTDIR}/RESULTS/GM/log_n000001.log \
		output/test2/GM_log_n000001.log \
		> test02_gm_multi.diff)
	-@(${DIFFNUM} -b -r=1e-5 -a=1e-8 \
		${TESTDIR}/RESULTS/GM/magnetometers_n*.mag \
		output/test2/magnetometers.mag \
		>> test02_gm_multi.diff)
	-@(${DIFFNUM} -b -r=1e-5 -a=1e-8 \
		${TESTDIR}/RESULTS/GM/geoindex_n*.log \
		output/test2/GM_geoindex.log \
		>> test02_gm_multi.diff)
	-@(${DIFFNUM} -b -r=1e-5 -a=1e-4 \
		${TESTDIR}/RESULTS/RESTART/GM/x_geoindex.rst \
		output/test2/GM_restart_geoindex.txt \
		>> test02_gm_multi.diff)
	-@(${DIFFNUM} -b -r=1e-5 \
		${TESTDIR}/RESULTS/IE/IE_t000321_104500.log \
		output/test2/IE.log \
		> test02_ie.diff) 
	-@(${DIFFNUM} -b -r=1e-3 -a=1e-8 -p="grep -v 'AUXDATA SAVEDATE='" \
		${TESTDIR}/RESULTS/IM/t20.dat \
		output/test2/IM_2d__max_t00000020.dat.gz \
		> test02_im_rcm_multi.diff)
	ls -l test02_*.diff

### test_swpc_build: configure SWMF with GM+IE+IM+RB only. ###

# Configre with 4 components: GM+IE+IM+RB. Remove unused models.
# Run test_swpc in the reduced SWMF build SWPC_build.

BUILD=SWPC_build
test_swpc_build:
	rm -rf ${BUILD}
	Scripts/Configure.pl -d=${BUILD}
	cd ${BUILD}; rm -rf IM/HEIDI IM/CIMI IE/RIM util/HYPRE
	cd ${BUILD}; ./Config.pl -install \
	   	     	-compiler=${FORTRAN_COMPILER_NAME},${C_COMPILER_NAME}
	cd ${BUILD}; make test_swpc

### test_swpc: GM+IE+IM+RB ###

test_swpc:
	@echo "test_swpc_compile..." > test_swpc_rcm_gm.diff
	@echo "test_swpc_compile..." > test_swpc_rcm_ie.diff
	@echo "test_swpc_compile..." > test_swpc_rcm_rb.diff
	${MAKE} test_swpc_compile
	@echo "test_swpc_rundir..." >> test_swpc_rcm_gm.diff
	@echo "test_swpc_rundir..." >> test_swpc_rcm_ie.diff
	@echo "test_swpc_rundir..." >> test_swpc_rcm_rb.diff
	${MAKE} test_swpc_rundir
	@echo "test_swpc_run..." >> test_swpc_rcm_gm.diff
	@echo "test_swpc_run..." >> test_swpc_rcm_ie.diff
	@echo "test_swpc_run..." >> test_swpc_rcm_rb.diff
	${MAKE} test_swpc_run
	@echo "test_swpc_restart..." >> test_swpc_rcm_gm.diff
	@echo "test_swpc_restart..." >> test_swpc_rcm_ie.diff
	@echo "test_swpc_restart..." >> test_swpc_rcm_rb.diff
	${MAKE} test_swpc_restart
	@echo "test_swpc_check..." >> test_swpc_rcm_gm.diff
	@echo "test_swpc_check..." >> test_swpc_rcm_ie.diff
	@echo "test_swpc_check..." >> test_swpc_rcm_rb.diff
	${MAKE} test_swpc_check

test_swpc_compile:
	./Config.pl -v=Empty,GM/BATSRUS,IE/Ridley_serial,IM/RCM2,RB/RBE
	./Config.pl -o=GM:u=Default,e=Mhd,ng=2,g=8,8,8,IE:g=91,181
	${MAKE} SWMF
	make PIDL INTERPOLATE

test_swpc_rundir:
	rm -rf ${TESTDIR}
	${MAKE} rundir RUNDIR=${DIR}/${TESTDIR}
	cp Param/SWPC/*.in* Param/SWPC/*.dat ${TESTDIR}/
	-Scripts/TestParam.pl -F ${TESTDIR}/PARAM.in_SWPC_v2_init
	perl -pi -e 'if(/MaxIter|MaxBlock/){s/700/70/; s/1500/200/; s/5000/350/}; s/#BORIS/BORIS/; s/^\#(COMPONENTMAP.*production)/$$1/i; s/^(COMPONENTMAP.*nightly)/\#$$1/i' ${TESTDIR}/PARAM.in_SWPC_v2_init
	cp ${TESTDIR}/PARAM.in_SWPC_v2_init ${TESTDIR}/PARAM.in
	ln -s ${BINDIR}/INTERPOLATE.exe ${TESTDIR}/

test_swpc_run:
	cd ${TESTDIR}; ${MPIRUN} ./SWMF.exe > runlog
	cd ${TESTDIR}; ./Restart.pl -i SWMF_RESTART.20140410_000200

test_swpc_restart:
	-Scripts/TestParam.pl -F ${TESTDIR}/PARAM.in_SWPC_v2_restart
	perl -pi -e 's/#BORIS/BORIS/; s/^\#(COMPONENTMAP.*production)/$$1/i; s/^(COMPONENTMAP.*nightly)/\#$$1/i' ${TESTDIR}/PARAM.in_SWPC_v2_restart
	cd ${TESTDIR}; cp PARAM.in_SWPC_v2_restart PARAM.in
	cd ${TESTDIR}; ${MPIRUN} ./SWMF.exe > runlog_restart
	cd ${TESTDIR}; ./PostProc.pl -noptec
	cd ${TESTDIR}; ./INTERPOLATE.exe < INTERPOLATE.in

test_swpc_check:
	-@(${DIFFNUM} -b -r=1e-5 -a=2e-8 \
		${TESTDIR}/GM/IO2/log_e20140410-000200.log \
		Param/SWPC/TestOutput/log_e20140410-000200.log \
		> test_swpc_rcm_gm.diff)
	-@(${DIFFNUM} -b -r=2e-4 -a=1e-5 \
		${TESTDIR}/GM/IO2/magnetometers_e20140410-000200.mag \
		Param/SWPC/TestOutput/magnetometers_e20140410-000200.mag \
		>> test_swpc_rcm_gm.diff)
	-@(${DIFFNUM} -b -r=2e-4 -a=1e-5 \
		${TESTDIR}/GM/IO2/ABK.txt \
		Param/SWPC/TestOutput/ABK.txt \
		>> test_swpc_rcm_gm.diff)
	-@(${DIFFNUM} -b -r=2e-4 -a=1e-4 \
		${TESTDIR}/GM/IO2/mag_grid_e20140410-000300.out \
		Param/SWPC/TestOutput/mag_grid_e20140410-000300.out.gz \
		>> test_swpc_rcm_gm.diff)
	-@(${DIFFNUM} -b -r=1e-5 -a=1e-8 \
		${TESTDIR}/GM/IO2/geoindex_e20140410-000200.log \
		Param/SWPC/TestOutput/geoindex_e20140410-000200.log \
		>> test_swpc_rcm_gm.diff)
	-@(${DIFFNUM} -b -r=1e-5 \
		${TESTDIR}/IE/ionosphere/IE_t140410_000200.log \
		Param/SWPC/TestOutput/IE_t140410_000200.log \
		> test_swpc_rcm_ie.diff)
	-@(${DIFFNUM} -b -r=2e-4 -a=6e-5 \
		${TESTDIR}/IE/ionosphere/it140410_000300_000.idl \
		Param/SWPC/TestOutput/it140410_000300_000.idl.gz \
		>> test_swpc_rcm_ie.diff)
	-@(${DIFFNUM} -b -r=1e-3 \
                ${TESTDIR}/RB/plots/20140410_000300_e.fls \
		Param/SWPC/TestOutput/RB_20140410_000300_e.fls.gz \
		> test_swpc_rcm_rb.diff)
	ls -l test_swpc_rcm_??.diff

### test_swpc_gpu: GM+IE+IM and GM runs on GPU or in GPU compatible mode ###

test_swpc_gpu:
	@echo "test_swpc_gpu_compile..." > test_swpc_gpu_gm.diff
	@echo "test_swpc_gpu_compile..." > test_swpc_gpu_ie.diff
	${MAKE} test_swpc_gpu_compile
	@echo "test_swpc_gpu_rundir..." >> test_swpc_gpu_gm.diff
	@echo "test_swpc_gpu_rundir..." >> test_swpc_gpu_ie.diff
	${MAKE} test_swpc_gpu_rundir
	@echo "test_swpc_gpu_run..." >> test_swpc_gpu_gm.diff
	@echo "test_swpc_gpu_run..." >> test_swpc_gpu_ie.diff
	${MAKE} test_swpc_gpu_run NP=2
	@echo "test_swpc_gpu_restart..." >> test_swpc_gpu_gm.diff
	@echo "test_swpc_gpu_restart..." >> test_swpc_gpu_ie.diff
	${MAKE} test_swpc_gpu_restart NP=2
	@echo "test_swpc_gpu_check..." >> test_swpc_gpu_gm.diff
	@echo "test_swpc_gpu_check..." >> test_swpc_gpu_ie.diff
	${MAKE} test_swpc_gpu_check

test_swpc_gpu_compile:
	./Config.pl -v=Empty,GM/BATSRUS,IE/Ridley_serial,IM/RCM2
	./Config.pl -o=GM:u=Default,e=Mhd,ng=2,g=8,8,8,IE:g=91,181
	./Config.pl -o=GM:opt=../../Param/SWPC/PARAM.in_SWPC_gpu_init
	./Config.pl -acc
	${MAKE} SWMF
	make PIDL

test_swpc_gpu_rundir:
	rm -rf ${TESTDIR}
	${MAKE} rundir RUNDIR=${DIR}/${TESTDIR}
	cp Param/SWPC/PARAM.in_*gpu* Param/SWPC/SATELLITES.in Param/SWPC/*.dat ${TESTDIR}/
	-Scripts/TestParam.pl -F ${TESTDIR}/PARAM.in_SWPC_gpu_init
	perl -pi -e 'if(/MaxIter|MaxBlock|Resolution/){s/700/70/; s/1500/200/; s/5000/350/; s/0\.125/2.0/}; s/#BORIS/BORIS/;' \
	     ${TESTDIR}/PARAM.in_SWPC_gpu_init
	cp ${TESTDIR}/PARAM.in_SWPC_gpu_init ${TESTDIR}/PARAM.in

test_swpc_gpu_run:
	cd ${TESTDIR}; ${MPIRUN} ./SWMF.exe | tee runlog
	cd ${TESTDIR}; ./Restart.pl -i SWMF_RESTART.20140410_000200

test_swpc_gpu_restart:
	-Scripts/TestParam.pl -F ${TESTDIR}/PARAM.in_SWPC_gpu_restart
	perl -pi -e 's/#BORIS/BORIS/' ${TESTDIR}/PARAM.in_SWPC_gpu_restart
	cd ${TESTDIR}; cp PARAM.in_SWPC_gpu_restart PARAM.in
	cd ${TESTDIR}; ${MPIRUN} ./SWMF.exe | tee runlog_restart
	cd ${TESTDIR}; ./PostProc.pl -noptec

test_swpc_gpu_check:
	-@(${DIFFNUM} -b -r=1e-5 -a=2e-8 \
		${TESTDIR}/GM/IO2/log_e20140410-000200.log \
		Param/SWPC/TestOutputGPU/log_e20140410-000200.log \
		> test_swpc_gpu_gm.diff)
	-@(${DIFFNUM} -b -r=2e-4 -a=1e-5 \
		${TESTDIR}/GM/IO2/magnetometers_e20140410-000200.mag \
		Param/SWPC/TestOutputGPU/magnetometers_e20140410-000200.mag \
		>> test_swpc_gpu_gm.diff)
	-@(${DIFFNUM} -b -r=2e-4 -a=1e-4 \
		${TESTDIR}/GM/IO2/mag_grid_e20140410-000300.out \
		Param/SWPC/TestOutputGPU/mag_grid_e20140410-000300.out.gz \
		>> test_swpc_gpu_gm.diff)
	-@(${DIFFNUM} -b -r=1e-5 -a=1e-8 \
		${TESTDIR}/GM/IO2/geoindex_e20140410-000200.log \
		Param/SWPC/TestOutputGPU/geoindex_e20140410-000200.log \
		>> test_swpc_gpu_gm.diff)
	-@(${DIFFNUM} -b -r=1e-5 \
		${TESTDIR}/IE/ionosphere/IE_t140410_000200.log \
		Param/SWPC/TestOutputGPU/IE_t140410_000200.log \
		> test_swpc_gpu_ie.diff)
	-@(${DIFFNUM} -b -r=1e-4 -a=6e-5 \
		${TESTDIR}/IE/ionosphere/it140410_000300_000.idl \
		Param/SWPC/TestOutputGPU/it140410_000300_000.idl.gz \
		>> test_swpc_gpu_ie.diff)
	ls -l test_swpc_gpu_??.diff


### test_swpc_large_gpu: GM+IE+IM and GM runs on GPU for performance check ###

test_swpc_large_gpu:
	@echo "test_swpc_large_gpu_compile..." > test_swpc_large_gpu_gm.diff
	@echo "test_swpc_large_gpu_compile..." > test_swpc_large_gpu_ie.diff
	${MAKE} test_swpc_large_gpu_compile
	@echo "test_swpc_large_gpu_rundir..." >> test_swpc_large_gpu_gm.diff
	@echo "test_swpc_large_gpu_rundir..." >> test_swpc_large_gpu_ie.diff
	${MAKE} test_swpc_large_gpu_rundir
	@echo "test_swpc_large_gpu_run..." >> test_swpc_large_gpu_gm.diff
	@echo "test_swpc_large_gpu_run..." >> test_swpc_large_gpu_ie.diff
	${MAKE} test_swpc_large_gpu_run NP=2
	@echo "test_swpc_large_gpu_check..." >> test_swpc_large_gpu_gm.diff
	@echo "test_swpc_large_gpu_check..." >> test_swpc_large_gpu_ie.diff
	${MAKE} test_swpc_large_gpu_check

test_swpc_large_gpu_compile:
	./Config.pl -v=Empty,GM/BATSRUS,IE/Ridley_serial,IM/RCM2
	./Config.pl -o=GM:u=Default,e=Mhd,ng=2,g=8,8,8,IE:g=91,181
	./Config.pl -o=GM:opt=../../Param/SWPC/PARAM.in_SWPC_large_gpu
	./Config.pl -acc
	${MAKE} SWMF
	make PIDL

test_swpc_large_gpu_rundir:
	rm -rf ${TESTDIR}
	${MAKE} rundir RUNDIR=${DIR}/${TESTDIR}
	cp Param/SWPC/PARAM.in_SWPC_large_gpu Param/SWPC/SATELLITES.in Param/SWPC/*.dat ${TESTDIR}/
	-Scripts/TestParam.pl -F ${TESTDIR}/PARAM.in_SWPC_large_gpu
	perl -pi -e 'if(/MaxIter|MaxBlock/){s/700/70/; s/1500/200/}' \
	     ${TESTDIR}/PARAM.in_SWPC_large_gpu
	cp ${TESTDIR}/PARAM.in_SWPC_large_gpu ${TESTDIR}/PARAM.in

test_swpc_large_gpu_run:
	cd ${TESTDIR}; ${MPIRUN} ./SWMF.exe | tee runlog
	cd ${TESTDIR}; ./PostProc.pl -noptec

test_swpc_large_gpu_check:
	-@(${DIFFNUM} -b -r=1e-5 -a=2e-8 \
		${TESTDIR}/GM/IO2/log_e20140410-000000.log \
		Param/SWPC/TestOutputLargeGPU/log_e20140410-000000.log \
		> test_swpc_large_gpu_gm.diff)
	-@(${DIFFNUM} -b -r=2e-4 -a=1e-5 \
		${TESTDIR}/GM/IO2/magnetometers_e20140410-000000.mag \
		Param/SWPC/TestOutputLargeGPU/magnetometers_e20140410-000000.mag \
		>> test_swpc_large_gpu_gm.diff)
	-@(${DIFFNUM} -b -r=2e-4 -a=1e-4 \
		${TESTDIR}/GM/IO2/mag_grid_e20140410-000000.out \
		Param/SWPC/TestOutputLargeGPU/mag_grid_e20140410-000000.out.gz \
		>> test_swpc_large_gpu_gm.diff)
	-@(${DIFFNUM} -b -r=1e-5 -a=1e-8 \
		${TESTDIR}/GM/IO2/geoindex_e20140410-000000.log \
		Param/SWPC/TestOutputLargeGPU/geoindex_e20140410-000000.log \
		>> test_swpc_large_gpu_gm.diff)
	-@(${DIFFNUM} -b -r=1e-5 \
		${TESTDIR}/IE/ionosphere/IE_t140410_000000.log \
		Param/SWPC/TestOutputLargeGPU/IE_t140410_000000.log \
		> test_swpc_large_gpu_ie.diff)
	-@(${DIFFNUM} -b -r=1e-4 -a=6e-5 \
		${TESTDIR}/IE/ionosphere/it140410_000100_000.idl \
		Param/SWPC/TestOutputLargeGPU/it140410_000100_000.idl.gz \
		>> test_swpc_large_gpu_ie.diff)
	ls -l test_swpc_large_gpu_??.diff

### test_swpc_multiion: GM+IE+IM with multiion MHD ###

test_swpc_multiion:
	@echo "test_swpc_multiion_compile..." > test_swpc_multiion_gm.diff
	@echo "test_swpc_multiion_compile..." > test_swpc_multiion_ie.diff
	${MAKE} test_swpc_multiion_compile
	@echo "test_swpc_multiion_rundir..." >> test_swpc_multiion_gm.diff
	@echo "test_swpc_multiion_rundir..." >> test_swpc_multiion_ie.diff
	${MAKE} test_swpc_multiion_rundir
	@echo "test_swpc_multiion_run..." >> test_swpc_multiion_gm.diff
	@echo "test_swpc_multiion_run..." >> test_swpc_multiion_ie.diff
	${MAKE} test_swpc_multiion_run
	@echo "test_swpc_multiion_restart..." >> test_swpc_multiion_gm.diff
	@echo "test_swpc_multiion_restart..." >> test_swpc_multiion_ie.diff
	${MAKE} test_swpc_multiion_restart
	@echo "test_swpc_multiion_check..." >> test_swpc_multiion_gm.diff
	@echo "test_swpc_multiion_check..." >> test_swpc_multiion_ie.diff
	${MAKE} test_swpc_multiion_check

test_swpc_multiion_compile:
	./Config.pl -v=Empty,GM/BATSRUS,IE/Ridley_serial,IM/RCM2
	./Config.pl -o=GM:u=Default,e=MultiIon,ng=2,g=8,8,8,IE:g=91,181
	${MAKE} SWMF
	make PIDL

test_swpc_multiion_rundir:
	rm -rf ${TESTDIR}
	${MAKE} rundir RUNDIR=${DIR}/${TESTDIR}
	cp Param/SWPC/*.in* Param/SWPC/*.dat ${TESTDIR}/
	-Scripts/TestParam.pl -F ${TESTDIR}/PARAM.in_multiion_v2_init
	perl -pi -e 'if(/MaxIter|MaxBlock/){s/700/70/; s/1500/200/; s/5000/350/}; s/#BORIS/BORIS/' ${TESTDIR}/PARAM.in_multiion_v2_init
	cp ${TESTDIR}/PARAM.in_multiion_v2_init ${TESTDIR}/PARAM.in

test_swpc_multiion_run:
	cd ${TESTDIR}; ${MPIRUN} ./SWMF.exe > runlog
	cd ${TESTDIR}; ./Restart.pl -i SWMF_RESTART.20140410_000200

test_swpc_multiion_restart:
	-Scripts/TestParam.pl -F ${TESTDIR}/PARAM.in_multiion_v2_restart
	perl -pi -e 's/#BORIS/BORIS/' ${TESTDIR}/PARAM.in_multiion_v2_restart
	cd ${TESTDIR}; cp PARAM.in_multiion_v2_restart PARAM.in
	cd ${TESTDIR}; ${MPIRUN} ./SWMF.exe > runlog_restart
	cd ${TESTDIR}; ./PostProc.pl -noptec

test_swpc_multiion_check:
	-@(${DIFFNUM} -b -r=1e-5 -a=1e-8 \
		${TESTDIR}/GM/IO2/log_e20140410-000200.log \
		Param/SWPC/TestOutputMultiIon/log_e20140410-000200.log \
		> test_swpc_multiion_gm.diff)
	-@(${DIFFNUM} -b -r=1e-3 -a=2e-5 \
		${TESTDIR}/GM/IO2/magnetometers_e20140410-000200.mag \
		Param/SWPC/TestOutputMultiIon/magnetometers*0200.mag \
		>> test_swpc_multiion_gm.diff)
	-@(${DIFFNUM} -b -r=2e-4 -a=1e-5 \
		${TESTDIR}/GM/IO2/mag_grid_e20140410-000300.out \
		Param/SWPC/TestOutputMultiIon/mag_grid_e20140410-000300.out \
		>> test_swpc_multiion_gm.diff)
	-@(${DIFFNUM} -b -r=1e-5 -a=1e-8 \
		${TESTDIR}/GM/IO2/geoindex_e20140410-000200.log \
		Param/SWPC/TestOutputMultiIon/geoindex_e20140410-000200.log \
		>> test_swpc_multiion_gm.diff)
	-@(${DIFFNUM} -b -r=1e-5 \
		${TESTDIR}/IE/ionosphere/IE_t140410_000200.log \
		Param/SWPC/TestOutputMultiIon/IE_t140410_000200.log \
		> test_swpc_multiion_ie.diff)
	-@(${DIFFNUM} -b -r=2e-4 -a=5e-5 \
		${TESTDIR}/IE/ionosphere/it140410_000300_000.idl \
		Param/SWPC/TestOutputMultiIon/it140410_000300_000.idl.gz \
		>> test_swpc_multiion_ie.diff)
	ls -l test_swpc_multiion_*.diff

### test_swpc_multispecies: GM+IE+IM with multispecies MHD ###

test_swpc_multispecies:
	@echo "test_swpc_multispecies_compile..." > test_swpc_multispecies_gm.diff
	@echo "test_swpc_multispecies_compile..." > test_swpc_multispecies_ie.diff
	${MAKE} test_swpc_multispecies_compile
	@echo "test_swpc_multispecies_rundir..." >> test_swpc_multispecies_gm.diff
	@echo "test_swpc_multispecies_rundir..." >> test_swpc_multispecies_ie.diff
	${MAKE} test_swpc_multispecies_rundir
	@echo "test_swpc_multispecies_run..." >> test_swpc_multispecies_gm.diff
	@echo "test_swpc_multispecies_run..." >> test_swpc_multispecies_ie.diff
	${MAKE} test_swpc_multispecies_run
	@echo "test_swpc_multispecies_restart..." >> test_swpc_multispecies_gm.diff
	@echo "test_swpc_multispecies_restart..." >> test_swpc_multispecies_ie.diff
	${MAKE} test_swpc_multispecies_restart
	@echo "test_swpc_multispecies_check..." >> test_swpc_multispecies_gm.diff
	@echo "test_swpc_multispecies_check..." >> test_swpc_multispecies_ie.diff
	${MAKE} test_swpc_multispecies_check

test_swpc_multispecies_compile:
	./Config.pl -v=Empty,GM/BATSRUS,IE/Ridley_serial,IM/RCM2
	./Config.pl -o=GM:u=Default,e=MhdHpOp,ng=2,g=8,8,8,IE:g=91,181
	${MAKE} SWMF
	make PIDL

test_swpc_multispecies_rundir:
	rm -rf ${TESTDIR}
	${MAKE} rundir RUNDIR=${DIR}/${TESTDIR}
	cp Param/SWPC/*.in* Param/SWPC/*.dat ${TESTDIR}/
	-Scripts/TestParam.pl -F ${TESTDIR}/PARAM.in_multispecies_v2_init
	perl -pi -e 'if(/MaxIter|MaxBlock/){s/700/70/; s/1500/200/; s/5000/350/}; s/#BORIS/BORIS/' ${TESTDIR}/PARAM.in_multispecies_v2_init
	cp ${TESTDIR}/PARAM.in_multispecies_v2_init ${TESTDIR}/PARAM.in

test_swpc_multispecies_run:
	cd ${TESTDIR}; ${MPIRUN} ./SWMF.exe > runlog
	cd ${TESTDIR}; ./Restart.pl -i SWMF_RESTART.20140410_000200

test_swpc_multispecies_restart:
	-Scripts/TestParam.pl -F ${TESTDIR}/PARAM.in_multispecies_v2_restart
	perl -pi -e 's/#BORIS/BORIS/' ${TESTDIR}/PARAM.in_multispecies_v2_restart
	cd ${TESTDIR}; cp PARAM.in_multispecies_v2_restart PARAM.in
	cd ${TESTDIR}; ${MPIRUN} ./SWMF.exe > runlog_restart
	cd ${TESTDIR}; ./PostProc.pl -noptec

test_swpc_multispecies_check:
	-@(${DIFFNUM} -b -r=1e-5 -a=1e-8 \
		${TESTDIR}/GM/IO2/log_e20140410-000200.log \
		Param/SWPC/TestOutputMultiSpecies/log_e20140410-000200.log \
		> test_swpc_multispecies_gm.diff)
	-@(${DIFFNUM} -b -r=1e-3 -a=1e-5 \
		${TESTDIR}/GM/IO2/magnetometers_e20140410-000200.mag \
		Param/SWPC/TestOutputMultiSpecies/magnetometers_e20140410-000200.mag \
		>> test_swpc_multispecies_gm.diff)
	-@(${DIFFNUM} -b -r=2e-4 -a=1e-5 \
		${TESTDIR}/GM/IO2/mag_grid_e20140410-000300.out \
		Param/SWPC/TestOutputMultiSpecies/mag_grid_e20140410-000300.out \
		>> test_swpc_multispecies_gm.diff)
	-@(${DIFFNUM} -b -r=1e-5 -a=1e-8 \
		${TESTDIR}/GM/IO2/geoindex_e20140410-000200.log \
		Param/SWPC/TestOutputMultiSpecies/geoindex_e20140410-000200.log \
		>> test_swpc_multispecies_gm.diff)
	-@(${DIFFNUM} -b -r=1e-5 \
		${TESTDIR}/IE/ionosphere/IE_t140410_000200.log \
		Param/SWPC/TestOutputMultiSpecies/IE_t140410_000200.log \
		> test_swpc_multispecies_ie.diff)
	-@(${DIFFNUM} -b -r=1e-4 -a=5e-5 \
		${TESTDIR}/IE/ionosphere/it140410_000300_000.idl \
		Param/SWPC/TestOutputMultiSpecies/it140410_000300_000.idl.gz \
		>> test_swpc_multispecies_ie.diff)
	ls -l test_swpc_multispecies_*.diff

### test_swpc_cimi: GM+IE+IM/CIMI2 ###

test_swpc_cimi:
	@echo "test_swpc_cimi_compile..." > test_swpc_cimi_gm.diff
	@echo "test_swpc_cimi_compile..." > test_swpc_cimi_ie.diff
	${MAKE} test_swpc_cimi_compile
	@echo "test_swpc_cimi_rundir..." >> test_swpc_cimi_gm.diff
	@echo "test_swpc_cimi_rundir..." >> test_swpc_cimi_ie.diff
	${MAKE} test_swpc_cimi_rundir
	@echo "test_swpc_cimi_run..." >> test_swpc_cimi_gm.diff
	@echo "test_swpc_cimi_run..." >> test_swpc_cimi_ie.diff
	${MAKE} test_swpc_cimi_run
	@echo "test_swpc_cimi_restart..." >> test_swpc_cimi_gm.diff
	@echo "test_swpc_cimi_restart..." >> test_swpc_cimi_ie.diff
	${MAKE} test_swpc_cimi_restart
	@echo "test_swpc_cimi_check..." >> test_swpc_cimi_gm.diff
	@echo "test_swpc_cimi_check..." >> test_swpc_cimi_ie.diff
	${MAKE} test_swpc_cimi_check

test_swpc_cimi_compile:
	./Config.pl -v=Empty,GM/BATSRUS,IE/Ridley_serial,IM/CIMI2
	./Config.pl -o=GM:u=Default,e=MhdAnisoP,ng=2,g=8,8,8,IE:g=91,181
	./Config.pl -o=IM:EarthHO,GridExpanded
	${MAKE} SWMF
	make PIDL

test_swpc_cimi_rundir:
	rm -rf ${TESTDIR}
	${MAKE} rundir RUNDIR=${DIR}/${TESTDIR}
	cp Param/SWPC/*.in* Param/SWPC/*.dat ${TESTDIR}/
	-Scripts/TestParam.pl -F ${TESTDIR}/PARAM.in_cimi_v2_init
	perl -pi -e 'if(/MaxIter|MaxBlock/){s/700/70/; s/1500/200/; s/5000/350/}; s/#BORIS/BORIS/; s/^\#(COMPONENTMAP.*production)/$$1/i; s/^(COMPONENTMAP.*nightly)/\#$$1/i' ${TESTDIR}/PARAM.in_cimi_v2_init
	cp ${TESTDIR}/PARAM.in_cimi_v2_init ${TESTDIR}/PARAM.in

test_swpc_cimi_run:
	cd ${TESTDIR}; ${MPIRUN} ./SWMF.exe > runlog
	cd ${TESTDIR}; ./Restart.pl -i SWMF_RESTART.20140410_000200

test_swpc_cimi_restart:
	-Scripts/TestParam.pl -F ${TESTDIR}/PARAM.in_cimi_v2_restart
	perl -pi -e 's/#BORIS/BORIS/; s/^\#(COMPONENTMAP.*production)/$$1/i; s/^(COMPONENTMAP.*nightly)/\#$$1/i' ${TESTDIR}/PARAM.in_cimi_v2_restart
	cd ${TESTDIR}; cp PARAM.in_cimi_v2_restart PARAM.in
	cd ${TESTDIR}; ${MPIRUN} ./SWMF.exe > runlog_restart
	cd ${TESTDIR}; ./PostProc.pl -noptec

test_swpc_cimi_check:
	-@(${DIFFNUM} -b -r=1e-5 -a=1e-8 \
		${TESTDIR}/GM/IO2/log_e20140410-000200.log \
		Param/SWPC/TestOutputCimi/log_e20140410-000200.log \
		> test_swpc_cimi_gm.diff)
	-@(${DIFFNUM} -b -r=2e-4 -a=1e-5 \
		${TESTDIR}/GM/IO2/magnetometers_e20140410-000200.mag \
		Param/SWPC/TestOutputCimi/magnetometers_e20140410-000200.mag \
		>> test_swpc_cimi_gm.diff)
	-@(${DIFFNUM} -b -r=2e-4 -a=1e-5 \
		${TESTDIR}/GM/IO2/mag_grid_e20140410-000300.out \
		Param/SWPC/TestOutputCimi/mag_grid_e20140410-000300.out \
		>> test_swpc_cimi_gm.diff)
	-@(${DIFFNUM} -b -r=1e-5 -a=1e-8 \
		${TESTDIR}/GM/IO2/geoindex_e20140410-000200.log \
		Param/SWPC/TestOutputCimi/geoindex_e20140410-000200.log \
		>> test_swpc_cimi_gm.diff)
	-@(${DIFFNUM} -b -r=1e-5 \
		${TESTDIR}/IE/ionosphere/IE_t140410_000200.log \
		Param/SWPC/TestOutputCimi/IE_t140410_000200.log \
		> test_swpc_cimi_ie.diff)
	-@(${DIFFNUM} -b -r=2e-4 -a=5e-5 \
		${TESTDIR}/IE/ionosphere/it140410_000300_000.idl \
		Param/SWPC/TestOutputCimi/it140410_000300_000.idl.gz \
		>> test_swpc_cimi_ie.diff)
	ls -l test_swpc_cimi_*.diff

### test_swpc_pwom: GM+IE+IM/RCM2+PW/PWOM ###

test_swpc_pwom:
	@echo "test_swpc_pwom_compile..." > test_swpc_pwom_gm.diff
	@echo "test_swpc_pwom_compile..." > test_swpc_pwom_ie.diff
	${MAKE} test_swpc_pwom_compile
	@echo "test_swpc_pwom_rundir..." >> test_swpc_pwom_gm.diff
	@echo "test_swpc_pwom_rundir..." >> test_swpc_pwom_ie.diff
	${MAKE} test_swpc_pwom_rundir
	@echo "test_swpc_pwom_run..." >> test_swpc_pwom_gm.diff
	@echo "test_swpc_pwom_run..." >> test_swpc_pwom_ie.diff
	${MAKE} test_swpc_pwom_run
	@echo "test_swpc_pwom_restart..." >> test_swpc_pwom_gm.diff
	@echo "test_swpc_pwom_restart..." >> test_swpc_pwom_ie.diff
	${MAKE} test_swpc_pwom_restart
	@echo "test_swpc_pwom_check..." >> test_swpc_pwom_gm.diff
	@echo "test_swpc_pwom_check..." >> test_swpc_pwom_ie.diff
	${MAKE} test_swpc_pwom_check

test_swpc_pwom_compile:
	./Config.pl -v=Empty,GM/BATSRUS,IE/Ridley_serial,IM/RCM2,PW/PWOM
	./Config.pl -o=GM:u=Default,e=Mhd,ng=2,g=8,8,8,IE:g=91,181,PW:Earth
	${MAKE} SWMF
	make PIDL

test_swpc_pwom_rundir:
	rm -rf ${TESTDIR}
	${MAKE} rundir RUNDIR=${DIR}/${TESTDIR}
	cp Param/SWPC/*.in* Param/SWPC/*.dat ${TESTDIR}/
	-Scripts/TestParam.pl -F ${TESTDIR}/PARAM.in_pwom_init
	perl -pi -e 'if(/MaxIter|MaxBlock|nTotalLine/){s/700/70/; s/1500/200/; s/5000/350/; s/252/4/;}; s/#BORIS/BORIS/; s/^\#(COMPONENTMAP.*production)/$$1/i; s/^(COMPONENTMAP.*nightly)/\#$$1/i' ${TESTDIR}/PARAM.in_pwom_init
	cp ${TESTDIR}/PARAM.in_pwom_init ${TESTDIR}/PARAM.in

test_swpc_pwom_run:
	cd ${TESTDIR}; ${MPIRUN} ./SWMF.exe > runlog
	cd ${TESTDIR}; ./Restart.pl -i SWMF_RESTART.20140410_000200

test_swpc_pwom_restart:
	-Scripts/TestParam.pl -F ${TESTDIR}/PARAM.in_pwom_restart
	perl -pi -e 's/252/4/; s/#BORIS/BORIS/; s/^\#(COMPONENTMAP.*production)/$$1/i; s/^(COMPONENTMAP.*nightly)/\#$$1/i' ${TESTDIR}/PARAM.in_pwom_restart
	cd ${TESTDIR}; cp PARAM.in_pwom_restart PARAM.in
	cd ${TESTDIR}; ${MPIRUN} ./SWMF.exe > runlog_restart
	cd ${TESTDIR}; ./PostProc.pl -noptec

test_swpc_pwom_check:
	-@(${DIFFNUM} -b -r=1e-5 -a=1e-8 \
		${TESTDIR}/GM/IO2/log_e20140410-000200.log \
		Param/SWPC/TestOutputPwom/log_e20140410-000200.log \
		> test_swpc_pwom_gm.diff)
	-@(${DIFFNUM} -b -r=2e-4 -a=1e-5 \
		${TESTDIR}/GM/IO2/magnetometers_e20140410-000200.mag \
		Param/SWPC/TestOutputPwom/magnetometers_e20140410-000200.mag \
		>> test_swpc_pwom_gm.diff)
	-@(${DIFFNUM} -b -r=2e-4 -a=1e-4 \
		${TESTDIR}/GM/IO2/mag_grid_e20140410-000300.out \
		Param/SWPC/TestOutputPwom/mag_grid_e20140410-000300.out \
		>> test_swpc_pwom_gm.diff)
	-@(${DIFFNUM} -b -r=1e-5 -a=1e-8 \
		${TESTDIR}/GM/IO2/geoindex_e20140410-000200.log \
		Param/SWPC/TestOutputPwom/geoindex_e20140410-000200.log \
		>> test_swpc_pwom_gm.diff)
	-@(${DIFFNUM} -b -r=1e-5 \
		${TESTDIR}/IE/ionosphere/IE_t140410_000200.log \
		Param/SWPC/TestOutputPwom/IE_t140410_000200.log \
		> test_swpc_pwom_ie.diff)
	-@(${DIFFNUM} -b -r=1e-4 -a=5e-5 \
		${TESTDIR}/IE/ionosphere/it140410_000300_000.idl \
		Param/SWPC/TestOutputPwom/it140410_000300_000.idl.gz \
		>> test_swpc_pwom_ie.diff)
	ls -l test_swpc_pwom_*.diff

### test_swpc_pwom_species: multispecies GM+IE+IM/RCM2+PW/PWOM ###

test_swpc_pwom_species:
	@echo "test_swpc_pwom_species_compile..." \
	      > test_swpc_pwom_species_gm.diff
	@echo "test_swpc_pwom_species_compile..." \
	      > test_swpc_pwom_species_ie.diff
	${MAKE} test_swpc_pwom_species_compile
	@echo "test_swpc_pwom_species_rundir..." \
	      >> test_swpc_pwom_species_gm.diff
	@echo "test_swpc_pwom_species_rundir..." \
	      >> test_swpc_pwom_species_ie.diff
	${MAKE} test_swpc_pwom_species_rundir
	@echo "test_swpc_pwom_species_run..." \
	      >> test_swpc_pwom_species_gm.diff
	@echo "test_swpc_pwom_species_run..." \
	      >> test_swpc_pwom_species_ie.diff
	${MAKE} test_swpc_pwom_species_run
	@echo "test_swpc_pwom_species_restart..." \
	      >> test_swpc_pwom_species_gm.diff
	@echo "test_swpc_pwom_species_restart..." \
	      >> test_swpc_pwom_species_ie.diff
	${MAKE} test_swpc_pwom_species_restart
	@echo "test_swpc_pwom_species_check..." \
	      >> test_swpc_pwom_species_gm.diff
	@echo "test_swpc_pwom_species_check..." \
	      >> test_swpc_pwom_species_ie.diff
	${MAKE} test_swpc_pwom_species_check

test_swpc_pwom_species_compile:
	./Config.pl -v=Empty,GM/BATSRUS,IE/Ridley_serial,IM/RCM2,PW/PWOM
	./Config.pl -o=GM:u=Default,e=MhdHpOp,ng=2,g=8,8,8,IE:g=91,181,PW:Earth
	${MAKE} SWMF
	make PIDL

test_swpc_pwom_species_rundir:
	rm -rf ${TESTDIR}
	${MAKE} rundir RUNDIR=${DIR}/${TESTDIR}
	cp Param/SWPC/*.in* Param/SWPC/*.dat ${TESTDIR}/
	-Scripts/TestParam.pl -F ${TESTDIR}/PARAM.in_pwom_species_init
	perl -pi -e 'if(/MaxIter|MaxBlock|nTotalLine/){s/700/70/; s/1500/200/; s/5000/350/; s/252/4/;}; s/#BORIS/BORIS/; s/^\#(COMPONENTMAP.*production)/$$1/i; s/^(COMPONENTMAP.*nightly)/\#$$1/i' ${TESTDIR}/PARAM.in_pwom_species_init
	cp ${TESTDIR}/PARAM.in_pwom_species_init ${TESTDIR}/PARAM.in

test_swpc_pwom_species_run:
	cd ${TESTDIR}; ${MPIRUN} ./SWMF.exe > runlog
	cd ${TESTDIR}; ./Restart.pl -i SWMF_RESTART.20140410_000200

test_swpc_pwom_species_restart:
	-Scripts/TestParam.pl -F ${TESTDIR}/PARAM.in_pwom_species_restart
	perl -pi -e 's/252/4/; s/#BORIS/BORIS/; s/^\#(COMPONENTMAP.*production)/$$1/i; s/^(COMPONENTMAP.*nightly)/\#$$1/i' ${TESTDIR}/PARAM.in_pwom_species_restart
	cd ${TESTDIR}; cp PARAM.in_pwom_species_restart PARAM.in
	cd ${TESTDIR}; ${MPIRUN} ./SWMF.exe > runlog_restart
	cd ${TESTDIR}; ./PostProc.pl -noptec

test_swpc_pwom_species_check:
	-@(${DIFFNUM} -b -r=1e-5 -a=1e-8 \
		${TESTDIR}/GM/IO2/log_e20140410-000200.log \
		Param/SWPC/TestOutputPwomSpecies/log_e20140410-000200.log \
		> test_swpc_pwom_species_gm.diff)
	-@(${DIFFNUM} -b -r=2e-4 -a=1e-5 \
		${TESTDIR}/GM/IO2/magnetometers_e20140410-000200.mag \
		Param/SWPC/TestOutputPwomSpecies/magnetometers_e20140410-000200.mag \
		>> test_swpc_pwom_species_gm.diff)
	-@(${DIFFNUM} -b -r=2e-4 -a=1e-4 \
		${TESTDIR}/GM/IO2/mag_grid_e20140410-000300.out \
		Param/SWPC/TestOutputPwomSpecies/mag_grid_e20140410-000300.out \
		>> test_swpc_pwom_species_gm.diff)
	-@(${DIFFNUM} -b -r=1e-5 -a=1e-8 \
		${TESTDIR}/GM/IO2/geoindex_e20140410-000200.log \
		Param/SWPC/TestOutputPwomSpecies/geoindex_e20140410-000200.log \
		>> test_swpc_pwom_species_gm.diff)
	-@(${DIFFNUM} -b -r=1e-5 \
		${TESTDIR}/IE/ionosphere/IE_t140410_000200.log \
		Param/SWPC/TestOutputPwomSpecies/IE_t140410_000200.log \
		> test_swpc_pwom_species_ie.diff)
	-@(${DIFFNUM} -b -r=1e-4 -a=5e-5 \
		${TESTDIR}/IE/ionosphere/it140410_000300_000.idl \
		Param/SWPC/TestOutputPwomSpecies/it140410_000300_000.idl.gz \
		>> test_swpc_pwom_species_ie.diff)
	ls -l test_swpc_pwom_species_*.diff

### test_swpc_amps_build: Configure SWMF with GM+IE+IM+PT+RB ###

test_swpc_amps_build:
	rm -rf ${BUILD}
	Scripts/Configure.pl -d=${BUILD} -on=PT
	cd ${BUILD}; rm -rf IM/HEIDI IM/CIMI IE/RIM util/HYPRE
	cd ${BUILD}; ./Config.pl -install \
	   	     	-compiler=${FORTRAN_COMPILER_NAME},${C_COMPILER_NAME}
	cd ${BUILD}; make test_swpc_amps

test_swpc_amps:
	@echo "test_swpc_amps_compile..." > test_swpc_amps_gm.diff
	@echo "test_swpc_amps_compile..." > test_swpc_amps_ie.diff
	@echo "test_swpc_amps_compile..." > test_swpc_amps_rb.diff
	${MAKE} test_swpc_amps_compile
	@echo "test_swpc_amps_rundir..." >> test_swpc_amps_gm.diff
	@echo "test_swpc_amps_rundir..." >> test_swpc_amps_ie.diff
	@echo "test_swpc_amps_rundir..." >> test_swpc_amps_rb.diff
	${MAKE} test_swpc_amps_rundir
	@echo "test_swpc_amps_run..." >> test_swpc_amps_gm.diff
	@echo "test_swpc_amps_run..." >> test_swpc_amps_ie.diff
	@echo "test_swpc_amps_run..." >> test_swpc_amps_rb.diff
	${MAKE} test_swpc_amps_run
	@echo "test_swpc_amps_restart..." >> test_swpc_amps_gm.diff
	@echo "test_swpc_amps_restart..." >> test_swpc_amps_ie.diff
	@echo "test_swpc_amps_restart..." >> test_swpc_amps_rb.diff
	${MAKE} test_swpc_amps_restart
	@echo "test_swpc_amps_check..." >> test_swpc_amps_gm.diff
	@echo "test_swpc_amps_check..." >> test_swpc_amps_ie.diff
	@echo "test_swpc_amps_check..." >> test_swpc_amps_rb.diff
	${MAKE} test_swpc_amps_check

test_swpc_amps_compile:
	./Config.pl -v=Empty,GM/BATSRUS,IE/Ridley_serial,IM/RCM2,RB/RBE,PT/AMPS
	./Config.pl -o=GM:u=Default,e=Mhd,ng=2,g=8,8,8,IE:g=91,181
	./Config.pl -o=PT:application=test/earth-test-swmf,spice-path=nospice,amps-test=on 
	${MAKE} SWMF
	make PIDL

test_swpc_amps_rundir:
	rm -rf ${TESTDIR}
	${MAKE} rundir RUNDIR=${DIR}/${TESTDIR}
	cp Param/SWPC/*.in* Param/SWPC/*.dat ${TESTDIR}/
	-Scripts/TestParam.pl -F ${TESTDIR}/PARAM.in_amps_init
	perl -pi -e 'if(/MaxIter/){s/700/70/; s/1500/200/}; s/#BORIS/BORIS/' \
	     ${TESTDIR}/PARAM.in_amps_init
	cp ${TESTDIR}/PARAM.in_amps_init ${TESTDIR}/PARAM.in

test_swpc_amps_run:
	cd ${TESTDIR}; ${MPIRUN} ./SWMF.exe > runlog
	cd ${TESTDIR}; ./Restart.pl -i SWMF_RESTART.20140410_000200

test_swpc_amps_restart:
	-Scripts/TestParam.pl -F ${TESTDIR}/PARAM.in_amps_restart
	perl -pi -e 's/#BORIS/BORIS/' ${TESTDIR}/PARAM.in_amps_restart
	cd ${TESTDIR}; cp PARAM.in_amps_restart PARAM.in
	cd ${TESTDIR}; ${MPIRUN} ./SWMF.exe > runlog_restart
	cd ${TESTDIR}; ./PostProc.pl -noptec

test_swpc_amps_check:
	-@(${DIFFNUM} -b -r=1e-5 -a=1e-8 \
		${TESTDIR}/GM/IO2/log_e20140410-000200.log \
		Param/SWPC/TestOutput/log_e20140410-000200.log \
		> test_swpc__amps_gm.diff)
	-@(${DIFFNUM} -b -r=2e-4 -a=1e-8 \
		${TESTDIR}/GM/IO2/magnetometers_e20140410-000200.mag \
		Param/SWPC/TestOutput/magnetometers_e20140410-000200.mag \
		>> test_swpc_amps_gm.diff)
	-@(${DIFFNUM} -b -r=2e-4 -a=1e-8 \
		${TESTDIR}/GM/IO2/mag_grid_e20140410-000300.out \
		Param/SWPC/TestOutput/mag_grid_e20140410-000300.out \
		>> test_swpc_amps_gm.diff)
	-@(${DIFFNUM} -b -r=1e-5 -a=1e-8 \
		${TESTDIR}/GM/IO2/geoindex_e20140410-000200.log \
		Param/SWPC/TestOutput/geoindex_e20140410-000200.log \
		>> test_swpc_amps_gm.diff)
	-@(${DIFFNUM} -b -r=1e-5 \
		${TESTDIR}/IE/ionosphere/IE_t140410_000200.log \
		Param/SWPC/TestOutput/IE_t140410_000200.log \
		> test_swpc_amps_ie.diff)
	-@(${DIFFNUM} -b -r=1e-4 -a=5e-5 \
		${TESTDIR}/IE/ionosphere/it140410_000300_000.idl \
		Param/SWPC/TestOutput/it140410_000300_000.idl.gz \
		>> test_swpc_amps_ie.diff)
	-@(${DIFFNUM} -b -r=1e-3 \
                ${TESTDIR}/RB/plots/20140410_000300_e.fls \
		Param/SWPC/TestOutput/RB_20140410_000300_e.fls.gz \
		> test_swpc_amps_rb.diff)
	ls -l test_swpc_amps_??.diff

### test4: GM+IE+HEIDI ###

test4:
	@echo "test4_compile..." > test04_gm.diff
	@echo "test4_compile..." > test04_ie.diff
	${MAKE} test4_compile
	@echo "test4_rundir..." >> test04_gm.diff
	@echo "test4_rundir..." >> test04_ie.diff
	${MAKE} test4_rundir
	@echo "test4_run..." >> test04_gm.diff
	@echo "test4_run..." >> test04_ie.diff
	${MAKE} test4_run
	@echo "test4_check..." >> test04_gm.diff
	@echo "test4_check..." >> test04_ie.diff
	${MAKE} test4_check

test4_compile:
	./Config.pl -v=Empty,GM/BATSRUS,IE/Ridley_serial,IM/HEIDI
	./Config.pl -o=GM:u=Default,e=Mhd,ng=2,g=8,8,8,IE:g=91,181
	${MAKE} SWMF
	make PIDL

test4_rundir:
	rm -rf ${TESTDIR}
	${MAKE} rundir RUNDIR=${DIR}/${TESTDIR}
	cp Param/PARAM.in.test.GMIEHEIDI ${TESTDIR}/PARAM.in
	-Scripts/TestParam.pl -F ${TESTDIR}/PARAM.in

test4_run:
	cd ${TESTDIR}; ${MPIRUN} ./SWMF.exe > runlog
	cd ${TESTDIR}; ./PostProc.pl -M RESULTS

test4_check:
	-@(${DIFFNUM} -b -r=1e-3 -a=1e-8 \
		${TESTDIR}/RESULTS/GM/log_n000000.log \
		output/test4/GM_log_n000000.log \
		> test04_gm.diff)
	-@(${DIFFNUM} -b -t -r=1e-5 \
		${TESTDIR}/RESULTS/IE/IE*.log \
		output/test4/IE.log \
		> test04_ie.diff)
	ls -l test04_*.diff

### test5: Eruptive Event generator coupled to the corona ###
test5:
	@echo "test5_compile..." > test05_ee.diff
	${MAKE} test5_compile
	@echo "test_rundir..." >> test05_ee.diff
	${MAKE} test5_rundir
	@echo "test5_run..." >> test05_ee.diff
	${MAKE} test5_run
	@echo "test5_check..." >> test05_ee.diff
	${MAKE} test5_check

test5_compile:
	./Config.pl -v=Empty,EE/BATSRUS,SC/BATSRUS
	./Config.pl -o=EE:u=Swarm,e=MhdEos,ng=2,g=10,10,10
	./Config.pl -o=SC:u=Awsom,e=Awsom,ng=2,g=4,4,4
	${MAKE} SWMF
	make PIDL

test5_rundir:
	rm -rf ${TESTDIR}
	${MAKE} rundir RUNDIR=${DIR}/${TESTDIR}
	cd ${TESTDIR}; cp ../GM/BATSRUS/data/FLUXEMERGENCE/*Sph.dat .
	cd ${TESTDIR}; cp ../GM/BATSRUS/data/FLUXEMERGENCE/EOS.* .
	cd ${TESTDIR}; gunzip *.gz

test5_run:
	cd ${TESTDIR}; cp Param/PARAM.in.test.EE.3D PARAM.in
	-Scripts/TestParam.pl -F ${TESTDIR}/PARAM.in
	cd ${TESTDIR}; ${MPIRUN} ./SWMF.exe > runlog_3d
	cd ${TESTDIR}; ./Restart.pl
	cd ${TESTDIR}; ./PostProc.pl -M RESULTS/run_3d
	cd ${TESTDIR}; cp Param/PARAM.in.test.EE PARAM.in
	-Scripts/TestParam.pl -F ${TESTDIR}/PARAM.in
	cd ${TESTDIR}; ${MPIRUN} ./SWMF.exe > runlog_restart
	cd ${TESTDIR}; ./PostProc.pl -M -cat RESULTS/run_restart; \
	               mv runlog_restart RESULTS/run_restart
	cd ${TESTDIR}; cp Param/PARAM.in.test.SC PARAM.in
	-Scripts/TestParam.pl -F ${TESTDIR}/PARAM.in
	cd ${TESTDIR}; ${MPIRUN} ./SWMF.exe > runlog
	cd ${TESTDIR}; ./Restart.pl
	cd ${TESTDIR}; ./PostProc.pl -M -cat RESULTS/run_SC
	cd ${TESTDIR}; cp Param/PARAM.in.test.EESC PARAM.in
	-Scripts/TestParam.pl -F ${TESTDIR}/PARAM.in
	cd ${TESTDIR}; ${MPIRUN} ./SWMF.exe > runlog
	cd ${TESTDIR}; ./PostProc.pl -M -cat RESULTS/coupled

test5_check:
	-(${DIFFNUM} -b -r=1e-5 -a=2e-15 \
		${TESTDIR}/RESULTS/coupled/EE/log_n000010.log \
		output/test5/EE.log \
		> test05_ee.diff)
	-(${DIFFNUM} -b -r=1e-5 \
		${TESTDIR}/RESULTS/coupled/SC/log_n000070.log \
		output/test5/SC.log \
		> test05_sc.diff)
	ls -l test05_*.diff


### test7: IH+OH test ###

test7:
	@echo "test7_compile..." > test07_ih.diff
	@echo "test7_compile..." > test07_oh.diff
	${MAKE} test7_compile
	@echo "test7_rundir..." >> test07_ih.diff
	@echo "test7_rundir..." >> test07_oh.diff
	${MAKE} test7_rundir
	@echo "test7_run..." >> test07_ih.diff
	@echo "test7_run..." >> test07_oh.diff
	${MAKE} test7_run
	@echo "test7_check..." >> test07_ih.diff
	@echo "test7_check..." >> test07_oh.diff
	${MAKE} test7_check

test7_compile:
	./Config.pl -v=Empty,IH/BATSRUS,OH/BATSRUS
	./Config.pl -o=IH:u=Waves,e=Mhd,ng=2,g=4,4,4 -o=OH:u=Waves,e=Mhd,ng=2,g=4,4,4
	${MAKE} SWMF
	make PIDL

test7_rundir:
	rm -rf ${TESTDIR}
	${MAKE} rundir RUNDIR=${DIR}/${TESTDIR}
	cp Param/PARAM.in.test.IHOH.CoupleSphToXyz ${TESTDIR}/PARAM.in
	-Scripts/TestParam.pl -F ${TESTDIR}/PARAM.in

test7_run:
	cd ${TESTDIR}; ${MPIRUN} ./SWMF.exe > runlog
	cd ${TESTDIR}; ./PostProc.pl -M RESULTS

test7_check:
	-(${DIFFNUM} -b -r=1e-5 -a=2e-29\
		${TESTDIR}/RESULTS/IH/log_n000000.log \
		output/test7/Sph_IH_log_n000000.log \
		> test07_ih.diff)
	-(${DIFFNUM} -b -r=1e-5 -a=1e-29\
		${TESTDIR}/RESULTS/OH/log_n000000.log \
		output/test7/Xyz_OH_log_n000000.log \
		> test07_oh.diff)
	ls -l test07_*.diff


### test8: SC + IH with Threaded-Field-Line boundary condition test ###

test8:
	@echo "test8_compile..." > test08_sc_thread.diff
	@echo "test8_compile..." > test08_ih.diff
	${MAKE} test8_compile
	@echo "test8_rundir..." >> test08_sc_thread.diff
	@echo "test8_rundir..." >> test08_ih.diff
	${MAKE} test8_rundir
	@echo "test8_run..." >> test08_sc_thread.diff
	@echo "test8_run..." >> test08_ih.diff
	${MAKE} test8_run
	@echo "test8_restart..." >> test08_sc_thread.diff
	@echo "test8_restart..." >> test08_ih.diff
	${MAKE} test8_restart
	@echo "test8_check..." >> test08_sc_thread.diff
	@echo "test8_check..." >> test08_ih.diff
	${MAKE} test8_check

test8_compile:
	./Config.pl -v=Empty,SC/BATSRUS,IH/BATSRUS
	./Config.pl -o=SC:u=Awsom,e=Awsom,ng=2,g=4,4,4
	./Config.pl -o=IH:u=Awsom,e=Awsom,ng=2,g=4,4,4
	${MAKE} SWMF
	make PIDL

test8_rundir:
	rm -rf ${TESTDIR}
	${MAKE} rundir RUNDIR=${DIR}/${TESTDIR}
	cp Param/PARAM.in.test.SCIH_threadbc ${TESTDIR}/PARAM.in
	-Scripts/TestParam.pl -F ${TESTDIR}/PARAM.in

test8_run:
	cd ${TESTDIR}; ${MPIRUN} ./SWMF.exe > runlog
	cd ${TESTDIR}; ./Restart.pl

test8_restart:
	cd ${TESTDIR}; mv PARAM.in PARAM.in_start; rm -f PARAM.in_orig_
	cd ${TESTDIR}; cp Param/PARAM.in.test.restart.SCIH_threadbc PARAM.in
	-Scripts/TestParam.pl -F ${TESTDIR}/PARAM.in
	cd ${TESTDIR}; ${MPIRUN} ./SWMF.exe > runlog_restart
	cd ${TESTDIR}; ./PostProc.pl -M -cat RESULTS

test8_check:
	-(${DIFFNUM} -b -r=1e-4 -a=1e-26 \
		${TESTDIR}/RESULTS/SC/log_n000000.log \
		output/test8/SC_log_n000000.log \
		> test08_sc_thread.diff)
	-(${DIFFNUM} -b -r=1e-5 -a=1e-26 \
		${TESTDIR}/RESULTS/SC/rfr_rwi_4_*.outs \
		output/test8/SC_rfr.ref.gz \
		>> test08_sc_thread.diff)
	-(${DIFFNUM} -b -r=1e-5 \
		${TESTDIR}/RESULTS/IH/log_n000000.log \
		output/test8/IH_log_n000000.log \
		> test08_ih.diff)
	-(${DIFFNUM} -b -r=1e-5 \
		${TESTDIR}/RESULTS/SC/los_sdo_aia*.out \
		SC/BATSRUS/data/output/test8/SC_los_sdo_aia.out.gz \
		> test08_sc_los.diff)
	ls -l test08_*.diff

### test9: SC + IH. Test 2nd and 5th order schemes ### 

test9:
	@echo "test9_compile..." > test09_sc_chromo.diff
	@echo "test9_compile..." > test09_ih.diff
	${MAKE} test9_compile
	@echo "test9_rundir..." >> test09_sc_chromo.diff
	@echo "test9_rundir..." >> test09_ih.diff
	${MAKE} test9_rundir
	@echo "test9_run..." >> test09_sc_chromo.diff
	@echo "test9_run..." >> test09_ih.diff
	${MAKE} test9_run
	@echo "test9_restart..." >> test09_sc_chromo.diff
	@echo "test9_restart..." >> test09_ih.diff
	${MAKE} test9_restart
	@echo "test9_check..." >> test09_sc_chromo.diff
	@echo "test9_check..." >> test09_ih.diff
	${MAKE} test9_check

test9_compile:
	./Config.pl -v=Empty,SC/BATSRUS,IH/BATSRUS
	./Config.pl -o=SC:u=Awsom,e=AwsomAnisoPi,ng=2,g=6,8,8
	./Config.pl -o=IH:u=Awsom,e=AwsomAnisoPi,ng=3,g=8,8,8
	${MAKE} SWMF
	make PIDL
	make FDIPS
	cd GM/BATSRUS/srcPostProc; make SPECTRUM 
	
test9_rundir:
	rm -rf ${TESTDIR}
	${MAKE} rundir RUNDIR=${DIR}/${TESTDIR}
	cp Param/PARAM.in.test.start.SCIH ${TESTDIR}/PARAM.in
	-Scripts/TestParam.pl -F ${TESTDIR}/PARAM.in
	cp GM/BATSRUS/Param/CORONA/RadCoolCorona_8.0.dat ${TESTDIR}/RadCoolCorona_8.0.dat
	cp Param/map_1.out ${TESTDIR}/SC/
	cp -f ${MAGNETOGRAMDIR}/FDIPS.in.orig ${TESTDIR}/SC/FDIPS.in
	perl -i -pe 's/dipole11.out/map_1.out/' ${TESTDIR}/SC/FDIPS.in
	perl -i -pe 's/#(CHANGEPOLARFIELD)/$$1/; s/20/90/ if /n(Theta|Phi)/' ${TESTDIR}/SC/FDIPS.in
	cd ${TESTDIR}/SC; ${PARALLEL} ${NPFLAG} 4 ./FDIPS.exe > runlog_fdips
	gunzip -c GM/BATSRUS/data/SPECTRUM/SPECTRUM_chianti_tbl.dat.gz \
        > ${TESTDIR}/SPECTRUM_chianti_tbl.dat
	cp -f Param/SPECTRUM.in.test ${TESTDIR}/SPECTRUM.in	
	ln -s ${BINDIR}/SPECTRUM.exe ${TESTDIR}/SPECTRUM.exe
		
test9_run:
	cd ${TESTDIR}; ${MPIRUN} ./SWMF.exe > runlog
	cd ${TESTDIR}; ./PostProc.pl -m RESULTS/run_start
	cd ${TESTDIR}; ./Restart.pl -i RESULTS/run_start/RESTART
	cd ${TESTDIR}; 	${MPIRUN} ./SPECTRUM.exe > spectrum.log
	
test9_restart:
	cd ${TESTDIR}; mv PARAM.in PARAM.in_start; rm -f PARAM.in_orig_
	cd ${TESTDIR}; cp Param/PARAM.in.test.restart.SCIH PARAM.in
	-Scripts/TestParam.pl -F ${TESTDIR}/PARAM.in
	cd ${TESTDIR}; ${MPIRUN} ./SWMF.exe > runlog
	cd ${TESTDIR}; ./PostProc.pl -m RESULTS/run_restart

test9_check:
	-(${DIFFNUM} -b -r=1e-5 -a=1e-27 \
		${TESTDIR}/RESULTS/run_start/SC/log_n000001.log \
		output/test9/SC_log_n000001.log \
		> test09_sc_chromo.diff)
	-(${DIFFNUM} -b -r=1e-5 \
		${TESTDIR}/RESULTS/run_restart/IH/log_n000005.log \
		output/test9/IH_log_n000005.log \
		> test09_ih.diff)
	-(${DIFFNUM} -b -r=2e-5 -a=1e-27 \
		${TESTDIR}/RESULTS/run_restart/IH/trj_earth_t00000010_n00000007.sat \
		output/test9/IH_trj_earth_t00000010_n00000007.sat \
		> test09_ih_trj_earth.diff)
	-(${DIFFNUM} -b -r=1e-4 -a=1e-27 \
                ${TESTDIR}/RESULTS/run_start/SC/los_sta_euvi_5_n0000005.out \
                output/test9/SC_los_sta_euvi_5_n0000005.out.gz \
                > test09_sc_los.diff)
	-(${DIFFNUM} -b -r=1e-4 -a=1e-27 \
                ${TESTDIR}/RESULTS/run_start/SC/los_stb_euvi_6_n0000005.out \
                output/test9/SC_los_stb_euvi_6_n0000005.out.gz \
                > test09_sc_los.diff)
	-(${DIFFNUM} -b -r=1e-5 -a=1e-12 \
                ${TESTDIR}/RESULTS/spectrum.out \
                output/test9/spectrum.out.gz \
		> test09_spectrum.diff)

	ls -l test09_*.diff


### test11: Dipole field with solar wind coupled with an exosphere model, to test the GM/PT coupler ###

test11:
	@echo "test11_compile..." > test11_gm.diff
	@echo "test11_compile..." > test11_pt.diff
	${MAKE} test11_compile
	@echo "test11_rundir..." >> test11_gm.diff
	@echo "test11_rundir..." >> test11_pt.diff
	${MAKE} test11_rundir
	@echo "test11_run..." >> test11_gm.diff
	@echo "test11_run..." >> test11_pt.diff
	${MAKE} test11_run
	@echo "test11_check..." >> test11_gm.diff
	@echo "test11_check..." >> test11_pt.diff
	${MAKE} test11_check

test11_compile:
	./Config.pl -v=Empty,GM/BATSRUS,PT/AMPS
	./Config.pl -o=GM:u=Default,e=Mhd,ng=2,g=8,8,8
	./Config.pl -o=PT:spice-path=nospice,application=coupler-test
	${MAKE} SWMF
	make PIDL

test11_rundir:
	rm -rf ${TESTDIR}
	${MAKE} rundir RUNDIR=${DIR}/${TESTDIR}
	cp Param/PARAM.in.test.GMPT ${TESTDIR}/PARAM.in
	-Scripts/TestParam.pl -F ${TESTDIR}/PARAM.in

test11_run:
	cd ${TESTDIR}; ${MPIRUN} ./SWMF.exe > runlog
	cd ${TESTDIR}; ./PostProc.pl -M RESULTS

# Note: skip columns that depend on number of processors
#       multiply density column by 1e20 to make it order unity
test11_check:
	-(${DIFFNUM} -b -a=1e-27 \
	        ${TESTDIR}/RESULTS/GM/log_n000001.log \
		output/test11/log_n000001.log \
		> test11_gm.diff)
	gunzip -c output/test11/pic.NA.s=0.out=19.dat.gz \
	       	> ${TESTDIR}/RESULTS/PT/pic.NA.ref
	cd ${TESTDIR}/RESULTS/PT; \
	        cp -f pic.NA.s=0.out=19.dat pic.NA.dat; \
	        perl -ai -ne 'print "@F[0..3] @F[6..14] ",1e20*$$F[15]," @F[16..$$#F]\n"' \
			pic.NA.ref pic.NA.dat 
	-(${DIFFNUM} -a=1e-19 -r=1e-6 -b \
	        ${TESTDIR}/RESULTS/PT/pic.NA.dat \
		${TESTDIR}/RESULTS/PT/pic.NA.ref \
                > test11_pt.diff)
		ls -l test11*.diff

### test13: Europa's exosphere model with GM and PT ###

test13:
	@echo "test13_compile..." > test13_pt.diff
	${MAKE} test13_compile
	@echo "test13_rundir..." >> test13_pt.diff
	${MAKE} test13_rundir
	@echo "test13_run..." >> test13_pt.diff
	${MAKE} test13_run
	@echo "test13_check..." >> test13_pt.diff
	${MAKE} test13_check

test13_compile:
	./Config.pl -v=Empty,GM/BATSRUS,PT/AMPS
	./Config.pl -o=GM:u=EuropaSingleMHDPe,e=MhdPe,ng=2,g=8,8,8
	./Config.pl -o=PT:spice-path=nospice,spice-kernels=nospice,application=europa-test,amps-test=on 
	${MAKE} SWMF

test13_rundir:
	rm -rf ${TESTDIR}
	${MAKE} rundir RUNDIR=${DIR}/${TESTDIR}
	cp Param/PARAM.in.test.GMPT.Europa ${TESTDIR}/PARAM.in
	cp -r PT/AMPS/data/input/Europa/TestCase/RESTART_n003000 ${TESTDIR}/
	cd ${TESTDIR}; ./Restart.pl -i RESTART_n003000
	-Scripts/TestParam.pl -F ${TESTDIR}/PARAM.in

test13_run:
	cd ${TESTDIR}; ${MPIRUN} ./SWMF.exe > runlog

test13_check:
	-(${DIFFNUM} -b \
		${TESTDIR}/PT/plots/amps.dat \
		output/test13/amps.ref_np${NAMPS} \
		> test13_pt.diff)
		ls -l test13*.diff

### test14: OH/BATSRUS (single-ion) coupled with PT/AMPS ###

test14:
	@echo "test14_compile..." > test14_ohpt.diff
	${MAKE} test14_compile
	@echo "test14_rundir..." >> test14_ohpt.diff
	${MAKE} test14_rundir
	@echo "test14_run..." >> test14_ohpt.diff
	${MAKE} test14_run
	@echo "test14_restart..." >> test14_ohpt.diff
	${MAKE} test14_restart
	@echo "test14_check..." >> test14_ohpt.diff
	${MAKE} test14_check

test14_compile:
	./Config.pl -v=Empty,OH/BATSRUS,PT/AMPS
	./Config.pl -o=OH:u=OuterHelio,e=Swh,ng=2,g=4,4,4
	./Config.pl -o=PT:spice-path=nospice,spice-kernels=nospice,application=oh,amps-test=on,no-signals 
	${MAKE} SWMF
	make PIDL

test14_rundir:
	rm -rf ${TESTDIR}
	${MAKE} rundir RUNDIR=${DIR}/${TESTDIR}

test14_run:
	cp Param/PARAM.in.test.OHPT ${TESTDIR}/PARAM.in
	-Scripts/TestParam.pl -F ${TESTDIR}/PARAM.in
	cd ${TESTDIR}; ${MPIRUN} ./SWMF.exe > runlog_start
	cd ${TESTDIR}; rm -rf RESTART_n000120; ./Restart.pl
	cd ${TESTDIR}; mv PARAM.in PARAM.in.start; rm -f PARAM.in_orig

test14_restart:
	cd ${TESTDIR}; cp Param/PARAM.in.test.restart.OHPT PARAM.in
	-Scripts/TestParam.pl -F ${TESTDIR}/PARAM.in
	cd ${TESTDIR}; ${MPIRUN} ./SWMF.exe > runlog
	cd ${TESTDIR}; rm -rf RESULTS; ./PostProc.pl -M -cat RESULTS

test14_check:
	-(${DIFFNUM} -b -r=1e-5 -a=1e-29\
		${TESTDIR}/RESULTS/OH/log_n000000.log \
		output/test14/OH_log_n000000.log_np${NAMPS} \
		> test14_ohpt.diff)
	ls -l test14_ohpt.diff

### test14_1d: OH/BATSRUS (single-ion) coupled with PT/AMPS in a shocktube ###

test14_1d:
	@echo "test14_1d_compile..." > test14_1d_ohpt.diff
	${MAKE} test14_1d_compile
	@echo "test14_1d_rundir..." >> test14_1d_ohpt.diff
	${MAKE} test14_1d_rundir
	@echo "test14_1d_run..." >> test14_1d_ohpt.diff
	${MAKE} test14_1d_run
	@echo "test14_1d_check..." >> test14_1d_ohpt.diff
	${MAKE} test14_1d_check

test14_1d_compile:
	./Config.pl -v=Empty,OH/BATSRUS,PT/AMPS
	./Config.pl -o=OH:u=OuterHelio,e=Swh,ng=2,g=4,4,4
	./Config.pl -o=PT:spice-path=nospice,spice-kernels=nospice,application=test/oh.shocktube,amps-test=on
	${MAKE} SWMF
	make PIDL

test14_1d_rundir:
	rm -rf ${TESTDIR}
	${MAKE} rundir RUNDIR=${DIR}/${TESTDIR}
	cp Param/PARAM.in.test.OHPTshocktube ${TESTDIR}/PARAM.in
	-Scripts/TestParam.pl -F ${TESTDIR}/PARAM.in

test14_1d_run:
	cd ${TESTDIR}; ${MPIRUN} ./SWMF.exe > runlog
	cd ${TESTDIR}; rm -rf RESULTS; ./PostProc.pl -M RESULTS

test14_1d_check:
	-(${DIFFNUM} -b -r=1e-3 -a=1e-21\
                ${TESTDIR}/RESULTS/OH/log_n000000.log \
                output/test14_1d/OH_log_n000000.log_np${NAMPS} \
		> test14_1d_ohpt.diff)
	-(${DIFFNUM} -b -r=1e-3 \
                ${TESTDIR}/RESULTS/PT/amps.dat \
                output/test14_1d/amps.dat_np${NAMPS} \
		>> test14_1d_ohpt.diff)	
	ls -l test14_1d_ohpt.diff

### test14pui: OH/BATSRUS with PUI (multi-ion) coupled with PT/AMPS ###

test14pui:
	@echo "test14pui_compile..." > test14pui_ohpt.diff
	${MAKE} test14pui_compile
	@echo "test14pui_rundir..." >> test14pui_ohpt.diff
	${MAKE} test14pui_rundir
	@echo "test14pui_run..." >> test14pui_ohpt.diff
	${MAKE} test14pui_run
	@echo "test14pui_check..." >> test14pui_ohpt.diff
	${MAKE} test14pui_check

test14pui_compile:
	./Config.pl -v=Empty,OH/BATSRUS,PT/AMPS
	./Config.pl -o=OH:u=OuterHelio,e=SwhPui,ng=2,g=4,4,4
	./Config.pl -o=PT:spice-path=nospice,spice-kernels=nospice,application=oh,amps-test=on 
	${MAKE} SWMF
	make PIDL

test14pui_rundir:
	rm -rf ${TESTDIR}
	${MAKE} rundir RUNDIR=${DIR}/${TESTDIR}
	cp Param/PARAM.in.test.OHPTpui ${TESTDIR}/PARAM.in
	-Scripts/TestParam.pl -F ${TESTDIR}/PARAM.in

test14pui_run:
	cd ${TESTDIR}; ${MPIRUN} ./SWMF.exe > runlog
	cd ${TESTDIR}; rm -rf RESULTS; ./PostProc.pl -M RESULTS

test14pui_check:
	-(${DIFFNUM} -b -r=1e-5 -a=2e-21 \
		${TESTDIR}/RESULTS/OH/log_n000000.log \
		output/test14pui/OH_log_n000000_np${NAMPS}.log \
		> test14pui_ohpt.diff)
	ls -l test14pui_ohpt.diff

### test14pui1d: OH/BATSRUS with PUI (multi-ion) coupled with PT/AMPS in shock tube ###

test14pui1d:
	@echo "test14pui1d_compile..." > test14pui1d_ohpt.diff
	${MAKE} test14pui1d_compile
	@echo "test14pui1d_rundir..." >> test14pui1d_ohpt.diff
	${MAKE} test14pui1d_rundir
	@echo "test14pui1d_run..." >> test14pui1d_ohpt.diff
	${MAKE} test14pui1d_run
	@echo "test14pui1d_check..." >> test14pui1d_ohpt.diff
	${MAKE} test14pui1d_check

test14pui1d_compile:
	./Config.pl -v=Empty,OH/BATSRUS,PT/AMPS
	./Config.pl -o=OH:u=OuterHelio,e=SwhPui,ng=2,g=4,4,4
	./Config.pl -o=PT:spice-path=nospice,spice-kernels=nospice,application=test/oh.shocktube,amps-test=on
	${MAKE} SWMF
	make PIDL

test14pui1d_rundir:
	rm -rf ${TESTDIR}
	${MAKE} rundir RUNDIR=${DIR}/${TESTDIR}
	cp Param/PARAM.in.test.OHPTpuishocktube ${TESTDIR}/PARAM.in
	-Scripts/TestParam.pl -F ${TESTDIR}/PARAM.in

test14pui1d_run:
	cd ${TESTDIR}; ${MPIRUN} ./SWMF.exe > runlog
	cd ${TESTDIR}; rm -rf RESULTS; ./PostProc.pl -M RESULTS

test14pui1d_check:
	-(${DIFFNUM} -b -r=1e-3 -a=1e-21 \
                ${TESTDIR}/RESULTS/OH/log_n000000.log \
                output/test14pui1d/OH_log_n000000.log_np${NAMPS} \
                > test14pui1d_ohpt.diff)
	ls -l test14pui1d_ohpt.diff

### test15: SP/MFLAMPA model coupled with SC and IH models ###

test15_build:
	rm -rf SP_build
	Scripts/Configure.pl -on=SP,SC,IH -off=IM,IE,RB -d=SP_build
	cd SP_build; ./Config.pl -install \
	   	     	-compiler=${FORTRAN_COMPILER_NAME},${C_COMPILER_NAME}
	cd SP_build; make test15

test15:
	rm -f test15_scihsp.diff
	@echo "test15_compile..." > test15_scihsp.diff
	${MAKE} test15_compile
	@echo "test15_rundir..." >> test15_scihsp.diff
	${MAKE} test15_rundir
	@echo "test15_run..." >> test15_scihsp.diff
	${MAKE} test15_run
	@echo "test15_restart..." >> test15_scihsp.diff
	${MAKE} test15_restart
	@echo "test15_check..." >> test15_scihsp.diff
	${MAKE} test15_check

test15_compile:
	./Config.pl -v=Empty,SC/BATSRUS,IH/BATSRUS,SP/MFLAMPA
	./Config.pl -o=SC:u=Awsom,e=Awsom,ng=2,g=4,4,4
	./Config.pl -o=IH:u=Awsom,e=Awsom,ng=2,g=4,4,4
	./Config.pl -o=SP:g=2000
	${MAKE} SWMF
	make PIDL

test15_rundir:
	rm -rf ${TESTDIR}
	${MAKE} rundir RUNDIR=${DIR}/${TESTDIR}
	cp Param/PARAM.in.test.SCIHSP ${TESTDIR}/PARAM.in
	cp GM/BATSRUS/Param/CORONA/RadCoolCorona.dat ${TESTDIR}/RadCoolCorona.dat
	cp -r SP/MFLAMPA/data/input/test15/RESTART_t0020.0000s ${TESTDIR}/
	cd ${TESTDIR}; ./Restart.pl -i RESTART_t0020.0000s
	-Scripts/TestParam.pl -F ${TESTDIR}/PARAM.in

test15_run:
	cd ${TESTDIR}; rm -rf RESTART_t0025.0000s
	cd ${TESTDIR}; ${MPIRUN} ./SWMF.exe > runlog
	cd ${TESTDIR}; ./Restart.pl

test15_restart:
	cd ${TESTDIR}; mv PARAM.in PARAM.in_start; rm -f PARAM.in_orig_
	cd ${TESTDIR}; cp Param/PARAM.in.test.restart.SCIHSP PARAM.in
	-Scripts/TestParam.pl -F ${TESTDIR}/PARAM.in
	cd ${TESTDIR}; ${MPIRUN} ./SWMF.exe > runlog_restart
	cd ${TESTDIR}; ./PostProc.pl -M RESULTS

test15_check:
	cat ${TESTDIR}/RESULTS/SP/MH_data_*_*_t00000030_n000010.out \
	    > ${TESTDIR}/RESULTS/SP/MH_data.outs
	-(${DIFFNUM} -a=1e-6 -r=1e-6 -b \
		${TESTDIR}/RESULTS/SP/MH_data.outs \
		output/test15/MH_data.ref.gz \
		> test15_scihsp.diff)
	ls -l test15_scihsp.diff

### test15_long: same as test15, but SC/IH start from scratch ###

test15_long:
	rm -f test15_long*.diff
	@echo "test15_compile..." >> test15_long_sp.diff
	${MAKE} test15_compile
	@echo "test15_long_rundir..." >> test15_long_sp.diff
	${MAKE} test15_long_rundir
	@echo "test15_long_run..." >> test15_long_sp.diff
	${MAKE} test15_long_run
	@echo "test15_long_restart..." >> test15_long_sp.diff
	${MAKE} test15_long_restart
	@echo "test15_long_check..." >> test15_long_sp.diff
	${MAKE} test15_long_check


test15_long_rundir:
	rm -rf ${TESTDIR}
	${MAKE} rundir RUNDIR=${DIR}/${TESTDIR}
	cp Param/PARAM.in.test.SCIHSP_long ${TESTDIR}/PARAM.in
	-Scripts/TestParam.pl -F ${TESTDIR}/PARAM.in

test15_long_run:
	cd ${TESTDIR}; rm -rf RESTART
	cd ${TESTDIR}; ${MPIRUN} ./SWMF.exe > runlog
	cd ${TESTDIR}; ./Restart.pl -o -W RESTART
	
test15_long_restart:
	cd ${TESTDIR}; rm -rf RESULTS runlog_restart
	cd ${TESTDIR}; ./Restart.pl -i RESTART
	cd ${TESTDIR}; mv PARAM.in PARAM.in_start; rm -f PARAM.in_orig_
	cd ${TESTDIR}; cp Param/PARAM.in.test.restart.SCIHSP_long PARAM.in
	-Scripts/TestParam.pl -F ${TESTDIR}/PARAM.in
	cd ${TESTDIR}; ${MPIRUN} ./SWMF.exe > runlog_restart
	cd ${TESTDIR}; ./PostProc.pl -M -cat RESULTS

test15_long_check:
	-(${DIFFNUM} -b -r=1e-5 -a=1e-26 \
		${TESTDIR}/RESULTS/SC/log_n000000.log \
		output/test15_long/SC_log_n000000.log \
		> test15_long_sc.diff)
	-(${DIFFNUM} -b -r=1e-5 -a=1e-26 \
		${TESTDIR}/RESULTS/SC/rfr_rwi_4_*.outs \
		output/test8/SC_rfr.ref.gz \
		> test15_long_sc_rfr.diff)
	-(${DIFFNUM} -b -r=1e-5 \
		${TESTDIR}/RESULTS/IH/log_n000000.log \
		output/test15_long/IH_log_n000000.log \
		> test15_long_ih.diff)
	-(${DIFFNUM} -b -r=1e-5 \
		${TESTDIR}/RESULTS/SC/los_sdo_aia*.out \
		SC/BATSRUS/data/output/test8/SC_los_sdo_aia.out.gz \
		> test15_long_sc_los.diff)
	cat ${TESTDIR}/RESULTS/SP/MH_data_*_*_t00000150_n000010.out \
	    > ${TESTDIR}/RESULTS/SP/MH_data.outs
	-(${DIFFNUM} -a=1e-6 -r=1e-6 -b \
		${TESTDIR}/RESULTS/SP/MH_data.outs \
		output/test15_long/MH_data.ref.gz \
		> test15_long_sp.diff)
	ls -l test15_long*.diff

### test15_single: same as test15_long, but with no restart ###

test15_single:
	rm -f test15_long*.diff
	@echo "test15_compile..." >> test15_long_sp.diff
	${MAKE} test15_compile
	@echo "test15_single_rundir..." >> test15_long_sp.diff
	${MAKE} test15_single_rundir
	@echo "test15_single_run..." >> test15_long_sp.diff
	${MAKE} test15_single_run
	@echo "test15_long_check..." >> test15_long_sp.diff
	${MAKE} test15_long_check

test15_single_rundir:
	rm -rf ${TESTDIR}
	${MAKE} rundir RUNDIR=${DIR}/${TESTDIR}
	share/Scripts/ParamConvert.pl \
	Param/PARAM.in.test.SCIHSP_single ${TESTDIR}/PARAM.in
	-Scripts/TestParam.pl -F ${TESTDIR}/PARAM.in

test15_single_run:
	cd ${TESTDIR}; ${MPIRUN} ./SWMF.exe > runlog
	cd ${TESTDIR}; ./PostProc.pl -M RESULTS


### test15_pt: same as test15, but PT/AMPS coupled with SC and IH ###

test15_pt:
	rm -f test15_scihpt.diff
	@echo "test15_pt_compile..." > test15_scihpt.diff
	${MAKE} test15_pt_compile
	@echo "test15_pt_rundir..." >> test15_scihpt.diff
	${MAKE} test15_pt_rundir
	@echo "test15_pt_run..." >> test15_scihpt.diff
	${MAKE} test15_pt_run
	@echo "test15_pt_check..." >> test15_scihpt.diff
	${MAKE} test15_pt_check

test15_pt_compile:
	./Config.pl -v=Empty,SC/BATSRUS,IH/BATSRUS,PT/AMPS
	./Config.pl -o=SC:u=Awsom,e=Awsom,ng=2,g=4,4,4
	./Config.pl -o=IH:u=Awsom,e=Awsom,ng=2,g=4,4,4
	./Config.pl -o=PT:spice-path=nospice,spice-kernels=nospice,\
	application=test/sep_focused_transport-ih-mflampa-fl,amps-test=on
	${MAKE} SWMF
	make PIDL

test15_pt_rundir:
	rm -rf ${TESTDIR}
	${MAKE} rundir RUNDIR=${DIR}/${TESTDIR}
	cp Param/PARAM.in.test.SCIHPT ${TESTDIR}/PARAM.in
	cp GM/BATSRUS/Param/CORONA/RadCoolCorona.dat ${TESTDIR}/RadCoolCorona.dat
	cp -r SP/MFLAMPA/data/input/test15/RESTART_t0020.0000s ${TESTDIR}/
	cd ${TESTDIR}; ./Restart.pl -i RESTART_t0020.0000s
	-Scripts/TestParam.pl -F ${TESTDIR}/PARAM.in

test15_pt_run:
	cd ${TESTDIR}; ${MPIRUN} ./SWMF.exe > runlog
	cd ${TESTDIR}; ./PostProc.pl -M RESULTS

test15_pt_check:
	cat ${TESTDIR}/RESULTS/PT/MH_data_*_*_t00000030.out \
	    > ${TESTDIR}/RESULTS/PT/MH_data.outs
	-(${DIFFNUM} -a=1e-6 -r=1e-6 -b \
		${TESTDIR}/RESULTS/PT/MH_data.outs \
		output/test15_pt/MH_data.ref.gz \
		> test15_scihpt.diff)
	ls -l test15_scihpt.diff

# test for MFLAMPA inside the SWMF

test_sp:
	@echo "test_sp_compile..." > test_sp.diff
	${MAKE} test_sp_compile
	@echo "test_sp_rundir..." >> test_sp.diff
	${MAKE} test_sp_rundir
	@echo "test_sp_run..." >> test_sp.diff
	${MAKE} test_sp_run
	@echo "test_sp_check..." >> test_sp.diff
	${MAKE} test_sp_check

test_sp_compile:
	./Config.pl -v=Empty,SP/MFLAMPA
	./Config.pl -o=SP:g=2000
	${MAKE} SWMF
	make PIDL

test_sp_rundir:
	rm -rf ${TESTDIR}
	${MAKE} rundir RUNDIR=${DIR}/${TESTDIR}
	mkdir ${TESTDIR}/MH_data/
	cp Param/PARAM.in.test.SP ${TESTDIR}/PARAM.in
	cd ${TESTDIR}; tar xzf ../SP/MFLAMPA/data/input/test15/MH_data.tgz
	-Scripts/TestParam.pl -F ${TESTDIR}/PARAM.in

test_sp_run:
	cd ${TESTDIR}; ${MPIRUN} ./SWMF.exe > runlog
	cd ${TESTDIR}; ./PostProc.pl -M RESULTS

test_sp_check:
	cat ${TESTDIR}/RESULTS/SP/MH_data_*_*_t00000030_n000010.out \
	    > ${TESTDIR}/RESULTS/SP/MH_data.outs
	-(${DIFFNUM} -a=1e-6 -r=1e-6 -b \
		${TESTDIR}/RESULTS/SP/MH_data.outs \
		SP/MFLAMPA/output/MH_data.outs.ref.gz \
		> test_sp.diff)
	ls -l test_sp.diff


### test16: GM/BATSRUS+PC/FLEKS. BATSRUS only provides the initial condition.
###         FLEKS runs with periodic boundary conditions. 

test16:
	@echo "test16_compile..." > test16.diff
	$(MAKE) test16_compile
	@echo "test16_rundir..." > test16.diff
	$(MAKE) test16_rundir
	@echo "test16_run..." > test16.diff
	$(MAKE) test16_run
	@echo "test16_check..." > test16.diff
	$(MAKE) test16_check

test16_compile:
	./Config.pl -v=Empty,PC/FLEKS,GM/BATSRUS -amrex
	./Config.pl -o=GM:u=Waves,e=MultiIonPe,ng=2,g=8,8,1
	$(MAKE) SWMF
	$(MAKE) PIDL
	
test16_rundir:
	rm -rf ${TESTDIR}
	$(MAKE) rundir RUNDIR=${DIR}/${TESTDIR} COMPDIR=GM,PC
	cp Param/PARAM.in.test.GMPC.FLEKS.periodic  ${TESTDIR}/PARAM.in
	-Scripts/TestParam.pl -F ${TESTDIR}/PARAM.in

test16_run:
	cd ${TESTDIR}; ${MPIRUN} ./SWMF.exe > runlog
	cd ${TESTDIR}; ./PostProc.pl RESULTS

test16_check:
	-@(${DIFFNUM} -b -r=1e-7 \
	   ${TESTDIR}/RESULTS/PC/z=0_fluid_region0_0_t00000002_n00000007.out \
	   output/test16/z=0.ref.gz \
		> test16.diff)
	-@(${DIFFNUM} -b -r=1e-12 \
	   ${TESTDIR}/RESULTS/PC/log_n00000000.log \
	   output/test16/PC.log \
		>> test16.diff)		
	ls -l test16.diff


### test17: GM/BATSRUS+PC/FLEKS. 2D Alfven wave ###

test17:
	@echo "test17_compile..." > test17.diff
	$(MAKE) test17_compile
	@echo "test17_rundir..." > test17.diff
	$(MAKE) test17_rundir
	@echo "test17_run..." > test17.diff
	$(MAKE) test17_run
	@echo "test17_check..." > test17.diff
	$(MAKE) test17_check

test17_compile:
	./Config.pl -v=Empty,PC/FLEKS,GM/BATSRUS -amrex
	./Config.pl -o=GM:u=Waves,e=MhdAnisoP,ng=2,g=8,8,1
	$(MAKE) SWMF
	$(MAKE) PIDL
	
test17_rundir:
	rm -rf ${TESTDIR}
	$(MAKE) rundir RUNDIR=${DIR}/${TESTDIR} COMPDIR=GM,PC
	cp Param/PARAM.in.test.GMPC.FLEKS.2D.2steps  ${TESTDIR}/PARAM.in
	-Scripts/TestParam.pl -F ${TESTDIR}/PARAM.in

test17_run:
	cd ${TESTDIR}; ${MPIRUN} ./SWMF.exe > runlog
	cd ${TESTDIR}; ./PostProc.pl RESULTS/2step
	cd ${TESTDIR}; ./Restart.pl -v -i RESULTS/2step/RESTART
	cp Param/PARAM.in.test.GMPC.FLEKS.2D.restart ${TESTDIR}/PARAM.in
	-Scripts/TestParam.pl -F ${TESTDIR}/PARAM.in
	cd ${TESTDIR}; ${MPIRUN} ./SWMF.exe > runlog
	cd ${TESTDIR}; ./PostProc.pl RESULTS/2restart
	cp Param/PARAM.in.test.GMPC.FLEKS.2D.4steps  ${TESTDIR}/PARAM.in
	-Scripts/TestParam.pl -F ${TESTDIR}/PARAM.in
	cd ${TESTDIR}; ${MPIRUN} ./SWMF.exe > runlog
	cd ${TESTDIR}; ./PostProc.pl RESULTS/4steps
	cp Param/PARAM.in.test.GMPC.FLEKS.2D.adapt  ${TESTDIR}/PARAM.in
	-Scripts/TestParam.pl -F ${TESTDIR}/PARAM.in
	cd ${TESTDIR}; ${MPIRUN} ./SWMF.exe > runlog
	cd ${TESTDIR}; ./PostProc.pl RESULTS/adapt
	

test17_check:
	-@(${DIFFNUM} -b -r=1e-7 -a=1e-14 \
	   ${TESTDIR}/RESULTS/2restart/PC/z=0_fluid_region0_0_t00000005_n00000004.out \
	   output/test17/z=0_FLEKS0.ref.gz \
		> test17.diff)
	-@(${DIFFNUM} -b -r=1e-7 -a=1e-14 \
	   ${TESTDIR}/RESULTS/2restart/PC/z=0_fluid_region1_0_t00000004_n00000007.out \
	   output/test17/z=0_FLEKS1.ref.gz \
		>> test17.diff)
	-@(${DIFFNUM} -b -r=1e-7 -a=1e-14 \
	   ${TESTDIR}/RESULTS/4steps/PC/z=0_fluid_region0_0_t00000005_n00000004.out \
	   output/test17/z=0_FLEKS0.ref.gz \
		>> test17.diff)
	-@(${DIFFNUM} -b -r=1e-7 -a=1e-14 \
	   ${TESTDIR}/RESULTS/4steps/PC/z=0_fluid_region1_0_t00000004_n00000007.out \
	   output/test17/z=0_FLEKS1.ref.gz \
		>> test17.diff)					
	-@(${DIFFNUM} -b -r=1e-7 -a=1e-14 \
	   ${TESTDIR}/RESULTS/2restart/GM/z=0_var_1_t00000005_n00000004.out \
	   output/test17/z=0_GM.ref.gz \
		>> test17.diff)
	-@(${DIFFNUM} -b -r=1e-7 -a=1e-14 \
	   ${TESTDIR}/RESULTS/4steps/GM/z=0_var_1_t00000005_n00000004.out \
	   output/test17/z=0_GM.ref.gz \
		>> test17.diff)
	-@(${DIFFNUM} -b -r=1e-7 -a=1e-14 \
           ${TESTDIR}/RESULTS/adapt/PC/z=0_fluid_region0_0_t00000010_n00000008.out \
           output/test17/z=0_FLEKS_adapt.ref.gz \
                >> test17.diff)
	-@(${DIFFNUM} -b -r=1e-7 -a=1e-14 \
           ${TESTDIR}/RESULTS/adapt/GM/z=0_var_1_t00000010_n00000008.out \
           output/test17/z=0_GM_adapt.ref.gz \
                >> test17.diff)

	ls -l test17.diff


### test18: BATSRUS+FLEKS. 3D reconnection ###

test18:
	@echo "test18_compile..." > test18.diff
	$(MAKE) test18_compile
	@echo "test18_rundir..." > test18.diff
	$(MAKE) test18_rundir
	@echo "test18_run..." > test18.diff
	$(MAKE) test18_run
	@echo "test18_check..." > test18.diff
	$(MAKE) test18_check

test18_compile:
	./Config.pl -v=Empty,PC/FLEKS,GM/BATSRUS -amrex
	./Config.pl -o=GM:u=GemReconnect,e=MhdHyp,ng=2,g=4,4,4
	$(MAKE) SWMF
	$(MAKE) PIDL
	
test18_rundir:
	rm -rf ${TESTDIR}
	$(MAKE) rundir RUNDIR=${DIR}/${TESTDIR} COMPDIR=GM,PC
	cp Param/PARAM.in.test.GMPC.FLEKS.3D  ${TESTDIR}/PARAM.in
	-Scripts/TestParam.pl -F ${TESTDIR}/PARAM.in

test18_run:
	cd ${TESTDIR}; ${MPIRUN} ./SWMF.exe > runlog
	cd ${TESTDIR}; ./PostProc.pl RESULTS

test18_check:
	-@(${DIFFNUM} -b -r=1e-6 -a=1e-15 \
	     ${TESTDIR}/RESULTS/PC/3d_var_region0_0_t00000016_n00000016.out \
	     output/test18/3d.ref.gz \
	     > test18.diff)
	ls -l test18.diff

### test19: SC+IH+OH+SP ###

test19:
	rm -f test19*.diff
	@echo "test19_compile..." >> test19_sp.diff
	${MAKE} test19_compile
	@echo "test19_rundir..." >> test19_sp.diff
	${MAKE} test19_rundir
	@echo "test19_run..." >> test19_sp.diff
	${MAKE} test19_run
	@echo "test19_restart..." >> test19_sp.diff
	${MAKE} test19_restart
	@echo "test19_check..." >> test19_sp.diff
	${MAKE} test19_check

test19_compile:
	./Config.pl -v=Empty,SC/BATSRUS,IH/BATSRUS,OH/BATSRUS,SP/MFLAMPA
	./Config.pl -o=SC:u=Awsom,e=Awsom,ng=2,g=4,4,4
	./Config.pl -o=IH:u=Awsom,e=Awsom,ng=2,g=4,4,4
	./Config.pl -o=OH:u=Awsom,e=Awsom,ng=2,g=4,4,4
	./Config.pl -o=SP:g=2000
	${MAKE} SWMF
	make PIDL

test19_rundir:
	rm -rf ${TESTDIR}
	${MAKE} rundir RUNDIR=${DIR}/${TESTDIR}
	cp Param/PARAM.in.test.SCIHOHSP ${TESTDIR}/PARAM.in
	-Scripts/TestParam.pl -F ${TESTDIR}/PARAM.in

test19_run:
	cd ${TESTDIR}; rm -rf RESTART
	cd ${TESTDIR}; ${MPIRUN} ./SWMF.exe > runlog
	cd ${TESTDIR}; ./Restart.pl -o -W RESTART

test19_restart:
	cd ${TESTDIR}; rm -rf RESULTS runlog_restart
	cd ${TESTDIR}; ./Restart.pl -i RESTART
	cd ${TESTDIR}; mv PARAM.in PARAM.in_start; rm -f PARAM.in_orig_
	cd ${TESTDIR}; cp Param/PARAM.in.test.restart.SCIHOHSP PARAM.in
	-Scripts/TestParam.pl -F ${TESTDIR}/PARAM.in
	cd ${TESTDIR}; ${MPIRUN} ./SWMF.exe > runlog_restart
	cd ${TESTDIR}; ./PostProc.pl -M -cat RESULTS

test19_check:
	-(${DIFFNUM} -b -r=1e-5 -a=1e-26 \
		${TESTDIR}/RESULTS/SC/log_n000000.log \
		output/test15_long/SC_log_n000000.log \
		> test19_sc.diff)
	-(${DIFFNUM} -b -r=1e-5 -a=1e-26 \
		${TESTDIR}/RESULTS/SC/rfr_rwi_4_*.outs \
		output/test8/SC_rfr.ref.gz \
		> test19_sc_rfr.diff)
	-(${DIFFNUM} -b -r=1e-5 \
		${TESTDIR}/RESULTS/IH/log_n000000.log \
		output/test19/IH_log_n000000.log \
		> test19_ih.diff)
	-(${DIFFNUM} -b -r=1e-5 \
		${TESTDIR}/RESULTS/SC/los_sdo_aia*.out \
		SC/BATSRUS/data/output/test8/SC_los_sdo_aia.out.gz \
		> test19_sc_los.diff)
	cat ${TESTDIR}/RESULTS/SP/MH_data_*_*_t00000150_n000010.out \
	    > ${TESTDIR}/RESULTS/SP/MH_data.outs
	-(${DIFFNUM} -a=1e-6 -r=1e-6 -b \
		${TESTDIR}/RESULTS/SP/MH_data.outs \
		output/test19/MH_data.ref.gz \
		> test19_sp.diff)
	ls -l test19*.diff


### test_ramscb: GM+IE+IM/RAMSCB ###

test_ramscb:
	@echo "test_ramscb_compile..." > test_ramscb_gm.diff
	@echo "test_ramscb_compile..." > test_ramscb_ie.diff
	@echo "test_ramscb_compile..." > test_ramscb_im.diff
	${MAKE} test_ramscb_compile
	@echo "test_ramscb_rundir..." >> test_ramscb_gm.diff
	@echo "test_ramscb_rundir..." >> test_ramscb_ie.diff
	@echo "test_ramscb_rundir..." >> test_ramscb_im.diff
	${MAKE} test_ramscb_rundir
	@echo "test_ramscb_run..." >> test_ramscb_gm.diff
	@echo "test_ramscb_run..." >> test_ramscb_ie.diff
	@echo "test_ramscb_run..." >> test_ramscb_im.diff
	${MAKE} test_ramscb_run
	@echo "test_ramscb_check..." >> test_ramscb_gm.diff
	@echo "test_ramscb_check..." >> test_ramscb_ie.diff
	@echo "test_ramscb_check..." >> test_ramscb_im.diff
	${MAKE} test_ramscb_check

test_ramscb_compile:
	./Config.pl -v=Empty,GM/BATSRUS,IE/Ridley_serial,IM/RAM_SCB
	./Config.pl -o=GM:u=Default,e=Mhd,ng=2,g=8,8,8,IE:g=91,181
	${MAKE} SWMF
	make PIDL

test_ramscb_rundir:
	rm -rf ${TESTDIR}
	${MAKE} rundir RUNDIR=${DIR}/${TESTDIR}
	cp ${IMDIR}/input/PARAM.in.test.gmieim.simple ${TESTDIR}/PARAM.in
	-Scripts/TestParam.pl -F ${TESTDIR}/PARAM.in

test_ramscb_run:
	cd ${TESTDIR}; ${MPIRUN} ./SWMF.exe > runlog
	cd ${TESTDIR}; ./PostProc.pl -M RESULTS

test_ramscb_check:
	-@(${DIFFNUM} -b -r=1e-3 -a=1e-8 \
		${TESTDIR}/RESULTS/GM/log_n000001.log \
		${IMDIR}/output/test9/GM_log_n000001.log \
		> test_ramscb_gm.diff)
	-@(${DIFFNUM} -b -r=1e-3 -a=1e-8 \
		${TESTDIR}/RESULTS/IM/ram/log_n000001.log \
		${IMDIR}/output/test9/IM_log_n000001.log \
		> test_ramscb_im.diff)
	-@(${DIFFNUM} -b -r=1e-6 \
		${TESTDIR}/RESULTS/IE/IE*.log \
		${IMDIR}/output/test9/IE.log \
		> test_ramscb_ie.diff)
	ls -l test_ramscb*.diff

######## test PC/AMPS coupled with GM/BATSRUS #########

test12:
	-@(if [ -d PC/AMPS ]; then $(MAKE) test12_aniso_AMPS TESTDIR=run_test12_aniso_AMPS; fi)
	-@(if [ -d PC/AMPS ]; then $(MAKE) test12_aniso_AMPS_divE_corr TESTDIR=run_test12_aniso_AMPS_divE_corr; fi)
	-@(if [ -d PC/AMPS ]; then $(MAKE) test12_aniso_AMPS_dynamic   TESTDIR=run_test12_aniso_AMPS_dynamic; fi)
	-@(if [ -d PC/AMPS ]; then $(MAKE) test12_aniso_AMPS_openmp    TESTDIR=run_test12_aniso_AMPS_openmp; fi)
	-@(if [ -d PC/AMPS ]; then ./Config.pl -noopenmp -o=PC:openmp=off; fi)

### test12_aniso_AMPS_dynamic_fluxrope: test PC coupled with GM/BATSRUS ###
test12_aniso_AMPS_dynamic_fluxrope:
	@echo "test12_aniso_AMPS_dynamic_fluxrope_compile..." > test12_aniso_AMPS_dynamic_fluxrope_gmpc.diff
	$(MAKE) test12_aniso_AMPS_dynamic_fluxrope_compile
	@echo "test12_aniso_AMPS_dynamic_fluxrope_rundir..." > test12_aniso_AMPS_dynamic_fluxrope_gmpc.diff
	$(MAKE) test12_aniso_AMPS_dynamic_fluxrope_rundir
	@echo "test12_aniso_AMPS_dynamic_fluxrope_run..." > test12_aniso_AMPS_dynamic_fluxrope_gmpc.diff
	$(MAKE) test12_aniso_AMPS_dynamic_fluxrope_run
	@echo "test12_aniso_AMPS_dynamic_fluxrope_check..." > test12_aniso_AMPS_dynamic_fluxrope_gmpc.diff
	$(MAKE) test12_aniso_AMPS_dynamic_fluxrope_check

test12_aniso_AMPS_dynamic_fluxrope_compile:
	./Config.pl -v=Empty,GM/BATSRUS,PC/AMPS
	./Config.pl -o=GM:u=FluxRope,e=MhdHypPe,ng=2,g=8,8,4
	./Config.pl -o=PC:spice-path=nospice,spice-kernels=nospice,application=fluxrope-dynamic,amps-test=on,avx=256
	$(MAKE) SWMF
	$(MAKE) PIDL

test12_aniso_AMPS_dynamic_fluxrope_rundir:
	rm -rf ${TESTDIR}
	$(MAKE) rundir RUNDIR=${DIR}/${TESTDIR} COMPDIR=GM,PC
	cp Param/PARAM.in.test.GMPC.aniso.AMPS.fluxrope  ${TESTDIR}/PARAM.in
	-Scripts/TestParam.pl -F ${TESTDIR}/PARAM.in

test12_aniso_AMPS_dynamic_fluxrope_run:
	cd ${TESTDIR}; mpirun -np 8 ./SWMF.exe > runlog
	cd ${TESTDIR}; ./PostProc.pl -m RESULTS/fluxrope

test12_aniso_AMPS_dynamic_fluxrope_check:
	-@(${DIFFNUM} -b -a=2.0e-3 \
	   ${TESTDIR}/RESULTS/fluxrope/PC/z\=0_var_region0_0_t00000000_n00000010.out  \
	   output/test12/3d_aniso_fluxrope.out.gz \
		> test12_aniso_AMPS_dynamic_fluxrope_gmpc.diff)
	ls -l test12_aniso_AMPS_dynamic_fluxrope_gmpc.diff

#############################################################################

test12_aniso_AMPS:
	@echo "test12_aniso_AMPS_compile..." > test12_aniso_AMPS_gmpc.diff
	$(MAKE) test12_aniso_AMPS_compile
	@echo "test12_aniso_AMPS_rundir..." > test12_aniso_AMPS_gmpc.diff
	$(MAKE) test12_aniso_AMPS_rundir
	@echo "test12_aniso_AMPS_run..." > test12_aniso_AMPS_gmpc.diff
	$(MAKE) test12_aniso_AMPS_run
	@echo "test12_aniso_AMPS_check..." > test12_aniso_AMPS_gmpc.diff
	$(MAKE) test12_aniso_AMPS_check

test12_aniso_AMPS_compile:
	./Config.pl -v=Empty,GM/BATSRUS,PC/AMPS
	./Config.pl -o=GM:u=Waves,e=MhdAnisoP,ng=2,g=8,8,8
	./Config.pl -o=PC:spice-path=nospice,spice-kernels=nospice,application=fluid-pic-coupler,amps-test=on 
	$(MAKE) SWMF
	$(MAKE) PIDL
	
test12_aniso_AMPS_rundir:
	rm -rf ${TESTDIR}
	$(MAKE) rundir RUNDIR=${DIR}/${TESTDIR} COMPDIR=GM,PC
	cp Param/PARAM.in.test.GMPC.aniso.AMPS  ${TESTDIR}/PARAM.in
	-Scripts/TestParam.pl -F ${TESTDIR}/PARAM.in

test12_aniso_AMPS_run:
	cd ${TESTDIR}; ${MPIRUN} ./SWMF.exe > runlog
	cd ${TESTDIR}; ./PostProc.pl -m RESULTS/5step
	cp Param/PARAM.in.test.GMPC.aniso.AMPS.2step  ${TESTDIR}/PARAM.in
	-Scripts/TestParam.pl -F ${TESTDIR}/PARAM.in
	cd ${TESTDIR}; ${MPIRUN} ./SWMF.exe > runlog
	cd ${TESTDIR}; ./PostProc.pl -m RESULTS/2step
	cd ${TESTDIR}; ./Restart.pl -v -i RESULTS/2step/RESTART
	cp Param/PARAM.in.test.GMPC.aniso.AMPS.2step.restart ${TESTDIR}/PARAM.in
	-Scripts/TestParam.pl -F ${TESTDIR}/PARAM.in
	cd ${TESTDIR}; ${MPIRUN} ./SWMF.exe > runlog
	cd ${TESTDIR}; ./PostProc.pl -m RESULTS/2restart	

test12_aniso_AMPS_check:
	-@(${DIFFNUM} -b -a=2.0e-3 \
	   ${TESTDIR}/RESULTS/5step/PC/3d_fluid_region0_0_t00000005_n00000005.out \
	   output/test12/3d_aniso.out.gz \
		> test12_aniso_AMPS_gmpc.diff)
	-@(${DIFFNUM} -b -a=2.0e-3 \
	   ${TESTDIR}/RESULTS/2restart/PC/3d_fluid_region0_0_t00000005_n00000005.out \
	   output/test12/3d_aniso.out.gz \
		>> test12_aniso_AMPS_gmpc.diff)
	ls -l test12_aniso_AMPS_gmpc.diff

#########################################################

test12_aniso_AMPS_divE_corr:
	@echo "test12_aniso_AMPS_divE_corr_compile..." > test12_aniso_AMPS_divE_corr_gmpc.diff
	$(MAKE) test12_aniso_AMPS_divE_corr_compile
	@echo "test12_aniso_AMPS_divE_corr_rundir..." > test12_aniso_AMPS_divE_corr_gmpc.diff
	$(MAKE) test12_aniso_AMPS_divE_corr_rundir
	@echo "test12_aniso_AMPS_divE_corr_run..." > test12_aniso_AMPS_divE_corr_gmpc.diff
	$(MAKE) test12_aniso_AMPS_divE_corr_run
	@echo "test12_aniso_AMPS_divE_corr_check..." > test12_aniso_AMPS_divE_corr_gmpc.diff
	$(MAKE) test12_aniso_AMPS_divE_corr_check

test12_aniso_AMPS_divE_corr_compile:
	./Config.pl -v=Empty,GM/BATSRUS,PC/AMPS
	./Config.pl -o=GM:u=Waves,e=MhdAnisoP,ng=2,g=8,8,8
	./Config.pl -o=PC:spice-path=nospice,spice-kernels=nospice,application=fluid-pic-coupler-with-correction,amps-test=on 
	$(MAKE) SWMF
	$(MAKE) PIDL
	
test12_aniso_AMPS_divE_corr_rundir:
	rm -rf ${TESTDIR}
	$(MAKE) rundir RUNDIR=${DIR}/${TESTDIR} COMPDIR=GM,PC
	cp Param/PARAM.in.test.GMPC.aniso.AMPS.corr  ${TESTDIR}/PARAM.in
	-Scripts/TestParam.pl -F ${TESTDIR}/PARAM.in

test12_aniso_AMPS_divE_corr_run:
	cd ${TESTDIR}; ${MPIRUN} ./SWMF.exe > runlog
	cd ${TESTDIR}; ./PostProc.pl -m RESULTS/5step_corr

test12_aniso_AMPS_divE_corr_check:
	-@(${DIFFNUM} -b -a=2.0e-3 \
	   ${TESTDIR}/RESULTS/5step_corr/PC/3d_fluid_region0_0_t00000005_n00000005.out \
	   output/test12/3d_aniso_divE.out.gz \
		> test12_aniso_AMPS_divE_corr_gmpc.diff)
	ls -l test12_aniso_AMPS_divE_corr_gmpc.diff

#########################################################

test12_aniso_AMPS_openmp:
	@echo "test12_aniso_AMPS_openmp_compile..." > test12_aniso_AMPS_openmp_gmpc.diff
	$(MAKE) test12_aniso_AMPS_openmp_compile
	@echo "test12_aniso_AMPS_openmp_rundir..." > test12_aniso_AMPS_openmp_gmpc.diff
	$(MAKE) test12_aniso_AMPS_openmp_rundir
	@echo "test12_aniso_AMPS_openmp_run..." > test12_aniso_AMPS_openmp_gmpc.diff
	$(MAKE) test12_aniso_AMPS_openmp_run
	@echo "test12_aniso_AMPS_openmp_check..." > test12_aniso_AMPS_openmp_gmpc.diff
	$(MAKE) test12_aniso_AMPS_openmp_check

test12_aniso_AMPS_openmp_compile:
	./Config.pl -v=Empty,GM/BATSRUS,PC/AMPS
	./Config.pl -o=GM:u=Waves,e=MhdAnisoP,ng=2,g=8,8,8
	./Config.pl -openmp -o=PC:spice-path=nospice,spice-kernels=nospice,application=fluid-pic-coupler,amps-test=on,openmp=on 
	$(MAKE) SWMF
	$(MAKE) PIDL
	
test12_aniso_AMPS_openmp_rundir:
	rm -rf ${TESTDIR}
	$(MAKE) rundir RUNDIR=${DIR}/${TESTDIR} COMPDIR=GM,PC
	cp Param/PARAM.in.test.GMPC.aniso.AMPS  ${TESTDIR}/PARAM.in
	-Scripts/TestParam.pl -F ${TESTDIR}/PARAM.in

test12_aniso_AMPS_openmp_run:
	cd ${TESTDIR}; ${OMPIRUN} ./SWMF.exe > runlog
	cd ${TESTDIR}; ./PostProc.pl -m RESULTS

test12_aniso_AMPS_openmp_check:
	-@(${DIFFNUM} -b -a=2.0e-3 \
	   ${TESTDIR}/RESULTS/PC/3d_fluid_region0_0_t00000005_n00000005.out \
	   output/test12/3d_aniso.out.gz \
		> test12_aniso_AMPS_openmp_gmpc.diff)
	ls -l test12_aniso_AMPS_openmp_gmpc.diff

#########################################################
test12_aniso_AMPS_dynamic:
	@echo "test12_aniso_AMPS_dynamic_compile..." > test12_aniso_AMPS_dynamic_gmpc.diff
	$(MAKE) test12_aniso_AMPS_dynamic_compile
	@echo "test12_aniso_AMPS_dynamic_rundir..." > test12_aniso_AMPS_dynamic_gmpc.diff
	$(MAKE) test12_aniso_AMPS_dynamic_rundir
	@echo "test12_aniso_AMPS_dynamic_run..." > test12_aniso_AMPS_dynamic_gmpc.diff
	$(MAKE) test12_aniso_AMPS_dynamic_run
	@echo "test12_aniso_AMPS_dynamic_check..." > test12_aniso_AMPS_dynamic_gmpc.diff
	$(MAKE) test12_aniso_AMPS_dynamic_check

test12_aniso_AMPS_dynamic_compile:
	./Config.pl -v=Empty,GM/BATSRUS,PC/AMPS
	./Config.pl -o=GM:u=Waves,e=MhdAnisoP,ng=2,g=8,8,8
	./Config.pl -o=PC:spice-path=nospice,spice-kernels=nospice,application=fluid-pic-coupler-dynamic,amps-test=on 
	$(MAKE) SWMF
	$(MAKE) PIDL
	
test12_aniso_AMPS_dynamic_rundir:
	rm -rf ${TESTDIR}
	$(MAKE) rundir RUNDIR=${DIR}/${TESTDIR} COMPDIR=GM,PC
	cp Param/PARAM.in.test.GMPC.aniso.AMPS.dynamic  ${TESTDIR}/PARAM.in
	-Scripts/TestParam.pl -F ${TESTDIR}/PARAM.in

test12_aniso_AMPS_dynamic_run:
	cd ${TESTDIR}; ${MPIRUN} ./SWMF.exe > runlog
	cd ${TESTDIR}; ./PostProc.pl -m RESULTS

test12_aniso_AMPS_dynamic_check:
	-@(${DIFFNUM} -b -a=2.0e-3 \
	   ${TESTDIR}/RESULTS/PC/3d_fluid_region0_0_t00000010_n00000010.out \
	   output/test12/3d_AMPS_aniso_dynamic_gmpc.out.gz \
		> test12_aniso_AMPS_dynamic_gmpc.diff)
	ls -l test12_aniso_AMPS_dynamic_gmpc.diff

### test12_ipic3d: test PC/IPIC3D2 coupled with GM/BATSRUS ###

test12_ipic3d:
	-@(if [ -d PC/IPIC3D2 ]; then $(MAKE) test12_aniso TESTDIR=run_test12_aniso; fi)
	-@(if [ -d PC/IPIC3D2 ]; then $(MAKE) test12_iso   TESTDIR=run_test12_iso;   fi)
	-@(if [ -d PC/IPIC3D2 ]; then $(MAKE) test12_3d    TESTDIR=run_test12_3d;    fi)

#########################################################
test12_aniso:
	@echo "test12_aniso_compile..." > test12_aniso_gmpc.diff
	$(MAKE) test12_aniso_compile
	@echo "test12_aniso_rundir..." > test12_aniso_gmpc.diff
	$(MAKE) test12_aniso_rundir
	@echo "test12_aniso_run..." > test12_aniso_gmpc.diff
	$(MAKE) test12_aniso_run
	@echo "test12_aniso_check..." > test12_aniso_gmpc.diff
	$(MAKE) test12_aniso_check

test12_aniso_compile:
	./Config.pl -v=Empty,PC/IPIC3D2,GM/BATSRUS
	./Config.pl -o=GM:u=Waves,e=MhdAnisoP,ng=2,g=8,8,1
	$(MAKE) SWMF
	$(MAKE) PIDL
	
test12_aniso_rundir:
	rm -rf ${TESTDIR}
	$(MAKE) rundir RUNDIR=${DIR}/${TESTDIR} COMPDIR=GM,PC
	cp Param/PARAM.in.test.GMPC.start.aniso  ${TESTDIR}/PARAM.in
	-Scripts/TestParam.pl -F ${TESTDIR}/PARAM.in

test12_aniso_run:
	cd ${TESTDIR}; ${MPIRUN} ./SWMF.exe > runlog
	cd ${TESTDIR}; ./PostProc.pl -m RESULTS/4step
	cp Param/PARAM.in.test.GMPC.start.2step.aniso ${TESTDIR}/PARAM.in
	-Scripts/TestParam.pl -F ${TESTDIR}/PARAM.in
	cd ${TESTDIR}; ${MPIRUN} ./SWMF.exe > runlog
	cd ${TESTDIR}; ./PostProc.pl -m RESULTS/2step
	cd ${TESTDIR}; ./Restart.pl -v -i RESULTS/2step/RESTART
	cp Param/PARAM.in.test.GMPC.restart.2step.aniso ${TESTDIR}/PARAM.in
	-Scripts/TestParam.pl -F ${TESTDIR}/PARAM.in
	cd ${TESTDIR}; ${MPIRUN} ./SWMF.exe > runlog
	cd ${TESTDIR}; ./PostProc.pl -m RESULTS/2restart	

test12_aniso_check:
	-@(${DIFFNUM} -b -a=2.0e-3 \
	   ${TESTDIR}/RESULTS/4step/GM/z=0_var_1_t00000005_n00000004.out \
	   output/test12/z=0_var.out.ref.aniso.gz \
		> test12_aniso_gmpc.diff)
	-@(${DIFFNUM} -b -a=2.0e-3 \
	   ${TESTDIR}/RESULTS/2restart/GM/z=0_var_1_t00000005_n00000004.out \
	   output/test12/z=0_var.out.ref.aniso.gz \
		>> test12_aniso_gmpc.diff)
	ls -l test12_aniso_gmpc.diff

#### test12_3d: PC/IPIC3D2 + GM/BATSRUS for 3D reconetion ###

test12_3d:
	@echo "test12_3d_compile..." > test12_3d_gm.diff
	$(MAKE) test12_3d_compile
	@echo "test12_3d_rundir..." > test12_3d_gm.diff
	$(MAKE) test12_3d_rundir
	@echo "test12_3d_run..." > test12_3d_gm.diff
	$(MAKE) test12_3d_run
	@echo "test12_3d_check..." > test12_3d_gm.diff
	$(MAKE) test12_3d_check

test12_3d_compile:
	./Config.pl -v=Empty,PC/IPIC3D2,GM/BATSRUS
	./Config.pl -o=GM:u=GemReconnect,e=MhdHyp,ng=2,g=4,4,4
	$(MAKE) SWMF
	$(MAKE) PIDL
	
test12_3d_rundir:
	rm -rf ${TESTDIR}
	$(MAKE) rundir RUNDIR=${DIR}/${TESTDIR} COMPDIR=GM,PC
	cp Param/PARAM.in.test.GMPC.3D.start  ${TESTDIR}/PARAM.in
	-Scripts/TestParam.pl -F ${TESTDIR}/PARAM.in

test12_3d_run:
	cd ${TESTDIR}; ${MPIRUN} ./SWMF.exe > runlog
	cd ${TESTDIR}; ./PostProc.pl RESULTS/4step

test12_3d_check:
	-@(${DIFFNUM} -b -a=1e-9 \
	   ${TESTDIR}/RESULTS/4step/GM/x=0_var_1_t00000000_n00000004.out \
	   output/test12/x=0_var.out.ref.3d.gz \
		> test12_3d_gm.diff)
	-@(${DIFFNUM} -b -a=1e-9 \
	   ${TESTDIR}/RESULTS/4step/GM/y=0_var_2_t00000000_n00000004.out \
	   output/test12/y=0_var.out.ref.3d.gz \
		>> test12_3d_gm.diff)
	-@(${DIFFNUM} -b -a=1e-9 \
	   ${TESTDIR}/RESULTS/4step/GM/z=0_var_3_t00000000_n00000004.out \
	   output/test12/z=0_var.out.ref.3d.gz \
		>> test12_3d_gm.diff)
	-@(${DIFFNUM} -b -a=1e-7 -r=1e-7 \
	   ${TESTDIR}/RESULTS/4step/PC/cut_all_region0_0_t00000000_n00000004.out \
	   output/test12/cut_all.out.ref.3d.gz \
		> test12_3d_pc.diff)
	ls -l test12_3d_*.diff

#### test12_iso: PC/IPIC3D2 + GM/BATSRUS with isotropic MHD ###

test12_iso:
	@echo "test12_iso_compile..." > test12_iso_gm.diff
	$(MAKE) test12_iso_compile
	@echo "test12_iso_rundir..." > test12_iso_gm.diff
	$(MAKE) test12_iso_rundir
	@echo "test12_iso_run..." > test12_iso_gm.diff
	$(MAKE) test12_iso_run
	@echo "test12_iso_check..." > test12_iso_gm.diff
	$(MAKE) test12_iso_check

test12_iso_compile:
	./Config.pl -v=Empty,PC/IPIC3D2,GM/BATSRUS
	./Config.pl -o=GM:u=Waves,e=MultiIonPe,ng=2,g=8,8,1
	$(MAKE) SWMF
	$(MAKE) PIDL
	
test12_iso_rundir:
	rm -rf ${TESTDIR}
	$(MAKE) rundir RUNDIR=${DIR}/${TESTDIR} COMPDIR=GM,PC
	cp Param/PARAM.in.test.GMPC.start  ${TESTDIR}/PARAM.in
	-Scripts/TestParam.pl -F ${TESTDIR}/PARAM.in

test12_iso_run:
	cd ${TESTDIR}; ${MPIRUN} ./SWMF.exe > runlog
	cd ${TESTDIR}; ./PostProc.pl -m RESULTS

test12_iso_check:
	-@(${DIFFNUM} -t -a=1e-9 \
	   ${TESTDIR}/RESULTS/GM/z=0_all_1_t00000002_n00000004.out \
	   output/test12/z=0_all.out.ref.iso.gz \
		> test12_iso_gm.diff)
	-@(${DIFFNUM} -b -a=1e-8 \
	   ${TESTDIR}/RESULTS/PC/z=0_fluid_region0_0_t00000002_n00000004.out \
	   output/test12/z=0_fluid.out.ref.iso.gz \
		> test12_iso_pc.diff)
	-@(${DIFFNUM} -b -a=1e-9 \
	   ${TESTDIR}/RESULTS/PC/cut_particles1_region0_1_t00000002_n00000004.out \
	   output/test12/cut_particles1.out.ref.iso.gz \
		>> test12_iso_pc.diff)
	ls -l test12_iso*.diff

### test_chgl: test SC/BATSRUS + IH/BATSRUS with aligned field and stream lines ###

test_chgl:
	@echo "test_chgl_compile..." > test_chgl_sc.diff
	@echo "test_chgl_compile..." > test_chgl_ih.diff
	${MAKE} test_chgl_compile
	@echo "test_chgl_rundir..." >> test_chgl_sc.diff
	@echo "test_chgl_rundir..." >> test_chgl_ih.diff
	${MAKE} test_chgl_rundir
	@echo "test_chgl_run..." >> test_chgl_sc.diff
	@echo "test_chgl_run..." >> test_chgl_ih.diff
	${MAKE} test_chgl_run
	@echo "test_chgl_check..." >> test_chgl_sc.diff
	@echo "test_chgl_check..." >> test_chgl_ih.diff
	${MAKE} test_chgl_check

test_chgl_compile:
	./Config.pl -v=Empty,SC/BATSRUS,IH/BATSRUS
	./Config.pl -o=SC:u=Awsom,e=AwsomChGL,ng=2,g=6,4,4
	./Config.pl -o=IH:u=Awsom,e=AwsomChGL,ng=2,g=4,4,4
	${MAKE} SWMF
	make PIDL

test_chgl_rundir:
	rm -rf ${TESTDIR}
	$(MAKE) rundir RUNDIR=${DIR}/${TESTDIR} COMPDIR=SC,IH
	cp ${TESTDIR}/SC/Param/PARAM.in.ChGL_short  ${TESTDIR}/PARAM.in
	cp SC/BATSRUS/data/input/Gong_ChGL.dat ${TESTDIR}/SC/.
	-Scripts/TestParam.pl -F ${TESTDIR}/PARAM.in

test_chgl_run:
	cd ${TESTDIR}; ${MPIRUN} ./SWMF.exe > runlog
	cd ${TESTDIR}; ./PostProc.pl -M -cat RESULTS

test_chgl_check:
	-(${DIFFNUM} -b -r=1e-5 -a=1e-26 \
		${TESTDIR}/RESULTS/SC/log_n000000.log \
		SC/BATSRUS/data/output/test_chgl/SC_log_n000000.log \
		> test_chgl_sc.diff)
	-(${DIFFNUM} -b -r=1e-5 \
		${TESTDIR}/RESULTS/IH/log_n000000.log \
		SC/BATSRUS/data/output/test_chgl/IH_log_n000000.log \
		> test_chgl_ih.diff)
	ls -l test_chgl_*.diff