sum_jup.pro (Kronos IDL library)

;+
; Contains the sum_jup procedure
;
; :Author:
;    Philippe Zarka
;
; :History:
;    2003/12/08: Created
;
;    2003/12/08: Last Edit
;-
;
;+
; summary plots Kronos entre datedeb et datefin, 12h / page, echelle lin.
;
; :Uses:
;    aj_amj, dyn, extrait_mode, liste_comp, plot_spdyn
;
; :Params:
;    datedeb: in, required, type=sometype
;       A parameter named datedeb
;    datefin: in, required, type=sometype
;       A parameter named datefin
;
; :Keywords:
;    nblocks: in, optional, type=sometype
;       A keyword named nblocks
;-
pro SUM_JUP, datedeb,datefin, nblocks = nblocks


if not(keyword_set(nblocks)) then nblocks = 1

set_plot,'ps'

  device, yoff=27.7, xsize=27., xoff=0.8, ysize=19.5, /landscape    ; sum/gif
; device, yoff=29., xsize=28.3, xoff=0.8, ysize=19.5, /landscape        ; big
; device, yoff=27., xsize=24.,  xoff=2.,  ysize=17.,  /landscape
; device, yoff=0.7, ysize=28.3, xoff=0.8, xsize=19.5, /portrait         ; big
; device, yoff=2.3, ysize=25.,  xoff=1.5, xsize=18.,  /portrait

!p.font=0

path='/user4/rpws/'
nomfich='/poub/rpws/sum.lis'

for i=datedeb,datefin do begin

  nom=strtrim(long(i),2)+'_sumjup1'
  device,file=path+'sum/'+nom+'.ps'
  !p.multi=[0,1,2]
  ii=strtrim(aj_amj(i),2)

  PLOT_SPDYN, i,  0,12,  0,1,1, 1.,16.125, 1,0, 72,2, $
	'     '+strmid(ii,0,4)+'/'+strmid(ii,4,2)+'/'+strmid(ii,6,2), $
        nblocks=nblocks
  LISTE_COMP, i, 0, i, 12, nomfich, 0

  EXTRAIT_MODE, nomfich, 92, 'Df', aj,hd,hf
  plot, [0,12],[0,0],/noerase,xra=[0,12],xstyle=13,yra=[-0.1,0.1],ystyle=13,$
	pos=[0.08,0.045,0.88,0.055]
  if aj(0) ne -1 then $
	for j=0,n_elements(aj)-1 do oplot,[hd(j),hf(j)],[0,0],thick=8.
  xyouts,12.3,-0.05,'df-H1',charsize=0.7
  EXTRAIT_MODE, nomfich, 130, 'Df', aj,hd,hf
  plot, [0,12],[0,0],/noerase,xra=[0,12],xstyle=13,yra=[-0.1,0.1],ystyle=13,$
	pos=[0.08,0.060,0.88,0.070]
  if aj(0) ne -1 then $
	for j=0,n_elements(aj)-1 do oplot,[hd(j),hf(j)],[0,0],thick=8.
  xyouts,12.3,-0.05,'df-H2',charsize=0.7

  device,/close

  spawn,'pstogif -depth=8 '+path+'sum/'+nom+'.ps'
  read_gif,path+'sum/'+nom+'.gif',x, R, G, B
; tvlct, R, G, B
  write_gif,path+'gif/'+nom+'.gif',rotate(x,1), R, G, B
  spawn,'/rm '+path+'sum/'+nom+'.gif'
  spawn,'/rm '+path+'sum/'+nom+'.ppm'

  nom=strtrim(long(i),2)+'_sumjup2'
  device,file=path+'sum/'+nom+'.ps'
  !p.multi=[0,1,2]
  ii=strtrim(aj_amj(i),2)

  PLOT_SPDYN, i, 12,24,  0,1,1, 1.,16.125, 1,0, 72,2, $
	'     '+strmid(ii,0,4)+'/'+strmid(ii,4,2)+'/'+strmid(ii,6,2), $
        nblocks=nblocks
  LISTE_COMP, i, 12, i, 24, nomfich, 0

  EXTRAIT_MODE, nomfich, 92, 'Df', aj,hd,hf
  plot, [12,24],[0,0],/noerase,xra=[12,24],xstyle=13,yra=[-0.1,0.1],ystyle=13,$
	pos=[0.08,0.045,0.88,0.055]
  if aj(0) ne -1 then $
	for j=0,n_elements(aj)-1 do oplot,[hd(j),hf(j)],[0,0],thick=8.
  xyouts,24.3,-0.05,'df-H1',charsize=0.7
  EXTRAIT_MODE, nomfich, 130, 'Df', aj,hd,hf
  plot, [12,24],[0,0],/noerase,xra=[12,24],xstyle=13,yra=[-0.1,0.1],ystyle=13,$
	pos=[0.08,0.060,0.88,0.070]
  if aj(0) ne -1 then $
	for j=0,n_elements(aj)-1 do oplot,[hd(j),hf(j)],[0,0],thick=8.
  xyouts,24.3,-0.05,'df-H2',charsize=0.7

  device,/close

  spawn,'pstogif -depth=8 '+path+'sum/'+nom+'.ps'
  read_gif,path+'sum/'+nom+'.gif',x, R, G, B
; tvlct, R, G, B
  write_gif,path+'gif/'+nom+'.gif',rotate(x,1), R, G, B
  spawn,'/rm '+path+'sum/'+nom+'.gif'
  spawn,'/rm '+path+'sum/'+nom+'.ppm'

endfor

return
end