;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;

PRO TW,mode,COLOR=color,THICK=thick

; psf,'TW','m-TW.ps',font=20,scale=0.5
; psf,'TW,4','m-TW.ps',font=20,scale=0.5

str_large=1
@strings
@physconst
spectrum

IF N_ELEMENTS(mode) EQ 0 THEN mode=0

IF N_ELEMENTS(thick) EQ 0 THEN thick=!P.THICK
IF N_ELEMENTS(color) EQ 0 THEN color=0
IF color THEN BEGIN
    c=isocolors(4)
    l=REPLICATE(0,4)
ENDIF ELSE BEGIN
    c=[1,0,1,0]
    l=[0,0,2,2]
    t=[3.5,1,3.5,1]*thick
    grey=180
    tvlct,grey,grey,grey,1
ENDELSE

;ytitle='specific rotational energy / gravitational potential'
;ytitle='e!brot!n / !4U!3!bgrav!n'
ytitle='log(!i !nE!brot!n / |E!bpot!n|!i !n)'
xtitle='m / '+str_msun
xmargin=[6,0]
ymargin=[3.4,0]
xrange=[0,299.999D0]
yrange=[-6.2,0]
logy=1

path='/u/alex/kepler/zero/ed300z0'
dump2=''
legx=''
legy=['T/W','t/w']
nl=N_ELEMENTS(legy)
psel=[1,1]
CASE mode OF 
    0: dump='ed300z0a#presn2'
    1: dump='ed300z0a#presn3'
    2: dump='ed300z0a#presn4'
    3: BEGIN 
        dump=['ed300z0a#presn2','ed300z0a#presn3']
        legx=', '+['T!bc!n!i !n=!i !n10!a10!n!i !nK', $
              str_rho+' !i !n=!i !n10!a10!n!i !ng!i !ncm!a-3!n']
    END
    4: BEGIN 
        dump=['ed300z0a#presn2','ed300z0a#presn3']
        legx=str_rho+'!bc!n!i !n=!i !n'+ $
          ['5.5!i !nx!i !n10!a7!n!i   !n','1.0!i !nx!i !n10!a10!n']+ $
          '!i !ng!i !ncm!a-3!n'
        legy=['','']
        psel=[0,1]
        c=[1,0,1,0]
        l=[0,0,0,2]
        t=[3.5,1,1,1]*thick
        xrange=[0,199.999]
        yrange=[-2.5,0.]
        xmargin=[7.5,0]
        ytitle='log(!i !nt!i !n/!i !nw!i !n)'
    END
ENDCASE

nn=N_ELEMENTS(dump)

jj=0
IF psel(0) THEN BEGIN
    FOR j=0,nn-1 DO BEGIN
        k=0
        i=2*j+k
        name=path+'/'+dump(j)
        loaddump,name
        plotdump,'zm_sun','t/w_ext',color=c(i),line=l(i),THICK=t(i),$
          XTITLE=xtitle,YTITLE=ytitle,$ 
          XMARGIN=xmargin,YMARGIN=ymargin,$ 
          XRANGE=xrange,YRANGE=yrange,$ 
          LOGY=logy,ylog=0, $
          OVER=(jj NE 0)
        jj=jj+1
    ENDFOR    
ENDIF
    
IF psel(1) THEN BEGIN
    FOR j=0,nn-1 DO BEGIN
        k=1
        i=2*j+k
        cd,path,CURRENT=current
        startkepler,'xxx',dump(j),'s'
        plotkepler,'zm','angeek',XSCALE=XMSUN,YSCALE=2, $
          color=c(i),line=l(i),THICK=t(i),$
        XTITLE=xtitle,YTITLE=ytitle,$
          XMARGIN=xmargin,YMARGIN=ymargin,$
          XRANGE=xrange,YRANGE=yrange,$
          LOGY=logy,ylog=0, $
          OVER=(jj NE 0)
        endkepler
        jj=jj+1
    ENDFOR    
ENDIF

oplot,xrange,alog10(0.1376)*[1,1],LINE=1
xyouts,0.9*xrange(1),-0.8,'secular stability limit',ALIGNMENT=1

oplot,xrange,alog10(0.26)*[1,1],LINE=1
xyouts,0.9*xrange(1),-0.52,'dynamical stability limit',ALIGNMENT=1


LEGENDE,0,/LEFT,/TOP
LEGENDE_TEXT,'300!i !n'+str_msun

LEGENDE,6,/left,/bottom,ypos=0.00,LINEHEIGHT=1.5,CHARSIZE=0.9
FOR j=0,nn-1 DO BEGIN
    FOR k=0,nl-1 DO BEGIN
        IF psel(k) NE 0 THEN BEGIN
            i=2*j+k
            leg=legy(k)+legx(j)
            LEGENDE_LINE,l(i),leg,COLOR=c(i),THICK=t(i)
        ENDIF
    ENDFOR
ENDFOR

END

;-----------------------------------------------------------------------
;OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO
;-----------------------------------------------------------------------

PRO kepplot9,file,coremass,YSEL=ysel,PATH=path,COLOR=color,PS=ps

IF N_ELEMENTS(ps) EQ 0 THEN ps=0
IF ps NE 0 THEN BEGIN 
    CASE ps OF
        1: psf,"kepplot9,'ed300z0',PATH='/u/alex/a/kepler/zero/ed300z0'","~/ed300presn.ps",XSIZE=5.4,YSIZE=9,/INCH,RESCALE=2
        2: psf,"kepplot9,'ed250z0',PATH='/u/alex/kepler/zero/ed250z0'","ed250dense.ps",XSIZE=27.5,YSIZE=31.5,scale=0.5
        3: psf,"kepplot9,'ed100z-4',PATH='/u/alex/kepler/ultra/ed100z-4'","presn.ps",XSIZE=27.5,YSIZE=46,scale=0.5
        ELSE:
    ENDCASE
    RETURN
ENDIF

IF N_ELEMENTS(file) EQ 0 THEN file='E12'
IF N_ELEMENTS(coremass) EQ 0 THEN coremass=1
IF N_ELEMENTS(color) EQ 0 THEN color=0

str_large=1
@strings
@physconst

ny=4
nc=7

grey=180
cc=2
TVLCT,grey,grey,grey,cc

IF N_ELEMENTS(color) EQ 0 THEN color=0
IF color THEN BEGIN
    spectrum
    c0=isocolors(nc)
    ci=(INDGEN(nc)*((nc+1) / 2)) mod nc
    ci=c0(ci)
;    c=REPLICATE(ci,ny)
    c=[[ci],[ci],[ci],[ci]]
    l=REPLICATE(0,nc,ny)
    thick=!P.THICK
    t=REPLICATE(thick,nc,ny)
ENDIF ELSE BEGIN
    c=REPLICATE(0,nc,ny)
    l=INDGEN(nc)
    l(nc-1)=0
    l=[[l],[l],[l],[l]]
;    l=REPLICATE(l,ny)
    t=REPLICATE(1.*!P.THICK,nc,ny)
    c(6,*)=cc
    t(6,*)=4.0*!P.THICK
ENDELSE

yticklensave=!Y.TICKLEN
!Y.TICKLEN=0.012

xmargin=[6.25,6.5]
ymargin=3.4
xcharsize=[1.D-12,1.D-12,1.D-12,1.D-12]
yrange=[[-4.2,0],[4.5,10.5],[-1.6,2.6],[13.5,18.5]]
yrange2=[[0,0],[-9,11],[0.35,0.87],[-10.5,1.5]]
ypos=[0.85,0.05,1,1]
xpos=[0.95,0.05,0.95,0.95]
pypos=[0,1,0,1]
pxpos=[1,1,1,1]
ysel_def=[1,1,0,1]
xrange=[0,50]

dump2=''
vexpand=0
CASE file OF 
    'ed300z0': BEGIN xrange=[0,300.0] & yrange(0,1)=4.5001 & yrange2(*,1)=[-10.5,9] & yrange(*,3)=[16,19.8] & yrange2(*,3)=[-8.7,+0.5]  & ysel_def=[1,1,0,1] & dump='ed300z0b#presn2' & xpos(3)=0.25 & END 
;    'ed300z0': BEGIN xrange=[0,300.0] & yrange(0,1)=4.5001 & yrange2(*,1)=[-10.5,9] & yrange(*,3)=[16,19.8] & yrange2(*,3)=[-8.7,+0.5]  & ysel_def=[1,1,0,1] & dump='ed300z0b#presn2' & xpos(3)=0.25 & dump2='ed300z0a#presn3' & END 
    'ed250z0': BEGIN xrange=[0,250.0] & yrange(0,1)=4.5001 & yrange2(*,1)=[-10.5,9] & yrange(*,3)=[16,19.8] & yrange2(*,3)=[-8.7,+0.5]  & ysel_def=[1,1,0,0] & dump='ed250z0#dense' & xpos(3)=0.25 & vexpand=1 & xpos(1)=0.95 & ypos(1)=0.5 & END 
    'ed100z-4': BEGIN xrange=[0,10.0] & yrange(0,1)=4.5001 & yrange2(*,1)=[-10.5,9] & yrange(*,3)=[15,19.8] & & ysel_def=[1,1,1,1] & dump='ed100z-4#collaps' & xpos(0)=0.6 & xpos(3)=0.25 & vexpand=1 & xpos(1)=0.95 & ypos(1)=0.5 & END 
ELSE: 
ENDCASE

IF N_ELEMENTS(ysel) EQ 1 THEN BEGIN
    x=ysel-1
    ysel=REPLICATE(0,4)
    ysel(x)=1
ENDIF

IF N_ELEMENTS(ysel) NE 4 THEN ysel=ysel_def
x=WHERE(ysel NE 0)
IF x(0) EQ -1 THEN BEGIN
    PRINT,' [YSEL]  ERROR'
    RETURN
END
ny=TOTAL(ysel NE 0)

IF (ny EQ 1) AND (ysel(0) EQ 1) THEN xmargin(1)=0

!P.MULTI=[0,1,ny,0,0] 

IF ny LE 2 THEN BEGIN
    ycharsize=1
    xcharsize(ny-1)=1
ENDIF ELSE BEGIN
    ycharsize=2
    xcharsize(ny-1)=2
ENDELSE

xmargin=xmargin*xcharsize(ny-1)
yymargin=ymargin*ycharsize*[[1,0],[2,-1],[3,-2],[4,-3]]/(1.D0*ny)

IF N_ELEMENTS(path) EQ 0 THEN path='/afs/ipp/u/aeh/kepler/link/presn'
IF N_ELEMENTS(dump) EQ 0 THEN dump=file+'#presn'

cd,path,CURRENT=current
startkepler,'xxx',dump,'s'

legp=['A','B','C','D']

jj=0
;AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA
IF ysel(0) EQ 1 THEN BEGIN
j=0

leg=replicate('',nc)
legi=replicate(-1,nc)

plot,[0,1],/NODATA, $
          YRANGE=yrange(*,j),$
          XRANGE=xrange,/XSTYLE,/YSTYLE, $
          YTITLE='log(!i !nmass fraction!i !n)',$
          YMARGIN=yymargin(*,jj),XMARGIN=xmargin,$
          XCHARSIZE=Xcharsize(jj),YCHARSIZE=ycharsize,$
          XTITLE='m / '+str_msun


iso_legi=[0,1,2,3,4,5]
iso_col=[6,0,1,2,3,4]
iso='ion'+['nt1','h1','he4','n14','o16','si28']
iso_leg=['!a  !nn','!a 1!nH','!a 4!nHe','!a14!nN','!a16!nO','!a28!nSi']
niso=N_ELEMENTS(iso)

FOR k=0,niso-1 DO BEGIN
    IF iso(k) NE '' THEN BEGIN
        i=iso_col(k)
        plotkepler,'zm',iso(k),XSCALE=XMSUN,$
          YLOG=0,YRANGE=yrange(*,j),$
          COLOR=c(i,j),LINE=l(i,j),THICK=t(i,j),$
          XRANGE=xrange,logy=1,$
          ytitle='log(!i !nmass fraction!i !n)',$
          YMARGIN=yymargin(*,jj),XMARGIN=xmargin,$
          XCHARSIZE=Xcharsize(jj),YCHARSIZE=ycharsize,$
          XTITLE='m / '+str_msun,/OVER
        legi(iso_legi(k))=i
        leg(i)=iso_leg(k)
    ENDIF
ENDFOR

getkepler,'zm',zm
zm_sun=zm/XMSUN

i=iso_col(1)
getkepler,"ionpn1",ionpn1
pn1=alog10(ionpn1)
oplot,zm_sun,pn1,COLOR=c(i,j),LINE=l(i,j),THICK=t(i,j)

i=5
getkepler,"ion'fe'",ion_fe_
getkepler,'ionfe54',ionfe54
getkepler,'ionni56',ionni56
getkepler,'ionfe56',ionfe56
iron=alog10(ion_fe_+ionfe54+ionni56+ionfe56)
oplot,zm_sun,iron,COLOR=c(i,j),LINE=l(i,j),THICK=t(i,j)
legi(6)=i
leg(i)="'iron'"


nl=N_ELEMENTS(leg)
LEGENDE,4,LENGTH=5,/LEFT,/TOP,YPOS=ypos(j),XPOS=xpos(j), $
  LINEHEIGHT=1.5,CHARSIZE=0.9, $
  FRAMEBACKGROUND=!P.BACKGROUND,ORIENTATION=0, $
  FRAMELINES=nl,FRAMEBORDER=1,FRAMECOLOR=0, $
  DRAWFRAME=1,FRAMETOP=1
FOR li=0,nl-1,1 DO BEGIN
    i=legi(li)
    if i GE 0 THEN BEGIN
        LEGENDE_LINE,l(i,j),leg(i),THICK=t(i,j),COLOR=c(i,j)
    ENDIF
ENDFOR

IF ny NE 1 THEN BEGIN
    LEGENDE,1,/RIGHT,/BOTTOM,YPOS=pypos(jj),XPOS=pxpos(jj),CHARSIZE=1.41
    LEGENDE_TEXT,legp(jj)
ENDIF

jj=jj+1
ENDIF
;BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB
IF ysel(1) EQ 1 THEN BEGIN
j=1
leg=replicate('',nc)
legi=replicate(-1,nc)

IF jj NE 0 THEN !P.MULTI=[(ny-jj),1,ny,0,0]

plot,[0,1],/NODATA, $
  YRANGE=yrange(*,j),$
  XRANGE=xrange,/XSTYLE,YSTYLE=9, $
  YTITLE='log(!i !nT / K!i !n) ,     log(!i !nv!binfall!n / cm!i !ns!a-1!n!i !n)',$
  YMARGIN=yymargin(*,jj),XMARGIN=xmargin,$
  XCHARSIZE=Xcharsize(jj),YCHARSIZE=ycharsize,$
  XTITLE='m / '+str_msun

Y_tmp=!Y

AXIS,/YAXIS,/SAVE, $
  YRANGE=yrange2(*,j),$
  YSTYLE=1, $
  YTITLE='log(!i !n'+str_rho+' / '+str_gccm+'!i !n)',$
  YMARGIN=yymargin(*,jj),XMARGIN=xmargin,$
  XCHARSIZE=Xcharsize(jj),YCHARSIZE=ycharsize,$
  XTITLE='m / '+str_msun

i=6
plotkepler,'zm','dn',XSCALE=XMSUN,$
  YLOG=0,YRANGE=yrange2(*,j),$
  COLOR=c(i,j),LINE=l(i,j),THICK=t(i,j),$
  XRANGE=xrange,logy=1,$
  ytitle='',$
  YMARGIN=yymargin(*,jj),XMARGIN=xmargin,$
  XCHARSIZE=Xcharsize(jj),YCHARSIZE=ycharsize,$
  XTITLE='m / '+str_msun,/OVER
legi(1)=i
leg(i)=str_rho

!Y=Y_tmp

i=0
plotkepler,'zm','tn',XSCALE=XMSUN,$
  YLOG=0,YRANGE=yrange(*,j),$
  COLOR=c(i,j),LINE=l(i,j),THICK=t(i,j),$
  XRANGE=xrange,logy=1,$
  YMARGIN=yymargin(*,jj),XMARGIN=xmargin,$
  XCHARSIZE=Xcharsize(jj),YCHARSIZE=ycharsize,$
  /OVER
legi(0)=i
leg(i)='T'

i=2
plotkepler,'zm','un',XSCALE=XMSUN,$
  YLOG=0,YRANGE=yrange2(*,j), $
  YSCALE=-1., YMAX=-1., $
  COLOR=c(i,j),LINE=l(i,j),THICK=t(i,j),$
  XRANGE=xrange,logy=1,$
  ytitle='',$
  YMARGIN=yymargin(*,jj),XMARGIN=xmargin,$
  XCHARSIZE=Xcharsize(jj),YCHARSIZE=ycharsize,$
  XTITLE='m / '+str_msun,/OVER
legi(2)=i
leg(i)='v!binfall!n'

IF vexpand NE 0 THEN BEGIN
    i=1
    plotkepler,'zm','un',XSCALE=XMSUN,$
      YLOG=0,YRANGE=yrange2(*,j), $
      YSCALE=+1., YMIN=1., $
      COLOR=c(i,j),LINE=l(i,j),THICK=t(i,j),$
      XRANGE=xrange,logy=1,$
      ytitle='',$
      YMARGIN=yymargin(*,jj),XMARGIN=xmargin,$
      XCHARSIZE=Xcharsize(jj),YCHARSIZE=ycharsize,$
      XTITLE='m / '+str_msun,/OVER
    legi(4)=i
    leg(i)='v!bexpand!n'
ENDIF

IF dump2 NE '' THEN BEGIN
    i=3
    loaddump,dump2
    plotdump,'zm_sun','un',/logy,YMUL=-1, $
      LINE=l(i,j),THICK=t(i,j),COLOR=c(i,j), $
      /OVER
    legi(3)=i
    leg(i)='v!binfall!n (!i !n'+str_rho+'!bc!n!i !n=!i !n10!a10!n!i !ng!i !ncm!a-3!n!i !n)'
ENDIF

LEGENDE,6,LENGTH=5,/LEFT,/BOTTOM,LINEHEIGHT=1.5,CHARSIZE=1,XPOS=xpos(1),YPOS=ypos(1)
nl=N_ELEMENTS(leg)
FOR li=0,nl-1,1 DO BEGIN
    i=legi(li)
    if i GE 0 THEN BEGIN
        LEGENDE_LINE,l(i,j),leg(i),THICK=t(i,j),COLOR=c(i,j)
    ENDIF
ENDFOR

IF ny NE 1 THEN BEGIN
    LEGENDE,1,/RIGHT,/BOTTOM,YPOS=pypos(jj),XPOS=pxpos(jj),CHARSIZE=1.41
    LEGENDE_TEXT,legp(jj)
ENDIF

jj=jj+1
ENDIF
;CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
IF ysel(2) EQ 1 THEN BEGIN
j=2

IF jj NE 0 THEN !P.MULTI=[(ny-jj),1,ny,0,0]

plot,[0,1],/NODATA, YSTYLE=4, $
  XRANGE=xrange,/XSTYLE, $
  YMARGIN=yymargin(*,jj),XMARGIN=xmargin,$
  XCHARSIZE=Xcharsize(jj),YCHARSIZE=ycharsize,$
  XTITLE='m / '+str_msun

AXIS,YAXIS=1,/SAVE, $
  XRANGE=xrange,/XSTYLE, $
  YRANGE=yrange2(*,j),$
  YSTYLE=1, $
  YTITLE='electrons per baryon',$
  YMARGIN=yymargin(*,jj),XMARGIN=xmargin,$
  XCHARSIZE=Xcharsize(jj),YCHARSIZE=ycharsize,$
  XTITLE='m / '+str_msun

i=6
plotkepler,'zm','ionye',XSCALE=XMSUN,$
  YLOG=0,YRANGE=yrange2(*,j),$
  COLOR=c(i,j),LINE=l(i,j),THICK=t(i,j),$
  XRANGE=xrange,logy=0,$
  ytitle='',$
  YMARGIN=yymargin(*,jj),XMARGIN=xmargin,$
  XCHARSIZE=Xcharsize(jj),YCHARSIZE=ycharsize,$
  XTITLE='m / '+str_msun,/OVER

LEGENDE,2,LENGTH=5,/RIGHT,/TOP,LINEHEIGHT=1,CHARSIZE=0.9
leg=['Y!be!n']
LEGENDE_LINE,l(i,j),leg,THICK=t(i,j),COLOR=c(i,j)

AXIS,YAXIS=0,/SAVE, $
  YRANGE=yrange(*,j), YSTYLE=1, $
  YTITLE='log(!i !nentropy per baryon / k!bB!n!i !n)',$
  YMARGIN=yymargin(*,jj),XMARGIN=xmargin,$
  XCHARSIZE=Xcharsize(jj),YCHARSIZE=ycharsize,$
  XTITLE='m / '+str_msun

iso='sig'+['','i','e','p','r','ion']
leg=['total','ions','electrons','pairs','radiation','ionization']
FOR i=0,5 DO BEGIN
    plotkepler,'zm',iso(i),XSCALE=XMSUN,$
      YLOG=0,YRANGE=yrange(*,j),$
      COLOR=c(i,j),LINE=l(i,j),THICK=t(i,j),$
      XRANGE=xrange,logy=1,$
      ytitle='log(!i !nmass fraction!i !n)',$
      YMARGIN=yymargin(*,jj),XMARGIN=xmargin,$
      XCHARSIZE=Xcharsize(jj),YCHARSIZE=ycharsize,$
      XTITLE='m / '+str_msun,/OVER
ENDFOR

LEGENDE,6,LENGTH=5,/LEFT,/TOP,LINEHEIGHT=1,CHARSIZE=0.9
nl=6
FOR i=0,nl-1,1 DO BEGIN
    if leg(i) NE '' THEN BEGIN
        LEGENDE_LINE,l(i,j),leg(i),THICK=t(i,j),COLOR=c(i,j)
    ENDIF
ENDFOR

IF ny NE 1 THEN BEGIN
    LEGENDE,1,/RIGHT,/BOTTOM,YPOS=pypos(jj),XPOS=pxpos(jj),CHARSIZE=1.41
    LEGENDE_TEXT,legp(jj)
ENDIF

jj=jj+1
ENDIF
;DDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDD
IF ysel(3) EQ 1 THEN BEGIN
IF jj NE 0 THEN !P.MULTI=[(ny-jj),1,ny,0,0]
j=3

loaddump,dump

leg=replicate('',nc)
legi=replicate(-1,nc)

plot,[0,1],/NODATA, $
  YRANGE=yrange(*,j),$
  XRANGE=xrange,/XSTYLE,YSTYLE=9, $
  YTITLE='log(!i !nj / cm!a2!n!i !ns!a-1!n!i !n)',$
  YMARGIN=yymargin(*,jj),XMARGIN=xmargin,$
  XCHARSIZE=Xcharsize(jj),YCHARSIZE=ycharsize,$
  XTITLE='m / '+str_msun

Y_tmp=!Y

AXIS,/YAXIS,/SAVE, $
  YRANGE=yrange2(*,j),$
  YSTYLE=1, $
  YTITLE='log(!i !n'+str_omega+' / rad!i !ns!a-1!n!i !n)',$
  YMARGIN=yymargin(*,jj),XMARGIN=xmargin,$
  XCHARSIZE=Xcharsize(jj),YCHARSIZE=ycharsize,$
  XTITLE='m / '+str_msun

i=6
plotkepler,'zm','angw',XSCALE=XMSUN,$
  YLOG=0,YRANGE=yrange2(*,j),$
  COLOR=c(i,j),LINE=l(i,j),THICK=t(i,j),$
  XRANGE=xrange,logy=1,$
  ytitle='',$
  YMARGIN=yymargin(*,jj),XMARGIN=xmargin,$
  XCHARSIZE=Xcharsize(jj),YCHARSIZE=ycharsize,$
  XTITLE='m / '+str_msun,/OVER
legi(0)=i
leg(i)=str_omega

!Y=Y_tmp

i=0
plotkepler,'zm','angj',XSCALE=XMSUN, YSCALE=2.d0/3.d0, $
  YLOG=0,YRANGE=yrange(*,j),$
  COLOR=c(i,j),LINE=l(i,j),THICK=t(i,j),$
  XRANGE=xrange,logy=1,$
  ytitle='log(!i !nmass fraction!i !n)',$
  YMARGIN=yymargin(*,jj),XMARGIN=xmargin,$
  XCHARSIZE=Xcharsize(jj),YCHARSIZE=ycharsize,$
  XTITLE='m / '+str_msun,/OVER
legi(4)=i
leg(i)='j!beq!n'

i=1
plotdump,'zm_sun','anglst',/logy, $
  LINE=l(i,j),THICK=t(i,j),COLOR=c(i,j), $
  /OVER
legi(3)=i
leg(i)='j!bLSO!n'

i=2
plotdump,'zm_sun','anglsts',/logy, $
  LINE=l(i,j),THICK=t(i,j),COLOR=c(i,j), $
  /OVER
legi(1)=i
leg(i)='j!bLSO,Schwarzschild!n'

i=3
plotdump,'zm_sun','anglsts',/logy, YMULT=1.d0/3.d0, $
  LINE=l(i,j),THICK=t(i,j),COLOR=c(i,j), $
  /OVER
legi(2)=i
leg(i)='j!bLSO,Kerr!n'


LEGENDE,6,LENGTH=5,/LEFT,/BOTTOM,LINEHEIGHT=1.35,CHARSIZE=1,XPOS=xpos(j)
nl=N_ELEMENTS(leg)
FOR li=0,nl-1,1 DO BEGIN
    i=legi(li)
    if i GE 0 THEN BEGIN
        LEGENDE_LINE,l(i,j),leg(i),THICK=t(i,j),COLOR=c(i,j)
    ENDIF
ENDFOR

IF ny NE 1 THEN BEGIN
    LEGENDE,1,/RIGHT,/BOTTOM,YPOS=pypos(jj),XPOS=pxpos(jj),CHARSIZE=1.41
    LEGENDE_TEXT,legp(jj)
ENDIF

ENDIF

;========================================================================
        
endkepler
cd,current

!P.multi=0
!Y.TICKLEN=yticklensave

END

;-----------------------------------------------------------------------
;OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO
;-----------------------------------------------------------------------

pro homplot,mode,xc

str_large=1
@strings
@physconst

IF N_ELEMENTS(mode) EQ 0 THEN mode=0
IF N_ELEMENTS(xc) EQ 0 THEN xc=0

dir = '~/kepler/zero/ed300z0/'

dumps=dir+'ed300z0a#'+['19000','presn2','20000','presn3','presn4']
l=[4,2,3,0,1]
xcoord='rn'
xr=[0,2.d10]*0.99999D0
yr=[-2.d9,0]
xmargin=[12,0]
ymargin=[3,0]

case mode of
    1: xr=[0,5.d9]*0.99999D0
    2: xr=[0,2.d9]*0.99999D0
else:
endcase

IF xc EQ 1 THEN BEGIN
    xcoord='zm_sun'
    xr=[0,300]*0.99999D0
ENDIF

n=N_ELEMENTS(dumps)

LEGENDE,31,/RIGHT,/BOTTOM,LINEHEIGHT=1.25

FOR i=0,n-1 DO BEGIN
    loaddump,dumps(i)
    plotdump,xcoord,'un',XRANGE=xr,YRANGE=yr,LINE=l(i),OVER=(i NE 0), $
      XMARGIN=xmargin,YMARGIN=ymargin

    getdump,'dn',dn
    getdump,'tn',tn
    leg=str_rho+'!bc!n!i !n=!i !n'+STRING(dn(0),FORMAT='(E7.1)')+ $
      '!i !ng!i !ncm!a-3!n, '+ $
      'T!bc!n!i !n=!i !n'+STRING(tn(0),FORMAT='(E7.1)')+'!i !nK'
    LEGENDE_LINE,l(i),leg
ENDFOR

end