fixing issues with NDBC hdf5

NDBC_hdf5/get_NDBC_air.m

@@ -45,18 +45,25 @@
 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');
+    try
+        wspd = ncread(fname,['/payload_2/anemometer_1/wind_speed']);
+        payload = 'payload_2';
+    catch
+        payload = 'none';aa = 0;
+        return
+    end
+
+end
+try
+    payatt = ncreadatt(fname,['/',payload],'description');
     np = length(payatt);sp = ['%',num2str(np),'s\n'];
-    fprintf(1,'Using payload_2\n');
+    fprintf(1,'Using %9s \n',payload);
     fprintf(1,sp,payatt);
+catch
+    fprintf(1,'No Description of %9s \n',payload);
 end
+
 wspd_q = ncread(fname,['/',payload,'/anemometer_1/wind_speed_qc']);
 ii = wspd_q ~= 0;
 wspd(ii) = -99.99;

NDBC_hdf5/get_NDBC_spec.m

@@ -43,8 +43,8 @@
 catch
     fprintf(1,'Wave Spectra information not available for station %5s\n', ...
         fname(end-21:end-17));
-    times = double(h5read(fname,'/time'));
-    c11 = repmat(-999.00,[length(freq) length(times)]);
+    time_spec = double(h5read(fname,'/time'));
+    c11 = repmat(-999.00,[length(freq) length(time_spec)]);
     c11_info.FillValue = 0.0;
 end
 
@@ -87,6 +87,12 @@
         tm = 1./(fm + 1.0e-10);
         tp = 1./(fp + 1.0e-10);
     end
+
+    r1(r1(:,zz) == c11_info.FillValue,zz) = 0;
+    r2(r2(:,zz) == c11_info.FillValue,zz) = 0;
+    alpha1(alpha1(:,zz) == c11_info.FillValue,zz) = 0;
+    alpha2(alpha2(:,zz) == c11_info.FillValue,zz) = 0;
+
 %end
 
     if mp > 0

NDBC_hdf5/process_NODC_netcdf.m

@@ -41,20 +41,21 @@ function process_NODC_netcdf(year,mon,varargin)
     fname1 = dir('*.nc');
     np = size(fname1,1);
 else
-    try
-        fname1.name = [pname,'/',yearc,'/',mond,'/NDBC_',statc,'_',yearc, ...
-            monc,'_D1_v00.nc'];
-        binfo = ncinfo(fname1.name);
-    catch
-        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;
+%     try
+%         fname1.name = [pname,'/',yearc,'/',mond,'/NDBC_',statc,'_',yearc, ...
+%             monc,'_D1_v00.nc'];
+%         binfo = ncinfo(fname1.name);
+%     catch
+%         fname1.name= [pname,'/',yearc,'/',mond,'/NDBC_',statc,'_',yearc, ...
+%             monc,'_D3_v00.nc'];
+%         try
+%             binfo = ncinfo(fname1.name);
+%         catch
+%             fprintf(1,'Need specific name\n');
+%         end
+%     end
+    fname1 = dir(['*',statc,'*.nc']);
+    np = size(fname1,1);
 end
 
 onlform = ['%5s%6i%3i%3i%3i%3i%5i%5i%5i%5i%5i%5i%7.1f%7.2f%7.2f%7.2f', ...
@@ -65,8 +66,16 @@ function process_NODC_netcdf(year,mon,varargin)
 for zz = 1:np
     fname = fname1(zz).name;
     statc = fname(end-21:end-17);
-    nn = 1;
+    if zz > 1 & strcmp(statc,statcold)
+        nn = 2;
+    else
+        nn = 1;
+    end
     [aa payload] = get_NDBC_air(fname,pname);
+    if strcmp(payload,'none')
+        statcold = ' ';
+        continue
+    end
 %     if zz > 1 & files(zz).name(6:10) == stat
 %         nn = 2;
 %     end
@@ -127,4 +136,10 @@ function process_NODC_netcdf(year,mon,varargin)
         fclose(fid);
         fclose(fid2);
     end
-end    
+    statcold = statc;
+end    
+delete('*.nc');
+cd('ADCP');
+delete('*.nc');
+cd ../
+rmdir('ADCP');