sandbox/viscous_collapse.plot

    1
    2
    3
    4
    5
    6
    7
    8
    9
    10
    11
    12
    13
    14
    15
    16
    17
    18
    19
    20
    21
    22
    23
    24
    25
    26
    27
    28
    29
    30
    31
    32
    33
    34
    35
    36
    37
    38
    39
    40
    41
    42
    43
    44
    45
    
    /** gnuplot */
    
    !awk '{ if ($1 == "file:") file = $2; else print $0 > file; }' < viscous_collapse.log
    
    set term png
    
    set output viscous_collapse1.png
    
    reset
    set title 'Self Similar Solution for a viscous collapse, t=200, 400...1000'
    
    
    set xlabel 'x'
    set ylabel 'h(x t^(-1/3),t) t^(1/3)'
    
    s=2.5
    b=2.18812
    h(x) = 0.543217  * (1-x*x/b/b)**(1./3)
    bed(x) = 0
    t5=200**.2
    
    set key top right
    
    plot [-3:3] \
          '< grep ^p viscous_collapse.out' u 2:3 w l lc 3 t 'Numerical', \
          h(x) lw 1 lc 1 lt 1 t 'self similar',\
          'viscous_collapse.ref' u ($1/t5):($2*t5) t'order 1, HLL t=200 N=512' w l  lc 2
    
    
    # pause -1
    set output viscous_collapse2.png
    reset
    set xlabel 'x'
    set ylabel 'h(x,t)'
    set arrow from 0,.27 to 0,0.08 
    set label "t" at 0.2,0.1
    
    set arrow from 3,.05 to 10,0.05 
    set label "t" at 10,0.06 
    
    plot [][:.3]'eta-0' not  w l lc -1,'eta-1'not  w l lc -1,\
        'eta-2'not  w l lc -1,'eta-3'not  w l lc -1,'eta-4'not  w l lc -1,'eta-5'not  w l lc -1