<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) &
- euler_tend_u(k,iEdge) = euler_tend_u(k,iEdge) + 2.0*v_mom_eddy_visc2*( &
+ euler_tend_u(k,iEdge) = euler_tend_u(k,iEdge) + 3.0*v_mom_eddy_visc2*( &
(u(k+1,iEdge)-u(k ,iEdge))/(zp-z0) &
-(u(k ,iEdge)-u(k-1,iEdge))/(z0-zm) )/(0.5*(zp-zm))
else
</font>
</pre>