PRO edynage

;initialize variables

THTIO=0.0D0
EDYN=0.0D0
TCTIO=0.0D0
RCORE=0.0D0
EDYN=0.0D0
SEE=0.0D0
XI=0.0D0
XMTOT=0.0D0
XJTOT=0.0D0
XPTOT=0.0D0
XKTOT=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 file, which includes the
; parameters

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

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

; read the variables described in the gridfiles

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) = EDYN
ENDFOR


SET_PLOT, 'PS'
DEVICE, file='edynage.ps', /COLOR
LOADCT, 4
PLOT, datax,datay, xCHARSIZE=1.4, yCHARSIZE=1.4, color=1, $
  	xtitle='t/t!Drh!N', ytitle='e!Ddyn!N',linestyle=0, $
	XRANGE=[-2,12], YRANGE=[-0.02,0.4], xstyle=1, ystyle=1

XYOUTS, 2500,11500, '!7x!X!D0!N=0.9', CHARSIZE=1.1,/DEVICE
XYOUTS, 2500,9000, '!7x!X!D0!N=0.6', CHARSIZE=1.1,/DEVICE
XYOUTS, 2500,4500, '!7x!X!D0!N=0.3', CHARSIZE=1.1,/DEVICE
XYOUTS, 2500,2000, '!7x!X!D0!N=0.0', CHARSIZE=1.1,/DEVICE
XYOUTS, 12000,7000, 'Wo=6.0', CHARSIZE=1.1,/DEVICE

CLOSE, 1

; now overplot next data

; -------------------------------------------------
OPENR, 1, '../run060_030/time060_030_02' ;file to open
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) = EDYN
ENDFOR

OPLOT, datax, datay, linestyle=0, color=40
CLOSE, 1
; -------------------------------------------------
OPENR, 1, '../run060_060/time060_060_02' ;file to open
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) = EDYN
ENDFOR

OPLOT, datax, datay, linestyle=0, color=60
CLOSE, 1
; -------------------------------------------------
OPENR, 1, '../run060_090/time060_090_02' ;file to open
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) = EDYN
ENDFOR

OPLOT, datax, datay, linestyle=0, color=80
CLOSE, 1
; -------------------------------------------------


DEVICE, /CLOSE
SET_PLOT, 'X'

END