PRO df

Y = DBLARR(6100)
R = DBLARR(6100)
E = DBLARR(6100)
Jz=DBLARR(6100)
distf=DBLARR(6100)

datax = DBLARR(61)
datay = DBLARR(61)


XY=0.0D0
XR=0.0D0
XE=0.0D0
XJ= 0.0D0
GEJ= 0.0D0

; read the file containing variables 

OPENR, 1, '../run060_090/DF_060_090_00001'

FOR i=0,6099 DO BEGIN

    READF, 1, XY, XR, XE, XJ, GEJ
    Y(i) = XY
    R(i) = XR
    E(i) = XE
    Jz(i) = XJ	
    distf(i) = GEJ

ENDFOR

j=0

FOR i=0,60 DO BEGIN
	datax(j)=Jz(i)
	datay(j)=distf(i)
;	PRINT, datax(j), datay(j)
	j=j+1
ENDFOR

SET_PLOT, 'PS'
DEVICE, /COLOR, file='df.ps'
	
PLOT, datax,datay, $
  	xtitle='J!Dz!N', ytitle='DF(E,J!Dz!N)',$
	YRANGE=[0.0000001,1], /YLOG, YSTYLE = 1, $
	XRANGE=[-2,2]
XYOUTS, 3000,5000, '!7x!X!D0!N=0.9', CHARSIZE=1.1,/DEVICE
XYOUTS, 3000,6000, 'Wo=6.0', CHARSIZE=1.1,/DEVICE

; choose some points of the grid (total set: 61 a.m. values for each
; energy value). We choose 17 of 100 energy values.

j=0

FOR i=61,6099 DO BEGIN

JUMP:  	datax(j)=Jz(i)
	datay(j)=distf(i)
	PRINT, datax(j), datay(j),j,i
	j=j+1	
	IF (j EQ 61) THEN BEGIN
		OPLOT, datax, datay, linestyle=0
		PRINT, 'overploted'
		j=0
		i=i+306
		IF i GE 6100 THEN GOTO, JUMP2
		PRINT, i,j
		GOTO, JUMP
	ENDIF

ENDFOR

JUMP2: PRINT, 'END'

DEVICE, /CLOSE
SET_PLOT, 'X'

CLOSE, 1

END