load('NLDASWINDDAILY.mat','lat','lon','el');
motemp=matfile('TMXMN');
mowind=matfile('NLDASWINDDAILY');
moshum=matfile('SHUM_NLDAS');
morad=matfile('NLDASRADDAILY');
cd C:\NLDAS_DATA\Location_specific_data
for loni=1
    cd(['lat_' num2str(lat(loni,1))]);
    for lati=1 fprintf('**loni = %d & lati = %d**\n',loni,lati);
        templat=lat(loni,lati);
        templon=lon(loni,lati);
        clear data
        load(['NLDAS_' num2str(templat) '_' num2str(templon) '.mat']);
        tic
        parfor i=1:32 %fprintf('Now on year %d\n',i);
            %             data((1+(366*(i-1)):366*i),1)=1979+(i-1);
            %             data((1+(366*(i-1)):366*i),2)=1:366;
            data((1+(366*(i-1)):366*i),3)=squeeze(motemp.TMIN(loni,lati,:,i));
            data((1+(366*(i-1)):366*i),4)=squeeze(motemp.TMAX(loni,lati,:,i));
            data((1+(366*(i-1)):366*i),5)=squeeze(mowind.WIND(loni,lati,:,i));
            data((1+(366*(i-1)):366*i),6)=squeeze(moshum.DATA(loni,lati,:,i));
            data((1+(366*(i-1)):366*i),7)=squeeze(morad.RADIATION(loni,lati,:,i));
        end
        %         end
        %         ell=el(loni,lati);
        %         for i=1:366*32 %fprintf('Now on  %d\n',i);
        %             temp(i,8)=calc_es(temp(i,3),temp(i,4));
        %             [temprh,tempea,temptd]=EAcalcfromNLDAS(temp(i,6),temp(i,8),ell);
        %             temp(i,9)=temprh;
        %             temp(i,10)=tempea;
        %             temp(i,11)=temptd+273;
        %             clear temprh tempea temptd
        %         end
        %         data=temp;clear temp;clear ell;
        %         mkdir Location_specific_data; cd Location_specific_data
        save(['NLDAS_' num2str(templat) '_' num2str(templon) '.mat'], 'data')
        toc
        %         clear memory;
    end
    cd ..
end