function ratoplot( X, F, T, ampthresh, render_method )
  ms = 1 ;
  m = '.' ;
  pos_colon = strfind(render_method, '_');
  ms = str2num(render_method(pos_colon(1)+1:pos_colon(2)-1)); 
  m = render_method(pos_colon(2)+1:end) ;
  render_method = render_method(1:pos_colon(1)-1); 
  
  if strcmp(m,'dot')
    m = '.';
  elseif strcmp(m,'circle')
    m = 'o';
  end
  
  hold on
  
  if strcmp(render_method, 'fast') 
    plot(T, F, '.k', 'markersize', ms, 'marker', m);
    
  elseif strcmp(render_method, 'color')
    magX = log(abs(X));
    maxX = max(magX);
    minX = min(magX);
    midX = (maxX+minX)/2 ;
    top_quarter = (maxX+midX)/2 ;
    bot_quarter = (minX+midX)/2 ;
  
    plot_index = find(magX<=bot_quarter);
    plot(T(plot_index),F(plot_index),'.b', 'markersize', ms, 'marker', m );
    
    plot_index = find(magX<=midX & magX>bot_quarter);
    plot(T(plot_index),F(plot_index),'.g', 'markersize', ms, 'marker', m );
    
    plot_index = find(magX<=top_quarter & magX>midX);
    plot(T(plot_index),F(plot_index),'.y', 'markersize', ms, 'marker', m );
    
    plot_index = find(magX>top_quarter);
    plot(T(plot_index),F(plot_index),'.r', 'markersize', ms, 'marker', m );
 
  end
  
   hold off