% Make a regional temp-salt climatology for a given climatological
% date (set by clim_day) and region (set by ax)

% I have tried to make this generic by setting the region of interest
% at the start, so that in principle the code could be re-run without
% change for other lon/lat region of interest ... but no
% promises on this, a lot of stuff is still hardwired



% SOME CODE IS DISABLED BY PLACING WITHIN IF 0 -- END BLOCKS TO STRIP 
% THIS DOWN FOR TESTING OF PARALLEL MATLAB BY PARRY HUSBANDS



% -----------------------------------------------------------------
% define a region of interest
% limits of ROMS NJB grid
ax = [-75.4070 -71.6290 37.1280 41.0283]; 
subregion = 'NJB';

% Get all the NODC station data in the region of interest
if 0
  leo_clim_get_nodc_data % saves to 'MJJA_NODC_stdlev' <<< hardwired !!!!
end

% define a climatological day-of-year for the OI
clim_day = [NaN 7 1 0 0 0];

% define a target grid
if 0
  grd_dir  = '/home/e3/wilkin/wilkin/leo/';
  grd_file = 'njb1_grid_d.nc';
  g = [grd_dir grd_file];
end
if 0
  % model grid
  grd = roms_get_grid(g);
  Xg = grd.lon_rho;
  Yg = grd.lat_rho;
  disp([ ' Target coordinates from: ' grd_file])
else
  % simple cartesian grid for testing
  switch subregion
    case 'NJB'
      dlonlat = 0.1;
      [X,Y] = meshgrid(-75.5:0.1:-71.5,37.0:0.1:41.5);
    otherwise
      error(' Target grid undefined')
  end
  disp([ ' Target grid is uniform lon/lat at intervals of ' num2str(dlonlat)])
end

% OI NODC data to uniform lon/lat grid and NODC standard depths
nodc_clim_oi