program gridtest,11
use kvar; use kcnst; use kphysmod 
implicit none
integer, parameter :: ix=61
integer ii
real(mrl) F_tk
real, allocatable, dimension(:,:) :: za,xa,dummy
print *,'kickoff...'
call kickoff
open(unit=10,file='../play/grid.dat',status='unknown')
open(unit=20,file='../play/grids.dat',status='unknown')
print *,'will read...'
call read_namelist_file 
print *, 'will call allocate_arrays'
call allocate_arrays 
print *, 'did call allocate_arrays'
allocate(xa(ix,nf),za(ix,nf),dummy(ix,nf))
print *, "allocated"
dummy=0.
call repoint 
do ii=1,ix
	write(*,'(i3)',advance='no') ii
	hite=500+2000*(ii-1.)/(ix-1.)
	zbottom=500+500*cos(pi*(ii-1.)/(ix-1.))
	call nocturnal_grid
	zfdef=zf
	zfnext=zfdef
	call z_coords 
!	blgsq=1.
	if (gridadapt.ge.1) then
!	time=2.
		call blgrid_prep
		zf=zfnext
	end if
	do j=1,nf
		xa(ii,j)=ii-1
		za(ii,j)=zf(j)
    	dummy(ii,j)=1.9*F_tk(.01*(hite-za(ii,j)))-.95
	end do
end do
do n=1,nf
	write (10,*) n,zfdef(n),zf(n)
end do

write (20,*) "$in=",ix
write (20,*) "$jn=",nf
write (20,*) "$ixy=0"
write (20,'(e14.6)') dummy,xa,za
write (*,*) ' success'
end program gridtest