clear all
close all
clc

ny1 = 721;nx1 = 1440;
fil = ['ERA5_NRT_u10m_2023'];
d = dir(fil);
nt1 = (d.bytes/nx1/ny1/4);
lon =  0:0.25:359.75;
lat =  -90:0.25:90;

for  cc = 1:nt1
    Q = rdslice(fil,[nx1 ny1],cc,'real*4','b');
    figure;imagesc(lon,lat,Q');colorbar;axis xy;
    clear Q
    pause
end
