ROMS提取Forcing文件求助！！！

空空空空明

我在获取ROMS所需要的运行文件中的forcing文件（可能是叫应力文件？）的时候从NARR网站提取数据却总是提取不全，请问各位大佬这样要怎么办啊。我用的是COAWST模式中的ncei_2roms.m文件，正常应该是运行得到 romsforc_NARR.nc文件，但是实际我在matlab中查看nc文件中的数值时发现总是会有缺失，具体代码在下面。

提取NAM数据的网址时使用的这个'https://www.ncei.noaa.gov/thredds/dodsC/model-namanl-old/'或者如果用NARR和GFS能提取到也可以)

1. % (1) Select which variables to include in this netcdf forcing file.
   
2. % put a '1' if you want to include it, '0' otherwise.
   
3. % 1/0              Var description (Units)
   
4. get_lwrad = 1;    % gets downward and upward longwave and computes net flux of down-up (W/m2)
   
5.                   % this will also store lwrad_down so you can use LONGWAVE option.
   
6. get_swrad = 1;    % gets downward and updward shortwave and computes net down-up flux (W/m2)
   
7. get_rain  = 1;    % precipitation rate (kg/m2/s) at surface
   
8. get_Tair  = 1;    % surface air temperature (C) at 2m
   
9. get_Pair  = 1;    % pressure reduce to MSL (Pa)
   
10. get_Qair  = 1;    % relative_humidity (percent) at 2m
    
11. get_Wind  = 1;    % surface u- and v- winds (m/s) at 10m
    
12. 
    
13. % (2) Enter name of output ROMS forcing file
    
14. %ROMS_force_name = 'romsforc_GFS_Sandy2012.nc';
    
15. %ROMS_force_name = 'romsforc_GFS_Mangkhut2018.nc';
    
16. ROMS_force_name = 'romsforc_NAM_Mangkhut_0310_2.nc';
    
17. 
    
18. % (3) Enter start and end dates
    
19. time_start = datenum('15-Sep-2018');
    
20. time_end   = datenum('18-Sep-2018');
    
21. 
    
22. % (4) Select which data to obtain: NAM, NARR, both NAM+NARR -- or -- GFS.
    
23. get_NARR = 0;  % NARR-A grid 221 32km data, available 1979-2014
    
24. get_NAM  = 1;  % NAM grid 218 12km data
    
25. % --- or ---
    
26. get_GFS  = 0;  % GFS 0.5 degree
    
27. % GFS is 6 hr and NAM/NARR is 3 hr. I dont have time interpolation
    
28. % added in so you have to pick NAM/NARR or GFS.
    
29. 
    
30. % (5) Select to interpolate to a roms grid or a user defined grid.
    
31. % Set one of these to a 1, the other to a 0.
    
32. interpto_roms_grid = 1;
    
33. interpto_user_grid = 0;
    
34. 
    
35. % The code has been written so that longitude is in -deg E.
    
36. if interpto_roms_grid
    
37.   model_grid = 'F:\COAWST\Roms\Mangkhut\Mangkhut_roms_grid.nc';
    
38. elseif interpto_user_grid
    
39. % Provide lon_rho, lat_rho, and angle_rho.
    
40. % NAM / NARR grids are centered at~ -100 deg lon; GFS = 0:360 lon
    
41.   if (get_NARR); offset=-360; end
    
42.   if (get_NAM);  offset=-360; end
    
43.   if (get_GFS);  offset=0;    end
    
44. % You Probably want to make a finer resolution from 0.2 to 0.1.
    
45.   lon_rho = [100:0.2:130]+offset;
    
46.   lat_rho = [ 10:0.2:50 ];
    
47.   lon_rho = repmat(lon_rho,size(lat_rho,2),1)';
    
48.   lat_rho = repmat(lat_rho',1,size(lon_rho,1))';
    
49.   angle_rho = lon_rho*0;
    
50. else
    
51.   disp('pick a grid')
    
52. end
    
53. 
    
54. % (6) Select which method to use: ncread or nctoolbox
    
55. use_matlab = 0;
    
56. use_nctoolbox = 1;

复制代码