initial stages of new dtopo method. This includes modifications to the topo module creating new topo arrays at the spacing of the dtopo arrays. topo values are pulled from the finest topo files availabe for the region. The newer method of moving the topography has not been added to this commit...so aux arrays are still set from movetopo. This is a test of only the addition of new topo arrays from dtopo. Note: dtopo module has been added to the topo module...so all use dtopo have been removed. setprob.f must call dtopo before topo
fixed bug in topo module regarding forgotten (i)
added topo0work for saving original finest topography files covering dtopo region
added update_topo and fixed indice bug in topomodule. Features are complete at this point...however, code compiles and runs but does not generate waves given a dtopo uplift. entering debug phase 1
minor bug fix in b4step regarding logic in call to topo_update. still not moving topography correctly.
outline for new recursive rectintegral
version of cellgridintegrate to modify
rectintegral2 and cellgridintegrate2 with lots of debugging statements -- seems to be working
fixed rectintegral2 to check area before some recursive calls, seemed to fix but with 6 topo files.  Debug statements still in
test of new topo routines with 6 nested levels of topo grids
removed debug statements from rectintegral2 and cellgridintegrate2
call rectintegral starting with topo grid m not 1
added comment to unused file
Added check that topo arrays cover full domain using new recursive routine topoarea in topo_module.f90
stop with error message if grid cell doesn't overlap any topo grid
print more info if stopping with ntogo error
Skip setting topo in ghost cells outside physical domain to avoid errors if topo files don't extend beyond domain.  This should be tested more fully.
merged rjleveque/rectintegral into new_dtopo from @dlgeorge
merged changes from master
still debugging dtopo