<p><b>ringler@lanl.gov</b> 2011-02-08 13:38:50 -0700 (Tue, 08 Feb 2011)</p><p><br>
minor comments by reviewer for consideration (grep "TDR")<F9><br>
</p><hr noshade><pre><font color="gray">Modified: branches/ocean_projects/vert_adv_mrp/src/core_ocean/module_time_integration.F
===================================================================
--- branches/ocean_projects/vert_adv_mrp/src/core_ocean/module_time_integration.F 2011-02-08 19:53:07 UTC (rev 718)
+++ branches/ocean_projects/vert_adv_mrp/src/core_ocean/module_time_integration.F 2011-02-08 20:38:50 UTC (rev 719)
@@ -897,6 +897,12 @@
allocate(tracerTop(num_tracers,nVertLevels+1,nCells))
+!TDR: shouldn't we put tracer(:,1,:) into tracerTop(:,1,:)?
+!TDR: and put tracer(:,maxLevelCell(iCell),:) into tracerTop(:,maxLevelCell(iCell)+1,:)?
+!TDR: regardless of what we put, these layer interfaces should be no flux because w==0
+!TDR: still, this is a valid approximation at those interfaces and better than
+!TDR: what will be grabbed from uninitialized memory.
+
if (config_vert_tracer_adv.eq.'stencil'.and. &
config_vert_tracer_adv_order.eq.2) then
</font>
</pre>