<p><b>duda</b> 2010-06-10 14:40:08 -0600 (Thu, 10 Jun 2010)</p><p>BRANCH COMMIT<br>
<br>
Write CAM physics tendencies (theta_phys_tend, Ux_phys_tend, Uy_phys_tend,<br>
u_phys_tend, q*_phys_tend) to output files.<br>
<br>
<br>
M    src/core_hyd_atmos/Registry<br>
M    src/driver_cam_interface/module_mpas_cam_interface.F<br>
</p><hr noshade><pre><font color="gray">Modified: branches/mpas_cam_coupling/src/core_hyd_atmos/Registry
===================================================================
--- branches/mpas_cam_coupling/src/core_hyd_atmos/Registry        2010-06-10 15:25:39 UTC (rev 341)
+++ branches/mpas_cam_coupling/src/core_hyd_atmos/Registry        2010-06-10 20:40:08 UTC (rev 342)
@@ -161,11 +161,13 @@
 #var real    tracers_old ( nTracers nVertLevels nCells ) - tracers_old - -
 
 # Physics tendencies
-var real    theta_phys_tend ( nVertLevels nCells Time ) - theta_phys_tend - -
-var real    u_phys_tend ( nVertLevels nEdges Time ) - u_phys_tend - -
-var real    qv_phys_tend ( nVertLevels nCells Time ) - qv_phys_tend scalars_phys_tend moist_phys_tend
-var real    qc_phys_tend ( nVertLevels nCells Time ) - qc_phys_tend scalars_phys_tend moist_phys_tend
-var real    qr_phys_tend ( nVertLevels nCells Time ) - qr_phys_tend scalars_phys_tend moist_phys_tend
+var real    theta_phys_tend ( nVertLevels nCells Time ) o theta_phys_tend - -
+var real    Ux_phys_tend ( nVertLevels nCells Time ) o Ux_phys_tend - -      # original zonal tendency from CAM
+var real    Uy_phys_tend ( nVertLevels nCells Time ) o Uy_phys_tend - -      # original meridional tendency from CAM
+var real    u_phys_tend ( nVertLevels nEdges Time ) o u_phys_tend - -
+var real    qv_phys_tend ( nVertLevels nCells Time ) o qv_phys_tend scalars_phys_tend moist_phys_tend
+var real    qc_phys_tend ( nVertLevels nCells Time ) o qc_phys_tend scalars_phys_tend moist_phys_tend
+var real    qr_phys_tend ( nVertLevels nCells Time ) o qr_phys_tend scalars_phys_tend moist_phys_tend
 
 # Space needed for advection
 var real    deriv_two ( FIFTEEN TWO nEdges ) o deriv_two - -

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-06-10 15:25:39 UTC (rev 341)
+++ branches/mpas_cam_coupling/src/driver_cam_interface/module_mpas_cam_interface.F        2010-06-10 20:40:08 UTC (rev 342)
@@ -674,6 +674,7 @@
       integer, dimension(:,:), pointer :: edgesOnCell
       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
       type (block_type), pointer :: block
 
@@ -691,6 +692,8 @@
       scalars_phys_tend =&gt; block % time_levs(1) % state % scalars_phys_tend % array
       theta_tend =&gt; block % time_levs(1) % state % theta_phys_tend % array
       u_tend =&gt; block % time_levs(1) % state % u_phys_tend % array
+      Ux_phys_tend =&gt; block % time_levs(1) % state % Ux_phys_tend % array
+      Uy_phys_tend =&gt; block % time_levs(1) % state % Uy_phys_tend % array
 
       !
       !  Perform basic sanity check on expected and available field dimensions
@@ -742,6 +745,12 @@
                                                                                   edge_normal(3,iEdge) * north(3,iCell))
             end do
          end do
+
+         ! Store original zonal and meridional tendencies from CAM to be written to output files for testing
+         do k=1,block % mesh % nVertLevels
+            Ux_phys_tend(k,iCell) = Ux_tend_halo(k,iCell)
+            Uy_phys_tend(k,iCell) = Uy_tend_halo(k,iCell)
+         end do
       end do
 
       deallocate(Ux_tend_halo)

</font>
</pre>