User: szilagyi Date: 2014/01/03 02:14 AM

Modified: /branches/tapir/ param.ccl /branches/tapir/src/ NullEvol_InitialData.F90

Log: implement J=const (radially) as an option for initial data

File Changes:

Directory: /branches/tapir/src/ ===============================

File [modified]: NullEvol_InitialData.F90 Delta lines: +7 -0 =================================================================== --- branches/tapir/src/NullEvol_InitialData.F90 2014-01-03 07:20:33 UTC (rev 21) +++ branches/tapir/src/NullEvol_InitialData.F90 2014-01-03 08:14:54 UTC (rev 22) @@ -39,6 +39,13 @@ jcn = (0.,0.) jcs = (0.,0.)

+ else if (CCTK_EQUALS(initial_J_data,"constant_J")) then + + do i = 1, N_radial_pts + jcn(:,:,i) = j_wt(:,:,1) + jcs(:,:,i) = j_wt(:,:,2) + end do + else if (CCTK_EQUALS(initial_J_data,"smooth_J")) then

do i = 1, N_radial_pts

Directory: /branches/tapir/ ===========================

File [modified]: param.ccl Delta lines: +1 -0 =================================================================== --- branches/tapir/param.ccl 2014-01-03 07:20:33 UTC (rev 21) +++ branches/tapir/param.ccl 2014-01-03 08:14:54 UTC (rev 22) @@ -43,6 +43,7 @@ "extracted_J" :: "J as given by extraction" "polynomial_J" :: "J as third power polynomial to vanish J_xx at scri" "fitted_linearized_J" :: "linearized fitted J based on some previous run" + "constant_J" :: "J=const in the radial direction: J_ev=J_wt for all radial points" } "vanishing_J_scri"

private: