PRO dtime

;initialize variables

THTIO=0.0D0
TCTIO=0.0D0
RCORE=0.0D0
VROT=0.0D0
XSIGT=0.0D0
XRHO=0.0D0
EDYN=0.0D0
SEE=0.0D0
XI=0.0D0
XMTOT=0.0D0
XJTOT=0.0D0
XPTOT=0.0D0
XKTOT=0.0D0

; we change directories and read each file, which includes the
; parameters

OPENR, 1, '../run060_000/time060_000'

datax=DBLARR(48)
datay=DBLARR(48)

; read the variables

FOR i=0,47 DO BEGIN

  READF, 1, THTIO,TCTIO,RCORE,EDYN,SEE,XI,XMTOT,XJTOT,XPTOT,XKTOT, $
            EDYN,SEE,XI,XMTOT,XJTOT,XPTOT,XKTOT
    datax(i) = THTIO
    datay(i)= XRHO

ENDFOR

SET_PLOT, 'PS'
DEVICE, /COLOR, file='dtime.ps'
LOADCT, 4	

PLOT, datax,datay, $
  	xtitle='t/t!Drh!N', ytitle='density',$
	/YLOG, color=1, $
	XRANGE=[0,12], YRANGE=[0.01,1000]

XYOUTS, 6000,8000, '!7x!X!D0!N=0.9', CHARSIZE=1.1,/DEVICE
XYOUTS, 9000,8000, '!7x!X!D0!N=0.6', CHARSIZE=1.1,/DEVICE
XYOUTS, 12500,8000, '!7x!X!D0!N=0.3', CHARSIZE=1.1,/DEVICE
XYOUTS, 14500,8000, '!7x!X!D0!N=0.0', CHARSIZE=1.1,/DEVICE
XYOUTS, 14000,3000, 'Wo=6.0', CHARSIZE=1.1,/DEVICE
CLOSE,1
;------------------------------------------------------------
OPENR, 1, '../run060_030/time060_030'

datax=DBLARR(56)
datay=DBLARR(56)

FOR i=0,55 DO BEGIN

  READF, 1, THTIO,TCTIO,RCORE,EDYN,SEE,XI,XMTOT,XJTOT,XPTOT,XKTOT, $
            EDYN,SEE,XI,XMTOT,XJTOT,XPTOT,XKTOT
     datax(i)= THTIO
     datay(i)= XRHO

ENDFOR
OPLOT, datax, datay, color=100
CLOSE,1
;------------------------------------------------------------
OPENR, 1, '../run060_060/time060_060'

datax=DBLARR(59)
datay=DBLARR(59)

FOR i=0,58 DO BEGIN

   READF, 1, THTIO,TCTIO,RCORE,EDYN,SEE,XI,XMTOT,XJTOT,XPTOT,XKTOT, $
            EDYN,SEE,XI,XMTOT,XJTOT,XPTOT,XKTOT
     datax(i)= THTIO
     datay(i)= XRHO

ENDFOR
OPLOT, datax, datay, color=120
CLOSE,1
;------------------------------------------------------------
OPENR, 1, '../run060_090/time060_090'

datax=DBLARR(66)
datay=DBLARR(66)

FOR i=0,65 DO BEGIN

  READF, 1, THTIO,TCTIO,RCORE,EDYN,SEE,XI,XMTOT,XJTOT,XPTOT,XKTOT, $
            EDYN,SEE,XI,XMTOT,XJTOT,XPTOT,XKTOT
     datax(i)= THTIO
     datay(i)= XRHO

ENDFOR
OPLOT, datax, datay, color=140
CLOSE,1
;------------------------------------------------------------
DEVICE, /CLOSE
SET_PLOT, 'X'

END