Slow piece of code

Hello everyone. I have a working piece of code. But for large values it is really slow. Code is pasted below:

for i = 1:T

    [rowX1,colX1] = find(XI==X_trace1(i),1); 
    [rowY1,colY1] = find(YI==Y_trace1(i),1);

        if isempty(colX1) || isempty(rowY1)
           display('out of grid')
        else
lh1=m_quiver(XI(rowY1,colX1),YI(rowY1,colX1),UI(rowY1,colX1)*0.00027,VI(rowY1,colX1)*0.00027,0);              
hq1(i)=quiver(XI(rowY1,colX1),YI(rowY1,colX1),0.00027*UI(rowY1,colX1)/cos(YI(rowY1,colX1)*pi/180),VI(rowY1,colX1)*0.00027,'Autoscale','off');

        end
end

XI and YI (also UI and VI) are big matrices (2200 x 2200) produced using meshgrid. hq1 is the resultant array of this loop. For large value of T, this code takes long time. I would have no problem if this was the only part of my full code, But i have 10 similar pieces in full code. Is there any way to make this piece of code run faster. I will appreciate any kind of suggestion to improve this code.

thanks

[Walter Roberson]

Do I understand correctly that the output is hq1, a vector of T quiver plot handles? What do you do with all the different quiver plots?

[Aeronomy]

Hello,

I using m_quiver to plot quivers on map, and hq1 to find end point (lat and lon) of quivers using quiver. Actually I am using this code for tracing particle trajectories.

thanks

[Walter Roberson]

? quiver() returns the handle of the quivergroup which is a hggroup(). Are you later analyzing the contents of that quivergroup object ?

[Aeronomy]

This loop is a part of biger loop. After getting handle hq1, i am analyzing it to find lat and lon so that i can use those lat and lon during next big loop run.

[per isakson]

what does profile say?

[per isakson]

"The profile function helps you debug and optimize MATLAB® code files by tracking their execution time." Run

    doc profile