<html>#2413: Include "POWER" code to extrapolate waveforms to scri+ into ET
<table style='border-spacing: 1ex 0pt; '>
<tr><td style='text-align:right'> Reporter:</td><td>Roland Haas</td></tr>
<tr><td style='text-align:right'> Status:</td><td>open</td></tr>
<tr><td style='text-align:right'>Milestone:</td><td></td></tr>
<tr><td style='text-align:right'> Version:</td><td></td></tr>
<tr><td style='text-align:right'> Type:</td><td>task</td></tr>
<tr><td style='text-align:right'> Priority:</td><td>major</td></tr>
<tr><td style='text-align:right'>Component:</td><td>EinsteinToolkit thorn</td></tr>
</table>
<p>Comment (by Maria):</p>
<p>Roland, I added the parameters for outputting psi4.h5 to the nsnstohmns.par. This is what I added:</p>
<p>______________________________<br />
WeylScal4::fdOrder = 8<br />
WeylScal4::calc_scalars = "psis"<br />
WeylScal4::calc_invariants = "always"</p>
<p>Multipole::nradii = 7<br />
Multipole::radius[0] = 100<br />
Multipole::radius[1] = 115<br />
Multipole::radius[2] = 136<br />
Multipole::radius[3] = 167<br />
Multipole::radius[4] = 214<br />
Multipole::radius[5] = 300<br />
Multipole::radius[6] = 500<br />
Multipole::ntheta = 120<br />
Multipole::nphi = 240<br />
Multipole::variables = "WeylScal4::Psi4r{sw=-2 cmplx='WeylScal4::Psi4i' name='psi4'}"<br />
Multipole::out_every = 512<br />
Multipole::l_max = 4 <a href="https://bitbucket.org/einsteintoolkit/tickets/issues/8/unable-to-do-remote-submission-without" rel="nofollow" title="Unable to do remote submission without configuring your local machine" class="ap-connect-link"><s>#8</s></a><br />
Multipole::output_hdf5 = yes</p>
<p>Multipole::output_ascii = no</p>
<p>___________________________</p>
<p>One other change I made was to increase the boundary points and ghost points to 4.</p>
<p>I do have a problem with the nsnstohmns.par grid settings. The way is it now: N=396, dx = 18</p>
<p>I am getting interpolation error:</p>
<p>WARNING level 1 from host babiuct7910-1 process 0<br />
while executing schedule bin CCTK_ANALYSIS, routine Multipole::Multipole_Calc<br />
in thorn AEILocalInterp, file /home/babiuc/Cactus/arrangements/Numerical/AEILocalInterp/src/Hermite/../template.c:1109:<br />
-><br />
CCTK_InterpLocalUniform():<br />
interpolation point is either outside the grid,<br />
or inside but too close to the grid boundary!<br />
(this may be caused by a global interpolation with<br />
driver::ghost_size too small)<br />
0-origin interpolation point number pt=13469 of N_interp_points=13470<br />
interpolation point (x,y,z)=(237.748,116.683,424.1)<br />
grid x_
<p>--<br/>
Ticket URL: <a href='https://bitbucket.org/einsteintoolkit/tickets/issues/2413/include-power-code-to-extrapolate'>https://bitbucket.org/einsteintoolkit/tickets/issues/2413/include-power-code-to-extrapolate</a></p>
</html>