[Commits] [svn:einsteintoolkit] NullEvolve/branches/tapir/ (Rev. 22)

[bela at caltech.edu]

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: