<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 &
- + 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 &
+ + 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("calc. new prognostic vars")
</font>
</pre>