% roms_master_climatology_coawst_mw % % This routine : % - creates climatology, boundary, and initial condition files for ROMS: % coawst_clm.nc ; coawst_bdy.nc ; coawst_ini.nc % on a user-defined grid for a user-defined date. % % This is currently set up to use opendap calls to acquire data % from HYCOM + NCODA Global 1/12 Degree Analysis and interp to roms grid. % % based on efforts by: % written by Mingkui Li, May 2008 % Modified by Brandy Armstrong March 2009 % jcwarner April 20, 2009 % Ilgar Safak modified on June 27, 2012 such that now: % - HYCOM url is a user-definition % - "hc" is called from the structure "gn".(still needs to be tested with wet/dry). % - updatinit_coawst_mw.m modified to get desired time (T1) as a variable; % ocean_time=T1-datenum(1858,11,17,0,0,0) % Updates from Christie Hegermiller, Feb 2019 % %%%%%%%%%%%%%%%%%%%%% START OF USER INPUT %%%%%%%%%%%%%%%%%%%%%%%%%% % (1) ���뿪ʼ����(T1)�������Ի�ȡ�������� %T1 = datenum ( 2012 , 10 , 28 , 12 , 0 , 0 ); %��ʼ���� %T1 = datenum(2022,06,20,0,0,0); %��ʼ���� T1 = datenum(2018,09,13,0,0,0); %��ʼ���� %���������ļ���������Ƶ�� numdays = 5; dayFrequency = 1; % (2) Ϊ�����ʱ��T1����HYCOMĿ¼����ַ % see http://tds.hycom.org/thredds/catalog.html url = 'https://tds.hycom.org/thredds/dodsC/GLBa0.08/expt_91.2/2018'; % 2016-04 to 2018-11 %url = 'http://tds.hycom.org/thredds/dodsC/GLBv0.08/expt_93.0'; %url = 'https://tds.hycom.org/thredds/dodsC/GLBy0.08/expt_93.0' % (3) ���빤��Ŀ¼��wdr�� wdr = '/media/five/WinData/COAWST/matlab/Roms/Mangkhut_2'; % (4) ����ROMS����ģ�͵�·�������� modelgrid = 'Mangkhut_roms_grid2.nc' % (5) ��������ֱ�������-��Щ������Ҫ��ROMS����һ�¡� theta_s = 5.0; theta_b = 0.4; Tcline = 50.0; N = 16; Vtransform = 2; Vstretching = 4; %Vtransform = 1; %Vstretching = 1; %%%%%%%%%%%%%%%%%%%%% END OF USER INPUT %%%%%%%%%%%%%%%%%%%%%%%%%% eval(['cd ',wdr]) tic % Call to get HYCOM indices for the defined ROMS grid disp('getting roms grid, hycom grid, and overlapping indices') [gn, clm]=get_ijrg(url, modelgrid, theta_s, theta_b, Tcline, N, Vtransform, Vstretching); % Call to create the climatology (clm) file disp('going to create clm file') fn=updatclim_coawst_mw(T1, gn, clm, 'coawst_clm.nc', wdr, url) % Call to create the boundary (bdy) file disp('going to create bndry file') updatbdry_coawst_mw(fn, gn, 'coawst_bdy.nc', wdr) % Call to create the initial (ini) file disp('going to create init file') updatinit_coawst_mw(fn, gn, 'coawst_ini.nc', wdr, T1) toc %% Call to create the long climatology (clm) file if numdays>1 disp('going to create more days of clm and bnd files') if (ispc) eval(['!copy coawst_clm.nc coawst_clm_',datestr(T1,'yyyymmdd'),'.nc']) eval(['!copy coawst_bdy.nc coawst_bdy_',datestr(T1,'yyyymmdd'),'.nc']) else eval(['!cp coawst_clm.nc coawst_clm_',datestr(T1,'yyyymmdd'),'.nc']) eval(['!cp coawst_bdy.nc coawst_bdy_',datestr(T1,'yyyymmdd'),'.nc']) end for it=dayFrequency:dayFrequency:numdays-1 %1st day already created, NEED to set number of days at top! error =1; while error == 1; %da = 1; try fname=['coawst_clm_',datestr(T1+it,'yyyymmdd'),'.nc'] fn=updatclim_coawst_mw(T1+it,gn,clm,fname,wdr,url) fname=['coawst_bdy_',datestr(T1+it,'yyyymmdd'),'.nc']; updatbdry_coawst_mw(fn,gn,fname,wdr) error = 0; %da = da + 1; catch disp(['Unable to download HYCOM data at ' datestr(T1+it,'yyyymmdd')]); end end end %% get an organized list of dated files Dclm=dirsort('coawst_clm_*.nc'); Dbdy=dirsort('coawst_bdy_*.nc'); %names for merged climatology/boundary files fout='merged_coawst_clm.nc'; foutb='merged_coawst_bdy.nc'; %create netcdf files to merge climatology into create_roms_netcdf_clm_mwUL(fout,gn,length(Dclm));% converted to BI functions create_roms_netcdf_bndry_mwUL(foutb,gn,length(Dbdy));% converted to BI functions %% fill merged climatology files with data from each clm file % each file must contain only ONE time step %get variable names vinfo=ncinfo(fout); for nf=1:length(Dclm) fin=Dclm(nf).name; for nv=1:length({vinfo.Variables.Name}) if length({vinfo.Variables(nv).Dimensions.Name})==4; eval(['ncwrite(fout,''',vinfo.Variables(nv).Name,''',ncread(fin,''',vinfo.Variables(nv).Name,'''),[1 1 1 nf]);']); elseif length({vinfo.Variables(nv).Dimensions.Name})==3; eval(['ncwrite(fout,''',vinfo.Variables(nv).Name,''',ncread(fin,''',vinfo.Variables(nv).Name,'''),[1 1 nf]);']); elseif length({vinfo.Variables(nv).Dimensions.Name})==2; try eval(['ncwrite(fout,''',vinfo.Variables(nv).Name,''',ncread(fin,''',vinfo.Variables(nv).Name,'''),[1 nf]);']); catch display([vinfo.Variables(nv).Name ' is a dimension and has already been written to the file.']) end elseif length({vinfo.Variables(nv).Dimensions.Name})==1; try eval(['ncwrite(fout,''',vinfo.Variables(nv).Name,''',ncread(fin,''',vinfo.Variables(nv).Name,'''),[nf]);']); catch display([vinfo.Variables(nv).Name ' is a dimension and has already been written to the file.']) end end end end vinfo=ncinfo(foutb); for nf=1:length(Dbdy) for nv=1:length({vinfo.Variables.Name}) fin=Dbdy(nf).name; if length({vinfo.Variables(nv).Dimensions.Name})==4; eval(['ncwrite(foutb,''',vinfo.Variables(nv).Name,''',ncread(fin,''',vinfo.Variables(nv).Name,'''),[1 1 1 nf]);']); elseif length({vinfo.Variables(nv).Dimensions.Name})==3; eval(['ncwrite(foutb,''',vinfo.Variables(nv).Name,''',ncread(fin,''',vinfo.Variables(nv).Name,'''),[1 1 nf]);']); elseif length({vinfo.Variables(nv).Dimensions.Name})==2; try eval(['ncwrite(foutb,''',vinfo.Variables(nv).Name,''',ncread(fin,''',vinfo.Variables(nv).Name,'''),[1 nf]);']); catch display([vinfo.Variables(nv).Name ' is a dimension and has already been written to the file.']) end elseif length({vinfo.Variables(nv).Dimensions.Name})==1; try eval(['ncwrite(foutb,''',vinfo.Variables(nv).Name,''',ncread(fin,''',vinfo.Variables(nv).Name,'''),[nf]);']); catch display([vinfo.Variables(nv).Name ' is a dimension and has already been written to the file.']) end end end end end toc