subroutine readforcing
c
C	This subroutine reads in the forcing files: Horiz advection, PAR
c          vertical velocity, mixed layer depth, vertical diffusivity,
c          and temperature.
c
	include 'common_blocks.h'

	real swvel(500),day
c
	if(ntb.eq.1) then
	   open(unit=14,file=datdir//'AS_PAR.dat',status='old')
	   open(unit=15, file=datdir//'AS_mld.dat',status='old')
	   open(unit=16, file=datdir//'AS_w.dat',status='old')
	   open(unit=17,file=datdir//'AS_Kv.dat',status='old')
	   open(unit=18,file=datdir//'AS_T.dat',status='old')
	elseif(ntb.eq.2) then
	   open (unit=20,status='old',form='formatted',
     &       file=datdir//'EP_HA_Ni.dat')
	   open (unit=21,status='old',form='formatted',
     &       file=datdir//'EP_HA_Am.dat')
	   open (unit=22,status='old',form='formatted',
     &       file=datdir//'EP_HA_Fe.dat')
	   do ndz=1,nz
	      read(20,*)(HANi(ndz,ndt),ndt=1,15*30*24)
	      read(21,*)(HAAm(ndz,ndt),ndt=1,15*30*24)
	      read(22,*)(HAFe(ndz,ndt),ndt=1,15*30*24)
	   enddo	
	   open(unit=14,file=datdir//'EqPac_PAR.dat',status='old')
	   open(unit=15, file=datdir//'EqPac_mld.dat',status='old')
	   open(unit=16, file=datdir//'EqPac_w.dat',status='old')
	   open(unit=17,file=datdir//'EqPac_Kv.dat',status='old')
	   open(unit=18,file=datdir//'EqPac_T.dat',status='old')
	endif
c       
	read(16,*)(swvel(n),n=1,nz)
	read(15,*)day,adml
	ncz=0
	
	do jmain=1,isteps
	   read(14,*)cdays(jmain),qi(jmain)
	   qi(jmain) = qi(jmain)/0.4d0 !convert from PAR to total solar
c	   dml(jmain)=adml
	   dml(jmain)=dmin1(adml,110.d0)
	   read(17,*)day,(rkz(n,jmain),n=1,nz)
	   do n=1,nz
	      wvel(n,jmain)=-swvel(n)/86400.d0
	      if(n.gt.12)rkz(n,jmain)=1.e-4
	   enddo
	   do n=1,8
	      read(18,*)day,Tdat(jmain,n)
	   enddo
	   ncz=ncz+1
	   if(ncz.eq.24) then
	      ncz=0
	      read(15,*)day,adml
	      read(16,*)(swvel(n),n=1,nz)
	   endif
	enddo
	close(20)
	close(21)
	close(22)
	close(14)
	close(15)
	close(16)
	close(17)
	close(18)
c       
	return 
	end