program bc6

	implicit none
	
	! Declare variables
	
	integer :: ibcfcb
 	real :: hdry
	integer :: iwdflg
	real :: wetfct
	integer :: iwetit
	integer :: ihdwet
	
	integer :: ltype1, ltype2, ltype3, ltype4
	
	real :: trpy
	real,dimension(10) :: tran1, tran2, tran3, tran4
	real :: cnstntr
	integer :: iprn
	integer :: i
	real :: vcont
	
	! Define variables
	
	ibcfcb = 0 ! Flag for cell-by-cell flow written or not, =0 then not written
	hdry = 999.00 ! Head assigned to cells coverted to dry, similar to hnoflo in basic package
	iwdflg = 0 ! Flag for wetting capability, =0 then inactive
	wetfct = 0. ! Factor for when cell is converted from dry to wet, set to zero...
	iwetit = 0 ! Iteration variable, again, for this dry to wet cell thing, set to zero...
	ihdwet = 0 ! Flag for this wetting nightmare, set to zero!!!
	
	ltype1 = 10 ! Layer type and method for inter-block conductance, 0 = harmonic mean, 1 = arithmetic mean, 0= confined, 1 = unconfined
	ltype2 = 10 ! Layer type and method for inter-block conductance, 0 = harmonic mean, 1 = arithmetic mean, 0= confined, 1 = unconfined
	ltype3 = 10
	ltype4 = 10

	trpy = 1.0 ! Horizontal anisotopy for each layer, set to 1.0 for isotropic
	
	tran1(:)=1.0  ! Transmissivity along rows, low K, hom
	tran2(:)=10.0  ! Transmissivity along rows, high K, hom
	
	do i=1,10
		tran3(i)=i ! Transmissivity along rows, het K low to high
	end do

	do i=1,10
		tran4(i)=20-i ! Transmissivity along rows, het K high to low
	end do
	
	cnstntr = 1. ! Real-number constant that is a multiplier of the array being read
	iprn = 0 ! Flag that indicates if the array being read should be printed, if less than zero, array not printed

	vcont = 1.0e-10

	! Execute program
	
	! Output results
	
	write(*,1) ibcfcb, hdry, iwdflg, wetfct, iwetit, ihdwet
	1 format(i10,f10.2,i10,f10.2,2i10,"     ibcfcb, hdry, iwdflg, wetfct, iwetit, ihdwet")
	
	write(*,2) ltype1, ltype2, ltype3, ltype4 
	2 format(4i2)
	
	write(*,3) trpy
	3 format("CONSTANT",f10.2,"     trpy")
	
	write(*,4) cnstntr, iprn, (tran1(i),i=1,10), vcont
	4 format("INTERNAL",f10.2,"   (10f10.2)",i5,/,10f10.2,"     tran1",/,&
	"CONSTANT",es10.2,"     vcont1")
	
	write(*,5) cnstntr, iprn, (tran2(i),i=1,10), vcont
	5 format("INTERNAL",f10.2,"   (10f10.2)",i5,/,10f10.2,"     tran2",/,&
	"CONSTANT",es10.2,"     vcont2")
	
	write(*,6) cnstntr, iprn, (tran3(i),i=1,10), vcont
	6 format("INTERNAL",f10.2,"   (10f10.2)",i5,/,10f10.2,"     tran3",/,&
	"CONSTANT",es10.2,"     vcont3")
	
	write(*,7) cnstntr, iprn, (tran4(i),i=1,10)
	7 format("INTERNAL",f10.2,"   (10f10.2)",i5,/,10f10.2,"     tran4")
	
end program bc6