diff options
| author | schnetter <> | 2002-01-08 13:20:00 +0000 |
|---|---|---|
| committer | schnetter <> | 2002-01-08 13:20:00 +0000 |
| commit | 8b92cb0c94a56fd6e888569b43e151c5cb7008b7 (patch) | |
| tree | 8b4d9d65cf9ddc5d509964a019dbaa3da3386c04 /CarpetExtra/IDScalarWave | |
| parent | 6a7486e6ffe1f46c93f9b7accb4cb409928be663 (diff) | |
Added new 1/r initial data for WaveToy.
Added new 1/r initial data for WaveToy.
Added excision to WaveToy for 1/r data.
darcs-hash:20020108132057-07bb3-541b25e99f5d694d4baa82ce3d09ced4383b27c9.gz
Diffstat (limited to 'CarpetExtra/IDScalarWave')
| -rw-r--r-- | CarpetExtra/IDScalarWave/param.ccl | 3 | ||||
| -rw-r--r-- | CarpetExtra/IDScalarWave/src/InitialData.F77 | 16 |
2 files changed, 17 insertions, 2 deletions
diff --git a/CarpetExtra/IDScalarWave/param.ccl b/CarpetExtra/IDScalarWave/param.ccl index 9b9d7fc1a..75244b562 100644 --- a/CarpetExtra/IDScalarWave/param.ccl +++ b/CarpetExtra/IDScalarWave/param.ccl @@ -1,5 +1,5 @@ # Parameter definitions for thorn IDScalarWave -# $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/IDScalarWave/param.ccl,v 1.2 2001/03/05 14:31:46 eschnett Exp $ +# $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/IDScalarWave/param.ccl,v 1.3 2002/01/08 14:20:57 schnetter Exp $ shares: grid @@ -12,6 +12,7 @@ KEYWORD initial_data "Type of initial data" "plane" :: "Plane wave" "gaussian" :: "Gaussian wave" "box" :: "Box wave" + "1/r" :: "1/r data" "none" :: "No initial data, zero phi" } "gaussian" diff --git a/CarpetExtra/IDScalarWave/src/InitialData.F77 b/CarpetExtra/IDScalarWave/src/InitialData.F77 index 739c6d2d2..5f330c924 100644 --- a/CarpetExtra/IDScalarWave/src/InitialData.F77 +++ b/CarpetExtra/IDScalarWave/src/InitialData.F77 @@ -1,5 +1,5 @@ c -*-Fortran-*- -c $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/IDScalarWave/src/InitialData.F77,v 1.8 2001/12/14 18:00:01 schnetter Exp $ +c $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/IDScalarWave/src/InitialData.F77,v 1.9 2002/01/08 14:20:57 schnetter Exp $ /*@@ @file InitialData.F77 @@ -132,6 +132,20 @@ c Use kx,ky,kz as number of modes in each direction. end do end do + else if (CCTK_EQUALS(initial_data, "1/r")) then + + do k=1,cctk_lsh(3) + do j=1,cctk_lsh(2) + do i=1,cctk_lsh(1) + + phi(i,j,k) = 1 / spher3d_r(i,j,k) + phi_p(i,j,k) = 1 / spher3d_r(i,j,k) + phi_p_p(i,j,k) = 1 / spher3d_r(i,j,k) + + end do + end do + end do + else do k=1,cctk_lsh(3) |