Matlab - loop para netcdf

% Specify the folder where the files are.

myFolder = 'G:\SWICCA\DATA_validation\ESA\data\daily_files\COMBINED\v02.2\2010\files';

% Check to make sure that folder actually exists.  Warn user if it doesn't.

if ~isdir(myFolder)

  errorMessage = sprintf('Error: The following folder does not exist:\n%s', myFolder);

  uiwait(warndlg(errorMessage));

  return;

end

% Read the names of nc-files in the folder 

ncfiles = dir('*.nc') ;

Nfiles = length(ncfiles) ;

% Loop for each nc-file

for i = 1:Nfiles

    ncid=netcdf.open('ncfiles', 'NC_NOWRITE');

    varname= netcdf.inqVar(ncid(i),4);

    varid = netcdf.inqVarID(ncid(i),varname(i));

    data(i) = netcdf.getVar(ncid(i),varid(i));

    netcdf.close(ncid(i));

end

lon1=netcdf.getVar(ncid,0,0,1440);

lat1=netcdf.getVar(ncid,1,0,720);

%crop to interest area

[longrid,latgrid]=meshgrid(lon1,lat1);

data2=data((longrid>-1)&(longrid<0)&(latgrid>39)&(latgrid<40));

ncid1=netcdf.open('ESACCI-SOILMOISTURE-L3S-SSMV-COMBINED-20100101000000-fv02.2.nc', 'NC_NOWRITE');

varname = netcdf.inqVar(ncid1,4);

varid = netcdf.inqVarID(ncid1,varname);

data = netcdf.getVar(ncid1,varid);

lon1=netcdf.getVar(ncid1,0,0,1440);

lat1=netcdf.getVar(ncid1,1,0,720);

%crop to interest area

[longrid,latgrid]=meshgrid(lon1,lat1);

data2=data((longrid>-1)&(longrid<0)&(latgrid>39)&(latgrid<40));

function [fecha,LAT,LON,VX,VY,VEL] = leer_vientos(fname,latmin,latmax,lonmin,lonmax);

% Área de extracción de información

% latmax = xx;     % Extremo norte

% latmin = xx;   % Extremo sur

% lonmax = xx;     % Extremo este

% lonmin = xx;   % Extremo oeste

% Lee la fecha del archivo

fecha = double(ncread(fname,'time')/24 + datenum(1900,1,0));

% Lee los datos de latitud

LAT = ncread(fname,'latitude');

idla = find(and(LAT>=latmin,LAT<=latmax));

LA = LAT(idla);

% Lee los datos de longitud

LON = ncread(fname,'longitude');

idlo = find(and(LON>=lonmin,LON<=lonmax));

LO = LON(idlo);

[LAT,LON] = meshgrid(double(LA),double(LO));

% Lee los datos de viento (X)

vx = ncread(fname,'eastward_wind');

VX = vx(idlo,idla);

% Lee los datos de viento (Y)

vy = ncread(fname,'northward_wind');

VY = vy(idlo,idla);

VEL = double(sqrt(VX .^2 + VY.^2));