$Id: README,v 1.25 2022/10/05 01:20:48 brandenb Exp $

tex directories:
  tex/tina/GW
  tex/roper/GW
  tex/roper/helical
  tex/roper/PencilGW
  tex/sayan/GW/fig
  tex/tina/CircPol
  tex/rei/GW

GitHub:
  ~/GitHub/Ramkishor/GWspectra_from_mhd_turbulence/

Run for kinematic flows:
  kin16a

Analysis tool:
  IDL> .r ptaver

M1152e3_exp6
M1152c3_exp5

date: 2019-09-03 revision 1.38 was probably ok

F288e_sig1_t11_M4
freq,k[1],OmegaGW[1]=   2.6210769e-06       100.00000   1.3728826e-11

F1152d2_sig1_t11_M2c
freq,k[1],OmegaGW[1]=   2.6210769e-06      100.000   6.6163852e-13
freq,k[1],OmegaGW[1]=   2.6210769e-06      100.000   6.6163852e-13


/cfs/klemming/nobackup/b/brandenb/pencil-code/sayan/GW/512_1e2_2e4_1e8_5em6a
crashed, Thu Aug  6 06:24:14 CEST 2020

/cfs/klemming/nobackup/b/brandenb/pencil-code/sayan/GW/512_1e2_5e4_1e8_5em6a
crashed, Thu Aug  6 06:25:33 CEST 2020

/cfs/klemming/nobackup/b/brandenb/pencil-code/sayan/GW/512_1e2_1e4_1e8_5em6b
OK, was to demonstrate that intial seed doesn't matter; Thu Aug  6 06:28:05 CEST 2020

/cfs/klemming/nobackup/b/brandenb/pencil-code/sayan/GW/512_2e2_2e4_1e8_5em6a
OK, Thu Aug  6 06:30:55 CEST 2020

/cfs/klemming/nobackup/b/brandenb/pencil-code/sayan/GW/512_5e2_5e4_1e8_5em6a
crashed, Thu Aug  6 06:43:13 CEST 2020

SVT decomposition
-rwxrwxrwx 1 root root 29K 2020-11-05 20:56 pSCL2_M1152e_exp6k4_M4b_rerun.ps
-rwxrwxrwx 1 root root 30K 2020-11-05 20:56 pSCL2_F1152a_sig1_t11d_double.ps
-rwxrwxrwx 1 root root 37K 2020-11-05 20:56 pSCL2_E1152e_t11_M4d_double.ps

Run B10: 512_1e2_1e4_4e8_1em3a
Run B1: 512_1e2_1e4_4e8_1em6a
        512_1e2_1e4_4e8_1em6aD

Runs from Fig.1 if Kahniashvili+21
  F1152a_sig0_t11_M4_ramp01b
  F1152a_sig0_t11_M4_ramp02a
  F1152a_sig0_t11_M4_ramp05a
  F1152a_sig0_t11_M4_ramp1a
  F1152a_sig0_t11_M4_ramp2a

file='pe_2em2_kf6.txt'  ;(kf6_red)
file='pe_5em3_kf6.txt'  ;(kf6_blue)
file='pe_5em5.txt'  ;(kf600_red)
file='pe_ramp0_5em5.txt' ;(kf600_red_dotted)
file='pe_1em4.txt' ;(kf600_orange)
file='pe_2em6.txt' ;(kf600_blue)
file='pe_2em7.txt' ;(kf600_black)

-rw-r--r-- 1 brandenb brandenb 1,4K sep 21 04:47 wtab_wkf6.pro
cat pe_2em2_kf6.txt pe_5em3_kf6.txt >> pe_wkf6.txt

-rw-r--r-- 1 brandenb brandenb 1,4K sep  9 22:11 wtab_line.pro
cat pe_5em5.txt pe_ramp0_5em5.txt pe_1em4.txt pe_2em6.txt pe_2em7.txt >> pe_line.txt

pc_diffruna M512sig1_k6_ramp1b_f007 M512sig1_k6_ramp0b2_f007

                           k                   color
ik+dik,col,lct[dik]=       3     155       0   orange
ik+dik,col,lct[dik]=       4     122       0   red
ik+dik,col,lct[dik]=       5     122       1   green
ik+dik,col,lct[dik]=       6      55       2   blue
ik+dik,col,lct[dik]=       7     255       0   black


To do k-omega diagram, start with:
  ../idl/rslice_stress.pro
This yields fff1, 2, 3, and 4, and saves fff=fff1+i*fff2 as save file. Next, do
  .r ../idl/calc_ko.pro

vortical example; use t=2-3
  K512sig0_k6_ramp1a_cont

irrotational:
  K512irrok6c_rep1_cont