subroutine readdata
c---------------------------------------------------------------- 
c This subroutine reads in data to be included in cost function
c   computation
c----------------------------------------------------------------- 
c 
      include 'common_blocks.h'
c
ccccccccccccccc   READ IN DATA TO COMPUTE COST FN  ccccccccccccccccccc
c
      if(ntb.eq.1) then
         STlag=(800.d0-dz*nz)/wtrap
         open(unit=36,file=datdir//'AS_P.dat',status='old')
         open(unit=37,file=datdir//'AS_N.dat',status='old')
         open(unit=38,file=datdir//'AS_ST.dat',status='old')
         open(unit=39,file=datdir//'AS_Z.dat',status='old')
         open(unit=40,file=datdir//'AS_PrPr.dat',status='old')
         read(36,*)nPPdat
         read(37,*)nNdat
         read(38,*)nSTdat
         read(39,*)nZdat
         read(40,*)nPrPrdat
         do ndp=1,nPPdat
            read(36,*)dayChl(ndp),zChl(ndp),Chl(ndp)
         enddo
         do ndp=1,nNdat
            read(37,*)dayN(ndp),zN(ndp),aNut(ndp)            
         enddo
         do ndp=1,nSTdat        
            read(38,*)dayST(ndp),zST(ndp),timeST(ndp),ST(ndp)
            dayST(ndp)=dayST(ndp)-STlag
         enddo
         do ndp=1,nZdat        
            read(39,*)dayZ(ndp),zZ(ndp),aZoo(ndp)
         enddo
         do ndp=1,nPrPrdat  
            read(40,*)dayPrPr(ndp),zPrPr(ndp),aPrPr(ndp)
         enddo
      elseif(ntb.eq.2) then
         STlag=(880.d0-dz*nz)/wtrap
         open(unit=36,file=datdir//'EqPac_P.dat',status='old')
         open(unit=37,file=datdir//'EqPac_N.dat',status='old')
         open(unit=38,file=datdir//'EqPac_ST.dat',status='old')
         open(unit=39,file=datdir//'EqPac_Z.dat',status='old')
         open(unit=40,file=datdir//'EqPac_PrPr.dat',status='old')
         read(36,*)nPPdat
         read(37,*)nNdat
         read(38,*)nSTdat
         read(39,*)nZdat
         read(40,*)nPrPrdat
         do ndp=1,nPPdat
            read(36,*)dayChl(ndp),zChl(ndp),Chl(ndp)
         enddo
         do ndp=1,nNdat
            read(37,*)dayN(ndp),zN(ndp),aNut(ndp)            
         enddo
         do ndp=1,nSTdat        
            read(38,*)dayST(ndp),zST(ndp),timeST(ndp),ST(ndp)
            dayST(ndp)=dayST(ndp)-STlag
         enddo
         do ndp=1,nZdat        
            read(39,*)dayZ(ndp),zZ(ndp),smallZ,aZoo(ndp)
         enddo
         do ndp=1,nPrPrdat  
            read(40,*)dayPrPr(ndp),zPrPr(ndp),aPrPr(ndp)
         enddo
      endif
      close(36)
      close(37)
      close(38)
      close(39)
      close(40)
c
      return 
      end