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-2D-IDL-'+fname+'.dat'

para=fltarr(30,I1,J1)

header=''
openr,unit,dirname+'rcm-output-2D-'+fname+'.dat',/get_lun


for i=0,9 do begin ; if rcm_tecplot.out is complied using gfortran
;for i=0,10 do begin ; if rcm_tecplot.out is complied using ifort

readf,unit,header
endfor


 for j=0,J1-1 do begin
for i=0,I1-1 do begin 
temp=fltarr(5,6)
readf,unit,temp
para(0,i,j)=temp(0,0)
para(1,i,j)=temp(1,0)
para(2,i,j)=temp(2,0)
para(3,i,j)=temp(3,0)
para(4,i,j)=temp(4,0)
para(5,i,j)=temp(0,1)
para(6,i,j)=temp(1,1)
para(7,i,j)=temp(2,1)
para(8,i,j)=temp(3,1)
para(9,i,j)=temp(4,1)
para(10,i,j)=temp(0,2)
para(11,i,j)=temp(1,2)
para(12,i,j)=temp(2,2)
para(13,i,j)=temp(3,2)
para(14,i,j)=temp(4,2)
para(15,i,j)=temp(0,3)
para(16,i,j)=temp(1,3)
para(17,i,j)=temp(2,3)
para(18,i,j)=temp(3,3)
para(19,i,j)=temp(4,3)
para(20,i,j)=temp(0,4)
para(21,i,j)=temp(1,4)
para(22,i,j)=temp(2,4)
para(23,i,j)=temp(3,4)
para(24,i,j)=temp(4,4)
para(25,i,j)=temp(0,5)
para(26,i,j)=temp(1,5)
para(27,i,j)=temp(2,5)
para(28,i,j)=temp(3,5)
para(29,i,j)=temp(4,5)

endfor
endfor

free_lun,unit

OPENW, 1,dirname+'/rcm-output-2D-IDL-'+fname+'.dat'
form3='(f13.6)'
PRINTF, 1,para
close, 1 

end