In line 167 of your polarPcolor.m (caxis([quantile(Z(:),0.01),quantile(Z(:),0.99)]))), why did you use the 1% and 99% quantiles instead of the min and max? Just for my understanding, this function sets the min and max of the color scale, doesn't it (as in which color corresponds to which value)? I would assume that the quantile will skew your data if you have a wide range of values and a lot of data, which is the case in radar applications.

It appears to me that the optional argument 'typerose' may have a bug in it (or I may have misunderstood its intended use)..

As a quick example, using your example code:

R = linspace(3,10,100);
theta = linspace(0,180,360);
Z = linspace(0,10,360)'*linspace(0,10,100);
figure;
polarPcolor(R,theta,Z,'Ncircles',3, 'typeRose', 'meteo');

gives rise to the "default" orientation for your polarPcolor, with zero degrees up towards north, and theta going clockwise around the polar plot to the maximum value of Z directly at the south.. Which is all good, and as expected based on R, theta and Z.

However if one wants to change orientation to have zero deg theta at the east (similar to Matlab's own default value in polarplot), then with the following example code:

R = linspace(3,10,100);
theta = linspace(0,180,360);
Z = linspace(0,10,360)'*linspace(0,10,100);
figure;
polarPcolor(R,theta,Z,'Ncircles',3, 'typeRose', 'default');

we suddenly have shifted only the spokes, ticks etc. and revert orientation of them but the data itself is still plotted with the original orientation and from "north to south" (not "east to west" as the spokes/ticks etc imply).

See example plot:

As a side note, having "default" as an argument but not having this value be the actual "default" is perhaps making things more confusing than need be.