NDBC_hdf5 multiple stations

NDBC_hdf5/get_NDBC_air.m

@@ -1,4 +1,4 @@
-function [aa payload] = get_NDBC_air(fname)
+function [aa payload] = get_NDBC_air(fname,pdir)
 %
 %   read NODC netCDF4 file to get met information
 %   created 07/12 by TJ. Hesser
@@ -78,15 +78,26 @@
     re = 0;
 end
 if ~isnan(anemh)
-    dp = dir('C:\NDBC\alt_height.dat');
-    if exist('C:\NDBC\alt_height.dat','file')
-        fid2 = fopen(dp(1).name);
+    dp = [pdir,'alt_height.dat'];
+    if exist(dp,'file')
+        fid2 = fopen(dp);
         data = textscan(fid2,'%s%f');
-        ii = strcmp(data{1},fname(6:10));
+        ii = strcmp(data{1},fname(end-21:end-17));
         anemh = data{2}(ii);
+        fclose(fid2);
+        if isempty(anemh)
+            ques = ['What is the anemometer height for buoy ',fname(end-21:end-17),': '];
+            anemh = input(ques);
+            fid2 = fopen(dp,'a');
+            fprintf(fid2,'%5s %7.2f\n',fname(end-21:end-17),anemh);
+            fclose(fid2);
+        end
     else
         ques = ['What is the anemometer height for buoy ',fname(end-21:end-17),': '];
         anemh = input(ques);
+        fid2 = fopen(dp,'w');
+        fprintf(fid2,'%5s %7.2f\n',fname(end-21:end-17),anemh);
+        fclose(fid2);
     end
 end
 anemh = repmat(anemh,[length(year) 1]);

NDBC_hdf5/get_NDBC_air.m~

@@ -0,0 +1,174 @@
+function [aa payload] = get_NDBC_air(fname,pdir)
+%
+%   read NODC netCDF4 file to get met information
+%   created 07/12 by TJ. Hesser
+% 
+fprintf(1,'Analyizing Station %5s\n',fname(end-21:end-17));
+%Get attributes distributed
+fileatt = ncinfo(fname);
+fatt = {fileatt.Attributes.Name};
+
+% Get time out of the file and convert to year mon day hour min sec
+times = double(ncread(fname,'/time'));
+[year,mon,day,hour,min,sec] = datevec( ...
+        datenum(1970,1,1)+times./86400);
+
+% Get lat and lon out of file and convert to deg min sec  
+latm = strcmp(fatt,'nominal_latittude');
+lat = ncreadatt(fname,'/','nominal_latitude');
+[lad,lam,las] = ddeg2dms(abs(lat));
+lam = abs(lam);las = round(las);las = abs(las);
+lad = repmat(lad,[length(year) 1]);lam = repmat(lam,[length(year) 1]);
+las = repmat(las,[length(year) 1]);
+lon = ncreadatt(fname,'/','nominal_longitude');
+[lod,lom,los] = ddeg2dms(abs(lon));
+lom = abs(lom);los = abs(los);los = round(los);
+lod = repmat(lod,[length(year) 1]);lom = repmat(lom,[length(year) 1]);
+los = repmat(los,[length(year) 1]);
+if lon < 0
+    lod = -lod;
+end
+if lat < 0
+    lat = -lat;
+end
+
+
+% Get depth out of file
+try
+    dep = ncreadatt(fname,'/','sea_floor_depth_below_sea_level');
+catch
+    dep = -999.00;
+end
+dep = repmat(dep,[length(year) 1]);
+
+% Get windspeed, windgust, wind direction
+try
+    wspd = ncread(fname,['/payload_1/anemometer_1/wind_speed']);
+    payload = 'payload_1';
+    payatt = ncreadatt(fname,'/payload_1','description');
+    np = length(payatt);sp = ['%',num2str(np),'s\n'];
+    fprintf(1,'Using payload_1\n');
+    fprintf(1,sp,payatt);
+catch
+    wspd = ncread(fname,['/payload_2/anemometer_1/wind_speed']);
+    payload = 'payload_2';
+    payatt = ncreadatt(fname,'/payload_2','description');
+    np = length(payatt);sp = ['%',num2str(np),'s\n'];
+    fprintf(1,'Using payload_2\n');
+    fprintf(1,sp,payatt);
+end
+wspd_q = ncread(fname,['/',payload,'/anemometer_1/wind_speed_qc']);
+ii = wspd_q ~= 0;
+wspd(ii) = -99.99;
+if sum(ii) > 0
+  wspd2 = ncread(fname,['/',payload,'/anemometer_2/wind_speed']);
+  wspd_q2 = ncread(fname,['/',payload,'/anemometer_2/wind_speed_qc']);
+  jj = wspd_q2(ii) == 0;
+  wspd(ii(jj)) = wspd2(ii(jj));
+end
+wspd_info = ncinfo(fname,['/',payload,'/anemometer_1/wind_speed']);
+ii = wspd == wspd_info.FillValue;
+wspd(ii) = -99.99;
+
+% Get anemometer height
+try
+    anemh = ncreadatt(fname,['/',payload,'/anemometer_1','height_of_instrument']);
+catch
+    anemh = 'NaN';
+    re = 0;
+end
+if ~isnan(anemh)
+    if exist([pdir,alt_height.dat'],'file')
+        fid2 = fopen(dp(1).name);
+        data = textscan(fid2,'%s%f');
+        ii = strcmp(data{1},fname(6:10));
+        anemh = data{2}(ii);
+        if isempty(anemh
+    else
+        ques = ['What is the anemometer height for buoy ',fname(end-21:end-17),': '];
+        anemh = input(ques);
+    end
+end
+anemh = repmat(anemh,[length(year) 1]);
+
+wgus = ncread(fname,['/',payload,'/anemometer_1/wind_gust']);
+wgus_q = ncread(fname,['/',payload,'/anemometer_1/wind_gust_qc']);
+ii = wgus_q ~= 0;
+wgus(ii) = -99.99;
+if sum(ii) > 0
+  wgus2 = ncread(fname,['/',payload,'/anemometer_2/wind_gust']);
+  wgus_q2 = ncread(fname,['/',payload,'/anemometer_2/wind_gust_qc']);
+  jj = wgus_q2(ii) == 0;
+  wgus(ii(jj)) = wgus2(ii(jj));
+end
+wgus_info = ncinfo(fname,['/',payload,'/anemometer_1/wind_gust']);
+ii = wgus == wgus_info.FillValue;
+wgus(ii) = -99.99;
+
+wdir = ncread(fname,['/',payload,'/anemometer_1/wind_direction']);
+wdir_q = ncread(fname,['/',payload,'/anemometer_1/wind_direction_qc']);
+ii = wdir_q ~= 0;
+wdir(ii) = -999.0;
+if sum(ii) > 0
+  wdir2 = ncread(fname,['/',payload,'/anemometer_2/wind_direction']);
+  wdir_q2 = ncread(fname,['/',payload,'/anemometer_2/wind_direction_qc']);
+  jj = wdir_q2(ii) == 0;
+  wdir(ii(jj)) = wdir2(ii(jj));
+end
+wdir_info = ncinfo(fname,['/',payload,'/anemometer_1/wind_direction']);
+ii = wdir == wdir_info.FillValue;
+wdir(ii) = -999.0;
+wdir = single(wdir);
+
+% Get Air temperature and Sea tempurature
+atemp = ncread(fname,['/',payload,'/air_temperature_sensor_1/air_temperature']);
+atemp_q = ncread(fname,['/',payload,'/air_temperature_sensor_1/air_temperature_qc']);
+atemp_info = ncinfo(fname,['/',payload,'/air_temperature_sensor_1/air_temperature']);
+ii = atemp == atemp_info.FillValue;
+atemp = atemp - 272.15;
+atemp(ii) = -99.99;
+ii = atemp_q ~= 0;
+atemp(ii) = -99.99;
+try
+    if sum(ii) > 0
+    atemp2 = ncread(fname,['/',payload,'/air_temperature_sensor_2/air_temperature']);
+    atemp_q2 = ncread(fname,['/',payload,'/air_temperature_sensor_2/air_temperature']);
+    jj = atemp_q2(ii) == 0;
+    atemp(ii(jj)) = atemp2(ii(jj));
+    end
+catch
+end
+
+try
+    stemp = ncread(fname,['/',payload,'/ocean_temperature_sensor_1/sea_surface_temperature']);
+    stemp_q = ncread(fname,['/',payload,'/ocean_temperature_sensor_1/sea_surface_temperature_qc']);
+    stemp_info = ncinfo(fname,['/',payload,'/ocean_temperature_sensor_1/sea_surface_temperature']);
+    ii = stemp == stemp_info.FillValue;
+    stemp = stemp - 272.15;
+    stemp(ii) = -99.99;
+    ii = stemp_q ~= 0;
+    stemp(ii) = -99.99;
+catch
+    stemp = repmat(-99.99,[length(times) 1]);
+end
+% Get barometric pressure from file
+barp = ncread(fname,['/',payload,'/barometer_1/air_pressure_at_sea_level']);
+barp_q = ncread(fname,['/',payload,'/barometer_1/air_pressure_at_sea_level_qc']);
+barp_info = ncinfo(fname,['/',payload,'/barometer_1/air_pressure_at_sea_level']);
+ii = barp == barp_info.FillValue;
+barp = double(barp)./100;
+barp(ii) = -999.;
+ii = barp_q ~= 0;
+barp(ii) = -999.;
+if sum(ii) ~= 0
+    barp2 = ncread(fname,['/',payload,'/barometer_2/air_pressure_at_sea_level']);
+    barp_q2 = ncread(fname,['/',payload,'/barometer_2/air_pressure_at_sea_level_qc']);
+    jj = barp_q2(ii) == 0;
+    barp(ii(jj)) = barp2(ii(jj));
+end
+
+aa = [year,mon,day,hour,min,lad,lam,las,lod,lom,los,dep,anemh, ...
+    double(wspd),double(wgus),double(wdir),double(atemp),double(stemp), ...
+    double(barp)];
+
+

NDBC_hdf5/get_nodc_buoy.m

@@ -19,7 +19,11 @@
     mkdir(dir2);
 end
 cd(dir2);
-filefull = dir('*.nc');
+if ~strcmp(stat,'00000')
+    filefull = dir(['*',stat,'*.nc']);
+else
+    filefull = dir('*.nc');
+end
 if ~isempty(filefull)
     perr = 0;
     return
@@ -43,7 +47,7 @@
         num2str(year),'_',monc,'.spe*'],'.');
     return
 end
-if ~exist('stat','var')
+if strcmp(stat,'00000')
     mget(f,'*.nc');
 else
     mget(f,['*',stat,'*.nc']);

NDBC_hdf5/process_NODC_netcdf.m

@@ -1,6 +1,16 @@
-function process_NODC_netcdf(year,mon,statname,pname)
+function process_NODC_netcdf(year,mon,varargin)
 % retrieve and process NODC netCDF4 files
 % created 07/12 TJ. Hesser
+p = inputParser;
+p.addRequired('year');
+p.addRequired('mon');
+p.addOptional('statname','00000');
+p.addOptional('pname',pwd);
+parse(p,year,mon,varargin{:});
+
+statname = p.Results.statname;
+pname = p.Results.pname;
+
 if ~ischar(year)
     yearc = num2str(year);
 else
@@ -23,31 +33,40 @@ function process_NODC_netcdf(year,mon,statname,pname)
 else
     statc = statname;
 end
+
 get_nodc_buoy(year,mon,statc,pname)
 %pname = 'http://data.nodc.noaa.gov/thredds/dodsC/ndbc/cmanwx/';
 %pname = '';
-try
-    fname = [pname,'/',yearc,'/',mond,'/NDBC_',statc,'_',yearc, ...
-        monc,'_D1_v00.nc'];
-    binfo = ncinfo(fname);
-catch
-    fname = [pname,'/',yearc,'/',mond,'/NDBC_',statc,'_',yearc, ...
-        monc,'_D2_v00.nc'];
+if strcmp(statc,'00000')
+    fname1 = dir('*.nc');
+    np = size(fname1,1);
+else
     try
-        binfo = ncinfo(fname);
+        fname1.name = [pname,'/',yearc,'/',mond,'/NDBC_',statc,'_',yearc, ...
+            monc,'_D1_v00.nc'];
+        binfo = ncinfo(fname1.name);
     catch
-        fprintf(1,'Need specific name\n');
+        fname1.name= [pname,'/',yearc,'/',mond,'/NDBC_',statc,'_',yearc, ...
+            monc,'_D2_v00.nc'];
+        try
+            binfo = ncinfo(fname1.name);
+        catch
+            fprintf(1,'Need specific name\n');
+        end
     end
+    np = 1;
 end
 
-onlform = ['%5s%6i%3i%3i%3i%3i%5i%5i%5i%5i%5i%5i%7.1f%6.2f%6.2f%6.2f', ...
-    '%6.1f%7.2f%7.2f%7.1f%8.2f%8.2f%8.2f%8.2f%7.1f%7.1f%7.1f%7.1f', ...
+onlform = ['%5s%6i%3i%3i%3i%3i%5i%5i%5i%5i%5i%5i%7.1f%7.2f%7.2f%7.2f', ...
+    '%7.1f%7.2f%7.2f%7.1f%8.2f%8.2f%8.2f%8.2f%7.1f%7.1f%7.1f%7.1f', ...
     '%8.2f%8.2f%8.2f%8.2f%8.2f%10.4f%7.1f%7.1f%8.2f% 14.5E% 14.5E', ...
     '% 14.5E% 14.5E\n'];
 sp1form1 = ['%5s%6i%3i%3i%3i%3i'];
-%for zz = 1:size(files,1)
+for zz = 1:np
+    fname = fname1(zz).name;
+    statc = fname(end-21:end-17);
     nn = 1;
-    [aa payload] = get_NDBC_air(fname);
+    [aa payload] = get_NDBC_air(fname,pname);
 %     if zz > 1 & files(zz).name(6:10) == stat
 %         nn = 2;
 %     end
@@ -108,4 +127,4 @@ function process_NODC_netcdf(year,mon,statname,pname)
         fclose(fid);
         fclose(fid2);
     end
-%end    
+end