<p><b>mhecht@lanl.gov</b> 2010-04-21 09:42:44 -0600 (Wed, 21 Apr 2010)</p><p>Added initialization of tracers 1-4 under test case 6, for testing of<br>
port of advection schemes from hyd_atmos core to sw (and ocean)<br>
core. Tracers 1 and 3 initialized as 1.0, tracers 2 and 4 are<br>
initialized with the value of initial longitude.<br>
</p><hr noshade><pre><font color="gray">Modified: branches/ocean_projects/port_adv_mwh/src/core_sw/module_test_cases.F
===================================================================
--- branches/ocean_projects/port_adv_mwh/src/core_sw/module_test_cases.F        2010-04-21 14:25:49 UTC (rev 193)
+++ branches/ocean_projects/port_adv_mwh/src/core_sw/module_test_cases.F        2010-04-21 15:42:44 UTC (rev 194)
@@ -415,6 +415,13 @@
       real (kind=RKIND), parameter :: w = 7.848e-6
       real (kind=RKIND), parameter :: K = 7.848e-6
       real (kind=RKIND), parameter :: R = 4.0
+      ! --- mwh ---
+      real (kind=RKIND), parameter :: rr = pii/9.0
+      real (kind=RKIND), parameter :: theta_c = pii/6.0
+      real (kind=RKIND), parameter :: lambda_c = 3.0*pii/2.0
+ 
+      real (kind=RKIND) :: rtemp
+      ! --- end ---
 
       integer :: iCell, iEdge, iVtx
       real (kind=RKIND) :: u, v
@@ -472,6 +479,21 @@
 #endif
       end do
 
+      ! --- mwh ---
+      !
+      ! Initialize tracer fields
+      !
+      do iCell=1,grid % nCells
+         state % tracers % array(1,1,iCell) = 1.0
+         state % tracers % array(3,1,iCell) = 1.0
+         state % tracers % array(2,1,iCell) = grid % lonCell % array(iCell) &amp;
+              * 180.0/3.14159
+         state % tracers % array(4,1,iCell) = grid % lonCell % array(iCell) &amp;
+              * 180.0/3.14159
+      end do
+
+      ! --- end ---
+
    end subroutine sw_test_case_6
 
 

</font>
</pre>