pro test

 amp=fltarr(72,59,55)
 glat=fltarr(59)

 openw,20,'out_amp'

 openr,10,'zonwin_sep_all_Dam.dat'
 readf,10,amp

 for jlat=0,58 do begin
  glat(jlat) = -87 + jlat*3
  printf,20,jlat,glat(jlat),amp(2,jlat,2),amp(2,jlat,30),amp(2,jlat,54)
 endfor

 close,10

 end