code: def plot_site(fignum, SiteRec, data, key):
      
      print('Site mean data: ')
      print('   dec    inc n_lines n_planes kappa R alpha_95 comp coord')
      print(SiteRec['site_dec'], SiteRec['site_inc'], SiteRec['site_n_lines'], SiteRec['site_n_planes'], SiteRec['site_k'],
            SiteRec['site_r'], SiteRec['site_alpha95'], SiteRec['site_comp_name'], SiteRec['site_tilt_correction'])
      print('sample/specimen, dec, inc, n_specs/a95,| method codes ')
      for i in range(len(data)):
          print('%s: %s %s %s / %s | %s' % (data[i]['er_' + key + '_name'], data[i][key + '_dec'], data[i]
                                            [key + '_inc'], data[i][key + '_n'], data[i][key + '_alpha95'], data[i]['magic_method_codes']))
      plot_slnp(fignum, SiteRec, data, key)
      plot = input("s[a]ve plot, [q]uit or <return> to continue:   ")
      if plot == 'q':
          print("CUL8R")
          sys.exit()
      if plot == 'a':
          files = {}
          for key in list(EQ.keys()):
              files[key] = site + '_' + key + '.' + fmt
          save_plots(EQ, files)
docs: deprecated (used in ipmag)