; 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,*,*]

; Blue large peak
k=0
for i=0,1007 do begin $
  for j=0,1007 do begin &$
    v1b=max(abs(v1[0:50,i,j])) &$
    v1r=max(abs(v1[50:99,i,j])) &$
    if v1b gt .2 and v1r lt .05 then k=k+1 &$
  endfor &$
end
print,k

ib=complexarr(k)
k=-1
for i=0,1007 do begin $
  for j=0,1007 do begin &$
    v1b=max(abs(v1[0:50,i,j])) &$
    v1r=max(abs(v1[50:99,i,j])) &$
    if v1b gt .2 and v1r lt .05 then begin &$
      k=k+1 &$
      ib[k]=complex(i,j) &$
    endif &$
  endfor &$
end

; Red large peak
k=0
for i=0,1007 do begin $
  for j=0,1007 do begin &$
    v1b=max(abs(v1[0:50,i,j])) &$
    v1r=max(abs(v1[50:99,i,j])) &$
    if v1r gt .2 and v1b lt .05 then k=k+1 &$
  endfor &$
end
print,k

ir=complexarr(k)
k=-1
for i=0,1007 do begin $
  for j=0,1007 do begin &$
    v1b=max(abs(v1[0:50,i,j])) &$
    v1r=max(abs(v1[50:99,i,j])) &$
    if v1r gt .2 and v1b lt .05 then begin &$
      k=k+1 &$
      ir[k]=complex(i,j) &$
    endif &$
  endfor &$
end

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