use wgrib2api

	real, allocatable :: grid(:,:), lat(:,:), lon(:,:)
	integer (kind=8) :: date
	character (len=200) :: file, inv, var, invline
	character (len=400) :: gridline

	file = 'gep19.t00z.pgrb2af180'
	inv = '@mem:0'

!	make inv file, save in memory file #0
	iret = grb2_mk_inv(file, inv)
	if (iret.ne.0) stop 1

	var = ':TMP:500 mb:'
	date = 2009060500_8
	iret = grb2_inq(file,inv,var,ref_date=date,data2=grid,lat=lat,lon=lon, &
               desc=invline,grid_desc=gridline)
	if (iret.ne.1) stop 4
	write(*,*) 'T500(1,1) =',grid(1,1),' lat/lon=',lat(1,1), lon(1,1)
	write(*,*) 'date=',date, ' nx=', size(grid,1), ' ny=', size(grid,2)
	write(*,*) 'inventory=',trim(invline)
	write(*,*) 'grid=',trim(gridline)

	stop
	end