function oct_plot_nestdust(child_dust,thefield,thetime,skip)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
% Test the embedded dust forcing file.
%
%  Further Information:  
%  http://www.croco-ocean.org
%  
%  This file is part of CROCOTOOLS
%
%  CROCOTOOLS is free software; you can redistribute it and/or modify
%  it under the terms of the GNU General Public License as published
%  by the Free Software Foundation; either version 2 of the License,
%  or (at your option) any later version.
%
%  CROCOTOOLS is distributed in the hope that it will be useful, but
%  WITHOUT ANY WARRANTY; without even the implied warranty of
%  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
%  GNU General Public License for more details.
%
%  You should have received a copy of the GNU General Public License
%  along with this program; if not, write to the Free Software
%  Foundation, Inc., 59 Temple Place, Suite 330, Boston,
%  MA  02111-1307  USA
%
%  Copyright (c) 2004-2006 by Pierrick Penven 
%  e-mail:Pierrick.Penven@ird.fr  
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
npts=[0 0 0 0];
i=0;
for time=thetime
  i=i+1;
  
  subplot(2,length(thetime),i)


  ncid = netcdf.open(child_dust, 'NC_NOWRITE');
  parent_dust=ncid.parent_file(:);
  child_grd=ncid.grd_file(:);
  fieldc=netcdf.getVar(ncid, netcdf.inqVarID(ncid, thefield));
  fieldname=netcdf.getAtt(ncid, netcdf.inqVarID(ncid, thefield), 'long_name');
  result=close(ncid);

  ncid = netcdf.open(child_grd, 'NC_NOWRITE');
  parent_grd=ncid.parent_grid(:);
  refinecoeff=netcdf.getVar(ncid, netcdf.inqVarID(ncid, 'refine_coef'));
  lonc=netcdf.getVar(ncid, netcdf.inqVarID(ncid, 'lon_rho'));
  latc=netcdf.getVar(ncid, netcdf.inqVarID(ncid, 'lat_rho'));
  mask=netcdf.getVar(ncid, netcdf.inqVarID(ncid, 'mask_rho'));
  result=close(ncid);
  mask(mask==0)=NaN;
  pcolor(lonc,latc,mask.*fieldc)
  shading flat
  axis image
  caxis([min(min(fieldc)) max(max(fieldc))])
  colorbar
  axis([min(min(lonc)) max(max(lonc)) min(min(latc)) max(max(latc))])
  title(['\bf ',fieldname,' CHILD'])

  subplot(2,length(thetime),i+length(thetime))
  ncid = netcdf.open(parent_dust, 'NC_NOWRITE');
  field=netcdf.getVar(ncid, netcdf.inqVarID(ncid, thefield));
  fieldname=netcdf.getAtt(ncid, netcdf.inqVarID(ncid, thefield), 'long_name');
  result=close(ncid);
  
  ncid = netcdf.open(parent_grd, 'NC_NOWRITE');
  lon=netcdf.getVar(ncid, netcdf.inqVarID(ncid, 'lon_rho'));
  lat=netcdf.getVar(ncid, netcdf.inqVarID(ncid, 'lat_rho'));
  mask=netcdf.getVar(ncid, netcdf.inqVarID(ncid, 'mask_rho'));
  result=close(ncid);
  mask(mask==0)=NaN;
  pcolor(lon,lat,mask.*field)
  shading flat
  axis image
  caxis([min(min(fieldc)) max(max(fieldc))])
  colorbar
  axis([min(min(lonc)) max(max(lonc)) min(min(latc)) max(max(latc))])
  title(['\bf ',fieldname,' PARENT'])
end


return
