PRO vrot

W= DBLARR(10000)
Z= DBLARR(10000)
vrot= DBLARR(10000)

datax=DBLARR(100)
datay=DBLARR(100)
dataz=DBLARR(100,100)

XW= 0.0D0
XZ=0.0D0
XRHO=0.0D0
XROTT=0.0D0
XROTV = 0.0D0

; choose the datafile. This case: model (6.0,0.9), timestep 785
; (file: time060_090, column 2)
 
OPENR, 1, '../run060_090/ROTV_060_090_00785'

FOR i=0,9999 DO BEGIN

; read the variables

    READF, 1, XW, XZ, XRHO, XROTT, XROTV
    W(i) = XW
    Z(i) = XZ
    vrot(i) = XROTV
ENDFOR

k=0

FOR i=0,99 DO BEGIN
	FOR j=0,99 DO BEGIN
		datax(i)=W(k)
		datay(j)=Z(k)
		dataz(i,j)=vrot(k)
;		PRINT, datax(i), datay(j),dataz(i,j)
		k=k+1
	ENDFOR
ENDFOR

SET_PLOT, 'PS'
DEVICE, /COLOR, file='vrot.ps'

; plot contour. dataz is grid defined through datax and
; datay (vectors)

CONTOUR, dataz,datax,datay,nlevels=20,/ISOTROPIC,$
	xrange=[0,3],yrange=[0,3],/follow, $
	xtitle='!7q!X',ytitle='!8z!X'


DEVICE, /CLOSE
SET_PLOT, 'X'

CLOSE, 1

END