!
! f90 -o checkint checkint.o cutable.o -L$HOMES/lib -lcugen -lpcoco -lpfft \
! -lpbint -lpbend -lpmat -lpran -lpcomb -lpfig
!
!=============================================================================
program init_eqarea
!=============================================================================
! Test the migration method for generating the smooth cube grid
!=============================================================================
use cugen

implicit none
integer,parameter                :: n=40,norh=3,nq=2,ns=20,              &
                                    ntauh=10,ntau=ntauh*2,mtau=ntau,     &
                                    nouterit=40,ninnerit=80
!=============================================================================

!goto 3 !!!!!!!!! skip the generation process and just read the data in
call cucogen(norh,n,ntau,mtau,nouterit,ninnerit)
call cuwrite('cube_gen')
!3 call curead('abcdefgh')
call fini_map_table

end program init_eqarea