code: def surfplot(self, z, titletext):
      
      if self.latlon:
        plt.imshow(z, extent=(0, self.dx*z.shape[0], self.dy*z.shape[1], 0)) #,interpolation='nearest'
        plt.xlabel('longitude [deg E]', fontsize=12, fontweight='bold')
        plt.ylabel('latitude [deg N]', fontsize=12, fontweight='bold')
      else:
        plt.imshow(z, extent=(0, self.dx/1000.*z.shape[0], self.dy/1000.*z.shape[1], 0)) #,interpolation='nearest'
        plt.xlabel('x [km]', fontsize=12, fontweight='bold')
        plt.ylabel('y [km]', fontsize=12, fontweight='bold')
      plt.colorbar()
  
      plt.title(titletext,fontsize=16)
docs: Plot if you want to - for troubleshooting - 1 figure