From 63c13016173d70762f6f23c624f597484c906096 Mon Sep 17 00:00:00 2001 From: Erik Schnetter Date: Sat, 25 Feb 2006 20:49:00 +0000 Subject: IDScalarWaveMoL: Add initial data "level index" Add a new kind of initial data "level index", which fills each refinement levels with its refinement level index. This is good for testing. darcs-hash:20060225204907-dae7b-8b2c03bbb0b9293690ae8d18785b89f0e490acff.gz --- CarpetExtra/IDScalarWaveMoL/param.ccl | 1 + CarpetExtra/IDScalarWaveMoL/src/initialdata.F77 | 12 ++++++++++++ 2 files changed, 13 insertions(+) (limited to 'CarpetExtra') diff --git a/CarpetExtra/IDScalarWaveMoL/param.ccl b/CarpetExtra/IDScalarWaveMoL/param.ccl index c06c04c20..7dea15888 100644 --- a/CarpetExtra/IDScalarWaveMoL/param.ccl +++ b/CarpetExtra/IDScalarWaveMoL/param.ccl @@ -5,6 +5,7 @@ KEYWORD initial_data "Initial data" "plane wave" :: "Plane wave" "Gaussian pulse" :: "Gaussian pulse" "Gaussian" :: "Gaussian" + "level index" :: "Fill the refinement levels with their level index (for testing)" } "plane wave" diff --git a/CarpetExtra/IDScalarWaveMoL/src/initialdata.F77 b/CarpetExtra/IDScalarWaveMoL/src/initialdata.F77 index dc5ee0369..c86d0de1a 100644 --- a/CarpetExtra/IDScalarWaveMoL/src/initialdata.F77 +++ b/CarpetExtra/IDScalarWaveMoL/src/initialdata.F77 @@ -8,6 +8,9 @@ DECLARE_CCTK_ARGUMENTS DECLARE_CCTK_FUNCTIONS DECLARE_CCTK_PARAMETERS + CCTK_REAL one, two + parameter (one = 1) + parameter (two = 2) CCTK_REAL pi parameter (pi = 3.141592653589793238462643383279502884197169399375105820974944592307816406286208998628034825342117068d0) CCTK_REAL omega @@ -70,5 +73,14 @@ end do end do end do + else if (CCTK_EQUALS(initial_data, "level index")) then + do k=1,cctk_lsh(3) + do j=1,cctk_lsh(2) + do i=1,cctk_lsh(1) + phi(i,j,k) = log (one * cctk_levfac(1)) / log (two) + psi(i,j,k) = log (one * cctk_levfac(1)) / log (two) + end do + end do + end do end if end -- cgit v1.2.3