MCOR/moses_super_plot.pro at master · jacobdparker/MCOR · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
; Program to plot MOSES super exposures with a uniform scaling.  Used to make images for ParkerKankelborg.tex

restore,'moses_super.sav'


; Comment out old coyote graphics ... a thing of the past

;; cgwindow,'cgimage',super_plus,maxvalue=5,minvalue=0,ctindex=3,/keep_aspect_ratio,waspect=.5,wtitle='plus'

;; cgwindow,'cgimage',super_zero,maxvalue=5,minvalue=0,ctindex=3,/keep_aspect_ratio,waspect=.5,wtitle='zero'

;;  eit_pz = super_plus-super_zero
;;  pz = (eit_pz-mean(eit_pz))
;;  pz /= sqrt(total(total(pz^2,2),1))

;; cgwindow,'cgimage',pz,maxvalue=.0005,minvalue=-.0005,ctindex=3,/keep_aspect_ratio,waspect=.5
;; cgwindow,'cgpolygon', [0.35, 0.35, 0.65, 0.65, .35], [0.7, 0.97, 0.97, 0.7, 0.7], /NORMAL, COLOR='green',thick='10',/addcmd
;; cgwindow,'cgpolygon', [0.65, 0.65, 0.95, 0.95, .65], [0.25, 0.5, 0.5, 0.25, 0.25], /NORMAL, COLOR='green',thick='10',/addcmd
;; cgwindow,'cgpolygon', [0.02, 0.02, 0.07, 0.07, .02], [0.02, 0.2, 0.2, 0.02, 0.02], /NORMAL, COLOR='green',thick='10',/addcmd,wtitle='Plus - Zero'

  restore,'mosesI_pointing.sav'

  plus = image(super_plus,x,y,max_value = 5, min_value = 0, rgb_table = 3,buffer=1,axis_style=1,title = 'Time Averaged MOSES 304 m = +1 Image', xtitle ='X (arcsec)',ytitle = 'Y (arcsec)')
  zero = image(super_zero,x,y,max_value = 5, min_value = 0, rgb_table = 3,buffer=1,axis_style=1,title = 'Time Averaged MOSES 304 m = 0 Image', xtitle ='X (arcsec)',ytitle = 'Y (arcsec)')

  eit_pz = super_plus-super_zero
  pz = (eit_pz-mean(eit_pz))
  pz /= sqrt(total(total(pz^2,2),1))
  dif_plot = image(pz,x,y,max_value =.0005, min_value = -0.0005, rgb_table = 3,buffer=1,axis_style=1,title = 'Time Averaged MOSES Difference Image (m=+1 minus m=0)', xtitle ='X (arcsec)',ytitle = 'Y (arcsec)')

  colr = 'lime'

  ;; p = polygon([0.35, 0.35, 0.65, 0.65, .35], [0.65, 0.8, 0.8, 0.65, 0.65], /normal, target=dif_plot,fill_transparency=100,color=colr,thick=3)
  ;; p = polygon([0.65, 0.65, 0.95, 0.95, .65], [0.35, 0.5, 0.5, 0.35, 0.35], /normal, target=dif_plot,fill_transparency=100,color=colr,thick=3)
  ;; p = polygon([0.02, 0.02, 0.07, 0.07, .02], [0.2, 0.3, 0.3, 0.2, 0.2], /normal, target=dif_plot,fill_transparency=100,color=colr,thick=3)


  plus.save, "../MCOR_Paper/super_plus.eps"
  zero.save, "../MCOR_Paper/super_zero.eps"
  dif_plot.save, "../MCOR_Paper/super_pz.eps"


end