Todo List

Subprogram cumax (v1, v2, v3, u, v, w, umax, uxmax)

This routine is sloppy, could reuse state from elsewhere

Subprogram dudxyz (du, u, rm1, sm1, tm1, jm1, imsh, isd)

why this loop?

Subprogram filterq (v, f, nx, nz, w1, w2, ft, if3d, dmax)

Is this supposed to be nx*ny*nz?

Subprogram geom::init_geom ()

remove un{x,y,z} for regular geometry case

Subprogram hmh_gmres (res, h1, h2, wt, iter)

move these allocations to where they are needed

check if these inits are nessesary

Do we need this?

Subprogram hsmg_routines::h1mg_setup ()

Is it nessesary to set h1 and h2 here? They aren't inited until later

Subprogram poisson::spectral_solve (u, rhs)

replace this with coarse grid dssum

Subprogram soln::init_soln ()

Is qtl ever non-zero for .not. iflomach?

Subprogram vrdsmsh ()

clean up ta,tb,qmask

Subprogram xyzrst (xrm1, yrm1, zrm1, xsm1, ysm1, zsm1, XTM1, YTM1, ZTM1, IFAXIS)

why this loop?