<p><b>mpetersen@lanl.gov</b> 2011-02-01 14:00:50 -0700 (Tue, 01 Feb 2011)</p><p>BRANCH COMMIT in vert_adv_mrp.  Add some temporary code to run double gyre domain.<br>
</p><hr noshade><pre><font color="gray">Modified: branches/ocean_projects/vert_adv_mrp/src/core_ocean/module_mpas_core.F
===================================================================
--- branches/ocean_projects/vert_adv_mrp/src/core_ocean/module_mpas_core.F        2011-01-31 16:52:01 UTC (rev 710)
+++ branches/ocean_projects/vert_adv_mrp/src/core_ocean/module_mpas_core.F        2011-02-01 21:00:50 UTC (rev 711)
@@ -76,16 +76,47 @@
       type (block_type), intent(inout) :: block
       type (mesh_type), intent(inout) :: mesh
       real (kind=RKIND), intent(in) :: dt
-      integer :: i, iEdge, iCell
+      integer :: i, iEdge, iCell, k
    
    
       call compute_solve_diagnostics(dt, block % state % time_levs(1) % state, mesh)
-   
-      call rbfInterp_initialize(mesh)
-      call init_reconstruct(mesh)
-      call reconstruct(block % state % time_levs(1) % state, mesh)
+ 
+! mrp 110131 temp commented out:  
+!      call rbfInterp_initialize(mesh)
+!      call init_reconstruct(mesh)
+!      call reconstruct(block % state % time_levs(1) % state, mesh)
+! mrp 110131 temp commented out end
 
+      ! mrp 110131 testing only
+      block % mesh % u_src % array(:,:) = &amp;
+        block % mesh % u_src % array(:,:) *100 
+     print *, 'min,maxval(u_src) ',minval( block % mesh % u_src % array(:,:)),maxval( block % mesh % u_src % array(:,:))
 
+      do iCell=1,block % mesh % nCells
+         do k=1,mesh % nVertLevels
+            block % state % time_levs(1) % state % tracers % array( &amp;
+               1, k ,iCell) =  31-k
+            block % state % time_levs(1) % state % tracers % array( &amp;
+               2, k ,iCell) =  32 + k*0.32
+         enddo 
+            block % state % time_levs(1) % state % tracers % array( &amp;
+               3, :,iCell) =  1.
+            block % state % time_levs(1) % state % tracers % array( &amp;
+               4, :,iCell) =  1.
+
+            block % state % time_levs(1) % state % tracers % array( &amp;
+               3, 1:10,iCell) =  0.
+         do k=1,mesh % nVertLevels,2
+            block % state % time_levs(1) % state % tracers % array( &amp;
+               4, k ,iCell) =  0.
+         enddo 
+
+      end do
+
+      
+      ! mrp 110131 testing only end
+
+
       ! 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.
@@ -373,6 +404,10 @@
       nVertLevels = block % mesh % nVertLevels
       vertexDegree = block % mesh % vertexDegree
 
+! mrp 110131 for testing only
+         maxLevelCell(1:nCells) = nVertLevels
+! mrp 110131 for testing only end
+
       ! for z-grids, maxLevelCell should be in input state
       ! Isopycnal grid uses all vertical cells
       if (config_vert_grid_type.eq.'isopycnal') then

</font>
</pre>