C Program to read in UKMO files and add lat-long to them
C Expect format:
C   Header-line
C   lbnpt lbrow bzx bdx bzy bdy
C   data\ndata\n...
C
C Small note: if you are reading in gaussian-grid files 
C (eg ERA) then bzy and bdy will be zero. In which case
C you need to know the gaussian latitudes. They are
C (from 45 to the pole):
C     -45.420300      -46.541801      -47.663200      -48.784698      -49.906200      -51.027699      -52.149200      -53.270699
C     -54.392101      -55.513599      -56.635101      -57.756599      -58.877998      -59.999500      -61.120998      -62.242500
C     -63.363899      -64.485397      -65.606903      -66.728302      -67.849800      -68.971199      -70.092697      -71.214104
C     -72.335602      -73.457001      -74.578400      -75.699799      -76.821198      -77.942596      -79.064003      -80.185303
C     -81.306602      -82.427803      -83.548897      -84.669899      -85.790604      -86.910797      -88.029404      -89.141602

	character*79 filename,header
	integer lbnpt, lbrow, i, j 
	real bzx, bdx, bzy, bdy, x

	read(unit=*,fmt='(a)')filename
	open(unit=20,file=filename,status='old')
	read(unit=20,fmt='(a)')header
	read(unit=20,fmt=*)lbnpt, lbrow, bzx, bdx, bzy, bdy

	print *,header
	print *,'Lat Lon Value'
	do j=1,lbrow
	do i=1,lbnpt
	  read(unit=20,fmt=*)x
	  print *,bzy+j*bdy,bzx+i*bdx,x
        enddo
	enddo

	end