I have data netcdf from ERA hourly, i choose 1940-2024 for all of days ech month and also because of value of data i consider just 4 hours 0:00, 6:00, 12:00, 18:00 for every grid point and my parameters are u10 & v10. my goal is to have in a map average time series that compute total days we have wind speed more than 15 m/s, and we khnow that wind speed = sqrt (u^2+v^2),
my data time is for example for jan is as 1940-01-01-0:00, 1940-01-01-6:00, ….., 2024-01-31-18:00. i need atfirst get average hourly to be daily data for every grid point every year and than compute for every grid point total days have wind speed more than 15 m/s , and after that for every grid point compute average 85 yerly to show on the map , please guide me
my code is
clear
clc
year = 1940:2024;
filename1 = ‘uv1.nc’;filename2 = ‘uv2.nc’;filename3 = ‘uv3.nc’;filename4 = ‘uv4.nc’;filename5 = ‘uv5.nc’;filename6 = ‘uv6.nc’;
filename7 = ‘uv7.nc’;filename8 = ‘uv8.nc’;filename9 = ‘uv9.nc’;filename10 = ‘uv10.nc’;filename11 = ‘uv11.nc’;filename12 = ‘uv12.nc’;
u1 = ncread(filename1, ‘u10’);u2 = ncread(filename2, ‘u10’);u3 = ncread(filename3, ‘u10’);u4 = ncread(filename4, ‘u10’);
u5 = ncread(filename5, ‘u10’);u6 = ncread(filename6, ‘u10’);u7 = ncread(filename7, ‘u10’);u8 = ncread(filename8, ‘u10’);
u9 = ncread(filename9, ‘u10’);u10 = ncread(filename10, ‘u10’);u11 = ncread(filename11, ‘u10’);u12 = ncread(filename12, ‘u10’);
v1 = ncread(filename1, ‘v10’);v2 = ncread(filename2, ‘v10’);v3 = ncread(filename3, ‘v10’);v4 = ncread(filename4, ‘v10’);
v5 = ncread(filename5, ‘v10’);v6 = ncread(filename6, ‘v10’);v7 = ncread(filename7, ‘v10’);v8 = ncread(filename8, ‘v10’);
v9 = ncread(filename9, ‘v10’);v10 = ncread(filename10, ‘v10’);v11 = ncread(filename11, ‘v10’); v12 = ncread(filename12, ‘v10’);
time1 = double(ncread(filename1, ‘valid_time’));time2 = double(ncread(filename2, ‘valid_time’));time3 = double(ncread(filename3, ‘valid_time’));
time4 = double(ncread(filename4, ‘valid_time’));time5 = double(ncread(filename5, ‘valid_time’));time6 = double(ncread(filename6, ‘valid_time’));
time7 = double(ncread(filename7, ‘valid_time’));time8 = double(ncread(filename8, ‘valid_time’));time9 = double(ncread(filename9, ‘valid_time’));
time10 = double(ncread(filename10, ‘valid_time’));time11 = double(ncread(filename11, ‘valid_time’));time12 = double(ncread(filename12, ‘valid_time’));
lat = double(ncread(filename1,’latitude’));lon = double(ncread(filename1,’longitude’));[Lat,Lon] = meshgrid(lat,lon);
% % compute land-sea mask and wind speed for every month Jan to dec
land = island(Lat,Lon);
ws1 = sqrt(u1.^2 + v1.^2);ws2 = sqrt(u2.^2 + v2.^2);ws3 = sqrt(u3.^2 + v3.^2);ws4 = sqrt(u4.^2 + v4.^2);
ws5 = sqrt(u5.^2 + v5.^2);ws6 = sqrt(u6.^2 + v6.^2);ws7 = sqrt(u7.^2 + v7.^2);ws8 = sqrt(u8.^2 + v8.^2);
ws9 = sqrt(u9.^2 + v9.^2);ws10 = sqrt(u10.^2 + v10.^2);ws11 = sqrt(u11.^2 + v11.^2);ws12 = sqrt(u12.^2 + v12.^2);
ws1(land==0)=NaN;ws2(land==0)=NaN;ws3(land==0)=NaN;ws4(land==0)=NaN;ws5(land==0)=NaN;ws6(land==0)=NaN;
ws7(land==0)=NaN;ws8(land==0)=NaN;ws9(land==0)=NaN;ws10(land==0)=NaN;ws11(land==0)=NaN;ws12(land==0)=NaN;I have data netcdf from ERA hourly, i choose 1940-2024 for all of days ech month and also because of value of data i consider just 4 hours 0:00, 6:00, 12:00, 18:00 for every grid point and my parameters are u10 & v10. my goal is to have in a map average time series that compute total days we have wind speed more than 15 m/s, and we khnow that wind speed = sqrt (u^2+v^2),
my data time is for example for jan is as 1940-01-01-0:00, 1940-01-01-6:00, ….., 2024-01-31-18:00. i need atfirst get average hourly to be daily data for every grid point every year and than compute for every grid point total days have wind speed more than 15 m/s , and after that for every grid point compute average 85 yerly to show on the map , please guide me
my code is
clear
clc
year = 1940:2024;
filename1 = ‘uv1.nc’;filename2 = ‘uv2.nc’;filename3 = ‘uv3.nc’;filename4 = ‘uv4.nc’;filename5 = ‘uv5.nc’;filename6 = ‘uv6.nc’;
filename7 = ‘uv7.nc’;filename8 = ‘uv8.nc’;filename9 = ‘uv9.nc’;filename10 = ‘uv10.nc’;filename11 = ‘uv11.nc’;filename12 = ‘uv12.nc’;
u1 = ncread(filename1, ‘u10’);u2 = ncread(filename2, ‘u10’);u3 = ncread(filename3, ‘u10’);u4 = ncread(filename4, ‘u10’);
u5 = ncread(filename5, ‘u10’);u6 = ncread(filename6, ‘u10’);u7 = ncread(filename7, ‘u10’);u8 = ncread(filename8, ‘u10’);
u9 = ncread(filename9, ‘u10’);u10 = ncread(filename10, ‘u10’);u11 = ncread(filename11, ‘u10’);u12 = ncread(filename12, ‘u10’);
v1 = ncread(filename1, ‘v10’);v2 = ncread(filename2, ‘v10’);v3 = ncread(filename3, ‘v10’);v4 = ncread(filename4, ‘v10’);
v5 = ncread(filename5, ‘v10’);v6 = ncread(filename6, ‘v10’);v7 = ncread(filename7, ‘v10’);v8 = ncread(filename8, ‘v10’);
v9 = ncread(filename9, ‘v10’);v10 = ncread(filename10, ‘v10’);v11 = ncread(filename11, ‘v10’); v12 = ncread(filename12, ‘v10’);
time1 = double(ncread(filename1, ‘valid_time’));time2 = double(ncread(filename2, ‘valid_time’));time3 = double(ncread(filename3, ‘valid_time’));
time4 = double(ncread(filename4, ‘valid_time’));time5 = double(ncread(filename5, ‘valid_time’));time6 = double(ncread(filename6, ‘valid_time’));
time7 = double(ncread(filename7, ‘valid_time’));time8 = double(ncread(filename8, ‘valid_time’));time9 = double(ncread(filename9, ‘valid_time’));
time10 = double(ncread(filename10, ‘valid_time’));time11 = double(ncread(filename11, ‘valid_time’));time12 = double(ncread(filename12, ‘valid_time’));
lat = double(ncread(filename1,’latitude’));lon = double(ncread(filename1,’longitude’));[Lat,Lon] = meshgrid(lat,lon);
% % compute land-sea mask and wind speed for every month Jan to dec
land = island(Lat,Lon);
ws1 = sqrt(u1.^2 + v1.^2);ws2 = sqrt(u2.^2 + v2.^2);ws3 = sqrt(u3.^2 + v3.^2);ws4 = sqrt(u4.^2 + v4.^2);
ws5 = sqrt(u5.^2 + v5.^2);ws6 = sqrt(u6.^2 + v6.^2);ws7 = sqrt(u7.^2 + v7.^2);ws8 = sqrt(u8.^2 + v8.^2);
ws9 = sqrt(u9.^2 + v9.^2);ws10 = sqrt(u10.^2 + v10.^2);ws11 = sqrt(u11.^2 + v11.^2);ws12 = sqrt(u12.^2 + v12.^2);
ws1(land==0)=NaN;ws2(land==0)=NaN;ws3(land==0)=NaN;ws4(land==0)=NaN;ws5(land==0)=NaN;ws6(land==0)=NaN;
ws7(land==0)=NaN;ws8(land==0)=NaN;ws9(land==0)=NaN;ws10(land==0)=NaN;ws11(land==0)=NaN;ws12(land==0)=NaN; I have data netcdf from ERA hourly, i choose 1940-2024 for all of days ech month and also because of value of data i consider just 4 hours 0:00, 6:00, 12:00, 18:00 for every grid point and my parameters are u10 & v10. my goal is to have in a map average time series that compute total days we have wind speed more than 15 m/s, and we khnow that wind speed = sqrt (u^2+v^2),
my data time is for example for jan is as 1940-01-01-0:00, 1940-01-01-6:00, ….., 2024-01-31-18:00. i need atfirst get average hourly to be daily data for every grid point every year and than compute for every grid point total days have wind speed more than 15 m/s , and after that for every grid point compute average 85 yerly to show on the map , please guide me
my code is
clear
clc
year = 1940:2024;
filename1 = ‘uv1.nc’;filename2 = ‘uv2.nc’;filename3 = ‘uv3.nc’;filename4 = ‘uv4.nc’;filename5 = ‘uv5.nc’;filename6 = ‘uv6.nc’;
filename7 = ‘uv7.nc’;filename8 = ‘uv8.nc’;filename9 = ‘uv9.nc’;filename10 = ‘uv10.nc’;filename11 = ‘uv11.nc’;filename12 = ‘uv12.nc’;
u1 = ncread(filename1, ‘u10’);u2 = ncread(filename2, ‘u10’);u3 = ncread(filename3, ‘u10’);u4 = ncread(filename4, ‘u10’);
u5 = ncread(filename5, ‘u10’);u6 = ncread(filename6, ‘u10’);u7 = ncread(filename7, ‘u10’);u8 = ncread(filename8, ‘u10’);
u9 = ncread(filename9, ‘u10’);u10 = ncread(filename10, ‘u10’);u11 = ncread(filename11, ‘u10’);u12 = ncread(filename12, ‘u10’);
v1 = ncread(filename1, ‘v10’);v2 = ncread(filename2, ‘v10’);v3 = ncread(filename3, ‘v10’);v4 = ncread(filename4, ‘v10’);
v5 = ncread(filename5, ‘v10’);v6 = ncread(filename6, ‘v10’);v7 = ncread(filename7, ‘v10’);v8 = ncread(filename8, ‘v10’);
v9 = ncread(filename9, ‘v10’);v10 = ncread(filename10, ‘v10’);v11 = ncread(filename11, ‘v10’); v12 = ncread(filename12, ‘v10’);
time1 = double(ncread(filename1, ‘valid_time’));time2 = double(ncread(filename2, ‘valid_time’));time3 = double(ncread(filename3, ‘valid_time’));
time4 = double(ncread(filename4, ‘valid_time’));time5 = double(ncread(filename5, ‘valid_time’));time6 = double(ncread(filename6, ‘valid_time’));
time7 = double(ncread(filename7, ‘valid_time’));time8 = double(ncread(filename8, ‘valid_time’));time9 = double(ncread(filename9, ‘valid_time’));
time10 = double(ncread(filename10, ‘valid_time’));time11 = double(ncread(filename11, ‘valid_time’));time12 = double(ncread(filename12, ‘valid_time’));
lat = double(ncread(filename1,’latitude’));lon = double(ncread(filename1,’longitude’));[Lat,Lon] = meshgrid(lat,lon);
% % compute land-sea mask and wind speed for every month Jan to dec
land = island(Lat,Lon);
ws1 = sqrt(u1.^2 + v1.^2);ws2 = sqrt(u2.^2 + v2.^2);ws3 = sqrt(u3.^2 + v3.^2);ws4 = sqrt(u4.^2 + v4.^2);
ws5 = sqrt(u5.^2 + v5.^2);ws6 = sqrt(u6.^2 + v6.^2);ws7 = sqrt(u7.^2 + v7.^2);ws8 = sqrt(u8.^2 + v8.^2);
ws9 = sqrt(u9.^2 + v9.^2);ws10 = sqrt(u10.^2 + v10.^2);ws11 = sqrt(u11.^2 + v11.^2);ws12 = sqrt(u12.^2 + v12.^2);
ws1(land==0)=NaN;ws2(land==0)=NaN;ws3(land==0)=NaN;ws4(land==0)=NaN;ws5(land==0)=NaN;ws6(land==0)=NaN;
ws7(land==0)=NaN;ws8(land==0)=NaN;ws9(land==0)=NaN;ws10(land==0)=NaN;ws11(land==0)=NaN;ws12(land==0)=NaN; netcdf, climate, stormy MATLAB Answers — New Questions

​