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); -} +}}