dir = '/media/usbdisk/OCCAM/';
lon_pp = hdf5read(strcat(dir,'d6140to6145.h5m2'),'HDF4_DIMGROUP/LONGITUDE_U');
lat_pp = hdf5read(strcat(dir,'d6140to6145.h5m2'),'HDF4_DIMGROUP/LATITUDE_U');
[X Y] = meshgrid(lon_pp,lat_pp);
[lonA latA] = rotate_occam(X, Y);
lon = lonA;
%lon(lon>180) = lon(lon>180)-360;
lat = latA;
save('output_grid.mat', 'lon','lat');