Jason Sippel

Great code! After chatting with the author, I have found that the below addendum is useful for plotting a small subdomian of model output. This would be useful if, for example, you are tracking a tropical cyclone in a Lagrangian system with ilat,flat,ilon, and flon sampled from your model grid based on the number of grid points from the storm center.

If you don't do something similar to this, the physical size of the map frame on the screen may not match from one figure to the next.

Here's the code:

    tightmap
    xlim=get(gca,'xlim');
    ylim=get(gca,'ylim');
    
   cutfactor = 1.3;
   xcenter=mean(xlim);
   xoffset=diff(xlim)/2;
   ycenter=mean(ylim);
   yoffset=diff(ylim)/2;
   set(gca,'Xlim',[xcenter-cutfactor*xoffset xcenter+cutfactor*xoffset]);
   set(gca,'Ylim',[ycenter-cutfactor*yoffset ycenter+cutfactor*yoffset]);