cycling/filled_contours.m

   1 function filled_countours(ps,an)
   2 % plots filled contours for forecast and analysis times in seconds
   3
   4 tmax = max(ps.red.tign_g(:))-3600;
   5 tmin = min(ps.red.tign_g(:));
   6 clines = linspace(tmin,tmax,20);
   7
   8 if isfield(ps.red,'red')
   9     figure,contourf(ps.red.red.fxlong(1:10:end,1:10:end),ps.red.red.fxlat(1:10:end,1:10:end) ...
  10         ,ps.red.red.tign(1:10:end,1:10:end),clines,'k')
  11 else
  12     figure,contourf(ps.red.fxlong,ps.red.fxlat,ps.red.tign,clines,'k')
  13 end
  14 hold on
  15 scatter(ps.points(:,2),ps.points(:,1),'*r')
  16 xlabel('Lon'),ylabel('Lat'),title('Forecast and active fire detections')
  17 hold off
  18
  19
  20 figure,contourf(ps.red.fxlong,ps.red.fxlat,an,clines,'k')
  21 hold on
  22 scatter(ps.points(:,2),ps.points(:,1),'*r')
  23 xlabel('Lon'),ylabel('Lat'),title('Analysis and active fire detections')
  24 hold off
  25
  26 end