diff --git a/classROMSX.html b/classROMSX.html
index 118a55a0..d52e46db 100644
--- a/classROMSX.html
+++ b/classROMSX.html
@@ -4891,21 +4891,21 @@ <h2 class="memtitle"><span class="permalink"><a href="#a45591c69e52ac45b762eaa78
 <div class="line"><a name="l00235"></a><span class="lineno">  235</span>&#160;    <a class="code" href="classROMSX.html#afc892a41057b28364dab80adc896c7cb">vec_diff2_temp</a>[lev].reset(<span class="keyword">new</span> MultiFab(ba,dm,1,IntVect(<a class="code" href="IndexDefines_8H.html#a124da9e47fb1b1f3debaf49368f42e38">NGROW</a>,<a class="code" href="IndexDefines_8H.html#a124da9e47fb1b1f3debaf49368f42e38">NGROW</a>,0))); <span class="comment">// harmonic diffusivity temperature</span></div>
 <div class="line"><a name="l00236"></a><span class="lineno">  236</span>&#160; </div>
 <div class="line"><a name="l00237"></a><span class="lineno">  237</span>&#160;    <span class="comment">// maybe TODO: clean up component indexing in prestep?</span></div>
-<div class="line"><a name="l00238"></a><span class="lineno">  238</span>&#160;    <a class="code" href="classROMSX.html#ab1343869d59f89fd4773bc7fa230aa73">vec_ru</a>[lev].reset(<span class="keyword">new</span> MultiFab(ba,dm,2,IntVect(<a class="code" href="IndexDefines_8H.html#a124da9e47fb1b1f3debaf49368f42e38">NGROW</a>,<a class="code" href="IndexDefines_8H.html#a124da9e47fb1b1f3debaf49368f42e38">NGROW</a>,<a class="code" href="IndexDefines_8H.html#a124da9e47fb1b1f3debaf49368f42e38">NGROW</a>))); <span class="comment">// RHS u (incl horizontal and vertical advection)</span></div>
-<div class="line"><a name="l00239"></a><span class="lineno">  239</span>&#160;    <a class="code" href="classROMSX.html#ab5041f0f370ab1ec4f7f7e16127c9039">vec_rv</a>[lev].reset(<span class="keyword">new</span> MultiFab(ba,dm,2,IntVect(<a class="code" href="IndexDefines_8H.html#a124da9e47fb1b1f3debaf49368f42e38">NGROW</a>,<a class="code" href="IndexDefines_8H.html#a124da9e47fb1b1f3debaf49368f42e38">NGROW</a>,<a class="code" href="IndexDefines_8H.html#a124da9e47fb1b1f3debaf49368f42e38">NGROW</a>))); <span class="comment">// RHS v</span></div>
-<div class="line"><a name="l00240"></a><span class="lineno">  240</span>&#160;    <a class="code" href="classROMSX.html#ab5b7ff78a167e59ae727b22094b39079">vec_rufrc</a>[lev].reset(<span class="keyword">new</span> MultiFab(ba2d,dm,2,IntVect(<a class="code" href="IndexDefines_8H.html#a124da9e47fb1b1f3debaf49368f42e38">NGROW</a>,<a class="code" href="IndexDefines_8H.html#a124da9e47fb1b1f3debaf49368f42e38">NGROW</a>,0))); <span class="comment">//2d, (incl advection terms and surface/bottom stresses, integral over the whole columnn, k=0)</span></div>
-<div class="line"><a name="l00241"></a><span class="lineno">  241</span>&#160;    <a class="code" href="classROMSX.html#a52a5f265a370595f128046c31b9f31d5">vec_rvfrc</a>[lev].reset(<span class="keyword">new</span> MultiFab(ba2d,dm,2,IntVect(<a class="code" href="IndexDefines_8H.html#a124da9e47fb1b1f3debaf49368f42e38">NGROW</a>,<a class="code" href="IndexDefines_8H.html#a124da9e47fb1b1f3debaf49368f42e38">NGROW</a>,0))); <span class="comment">//2d, same as above but v</span></div>
-<div class="line"><a name="l00242"></a><span class="lineno">  242</span>&#160;    <a class="code" href="classROMSX.html#a2082155131b8b7e0dec29068a92cece2">vec_sustr</a>[lev].reset(<span class="keyword">new</span> MultiFab(ba2d,dm,1,IntVect(<a class="code" href="IndexDefines_8H.html#a124da9e47fb1b1f3debaf49368f42e38">NGROW</a>,<a class="code" href="IndexDefines_8H.html#a124da9e47fb1b1f3debaf49368f42e38">NGROW</a>,0))); <span class="comment">//2d, surface stress</span></div>
-<div class="line"><a name="l00243"></a><span class="lineno">  243</span>&#160;    <a class="code" href="classROMSX.html#a44a44e4d82bc38d55daa03f281437113">vec_svstr</a>[lev].reset(<span class="keyword">new</span> MultiFab(ba2d,dm,1,IntVect(<a class="code" href="IndexDefines_8H.html#a124da9e47fb1b1f3debaf49368f42e38">NGROW</a>,<a class="code" href="IndexDefines_8H.html#a124da9e47fb1b1f3debaf49368f42e38">NGROW</a>,0))); <span class="comment">//2d</span></div>
+<div class="line"><a name="l00238"></a><span class="lineno">  238</span>&#160;    <a class="code" href="classROMSX.html#ab1343869d59f89fd4773bc7fa230aa73">vec_ru</a>[lev].reset(<span class="keyword">new</span> MultiFab(convert(ba,IntVect(1,0,0)),dm,2,IntVect(<a class="code" href="IndexDefines_8H.html#a124da9e47fb1b1f3debaf49368f42e38">NGROW</a>,<a class="code" href="IndexDefines_8H.html#a124da9e47fb1b1f3debaf49368f42e38">NGROW</a>,<a class="code" href="IndexDefines_8H.html#a124da9e47fb1b1f3debaf49368f42e38">NGROW</a>))); <span class="comment">// RHS u (incl horizontal and vertical advection)</span></div>
+<div class="line"><a name="l00239"></a><span class="lineno">  239</span>&#160;    <a class="code" href="classROMSX.html#ab5041f0f370ab1ec4f7f7e16127c9039">vec_rv</a>[lev].reset(<span class="keyword">new</span> MultiFab(convert(ba,IntVect(0,1,0)),dm,2,IntVect(<a class="code" href="IndexDefines_8H.html#a124da9e47fb1b1f3debaf49368f42e38">NGROW</a>,<a class="code" href="IndexDefines_8H.html#a124da9e47fb1b1f3debaf49368f42e38">NGROW</a>,<a class="code" href="IndexDefines_8H.html#a124da9e47fb1b1f3debaf49368f42e38">NGROW</a>))); <span class="comment">// RHS v</span></div>
+<div class="line"><a name="l00240"></a><span class="lineno">  240</span>&#160;    <a class="code" href="classROMSX.html#ab5b7ff78a167e59ae727b22094b39079">vec_rufrc</a>[lev].reset(<span class="keyword">new</span> MultiFab(convert(ba2d,IntVect(1,0,0)),dm,2,IntVect(<a class="code" href="IndexDefines_8H.html#a124da9e47fb1b1f3debaf49368f42e38">NGROW</a>,<a class="code" href="IndexDefines_8H.html#a124da9e47fb1b1f3debaf49368f42e38">NGROW</a>,0))); <span class="comment">//2d, (incl advection terms and surface/bottom stresses, integral over the whole columnn, k=0)</span></div>
+<div class="line"><a name="l00241"></a><span class="lineno">  241</span>&#160;    <a class="code" href="classROMSX.html#a52a5f265a370595f128046c31b9f31d5">vec_rvfrc</a>[lev].reset(<span class="keyword">new</span> MultiFab(convert(ba2d,IntVect(0,1,0)),dm,2,IntVect(<a class="code" href="IndexDefines_8H.html#a124da9e47fb1b1f3debaf49368f42e38">NGROW</a>,<a class="code" href="IndexDefines_8H.html#a124da9e47fb1b1f3debaf49368f42e38">NGROW</a>,0))); <span class="comment">//2d, same as above but v</span></div>
+<div class="line"><a name="l00242"></a><span class="lineno">  242</span>&#160;    <a class="code" href="classROMSX.html#a2082155131b8b7e0dec29068a92cece2">vec_sustr</a>[lev].reset(<span class="keyword">new</span> MultiFab(convert(ba2d,IntVect(1,0,0)),dm,1,IntVect(<a class="code" href="IndexDefines_8H.html#a124da9e47fb1b1f3debaf49368f42e38">NGROW</a>,<a class="code" href="IndexDefines_8H.html#a124da9e47fb1b1f3debaf49368f42e38">NGROW</a>,0))); <span class="comment">//2d, surface stress</span></div>
+<div class="line"><a name="l00243"></a><span class="lineno">  243</span>&#160;    <a class="code" href="classROMSX.html#a44a44e4d82bc38d55daa03f281437113">vec_svstr</a>[lev].reset(<span class="keyword">new</span> MultiFab(convert(ba2d,IntVect(0,1,0)),dm,1,IntVect(<a class="code" href="IndexDefines_8H.html#a124da9e47fb1b1f3debaf49368f42e38">NGROW</a>,<a class="code" href="IndexDefines_8H.html#a124da9e47fb1b1f3debaf49368f42e38">NGROW</a>,0))); <span class="comment">//2d</span></div>
 <div class="line"><a name="l00244"></a><span class="lineno">  244</span>&#160;    <a class="code" href="classROMSX.html#ac49463f0fb789402ba8641f289906456">vec_rdrag</a>[lev].reset(<span class="keyword">new</span> MultiFab(ba2d,dm,1,IntVect(<a class="code" href="IndexDefines_8H.html#a124da9e47fb1b1f3debaf49368f42e38">NGROW</a>,<a class="code" href="IndexDefines_8H.html#a124da9e47fb1b1f3debaf49368f42e38">NGROW</a>,0))); <span class="comment">//2d, linear drag coefficient [m/s], defined at rho, somehow related to rdrg</span></div>
-<div class="line"><a name="l00245"></a><span class="lineno">  245</span>&#160;    <a class="code" href="classROMSX.html#a607e27d43f6eaebcf557c4ae17b8039a">vec_bustr</a>[lev].reset(<span class="keyword">new</span> MultiFab(ba2d,dm,1,IntVect(<a class="code" href="IndexDefines_8H.html#a124da9e47fb1b1f3debaf49368f42e38">NGROW</a>,<a class="code" href="IndexDefines_8H.html#a124da9e47fb1b1f3debaf49368f42e38">NGROW</a>,0))); <span class="comment">//2d, bottom stress</span></div>
-<div class="line"><a name="l00246"></a><span class="lineno">  246</span>&#160;    <a class="code" href="classROMSX.html#a5a981dafaf65d2d8d96f16e1e7467eb2">vec_bvstr</a>[lev].reset(<span class="keyword">new</span> MultiFab(ba2d,dm,1,IntVect(<a class="code" href="IndexDefines_8H.html#a124da9e47fb1b1f3debaf49368f42e38">NGROW</a>,<a class="code" href="IndexDefines_8H.html#a124da9e47fb1b1f3debaf49368f42e38">NGROW</a>,0)));</div>
+<div class="line"><a name="l00245"></a><span class="lineno">  245</span>&#160;    <a class="code" href="classROMSX.html#a607e27d43f6eaebcf557c4ae17b8039a">vec_bustr</a>[lev].reset(<span class="keyword">new</span> MultiFab(convert(ba2d,IntVect(1,0,0)),dm,1,IntVect(<a class="code" href="IndexDefines_8H.html#a124da9e47fb1b1f3debaf49368f42e38">NGROW</a>,<a class="code" href="IndexDefines_8H.html#a124da9e47fb1b1f3debaf49368f42e38">NGROW</a>,0))); <span class="comment">//2d, bottom stress</span></div>
+<div class="line"><a name="l00246"></a><span class="lineno">  246</span>&#160;    <a class="code" href="classROMSX.html#a5a981dafaf65d2d8d96f16e1e7467eb2">vec_bvstr</a>[lev].reset(<span class="keyword">new</span> MultiFab(convert(ba2d,IntVect(0,1,0)),dm,1,IntVect(<a class="code" href="IndexDefines_8H.html#a124da9e47fb1b1f3debaf49368f42e38">NGROW</a>,<a class="code" href="IndexDefines_8H.html#a124da9e47fb1b1f3debaf49368f42e38">NGROW</a>,0)));</div>
 <div class="line"><a name="l00247"></a><span class="lineno">  247</span>&#160; </div>
 <div class="line"><a name="l00248"></a><span class="lineno">  248</span>&#160;    <span class="comment">//all 2d -- all associated with the 2D advance</span></div>
-<div class="line"><a name="l00249"></a><span class="lineno">  249</span>&#160;    <a class="code" href="classROMSX.html#a18212e9657c0655cd29e0c6241c1224f">vec_DU_avg1</a>[lev].reset(<span class="keyword">new</span> MultiFab(ba2d,dm,1,IntVect(<a class="code" href="IndexDefines_8H.html#a124da9e47fb1b1f3debaf49368f42e38">NGROW</a>,<a class="code" href="IndexDefines_8H.html#a124da9e47fb1b1f3debaf49368f42e38">NGROW</a>,0))); <span class="comment">//2d DU: sum(height[incl free surface?] * u)</span></div>
-<div class="line"><a name="l00250"></a><span class="lineno">  250</span>&#160;    <a class="code" href="classROMSX.html#a883756f67229729f133b9bdf784c21ff">vec_DU_avg2</a>[lev].reset(<span class="keyword">new</span> MultiFab(ba2d,dm,1,IntVect(<a class="code" href="IndexDefines_8H.html#a124da9e47fb1b1f3debaf49368f42e38">NGROW</a>,<a class="code" href="IndexDefines_8H.html#a124da9e47fb1b1f3debaf49368f42e38">NGROW</a>,0))); <span class="comment">//2d like above, but correct(or)?</span></div>
-<div class="line"><a name="l00251"></a><span class="lineno">  251</span>&#160;    <a class="code" href="classROMSX.html#aa9738d90727df4f8aba2afb0ffed1e8f">vec_DV_avg1</a>[lev].reset(<span class="keyword">new</span> MultiFab(ba2d,dm,1,IntVect(<a class="code" href="IndexDefines_8H.html#a124da9e47fb1b1f3debaf49368f42e38">NGROW</a>,<a class="code" href="IndexDefines_8H.html#a124da9e47fb1b1f3debaf49368f42e38">NGROW</a>,0)));</div>
-<div class="line"><a name="l00252"></a><span class="lineno">  252</span>&#160;    <a class="code" href="classROMSX.html#a6847259dab9f6118da780bf623a9868a">vec_DV_avg2</a>[lev].reset(<span class="keyword">new</span> MultiFab(ba2d,dm,1,IntVect(<a class="code" href="IndexDefines_8H.html#a124da9e47fb1b1f3debaf49368f42e38">NGROW</a>,<a class="code" href="IndexDefines_8H.html#a124da9e47fb1b1f3debaf49368f42e38">NGROW</a>,0)));</div>
+<div class="line"><a name="l00249"></a><span class="lineno">  249</span>&#160;    <a class="code" href="classROMSX.html#a18212e9657c0655cd29e0c6241c1224f">vec_DU_avg1</a>[lev].reset(<span class="keyword">new</span> MultiFab(convert(ba2d,IntVect(1,0,0)),dm,1,IntVect(<a class="code" href="IndexDefines_8H.html#a124da9e47fb1b1f3debaf49368f42e38">NGROW</a>,<a class="code" href="IndexDefines_8H.html#a124da9e47fb1b1f3debaf49368f42e38">NGROW</a>,0))); <span class="comment">//2d DU: sum(height[incl free surface?] * u)</span></div>
+<div class="line"><a name="l00250"></a><span class="lineno">  250</span>&#160;    <a class="code" href="classROMSX.html#a883756f67229729f133b9bdf784c21ff">vec_DU_avg2</a>[lev].reset(<span class="keyword">new</span> MultiFab(convert(ba2d,IntVect(1,0,0)),dm,1,IntVect(<a class="code" href="IndexDefines_8H.html#a124da9e47fb1b1f3debaf49368f42e38">NGROW</a>,<a class="code" href="IndexDefines_8H.html#a124da9e47fb1b1f3debaf49368f42e38">NGROW</a>,0))); <span class="comment">//2d like above, but correct(or)?</span></div>
+<div class="line"><a name="l00251"></a><span class="lineno">  251</span>&#160;    <a class="code" href="classROMSX.html#aa9738d90727df4f8aba2afb0ffed1e8f">vec_DV_avg1</a>[lev].reset(<span class="keyword">new</span> MultiFab(convert(ba2d,IntVect(0,1,0)),dm,1,IntVect(<a class="code" href="IndexDefines_8H.html#a124da9e47fb1b1f3debaf49368f42e38">NGROW</a>,<a class="code" href="IndexDefines_8H.html#a124da9e47fb1b1f3debaf49368f42e38">NGROW</a>,0)));</div>
+<div class="line"><a name="l00252"></a><span class="lineno">  252</span>&#160;    <a class="code" href="classROMSX.html#a6847259dab9f6118da780bf623a9868a">vec_DV_avg2</a>[lev].reset(<span class="keyword">new</span> MultiFab(convert(ba2d,IntVect(0,1,0)),dm,1,IntVect(<a class="code" href="IndexDefines_8H.html#a124da9e47fb1b1f3debaf49368f42e38">NGROW</a>,<a class="code" href="IndexDefines_8H.html#a124da9e47fb1b1f3debaf49368f42e38">NGROW</a>,0)));</div>
 <div class="line"><a name="l00253"></a><span class="lineno">  253</span>&#160; </div>
 <div class="line"><a name="l00254"></a><span class="lineno">  254</span>&#160;    <a class="code" href="classROMSX.html#ab5d263be53bfdd397f6c64672ee77009">vec_rubar</a>[lev].reset(<span class="keyword">new</span> MultiFab(convert(ba2d,IntVect(1,0,0)),dm,4,IntVect(<a class="code" href="IndexDefines_8H.html#a124da9e47fb1b1f3debaf49368f42e38">NGROW</a>,<a class="code" href="IndexDefines_8H.html#a124da9e47fb1b1f3debaf49368f42e38">NGROW</a>,0))); <span class="comment">// 2d RHS ubar</span></div>
 <div class="line"><a name="l00255"></a><span class="lineno">  255</span>&#160;    <a class="code" href="classROMSX.html#a83d781ead821f8ed29f6f2b9ef70a9aa">vec_rvbar</a>[lev].reset(<span class="keyword">new</span> MultiFab(convert(ba2d,IntVect(0,1,0)),dm,4,IntVect(<a class="code" href="IndexDefines_8H.html#a124da9e47fb1b1f3debaf49368f42e38">NGROW</a>,<a class="code" href="IndexDefines_8H.html#a124da9e47fb1b1f3debaf49368f42e38">NGROW</a>,0)));</div>