<p><b>mhoffman@lanl.gov</b> 2012-06-15 11:22:19 -0600 (Fri, 15 Jun 2012)</p><p>BRANCH COMMIT -- land ice<br>
<br>
Making sure that when tracer advection is turned off, tracer values do not change at all, even during ice advance.<br>
</p><hr noshade><pre><font color="gray">Modified: branches/land_ice_projects/implement_core/src/core_land_ice/mpas_land_ice_time_integration_forwardeuler.F
===================================================================
--- branches/land_ice_projects/implement_core/src/core_land_ice/mpas_land_ice_time_integration_forwardeuler.F        2012-06-15 16:53:51 UTC (rev 1986)
+++ branches/land_ice_projects/implement_core/src/core_land_ice/mpas_land_ice_time_integration_forwardeuler.F        2012-06-15 17:22:19 UTC (rev 1987)
@@ -199,15 +199,17 @@
 
 
          ! Calculate new tracer values =================
-         do iTracer = 1, size(tracersNew, 1)
-             where (layerThicknessNew .gt. 0.0)
-                 tracersNew(iTracer,:,:) = (tracersOld(iTracer,:,:) * layerThicknessOld &amp;
-                     + tracer_tendency(iTracer,:,:) * dt / SecondsInYear)  /  (layerThicknessNew)
-             elsewhere
-                 ! May or may not want to assign tracer values to non-ice cells
-                 tracersNew(iTracer,:,:) = 0.0
-             end where
-         end do
+         if (config_tracer_advection .ne. 'None') then
+           do iTracer = 1, size(tracersNew, 1)
+               where (layerThicknessNew .gt. 0.0)
+                   tracersNew(iTracer,:,:) = (tracersOld(iTracer,:,:) * layerThicknessOld &amp;
+                       + tracer_tendency(iTracer,:,:) * dt / SecondsInYear)  /  (layerThicknessNew)
+               elsewhere
+                   ! May or may not want to assign tracer values to non-ice cells
+                   tracersNew(iTracer,:,:) = 0.0
+               end where
+           end do
+         endif
 
          call mpas_timer_stop(&quot;calc. new prognostic vars&quot;)
 

</font>
</pre>