diff options
author | Erik Schnetter <schnetter@aei.mpg.de> | 2005-05-01 17:15:00 +0000 |
---|---|---|
committer | Erik Schnetter <schnetter@aei.mpg.de> | 2005-05-01 17:15:00 +0000 |
commit | 0fbb3fc20f36bce10eb1f92921a3f947abef3c99 (patch) | |
tree | 0af3b28ff2cce5be1929d1af95633a26a7d56a77 /CarpetExtra/IDScalarWaveMoL/src | |
parent | a41c0f620092614b749e3e4f89696291f592881c (diff) |
IDScalarWaveMoL: Add initial data type "Gaussian"
darcs-hash:20050501171548-891bb-71cad901294c978432e36baf96614c8a506802bd.gz
Diffstat (limited to 'CarpetExtra/IDScalarWaveMoL/src')
-rw-r--r-- | CarpetExtra/IDScalarWaveMoL/src/initialdata.F77 | 27 |
1 files changed, 27 insertions, 0 deletions
diff --git a/CarpetExtra/IDScalarWaveMoL/src/initialdata.F77 b/CarpetExtra/IDScalarWaveMoL/src/initialdata.F77 index 17eed204b..dc5ee0369 100644 --- a/CarpetExtra/IDScalarWaveMoL/src/initialdata.F77 +++ b/CarpetExtra/IDScalarWaveMoL/src/initialdata.F77 @@ -43,5 +43,32 @@ end do end do end do + else if (CCTK_EQUALS(initial_data, "Gaussian")) then + do k=1,cctk_lsh(3) + do j=1,cctk_lsh(2) + do i=1,cctk_lsh(1) + tmp = sqrt( (x(i,j,k)-origin(1))**2 + $ + (y(i,j,k)-origin(2))**2 + $ + (z(i,j,k)-origin(3))**2) + if (tmp .gt. 1.0d-6 * sigma) then + phi(i,j,k) = amplitude/2 + $ * ( (cctk_time + tmp)/tmp * exp (-0.5d0 * ((cctk_time + tmp - radius) / sigma)**2) + $ - (cctk_time - tmp)/tmp * exp (-0.5d0 * ((cctk_time - tmp - radius) / sigma)**2)) + psi(i,j,k) = amplitude/2 + $ * (+ (1 + (cctk_time + tmp) * (cctk_time + tmp - radius) / sigma**2) / tmp * exp (-0.5d0 * ((cctk_time + tmp - radius) / sigma)**2) + $ - (1 + (cctk_time - tmp) * (cctk_time - tmp - radius) / sigma**2) / tmp * exp (-0.5d0 * ((cctk_time - tmp - radius) / sigma)**2)) + else + phi(i,j,k) = amplitude + $ * ( (radius * cctk_time - cctk_time**2 + sigma**2) / sigma**2 + $ + (radius**3 * cctk_time - cctk_time**4 + 6 * cctk_time**2 * sigma**2 - 3 * sigma**4 - 3 * radius**2 * (cctk_time - sigma)**2 + 3 * radius * (cctk_time**3 - 3 * cctk_time * sigma**2)) / (6 * sigma**6) * tmp**2) + $ * exp (-0.5d0 * (radius - cctk_time)**2 / sigma**2) + psi(i,j,k) = amplitude + $ * ( (radius**2 * cctk_time - 2 * radius * cctk_time**2 + cctk_time**3 + 2 * radius * sigma**2 - 3 * cctk_time * sigma**2) / sigma**4 + $ + (radius**4 * cctk_time + cctk_time**5 - 10 * cctk_time**3 * sigma**2 + 15 * cctk_time * sigma**4 - 4 * radius**4 * (cctk_time**2 - sigma**2) + 6 * radius**2 * (cctk_time**3 - 3 * cctk_time * sigma**2) - 4 * radius * (cctk_time**4 - 6 * cctk_time**2 * sigma**2 + 3 * sigma**4)) / (6 * sigma**8) * tmp**2) + $ * exp (-0.5d0 * (radius - cctk_time)**2 / sigma**2) + end if + end do + end do + end do end if end |