;# 6/26/16
;# should have done this sooner, maybe.

;# means764.pl plots U(z) from mean.0500 for gasp764p (t764_evolve) and gasp764m (t764_fixed)

   nz = 352
   zm = fltarr(nz)
   zp = fltarr(nz)
   um = fltarr(nz)
   up = fltarr(nz)

   ;opfile1 = './gasp764p-m/mean.0500'
   ;opfile2 = './gasp764m-m/mean.0500'
   opfile1 = './t764_evolve/mean.0500'
   opfile2 = './t764_fixed/mean.0500'

   close,1,2
   openr,1,opfile1
   openr,2,opfile2

   for i = 0, nz-1 do begin
     readf,1,zz1,uu1,aa1,bb1,cc1,dd1,ee1
     readf,2,zz2,uu2,aa2,bb2,cc2,dd2,ee2
     zp[i] = zz1
     zm[i] = zz2
     up[i] = uu1
     um[i] = uu2
   endfor
   close,1,2
   print,'z:',max(zp-zm),min(zp-zm)
   print,'u:',max(up-um),min(up-um)

   window,0,xsize=800,ysize=600
   plot,up,zp,charsize=2,xrange=[100,220]
   oplot,um,zm,line=2,thick=2

   end
