PRO ovsa_mesh, im, f, t, inew, fnew, tnew


  sz = SIZE(im)
  
  x = FLTARR(sz(1),sz(2))
  y = x
  z = x
 
  FOR i=0,sz(1)-1 DO x(i,*)=ALOG(f(i))
  FOR j=0,sz(2)-1 DO y(*,j)=t(j)

  FOR j=0,sz(2)-1  do begin
  hama=reform(im(*,j))
  if (n_elements(where(finite(hama))) lt 20) then hama(*)=!values.f_nan
  z(*,j)=hama
  endfor

  gs=[0.02,t(2)-t(1)]; (delta x, delta y)
  ps=[MIN(x), MIN(t), MAX(x), MAX(t)]; (min_xnew,min_ynew, max_xnew,max_xnew)
  n=LONG(sz(4))

  ok=WHERE(FINITE(z))
  x=x(ok)
  y=y(ok)
  z=alog(z(ok))
  
  TRIANGULATE,x,y,tr
  inew=exp(TRIGRID(x,y,z,tr,gs,ps,/missing))
  fnew=ps(0)+(ps(2)-ps(0))/(N_ELEMENTS(inew(*,0))-1)*FINDGEN(N_ELEMENTS(inew(*,0))) ;log value of fnew
  tnew=ps(1)+(ps(3)-ps(1))/(N_ELEMENTS(inew(0,*))-1)*FINDGEN(N_ELEMENTS(inew(0,*)))

END