dirname='/data/arc1/cat/RCM-MURI/RCM-2013017/run/'


I1=78  ; RCM grid
J1=51 ; RCM grid




hr=0
mn=1

;++++++++++++++++++++++++ generate the file name that match the output name ++++++++++++++++++++++++
tname='time '+string(format='(I2.2)',hr)+':'+string(format='(I2.2)',mn)
fname=string(format='(I6.6)',hr*10000L+mn*100L)

print,'rcm-output-3D-IDL-'+fname+'.dat'

para=fltarr(12,I1,J1)
eta=fltarr(I1,J1,200)
lamda=fltarr(I1,J1,200)
veff=fltarr(I1,J1)
header=''

openr,unit,dirname+'/rcm-output-3D-'+fname+'.dat',/get_lun

for i=0,13 do begin
readf,unit,header
endfor

for k=0,199 do begin

if k eq 0 then begin

 for j=0,J1-1 do begin
for i=0,I1-1 do begin 
temp=fltarr(14)
readf,unit,temp

para(0,i,j)=temp(0)
para(1,i,j)=temp(1)
para(2,i,j)=temp(2)
para(3,i,j)=temp(3)
para(4,i,j)=temp(4)
para(5,i,j)=temp(5)
para(6,i,j)=temp(6)
para(7,i,j)=temp(7)
para(8,i,j)=temp(8)
para(9,i,j)=temp(9)
para(10,i,j)=temp(10)
eta(i,j,k)=temp(11)
para(11,i,j)=temp(12)
lamda(i,j,k)=temp(13)

endfor
endfor

endif 

if k ne 0 then begin

for i=0,11 do begin ; if compling with gfortran
;for i=0,12 do begin ; if compling with ifort

readf,unit,header
endfor
 

for j=0,J1-1 do begin
for i=0,I1-1 do begin 
temp=fltarr(3)
readf,unit,temp
eta(i,j,k)=temp(0)
veff(i,j)=temp(1)
lamda(i,j,k)=temp(2)

endfor
endfor

endif

endfor
free_lun,unit

OPENW, 1,dirname+'/rcm-output-3D-IDL-'+fname+'.dat'

form3='(f13.6)'
PRINTF, 1,para,eta,lamda
close, 1 



end