<p><b>laura@ucar.edu</b> 2012-07-12 11:36:23 -0600 (Thu, 12 Jul 2012)</p><p>now use RvoRd for more accurate calculation of the potential temperature as functions of the modified temperature and watervapor in subroutines atm_compute_output_diagnostics, atm_compute_restart_diagnostics, and atm_init_coupled_diagnostics<br>
</p><hr noshade><pre><font color="gray">Modified: branches/atmos_physics/src/core_nhyd_atmos/mpas_atm_mpas_core.F
===================================================================
--- branches/atmos_physics/src/core_nhyd_atmos/mpas_atm_mpas_core.F        2012-07-12 17:33:19 UTC (rev 2017)
+++ branches/atmos_physics/src/core_nhyd_atmos/mpas_atm_mpas_core.F        2012-07-12 17:36:23 UTC (rev 2018)
@@ -366,6 +366,7 @@
    !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    
       use mpas_grid_types
+      use mpas_constants
    
       implicit none
    
@@ -378,7 +379,7 @@
 
       do iCell=1,grid%nCells
          do k=1,grid%nVertLevels
-            diag % theta % array(k,iCell) = state % theta_m % array(k,iCell) / (1.0 + 1.61 * state % scalars % array(state % index_qv,k,iCell))
+            diag % theta % array(k,iCell) = state % theta_m % array(k,iCell) / (1._RKIND + RvoRd * state % scalars % array(state % index_qv,k,iCell))
             diag % rho % array(k,iCell) = state % rho_zz % array(k,iCell) * grid % zz % array(k,iCell)
          end do
       end do
@@ -397,6 +398,7 @@
    !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    
       use mpas_grid_types
+      use mpas_constants
    
       implicit none
    
@@ -409,7 +411,7 @@
 
       do iCell=1,grid%nCells
          do k=1,grid%nVertLevels
-            diag % theta % array(k,iCell) = state % theta_m % array(k,iCell) / (1.0 + 1.61 * state % scalars % array(state % index_qv,k,iCell))
+            diag % theta % array(k,iCell) = state % theta_m % array(k,iCell) / (1._RKIND + RvoRd * state % scalars % array(state % index_qv,k,iCell))
             diag % rho % array(k,iCell) = state % rho_zz % array(k,iCell) * grid % zz % array(k,iCell)
          end do
       end do

Modified: branches/atmos_physics/src/core_nhyd_atmos/mpas_atm_time_integration.F
===================================================================
--- branches/atmos_physics/src/core_nhyd_atmos/mpas_atm_time_integration.F        2012-07-12 17:33:19 UTC (rev 2017)
+++ branches/atmos_physics/src/core_nhyd_atmos/mpas_atm_time_integration.F        2012-07-12 17:36:23 UTC (rev 2018)
@@ -3338,7 +3338,7 @@
 
       do iCell=1,grid%nCells
          do k=1,grid%nVertLevels
-            state % theta_m % array(k,iCell) = diag % theta % array(k,iCell) * (1.0 + 1.61 * state % scalars % array(state % index_qv,k,iCell))
+            state % theta_m % array(k,iCell) = diag % theta % array(k,iCell) * (1._RKIND + RvoRd * state % scalars % array(state % index_qv,k,iCell))
             state % rho_zz % array(k,iCell) = diag % rho % array(k,iCell) / grid % zz % array(k,iCell)
          end do
       end do

</font>
</pre>