<p><b>duda</b> 2010-07-15 12:01:57 -0600 (Thu, 15 Jul 2010)</p><p>BRANCH COMMIT<br>
<br>
In CAM_TO_MPAS(), halo exchange should be for scalars_phys_tend rather than scalars.<br>
Also, name the scalar tendency array scalars_tend just to be consistent with naming <br>
of u_tend and theta_tend.<br>
<br>
<br>
M module_mpas_cam_interface.F<br>
</p><hr noshade><pre><font color="gray">Modified: branches/mpas_cam_coupling/src/driver_cam_interface/module_mpas_cam_interface.F
===================================================================
--- branches/mpas_cam_coupling/src/driver_cam_interface/module_mpas_cam_interface.F        2010-07-14 20:17:36 UTC (rev 376)
+++ branches/mpas_cam_coupling/src/driver_cam_interface/module_mpas_cam_interface.F        2010-07-15 18:01:57 UTC (rev 377)
@@ -675,7 +675,7 @@
real (kind=RKIND), dimension(:,:), pointer :: east, north, edge_normal, pressure
real (kind=RKIND), dimension(:,:), pointer :: Ux_tend_halo, Uy_tend_halo, theta_tend, u_tend
real (kind=RKIND), dimension(:,:), pointer :: Ux_phys_tend, Uy_phys_tend
- real (kind=RKIND), dimension(:,:,:), pointer :: scalars, scalars_phys_tend
+ real (kind=RKIND), dimension(:,:,:), pointer :: scalars, scalars_tend
type (block_type), pointer :: block
write(0,*) 'Called CAM_TO_MPAS'
@@ -689,7 +689,7 @@
edge_normal => block % mesh % edgeNormalVectors % array
pressure => block % time_levs(1) % state % pressure % array
scalars => block % time_levs(1) % state % scalars % array
- scalars_phys_tend => block % time_levs(1) % state % scalars_phys_tend % array
+ scalars_tend => block % time_levs(1) % state % scalars_phys_tend % array
theta_tend => block % time_levs(1) % state % theta_phys_tend % array
u_tend => block % time_levs(1) % state % u_phys_tend % array
Ux_phys_tend => block % time_levs(1) % state % Ux_phys_tend % array
@@ -765,7 +765,7 @@
! properly assign variables from Tracer array provided from CAM physics
!
do iScalar=1,num_scalars
- scalars_phys_tend(iScalar,k,iCell) = (Tracer(iCell,k,iScalar) - scalars(iScalar,k,iCell)) / dt_physics
+ scalars_tend(iScalar,k,iCell) = (Tracer(iCell,k,iScalar) - scalars(iScalar,k,iCell)) / dt_physics
end do
end do
end do
@@ -780,7 +780,7 @@
call dmpar_exch_halo_field2dReal(domain % dminfo, block % time_levs(1) % state % u_phys_tend % array(:,:), &
block % mesh % nVertLevels, block % mesh % nEdges, &
block % parinfo % edgesToSend, block % parinfo % edgesToRecv)
- call dmpar_exch_halo_field3dReal(domain % dminfo, block % time_levs(1) % state % scalars % array(:,:,:), &
+ call dmpar_exch_halo_field3dReal(domain % dminfo, block % time_levs(1) % state % scalars_phys_tend % array(:,:,:), &
num_scalars, block % mesh % nVertLevels, block % mesh % nCells, &
block % parinfo % cellsToSend, block % parinfo % cellsToRecv)
</font>
</pre>