<p><b>ringler@lanl.gov</b> 2012-04-27 15:09:34 -0600 (Fri, 27 Apr 2012)</p><p><br>
add RBF reconstruction of wind stress field.<br>
<br>
this could should probably be removed before pushing onto trunk.<br>
</p><hr noshade><pre><font color="gray">Modified: branches/ocean_projects/monthly_forcing/src/core_ocean/Registry
===================================================================
--- branches/ocean_projects/monthly_forcing/src/core_ocean/Registry        2012-04-27 20:47:21 UTC (rev 1831)
+++ branches/ocean_projects/monthly_forcing/src/core_ocean/Registry        2012-04-27 21:09:34 UTC (rev 1832)
@@ -265,6 +265,11 @@
 var persistent real    uReconstructZ ( nVertLevels nCells Time ) 2 - uReconstructZ state - -
 var persistent real    uReconstructZonal ( nVertLevels nCells Time ) 2 o uReconstructZonal state - -
 var persistent real    uReconstructMeridional ( nVertLevels nCells Time ) 2 o uReconstructMeridional state - -
+var persistent real    uSrcReconstructX ( nVertLevels nCells Time ) 2 - uSrcReconstructX state - -
+var persistent real    uSrcReconstructY ( nVertLevels nCells Time ) 2 - uSrcReconstructY state - -
+var persistent real    uSrcReconstructZ ( nVertLevels nCells Time ) 2 - uSrcReconstructZ state - -
+var persistent real    uSrcReconstructZonal ( nVertLevels nCells Time ) 2 o uSrcReconstructZonal state - -
+var persistent real    uSrcReconstructMeridional ( nVertLevels nCells Time ) 2 o uSrcReconstructMeridional state - -
 var persistent real    MontPot ( nVertLevels nCells Time ) 2 - MontPot state - -
 var persistent real    pressure ( nVertLevels nCells Time ) 2 - pressure state - -
 var persistent real    wTop ( nVertLevelsP1 nCells Time ) 2 - wTop state - -

Modified: branches/ocean_projects/monthly_forcing/src/core_ocean/mpas_ocn_mpas_core.F
===================================================================
--- branches/ocean_projects/monthly_forcing/src/core_ocean/mpas_ocn_mpas_core.F        2012-04-27 20:47:21 UTC (rev 1831)
+++ branches/ocean_projects/monthly_forcing/src/core_ocean/mpas_ocn_mpas_core.F        2012-04-27 21:09:34 UTC (rev 1832)
@@ -278,6 +278,16 @@
                        block % state % time_levs(1) % state % uReconstructMeridional % array    &amp;
                       )
 
+!TDR
+      call mpas_reconstruct(mesh, mesh % u_src % array,                  &amp;
+                       block % state % time_levs(1) % state % uSrcReconstructX % array,            &amp;
+                       block % state % time_levs(1) % state % uSrcReconstructY % array,            &amp;
+                       block % state % time_levs(1) % state % uSrcReconstructZ % array,            &amp;
+                       block % state % time_levs(1) % state % uSrcReconstructZonal % array,        &amp;
+                       block % state % time_levs(1) % state % uSrcReconstructMeridional % array    &amp;
+                      )
+!TDR
+
       ! initialize velocities and tracers on land to be -1e34
       ! The reconstructed velocity on land will have values not exactly
       ! -1e34 due to the interpolation of reconstruction.

Modified: branches/ocean_projects/monthly_forcing/src/core_ocean/mpas_ocn_time_integration_rk4.F
===================================================================
--- branches/ocean_projects/monthly_forcing/src/core_ocean/mpas_ocn_time_integration_rk4.F        2012-04-27 20:47:21 UTC (rev 1831)
+++ branches/ocean_projects/monthly_forcing/src/core_ocean/mpas_ocn_time_integration_rk4.F        2012-04-27 21:09:34 UTC (rev 1832)
@@ -343,6 +343,16 @@
                           block % state % time_levs(2) % state % uReconstructMeridional % array    &amp;
                          )
 
+!TDR
+         call mpas_reconstruct(block % mesh, block % mesh % u_src % array,          &amp;
+                          block % state % time_levs(2) % state % uSrcReconstructX % array,            &amp;
+                          block % state % time_levs(2) % state % uSrcReconstructY % array,            &amp;
+                          block % state % time_levs(2) % state % uSrcReconstructZ % array,            &amp;
+                          block % state % time_levs(2) % state % uSrcReconstructZonal % array,        &amp;
+                          block % state % time_levs(2) % state % uSrcReconstructMeridional % array    &amp;
+                         )
+!TDR
+
          call ocn_time_average_accumulate(block % state % time_levs(2) % state, block % state % time_levs(1) % state)
 
          block =&gt; block % next

Modified: branches/ocean_projects/monthly_forcing/src/core_ocean/mpas_ocn_time_integration_split.F
===================================================================
--- branches/ocean_projects/monthly_forcing/src/core_ocean/mpas_ocn_time_integration_split.F        2012-04-27 20:47:21 UTC (rev 1831)
+++ branches/ocean_projects/monthly_forcing/src/core_ocean/mpas_ocn_time_integration_split.F        2012-04-27 21:09:34 UTC (rev 1832)
@@ -869,6 +869,7 @@
          end if
 
          call ocn_diagnostic_solve(dt, block % state % time_levs(2) % state, block % mesh)
+
          call mpas_reconstruct(block % mesh, block % state % time_levs(2) % state % u % array,          &amp;
             block % state % time_levs(2) % state % uReconstructX % array,            &amp;
             block % state % time_levs(2) % state % uReconstructY % array,            &amp;
@@ -876,6 +877,17 @@
             block % state % time_levs(2) % state % uReconstructZonal % array,        &amp;
             block % state % time_levs(2) % state % uReconstructMeridional % array)
 
+!TDR
+         call mpas_reconstruct(block % mesh, block % mesh % u_src % array,          &amp;
+                          block % state % time_levs(2) % state % uSrcReconstructX % array,            &amp;
+                          block % state % time_levs(2) % state % uSrcReconstructY % array,            &amp;
+                          block % state % time_levs(2) % state % uSrcReconstructZ % array,            &amp;
+                          block % state % time_levs(2) % state % uSrcReconstructZonal % array,        &amp;
+                          block % state % time_levs(2) % state % uSrcReconstructMeridional % array    &amp;
+                         )
+!TDR
+
+
          call ocn_time_average_accumulate(block % state % time_levs(2) % state, block % state % time_levs(1) % state)
 
 

</font>
</pre>