<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 => 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
+ Uy_phys_tend => 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>