*script to merge data from multiple netcdf files, e.g. hgt.1948.nc, hgt.1949.nc, etc
*edit start and end dates as needed
*view output with hgt500_merged.ctl (edit tdef for longer times)


'reinit'
iyear=1948

while (iyear <= 1949)
 'sdfopen hgt.'iyear'.nc'
 'q file'
 res=sublin(result,5)
 tsize=subwrd(res,12)
 say tsize
 tim = 1
 'set gxout fwrite'
 'set fwrite -ap merged_h500.grads'
 'set x 1 144'
 'set y 1 73'
while (tim <= tsize)
  'set t 'tim
  'd hgt(lev=500)'
  tim=tim+1
endwhile
 'disable fwrite'
 'close 1'
 iyear=iyear+1
endwhile  
'quit'