<p><b>ringler@lanl.gov</b> 2013-03-18 13:04:20 -0600 (Mon, 18 Mar 2013)</p><p><br>
fixing scaling factor<br>
</p><hr noshade><pre><font color="gray">Modified: branches/mpas_cam_coupling/src/core_hyd_atmos/module_time_integration.F
===================================================================
--- branches/mpas_cam_coupling/src/core_hyd_atmos/module_time_integration.F        2013-03-18 18:59:14 UTC (rev 2624)
+++ branches/mpas_cam_coupling/src/core_hyd_atmos/module_time_integration.F        2013-03-18 19:04:20 UTC (rev 2625)
@@ -746,9 +746,9 @@
 ! CSR
 
             ! tdr
-            ! (factor of 2.0 is to compensate for dt=dt/2 on the first RK step)
+            ! (factor of 3.0 is to compensate for dt=dt/3 on the first RK step)
                if(doEuler) then
-                  if(rk_step.eq.1) euler_tend_u(k,iEdge) = euler_tend_u(k,iEdge) + 2.0*u_diffusion
+                  if(rk_step.eq.1) euler_tend_u(k,iEdge) = euler_tend_u(k,iEdge) + 3.0*u_diffusion
                else
                   tend_u(k,iEdge) = tend_u(k,iEdge) + u_diffusion
                endif
@@ -834,9 +834,9 @@
                u_diffusion = meshScalingDel4(iEdge) * h_mom_eddy_visc4 * u_diffusion
 
                ! tdr
-               ! (factor of 2.0 is to compensate for dt=dt/2 on the first RK step)
+               ! (factor of 3.0 is to compensate for dt=dt/3 on the first RK step)
                if(doEuler) then
-                 if(rk_step.eq.1) euler_tend_u(k,iEdge) = euler_tend_u(k,iEdge) - 2.0*u_diffusion
+                 if(rk_step.eq.1) euler_tend_u(k,iEdge) = euler_tend_u(k,iEdge) - 3.0*u_diffusion
                else
                  tend_u(k,iEdge) = tend_u(k,iEdge) - u_diffusion
                endif
@@ -874,10 +874,10 @@
                zp = 0.5*(z3+z4)
     
                ! tdr
-               ! (factor of 2.0 is to compensate for dt=dt/2 on the first RK step)
+               ! (factor of 3.0 is to compensate for dt=dt/3 on the first RK step)
                if(doEuler) then
                    if(rk_step.eq.1)  &amp;
-                      euler_tend_u(k,iEdge) = euler_tend_u(k,iEdge) + 2.0*v_mom_eddy_visc2*(     &amp;
+                      euler_tend_u(k,iEdge) = euler_tend_u(k,iEdge) + 3.0*v_mom_eddy_visc2*(     &amp;
                                     (u(k+1,iEdge)-u(k  ,iEdge))/(zp-z0)                 &amp;
                                    -(u(k  ,iEdge)-u(k-1,iEdge))/(z0-zm) )/(0.5*(zp-zm))
                else

</font>
</pre>