diff options
Diffstat (limited to 'src')
-rw-r--r-- | src/InitialData.F77 | 33 |
1 files changed, 13 insertions, 20 deletions
diff --git a/src/InitialData.F77 b/src/InitialData.F77 index 85d1668..7d7463e 100644 --- a/src/InitialData.F77 +++ b/src/InitialData.F77 @@ -33,8 +33,7 @@ DECLARE_CCTK_FARGUMENTS DECLARE_CCTK_PARAMETERS - - INTEGER CCTK_Equals + DECLARE_CCTK_FUNCTIONS INTEGER i,j,k CCTK_REAL dt,omega, pi @@ -45,7 +44,7 @@ omega = sqrt(kx**2+ky**2+kz**2) - if (CCTK_Equals(initial_data,"plane").eq.1) then + if (CCTK_EQUALS(initial_data,"plane")) then do k=1,cctk_lsh(3) do j=1,cctk_lsh(2) @@ -60,22 +59,24 @@ end do end do - else if (CCTK_Equals(initial_data,"gaussian").eq.1) then + 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) + do k=1, cctk_lsh(3) + do j=1, cctk_lsh(2) + do i=1, cctk_lsh(1) - phi(i,j,k) = amplitude*exp( -(sqrt(x(i,j,k)**2 - & +y(i,j,k)**2+z(i,j,k)**2)-radius)**2/sigma**2) - phi_old(i,j,k) = amplitude*exp( -(sqrt(x(i,j,k)**2 - & +y(i,j,k)**2+z(i,j,k)**2)-radius-dt)**2/sigma**2) + phi(i,j,k) = amplitude*exp( + & -(sqrt(x(i,j,k)**2 + & +y(i,j,k)**2+z(i,j,k)**2)-radius)**2/sigma**2) + phi_old(i,j,k) = amplitude*exp( + & -(sqrt(x(i,j,k)**2 + & +y(i,j,k)**2+z(i,j,k)**2)-radius-dt)**2/sigma**2) end do end do end do - else if (CCTK_Equals(initial_data, "box").eq.1) then + else if (CCTK_EQUALS(initial_data, "box")) then c Use kx,ky,kz as number of modes in each direction. @@ -99,14 +100,6 @@ c Use kx,ky,kz as number of modes in each direction. end if -c Apply symmetry boundary conditions -c ---------------------------------- - call ApplySymmetry(cctkGH,"wavetoy::scalarevolve") - -c Synchronise -c ----------- - call CCTK_SyncGroup(cctkGH,"wavetoy::scalarevolve") - return end |