Changed to use timelevels - thanks for putting the rotation in PUGH Ed.

Note that the no-boundary condition tests assumed that phi_new had zero
at its boundaries.  This is, or course, no longer true with timelevel
rotation.

Tom


git-svn-id: http://svn.cactuscode.org/arrangements/CactusWave/IDScalarWaveC/trunk@6 071ba4cb-2270-47cd-a45a-3e542e115664

1 files changed, 6 insertions, 6 deletions

diff --git a/src/InitialData.c b/src/InitialData.c
index c03b04b..0d9b384 100644
--- a/src/InitialData.c
+++ b/src/InitialData.c
@@ -72,7 +72,7 @@ void IDScalarWaveC_InitialData(CCTK_ARGUMENTS)
           index =  CCTK_GFINDEX3D(cctkGH,i,j,k);
 
           phi[index]     = amplitude*cos(kx*x[index]+ky*y[index]+kz*z[index]+omega*cctk_time);
-          phi_old[index] = amplitude*cos(kx*x[index]+ky*y[index]+kz*z[index]+omega*(cctk_time-dt));
+          phi_p[index] = amplitude*cos(kx*x[index]+ky*y[index]+kz*z[index]+omega*(cctk_time-dt));
         }
       }
     }
@@ -94,7 +94,7 @@ void IDScalarWaveC_InitialData(CCTK_ARGUMENTS)
           R = sqrt(X*X + Y*Y + Z*Z);
 
           phi[index] = amplitude*exp( - sqr( (R - radius) / sigma ) );
-          phi_old[index] = amplitude*exp( - sqr( (R - radius - dt) / sigma ) );
+          phi_p[index] = amplitude*exp( - sqr( (R - radius - dt) / sigma ) );
 
         } 
       }
@@ -118,10 +118,10 @@ void IDScalarWaveC_InitialData(CCTK_ARGUMENTS)
                                  sin(kz*(z[index]-0.5)*pi)*
                                  cos(omega*cctk_time*pi);
 
-          phi_old[index] = amplitude*sin(kx*(x[index]-0.5)*pi)*
-                                         sin(ky*(y[index]-0.5)*pi)*
-                                         sin(kz*(z[index]-0.5)*pi)*
-                                         cos(omega*(cctk_time-dt)*pi);
+          phi_p[index] = amplitude*sin(kx*(x[index]-0.5)*pi)*
+                                       sin(ky*(y[index]-0.5)*pi)*
+                                       sin(kz*(z[index]-0.5)*pi)*
+                                       cos(omega*(cctk_time-dt)*pi);
         }
       }
     }