search

ja754969 / HYCOM_download

The codes for downloading data from https://www.hycom.org/
0 stars 0 forks source link

Downdload HYCOM data #2

Open ja754969 opened 1 year ago

ja754969 commented 1 year ago
ja754969 commented 1 year ago

1993

/datasets/GLBu0.08/expt_19.0/data/1993 https://tds.hycom.org/thredds/catalogs/GLBu0.08/expt_19.0.html

ja754969 commented 1 year ago

Download_HYCOMfromOPENDAP_ver3.m

clear,clc
latlim = [10 50];
lonlim = [110 270];
firstdate = datenum('2000010100','yyyymmddHH');
fn_i = ['https://tds.hycom.org/thredds/dodsC/GLBy0.08/expt_93.0/uv3z?lat[0:1:4250],lon[0:1:4499],time[0:1:0]'];  
first_date_download = datetime('2018/12/04 12:00:00','InputFormat','yyyy/MM/dd HH:mm:ss');
end_date_download = datetime('2021/05/22 09:00:00','InputFormat','yyyy/MM/dd HH:mm:ss');
%%
loop_count = 0;
while loop_count == 0
    try
    lon_i = ncread(fn_i,'lon');
    lat_i = ncread(fn_i,'lat');
    time_i = ncread(fn_i,'time');
    loop_count = loop_count+1;
    catch   %如果error,等待60秒並再執行一次
        disp('error')
        loop_count = loop_count - 1;
        pause(60)
    end
end
lonindex = find(lon_i>=lonlim(1)&lon_i<=lonlim(end))-1;
latindex = find(lat_i>=latlim(1)&lat_i<=latlim(end))-1;
%%
base_date_download = datetime(2000,01,01)+hours(time_i);

first_index = hours(first_date_download-base_date_download)/3;

end_index = hours(end_date_download-base_date_download)/3;
%%
time_index = first_index;
t = time_index; %起始時間index
while t <= end_index
    %https://tds.hycom.org/thredds/dodsC/GLBy0.08/expt_93.0/uv3z?lat[0:1:4250],lon[0:1:4499],time[0:1:10968],water_u[0:1:0][0:1:0][0:1:0][0:1:0],water_v[0:1:0][0:1:0][0:1:0][0:1:0]
    fn = ['https://tds.hycom.org/thredds/dodsC/GLBy0.08/expt_93.0/uv3z?lat[' ...
        num2str(latindex(1)) ':1:' num2str(latindex(end)) ...
        '],lon[' num2str(lonindex(1)) ':1:' num2str(lonindex(end)) ...
        '],time[' num2str(t) '],water_u[' num2str(t) ...
        '][0:1:0][' num2str(latindex(1)) ':1:' num2str(latindex(end)) ...
        '][' num2str(lonindex(1)) ':1:' num2str(lonindex(end)) ...
        '],water_v[' num2str(t) '][0:1:0][' num2str(latindex(1)) ':1:' ...
        num2str(latindex(end)) '][' num2str(lonindex(1)) ':1:' ...
        num2str(lonindex(end)) ']'];
    %     depth = ncread(fn,'depth');
    try
        lat = ncread(fn,'lat');
        lon = ncread(fn,'lon');
        u = ncread(fn,'water_u');
        v = ncread(fn,'water_v');
        time = ncread(fn,'time');
        date = datestr([firstdate + time/24],'yyyymmddHH')
        %   save(['F:\1_Tec_all\DATA\HYCOM\HYCOM_GLBy_surface_uv\' date],'lat','lon','u','v')
        saving_folder = './DATA/HYCOM/HYCOM_GLBy_surface_uv/';
        mkdir(saving_folder)
        save([saving_folder date],'lat','lon','u','v')
        %     clear depth lat lon u v time date
    catch   %如果error,等待60秒並再執行一次
        disp('error')
        t = t - 1;
        pause(60)

    end
    clear depth lat lon u v time date
    t = t + 1;
end