; Locations of Asymmetric V
 
;f=read_prf('mhd12_6km_342.6.prf')

v=reform(f.prf[*,3,*,*])
v1=v[0:99,*,*]
v2=v[100:199,*,*]

kb=0
kr=0
for i=0,1007 do begin $
  for j=0,1007 do begin &$
    v1x=max(v1[*,i,j],min=min) & vin=min &$
    i1=where(v1[*,i,j] eq v1x) &$
    i2=where(v1[*,i,j] eq v1n) &$
    i0=aver([i1,i2]) &$
    v1b=max(abs(v1[0:i0,i,j])) &$
    v1r=max(abs(v1[i0:99,i,j])) &$
    if v1b gt .2 and v1r lt .05 then begin kb=kb+1 & print,'b:',i0,v1b,v1r,kb & endif &$
    if v1r gt .2 and v1b lt .05 then begin kr=kr+1 & print,'r:',i0,v1b,v1r,kr & endif &$
  endfor &$
end
print,kb,kr

ib=complexarr(kb)
ir=complexarr(kr)
kb=-1
kr=-1
for i=0,1007 do begin $
  for j=0,1007 do begin &$
    v1x=max(v1[*,i,j],min=min) & v1n=min &$
    i1=where(v1[*,i,j] eq v1x) &$
    i2=where(v1[*,i,j] eq v1n) &$
    i0=aver([i1,i2]) &$
    v1b=max(abs(v1[0:i0,i,j])) &$
    v1r=max(abs(v1[i0:99,i,j])) &$
    if v1b gt .2 and v1r lt .05 then begin &$
      kb=kb+1 &$
      ib[kb]=complex(i,j) &$
    endif &$
    if v1r gt .2 and v1b lt .05 then begin &$
      kr=kr+1 &$
      ir[kr]=complex(i,j) &$
    endif &$
  endfor &$
end

save,/ver,file='Vrbasym.save',ib,ir