subroutine readinit
	include 'common_blocks.h'

	real*8 depth,binitC,binitZ,binitN,binitD,binitF,binitS
     &          ,binitP,binitA

	if(ntb.eq.1) then
	   open (unit=30,status='old',file=datdir//'AS_Chl.init')
	   open (unit=31,status='old',file=datdir//'AS_Z.init')
	   open (unit=32,status='old',file=datdir//'AS_N.init')
	   open (unit=33,status='old',file=datdir//'AS_D.init')
	   open (unit=34,status='old',file=datdir//'AS_Fe.init')
	   open (unit=35,status='old',file=datdir//'AS_SiO4.init')
	   open (unit=36,status='old',file=datdir//'AS_PO4.init')
	elseif(ntb.eq.2) then
	   open (unit=30,status='old',file=datdir//'EqPac_Chl.init')
	   open (unit=31,status='old',file=datdir//'EqPac_Z.init')
	   open (unit=32,status='old',file=datdir//'EqPac_N.init')
	   open (unit=33,status='old',file=datdir//'EqPac_D.init')
	   open (unit=34,status='old',file=datdir//'EqPac_Fe.init')
	   open (unit=35,status='old',file=datdir//'EqPac_SiO4.init')
	   open (unit=36,status='old',file=datdir//'EqPac_PO4.init')
	   open (unit=37,status='old',file=datdir//'EqPac_A.init')
	endif
	do n=1,nz+1
	   read(30,*)depth,binitC
	   read(31,*)depth,binitZ
	   read(32,*)depth,binitN
	   read(33,*)depth,binitD
	   read(34,*)depth,binitF
	   read(35,*)depth,binitS
	   read(36,*)depth,binitP
	   if(ntb.eq.2)read(37,*)depth,binitA
c
cccc  User add initial conditions here ccccccccccccccccccccccccc
c
	   bio(n,1,1)=binitC/n2chl
	   bio(n,1,2)=binitZ
	   bio(n,1,3)=binitD
	   bio(n,1,4)=binitN
c
cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
c
	enddo
	
	close(30)
	close(31)
	close(32)
	close(33)
	close(34)
	close(35)
	close(36)
	close(37)
c
        return
        end