<p><b>mpetersen@lanl.gov</b> 2012-05-23 08:50:23 -0600 (Wed, 23 May 2012)</p><p>BRANCH COMMIT: The transport_section branch is now complete. Because of an error in reading in 2D non-partitioned netcdf fields, this code currently reads in simple text for the variables sectionEdgeIndex, sectionEdgeSign as a temporary work-around. The text files are created by the matlab code in mpas_repo/branches/tools/transport_sections. The transport through each section, averaged over every timestep and summed over all edged and levels in each section, will then appear in the output file as the variable acc_uTransportSection.<br>
</p><hr noshade><pre><font color="gray">Modified: branches/ocean_projects/transport_sections/src/core_ocean/mpas_ocn_time_average.F
===================================================================
--- branches/ocean_projects/transport_sections/src/core_ocean/mpas_ocn_time_average.F 2012-05-23 14:44:20 UTC (rev 1931)
+++ branches/ocean_projects/transport_sections/src/core_ocean/mpas_ocn_time_average.F 2012-05-23 14:50:23 UTC (rev 1932)
@@ -140,6 +140,11 @@
acc_uTransportSection => state % acc_uTransportSection % array
old_acc_uTransportSection => old_state % acc_uTransportSection % array
+
+ ! Sum up transport over all edges in this processors partition,
+ ! for each section. Here we use h_edge for best precision.
+ ! To match with the matlab postprocessing code transport_sections.m,
+ ! replace h_edge(k,iEdge) with hZLevel(k).
do iSection = 1,nSections
tempSum = 0.0
do iLocalEdge = 1,nLocalEdgesInSection(iSection)
</font>
</pre>