plot_ns_maps.pro

pro plot_ns_maps, obsname=obsname,maindir=maindir,nsdir=nsdir

  ; Plot the individual maps and also a combined one
  ;
  ; Work in progress version - only a the most recent observations
  ;
  ; Options
  ; obsname       - Which NuSTAR observations (default ='201409')
  ; maindir       - Main directory of where the NuSTAR data is kept (default for IGH system but only need if no *.dat files)
  ; nsdir         - Specific directory where this NuSTAR obs is kepts (default for IGH system but only need if no *.dat files)
  ;
  ;
  ; 30-Aug-2016 IGH - Created
  ;~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

  if (n_elements(obsname) ne 1) then obsname='201607'
  if (n_elements(maindir) ne 1) then maindir='~/data/ns_data/

  ; What is the minimum energy we want for the image?
  min_eng=2.5

  ;-------------------------------------------

  if (obsname eq '201509') then nsdir='obs5/'
  if (obsname eq '201602') then nsdir='obs6/'
  if (obsname eq '201604') then nsdir='obs7/'
  if (obsname eq '201607') then nsdir='obs8/'
  ffa=file_search(maindir+nsdir,'*FPMA.fits')

  nf=n_elements(ffa)
  sr=2
  for i=0, nf-1 do begin
    fits2map,ffa[i],mm
    mms=mm
    mms.data=gauss_smooth(mms.data,sr)

    @post_outset
    !p.multi=0

    loadct,74,/silent
    reverse_ct
    tvlct,r,g,b,/get
    r[0]=0
    g[0]=0
    b[0]=0
    r[1]=255
    g[1]=255
    b[1]=255
    tvlct,r,g,b

    dnl=1e-3
    dmx=1e1
    plim=1500
    xr=[-1,1]*plim
    yr=[-1,1]*plim

    ftemp=strsplit(ffa[i],'/',/extract)
    ftemp=ftemp[n_elements(ftemp)-1]
    fname=(strsplit(ftemp,'.',/extract))[0]

    set_plot,'ps'
    device, /encapsulated, /color, / isolatin1,/inches, $
      bits=8, xsize=4.5, ysize=4.5,file='maps/'+fname+'.eps'

    plot_map,mms,/log,dmin=dnl,dmax=dmx,position=[0.17,0.13,0.87,0.83],$
      title=mms.id,grid_sp=30,xrange=xr,yrange=yr,$
      gcolor=0,bot=1,chars=0.8,xtitle='x [arcsec]',ytitle='y [arcsec]'

    plot_map_cb_igh,alog10([dnl,dmx]),position=[0.15,0.96,0.85,0.98],color=0,chars=0.8,$
      cb_title='NuSTAR [log!D10!N count s!U-1!N]',bottom=1,format='(f4.1)'

    device,/close
    set_plot, mydevice

  endfor



  stop
end