diff options
| author | Erik Schnetter <schnetter@cct.lsu.edu> | 2006-02-25 20:49:00 +0000 |
|---|---|---|
| committer | Erik Schnetter <schnetter@cct.lsu.edu> | 2006-02-25 20:49:00 +0000 |
| commit | 63c13016173d70762f6f23c624f597484c906096 (patch) | |
| tree | 556b48a4df39606985a88f09e1711d5c3d27b2a9 /CarpetExtra | |
| parent | 6e478cfc32f7b117779366ecbd2430fcf6020bbb (diff) | |
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
Diffstat (limited to 'CarpetExtra')
| -rw-r--r-- | CarpetExtra/IDScalarWaveMoL/param.ccl | 1 | ||||
| -rw-r--r-- | CarpetExtra/IDScalarWaveMoL/src/initialdata.F77 | 12 |
2 files changed, 13 insertions, 0 deletions
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 |