dumpage zero out smaller than e-99
changed pressure initialization...rise to minimum for failure and have second time interval for over-failure-pressure..ignore dilatancy for determining minimum failure pressure
fixed bug in b4step2 b/(h+b), not b/h+b
Merge branch 'master' into procR_prechanges
fixed dumpgage for values less than e-99 in topo loop
added 1d fill between from 2d data to python plotting. No changes to running code
returning in qad...no change in test case
cleaned up mistakes in b4step 2 regarding rising pressure
changing to new riemann solver
changes only to matlab plotting routines...no changes to any running code
bug catch for looping in setaux and qinit. added m=0 capability in qinit...should not effect results significantly
minor changes to movetopo that should fix small bug...no move topo for natrelease...should not make a difference
fixed bug of rho in basal fluid resistance...fixed factor of 2 in shear induced dilatancy in src routine
changes to admissableq regarding m=0 and high pressure
realized an incompatibility regarding m=0 and dividing by zero...fixed with new kperm...but for testing, setting mlo=1.e-16
changing to exponential kperm...still old viscous number
resorting mlo to 1.d-16
restoring to old kperm
back to new kperm
restoring factor of 2 in D