<p><b>duda</b> 2010-07-19 12:16:23 -0600 (Mon, 19 Jul 2010)</p><p>BRANCH COMMIT<br>
<br>
Convert to wet mixing ratios when passing scalar fields to CAM, and convert<br>
to dry mixing ratios when receiving scalar fields from CAM.<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-16 23:31:17 UTC (rev 381)
+++ branches/mpas_cam_coupling/src/driver_cam_interface/module_mpas_cam_interface.F        2010-07-19 18:16:23 UTC (rev 382)
@@ -630,7 +630,7 @@
! are what CAM expects, rather than simply what we have in our scalars array
!
do iScalar=1,num_scalars
- Tracer(iCell,k,iScalar) = scalars(iScalar,k,iCell)
+ Tracer(iCell,k,iScalar) = scalars(iScalar,k,iCell) / (1.0 + scalars(index_qv,k,iCell))
end do
end do
end do
@@ -770,7 +770,9 @@
! properly assign variables from Tracer array provided from CAM physics
!
do iScalar=1,num_scalars
- scalars_tend(iScalar,k,iCell) = (Tracer(iCell,k,iScalar) - scalars(iScalar,k,iCell)) / dt_physics
+ scalars_tend(iScalar,k,iCell) = (Tracer(iCell,k,iScalar) * (1.0 + Tracer(iCell,k,index_qv)) &
+ - scalars(iScalar,k,iCell)&
+ ) / dt_physics
! Couple scalar tendencies with h valid at time when tendencies were computed
scalars_tend(iScalar,k,iCell) = scalars_tend(iScalar,k,iCell) * h(k,iCell)
</font>
</pre>