Custom Query (986 matches)

Corrected several adjoint files that missed including module:

       USE mod_ncparam

Many thanks to Alan Berry for ​bringing this to my attention. I also corrected coupled of typos.

Several Matlab scripts used to process 4D-Var observations were updated. The observations URL are now passed as an argument to the processing functions, so the extraction scripts are more generic. Since mostly of the 4D-Var observations are on OpenDAP servers, the scripts now support the native NetCDF interface for Matlab version 2012a or higher and SNCTOOLS with Java support for older Matlab versions.

Users can download SNCTOOLS from ​SourceForge. If doing so, notice that the startup.m script has now:

my_home = getenv('HOME');

...

% Load NetCDF Toolbox for OpenDAP support for versions 2008b or higher. 
% However, this is not needed if version 2012a or higher since Matlab
% native NetCDF interface supports OpenDAP.  Users need to change the
% paths for SNCTOOLS and JAVA.

v = version('-release');
vyear = str2num(v(1:4));
load_toolbox = vyear >= 2008;
if ((vyear == 2008 && v(5:5) == 'a') || vyear >= 2012),
  load_toolbox = false;
end

if (load_toolbox),
  addpath (strcat(my_home, '/ocean/matlab/snctools'), '-end');
  javaaddpath (strcat(my_home, '/ocean/matlab/classes/toolsUI-4.1.jar'), '-end');
  javaaddpath (strcat(my_home, '/ocean/matlab/classes/netcdfAll-4.2.jar'), '-end');
  javaaddpath (strcat(my_home, '/ocean/matlab/snctools/classes'), '-end');
  setpref('SNCTOOLS','USE_JAVA', true);
end

Corrected a bug in Matlab script interp_boundary.m. In line 201 we need to have:

        if (~isempty(ind)),
          if (RemoveNaN),
            B.(field)(ind) = FN(I.XR.(edge)(ind), I.YR.(edge)(ind));  % <===
            Rmin = min(Rmin, min(B.(field)(ind)));
            Rmax = max(Rmax, max(B.(field)(ind)));

            ind = find(isnan(B.(field)));
            if (~isempty(ind)),
              Ncount = length(ind);
            end       
          else
            Ncount = length(ind);
          end
        end

instead of

            B.(field)(ind) = FN(I.XR.(edge)(ind), I.YR.(field)(ind));

We need to have edge instead of field in I.YR.(?) structure. This bug is only relevant when removing interpolated variable NaNs values with a nearest neighbor interpolant. This is the second interpolation to remove outliers.

This is an important correction and users need to update this file if using this function to prepare boundary conditions.