program driver
c----------------------------------------------------------------
c
c     Marjorie Friedrichs 03/03/03
c
c     Subroutines
c       readdata.f:      reads in data for computation of cost fn
c       readforcing.f:   reads in forcing files
c       readinit.f:      reads in initial conditions
c                                 ecosystem parameters
c       model.f:         calls physical and biological subroutines
c          biosub.f:     general biological processes
c              derivs.f: specific ecosystem model equations
c          horizadv.f:   applies horizontal advection
c          tridag.f:     tridiagonal solver for diffusion
c          sink.f:       applies detrital sinking
c       output.f:        outputs model equivalents of data
c       calccost.f:      outputs cost function
c
c-----------------------------------------------------------------
c     
      include 'common_blocks.h'
c
      integer nsv,startday
c
      homedir='/home/marjy/RegTestBeds/Workshop1/FinalTestbed/'
      datdir='/home/marjy/RegTestBeds/testbedfiles/'
c
      write(*,*)'Enter 1 for Arabian Sea testbed,' 
      write(*,*)'   or 2 for Equatorial Pacific testbed'
      read(5,*)ntb
      write(*,*)'Enter number of state variables that will undergo'
      write(*,*)'   the physical processes of diffusion, advection'
      write(*,*)'   and mixing'
      read(5,*)nsv
      nbio=nsv+2            
      ichl=nbio-1           ! store chlorophyll as a diagnostic
      ipp=nbio              ! store primary production as a diagnostic
c     
      write(*,*)'Enter value for remineralization rate [1/day]:'
      read(5,*)reminrate
      write(*,*)'Enter average sinking velocity of detrital'
      write(*,*)' nitrogen collected at 800m traps [m/day]:'
      read(5,*)wtrap
c
C     Specify the domain, vertical resolution, no. of vertical 
C     points and run duration
c
      dz = 10.d0                ! vertical resolution (meters)
      nz = 15                   ! number of points in the vertical domain
      delt=3600.
      startday=273.d0
      isteps=10660
c
      read(5,*)(wnsv(n),n=1,nsv)
      read(5,*)(bbcnsv(n),n=1,nsv)
      read(5,*)(HAnsv(n),n=1,nsv)
      read(5,*)(aeonsv(n),n=1,nsv)
      read(5,*)(Cnsv(n),n=1,nsv)
      read(5,*)(unitnsv(n),n=1,nsv)
c
      call readdata
      call readforcing
      call readinit
      call model
      call calccost
      call output
c
      stop
      end