set(gca,'YDir','normal'); ir=ir(~isnan(ir));jr=jr(~isnan(jr)); dawn=[3.2472:25.2472 25.9139:31.9139 32.5806:62.5806]; dawnind=round((dawn*24-t(1))*NstepsPerHr+1); dawnir=interp1(ir,1:numel(ir),dawnind); arg=angle(cwt(ir+Nsteps*(jr-1))); ...