[Commits] [svn:einsteintoolkit] Paper_EinsteinToolkit_2010/examples/tov/ (Rev. 76)

[knarf at cct.lsu.edu]

User: knarf
Date: 2011/04/28 12:47 PM

Modified:
 /examples/tov/
  tov.rpar

Log:
 add perturbed tov parameters (not in ET yet

File Changes:

Directory: /examples/tov/
=========================

File [modified]: tov.rpar
Delta lines: +13 -3
===================================================================
--- examples/tov/tov.rpar	2011-04-28 17:23:25 UTC (rev 75)
+++ examples/tov/tov.rpar	2011-04-28 17:47:46 UTC (rev 76)
@@ -1,9 +1,18 @@
 #!/usr/bin/perl -w
 use strict;
 
+foreach my $type (("", "pert")) {
 foreach my $res ((16, 8, 4, 2)) {
 
   my $out_every = 256/$res;
+  my $K_TOV    = 100;
+  my $K_change = "";
+  if ($type eq "pert") {
+    $K_TOV     = 99;
+    $K_change  = "TOVSolver::TOV_Transform_to_K = 100";
+  }
+
+
   my $par = <<PAR_EOF;
 ## mclachlan tov_static
 ActiveThorns = "Time MoL"
@@ -184,7 +193,8 @@
 
 TOVSolver::TOV_Rho_Central[0] = 1.28e-3
 TOVSolver::TOV_Gamma[0]       = 2.0
-TOVSolver::TOV_K[0]           = 100.0
+TOVSolver::TOV_K[0]           = $K_TOV
+$K_change
 
 # I/O
 
@@ -225,8 +235,8 @@
 "
 PAR_EOF
 
-  open(OUT, ">tov_${res}.par");
+  open(OUT, ">tov_${res}_${type}.par");
   print OUT "$par";
   close(OUT);
-}
+}}