program read_co2

c Read in the monthly NPP (360*180) from 1998-2006. The unit
c is g/(m^2 mon). The data is from Jim Collatz.

      real(kind=4) NPP(360,180)

      open(unit=11,file='npp.9806', form='FORMATTED')

      lun_in=11
      irec=1
      do nyear=1,9
      do mon=1,12

c read in and flip from N-S to S-N
 
      read(lun_in,7) NPP
7     format(360f10.3)

      enddo
      enddo

      stop
      end