To read and write integrand data, write:
  .r rbackreact.pro

kmax  Ht  run
0.20  21  dt01_dk5em3_nproc8c  
0.20  21  dt02_dk5em3_nproc8a

dt01_dk5em3_nproc8c
grandxy.sav


cvs up -d dk1em3_Q1p65em4_mu09em4a
cvs up -d dk1em3_Q1p9em4_mu1em4a
cvs up -d dk1em3_Q2p7em4_mu1p3em4a
cvs up -d dk1em3_Q3p2em4_mu1p5em4a
cvs up -d dk1em3_Q4p1em4_mu1p8em4a
cvs up -d dk1em3_Q5p1em4_mu2p1em4a
cvs up -d dk1em3_Q6em4_mu2p45em4a
cvs up -d dk1em3_Q7em4_mu2p76em4a

cvs up -d dk1em3_Q1p65em4_mu09em4b
cvs up -d dk1em3_Q1p9em4_mu1em4b
cvs up -d dk1em3_Q2p7em4_mu1p3em4b
cvs up -d dk1em3_Q3p2em4_mu1p5em4b
cvs up -d dk1em3_Q4p1em4_mu1p8em4b
cvs up -d dk1em3_Q5p1em4_mu2p1em4b
cvs up -d dk1em3_Q6em4_mu2p45em4b
cvs up -d dk1em3_Q7em4_mu2p76em4b

cvs up -d dk1em3_Q6p9em4_mu2p725em4
cvs up -d dk1em3_Q9em4_mu3p45em4
cvs up -d dk1em3_Q6p3em4_mu2p525em4
cvs up -d dk1em3_Q4p5em4_mu1p93em4
cvs up -d dk1em3_Q5p6em4_mu2p295em4
cvs up -d dk1em3_Q10em4_mu3p85em4
cvs up -d dk1em3_Q5p9em4_mu2p395em4
cvs up -d dk1em3_Q8em4_mu3p1em4
cvs up -d dk1em3_Q6p6em4_mu2p625em4
cvs up -d dk1em3_Q4p8em4_mu2p03em4
cvs up -d dk1em3_Q5p3em4_mu2p195em4
cvs up -d dk1em3_Q7em4_mu2p76em4b

   6.9& $6.9*10^{-4}$ & $2.725 *10^{-4}$ \\ [1ex] % dk1em3_Q6p9em4_mu2p725em4
   6.6& $6.6*10^{-4}$ & $2.625 *10^{-4}$ \\ [1ex] % dk1em3_Q6p6em4_mu2p625em4
   6.3& $6.3*10^{-4}$ & $2.525 *10^{-4}$ \\ [1ex] % dk1em3_Q6p3em4_mu2p525em4
   5.9& $5.9*10^{-4}$ & $2.395 *10^{-4}$ \\ [1ex] % dk1em3_Q5p9em4_mu2p395em4
   5.6& $5.6*10^{-4}$ & $2.295 *10^{-4}$ \\ [1ex] % dk1em3_Q5p6em4_mu2p295em4
   5.3& $5.3*10^{-4}$ & $2.195 *10^{-4}$ \\ [1ex] % dk1em3_Q5p3em4_mu2p195em4
   4.8& $4.8*10^{-4}$ & $2.03 *10^{-4}$ \\ [1ex] % dk1em3_Q4p8em4_mu2p03em4
   4.5& $4.5*10^{-4}$ & $1.93 *10^{-4}$ \\ [1ex] % dk1em3_Q4p5em4_mu1p93em4

dk5em3_Q4p1em4_mu1p8em4_ddot_rest

dk5em3_Q4p1em4_mu1p8em4_ddot_rest2
dt=.1


/cfs/klemming/home/b/brandenb/scr/oksana/axionSU2back/Q3p2em4_mu1p5em4_2048a_nmin0m12_nmax04_wshift
3.2 run with -12 to 4 and shift, expect the multiple peaks to extend further in time and to cause oscillations.
currently N=19, still no second peak,
if smaller k, range, BR disappears, will affect BG evolution, later amplifies modes which enter later

/cfs/klemming/home/b/brandenb/scr/oksana/axionSU2back/Q3p2em4_mu1p5em4_2048a_nmin0m6_nmax04_wshift
same, but with less k range
second bump appears at lnk=5; could this be because of better resolution?

/cfs/klemming/home/b/brandenb/scr/oksana/axionSU2back/Q6p0em4_mu2p45em4_2048a_nmin0m0_nmax012
6.0 run, with 12 and same dt
still only one bump

/cfs/klemming/home/b/brandenb/scr/oksana/axionSU2back/Q6p0em4_mu2p45em4_2048a_nmin0m0_nmax018
same, but with 18
NaNs

/cfs/klemming/home/b/brandenb/scr/oksana/axionSU2back/Q6p0em4_mu2p45em4_2048a_nmin0m0_nmax018_dt
shorter time, also still NaNs (time step should have been even shorter)

-rw-r--r-- 1 brandenb brandenb 435 aug  7 14:14 Q3p2em4_mu1p5em4_4096a_nmin0m12_nmax04_wshift/run.in
-rw-r--r-- 1 brandenb brandenb 435 aug  7 14:13 Q3p2em4_mu1p5em4_2048a_nmin0m6_nmax04_wshift/run.in
-rw-r--r-- 1 brandenb brandenb 435 aug  7 14:11 Q3p2em4_mu1p5em4_2048a_nmin0m12_nmax04_wshift/run.in
-rw-r--r-- 1 brandenb brandenb 435 aug  7 07:13 Q3p2em4_mu1p5em4_2048a_nmin0m0_nmax012_noshift2/run.in  <==
-rw-r--r-- 1 brandenb brandenb 435 aug  7 00:15 Q3p2em4_mu1p5em4_2048a_nmin0m0_nmax012_oshift/run.in
-rw-r--r-- 1 brandenb brandenb 435 aug  7 00:06 Q3p2em4_mu1p5em4_2048a_nmin0m0_nmax012_noshift/run.in
-rw-r--r-- 1 brandenb brandenb 435 aug  6 18:53 Q3p2em4_mu1p5em4_2048a_nmin0m2_nmax04_wshift/run.in
-rw-r--r-- 1 brandenb brandenb 435 aug  6 17:17 Q3p2em4_mu1p5em4_2048a_nmin0m0_nmax06_wshift/run.in
-rw-r--r-- 1 brandenb brandenb 435 aug  6 15:24 Q3p2em4_mu1p5em4_2048a_nmin0m0_nmax06/run.in


/cfs/klemming/home/b/brandenb/scr/oksana/axionSU2back/ Q6p0em4_mu2p45em4_2048a_nmin0m0_nmax06_rep
to check same as yesterday
result is the same.

/cfs/klemming/home/b/brandenb/scr/oksana/axionSU2back/ Q3p2em4_mu1p5em4_2048a_nmin0m0_nmax015_noshift2
larger  range, more proces
also still needs more time
2063884 start_ru          4  1024    17:30:53       None     R               (null)            2023-08-08T13:59:39            2023-08-08T13:59:38

/cfs/klemming/home/b/brandenb/scr/oksana/axionSU2back/ dk5em3_Q3p2em4_mu1p5em4_dt3b
linear k range, to reproduce January results
still running, need to wait, but now ready. No turn upward. Is now in Notes

/cfs/klemming/home/b/brandenb/scr/oksana/axionSU2back/ dk5em6_Q3p2em4_mu1p5em4
k0=5e-6, otherwise same
develops along horizon line. Also no upturn, is now in Notes

2063926 start_ru          4  1024    17:50:06       None     R               (null)            2023-08-08T14:18:52            2023-08-08T14:18:51
2063857 start_ru          4  1024    17:18:20       None     R               (null)            2023-08-08T13:47:06            2023-08-08T13:47:05
2063716 start_ru          1   256    15:59:41       None     R               (null)            2023-08-08T12:28:27            2023-08-08T12:28:26



drwxr-xr-x 5 brandenb brandenb 4,0K aug  9 09:22 dk5em3_Q3p2em4_mu1p5em4_dt2
drwxr-xr-x 5 brandenb brandenb 4,0K aug  9 07:11 Q3p2em4_mu1p5em4_2048a_nmin0m0_nmax015_noshift2
drwxr-xr-x 5 brandenb brandenb 4,0K aug  9 04:48 dk5em6_Q3p2em4_mu1p5em4
drwxr-xr-x 3 brandenb brandenb 4,0K aug  8 20:43 idl
drwxr-xr-x 5 brandenb brandenb 4,0K aug  8 20:37 Q6p0em4_mu2p45em4_2048a_nmin0m0_nmax06_rep
drwxr-xr-x 4 brandenb brandenb 4,0K aug  8 13:52 Q3p2em4_mu1p5em4_2048a_nmin0m0_nmax015_noshift
drwxr-xr-x 4 brandenb brandenb 4,0K aug  8 13:35 dk5em3_Q3p2em4_mu1p5em4_dt_rep


   1  axi_Q                                   -0.00013644684

   2  axi_Qdot                                0.0014536784

   3  axi_chi                                 0.0051878765

   4  axi_chidot                              0.82380156

   5  axi_psi                                 3.6378930e+13

   6  axi_psidot                              3.6378957e+18

   7  axi_TR                                  -1.6617992e+11

   8  axi_TRdot                               -1.7979153e+16

   9  axi_impsi                               1.9482402e+13

  10  axi_impsidot                            1.9482417e+18

  11  axi_imTR                                -8.8995680e+10

  12  axi_imTRdot                             -9.6285257e+15

