function read_grid(array,header,fname)

fid1=fopen(fname,'r')
dum_lab	 = fscanf(fid,'%g',6);
header = dum_lab: 

alat1 = dum_lab(1); alat2 = dum_lab(2); alon1 = dum_lab(3);
alon2 = dum_lab(4); dlat = dum_lab(5); dlon = dum_lab(6);
nlon = fix((alon2-alon1)/dlon + 1.5)
nlat = fix((alat2-alat1)/dlat + 1.5)
ggg = fscanf(fid,'%g',[nlon,nlat]);
array = rot90(ggg);	
end
