diff options
author | simulation <simulation@nlo-ext3.ioq.uni-jena.de> | 2016-02-09 09:54:46 +0100 |
---|---|---|
committer | simulation <simulation@nlo-ext3.ioq.uni-jena.de> | 2016-02-09 09:54:46 +0100 |
commit | 6be0f285adb27f60f058c2937a7885899e17cb88 (patch) | |
tree | be50a15314fe28c5b7b8f721233a34e71ca5b3b4 /idl/idl2avi.body | |
download | lpic-master.tar.xz |
Diffstat (limited to 'idl/idl2avi.body')
-rw-r--r-- | idl/idl2avi.body | 84 |
1 files changed, 84 insertions, 0 deletions
diff --git a/idl/idl2avi.body b/idl/idl2avi.body new file mode 100644 index 0000000..75f1a9e --- /dev/null +++ b/idl/idl2avi.body @@ -0,0 +1,84 @@ +; ---------------------------------------------------------------------------------------- +; body for idlmovie.pro +; ---------------------------------------------------------------------------------------- + + framepre = "frame_" + framepost = ".eps" + + s1 = dimx + s2 = dimv + xmin = xmax0/dimx0 * cutx - xoffset + xmax = xmin + xmax0/dimx0 * dimx + ymax = vmax0 - vmax0/dimv0 * cutv + x = xmin + (xmax-xmin)/s1 * findgen(s1) + y = -ymax + 2.0*ymax/s2 * findgen(s2) + + print, cutv + print, ymax + +; set_plot,"x" +; window, 1, xpos=540, ypos=28, xsize=dimx+200, ysize=dimv+180, title="phasespace(x,v)" + xoff=120 + yoff=100 + + loadct,13 + tvlct, r, g, b, /get + r(0) = 255 + g(0) = 255 + b(0) = 255 + r(1) = 0 + g(1) = 0 + b(1) = 0 + tvlct, r, g, b + + print, file0 + + if (special_file lt -10) then begin + i=-1 + repeat begin + i = i+1 + f = floor( 1000.0*(file_begin+increment*i)+0.5 ) / 1000.0 + if (f lt 10) then s = string( format='(F5.3)',f ) + if (f ge 10) then s = string( format='(F6.3)',f ) + if (f lt 10) then s2 = '0'+s + if (f gt 10) then s2 = s + file = file0 + s +; print, file + openu,1,file + + set_plot,"ps" + device, /encapsulated, file=framepre+file0+s2+framepost, xsize=40, ysize=40, bits=8 + + plot, x, y, xstyle=1, ystyle=1, xrange=[xmin,xmax], yrange=[-ymax,ymax],$ + xtitle="!3x/!7k!3", ytitle="!3v/c!3",ticklen=-0.02,$ + charsize=1.5, color=1, position=[xoff-1,yoff-1,xoff+s1+1,yoff+s2+1],$ + /noerase, /nodata, /device + + a = assoc(1,bytarr(dimx0,dimv0,/nozero)) + phase = extrac(a(0)+1,cutx,cutv,dimx,dimv) + tvscl,alog10(phase),xoff,yoff + device,/close + close,1 + endrep until ( f gt file_end-increment ) + endif else begin + if (special_file lt 10) then s = string( format='(F5.3)',special_file ) + if (special_file lt 0) then s = string( format='(F6.3)',special_file ) + if (special_file gt 10) then s = string( format='(F6.3)',special_file ) + file = file0 + s + print, file + openu,1,file + a = assoc(1,bytarr(dimx0,dimv0,/nozero)) + phase = extrac(a(0)+1,cutx,cutv,dimx,dimv) + tvscl,alog10(phase),xoff,yoff + close,1 + endelse +end + + + + + + + + + |