program itgc
c
c     (c) Copyright 1991 to 1995 by Michael A. Beer, William D. Dorland, 
c     and Gregory W. Hammett. ALL RIGHTS RESERVED.
c
      implicit none
      include 'netcdf.inc'
      include 'itg.par'
      include 'itg.cmn'

      real time(mz),dt(mz)
      integer icontrol(2)
      integer ntim,ntimr
      real tim,dt1

c netcdf variables
      integer ncid,status
      character filename*80

      icontrol(1)=0
      icontrol(2)=0
      call itg(icontrol,time,dt)
c
c I find that the info written by prtime is useless. BD
c
c      call prtime

c Save the results to the *.res file:
      ntim=0
      ntimr=0
      if(binary) then
         open(unit=7,file=runname(1:lrunname)//'.bres',
     .        form='unformatted')
         call bwpunch(7,ntim,ntimr,tim,dt1,time,dt)
         close(7)
      else         
         filename=runname(1:lrunname)//'.nc'
         status = nf_create(filename, 0, ncid) ! overwrites old
         call wpunchnc(ncid,ntim,ntimr,tim,dt1,time,dt)
         status = nf_close(ncid)
c         open(unit=7,file=runname(1:lrunname)//'.res')
      endif

      stop         ! on unicos cray, also prints out CPU time summary.
      end