<p><b>laura@ucar.edu</b> 2012-08-06 09:12:54 -0600 (Mon, 06 Aug 2012)</p><p>changed RvoRd to lower case rvord<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-08-05 21:40:48 UTC (rev 2084)
+++ branches/atmos_physics/src/core_nhyd_atmos/mpas_atm_mpas_core.F        2012-08-06 15:12:54 UTC (rev 2085)
@@ -391,7 +391,7 @@
do iCell=1,grid%nCells
do k=1,grid%nVertLevels
- diag % theta % array(k,iCell) = state % theta_m % array(k,iCell) / (1._RKIND + RvoRd * 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
@@ -423,7 +423,7 @@
do iCell=1,grid%nCells
do k=1,grid%nVertLevels
- diag % theta % array(k,iCell) = state % theta_m % array(k,iCell) / (1._RKIND + RvoRd * 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-08-05 21:40:48 UTC (rev 2084)
+++ branches/atmos_physics/src/core_nhyd_atmos/mpas_atm_time_integration.F        2012-08-06 15:12:54 UTC (rev 2085)
@@ -3354,7 +3354,7 @@
do iCell=1,grid%nCells
do k=1,grid%nVertLevels
- state % theta_m % array(k,iCell) = diag % theta % array(k,iCell) * (1._RKIND + RvoRd * 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
Modified: branches/atmos_physics/src/framework/mpas_constants.F
===================================================================
--- branches/atmos_physics/src/framework/mpas_constants.F        2012-08-05 21:40:48 UTC (rev 2084)
+++ branches/atmos_physics/src/framework/mpas_constants.F        2012-08-06 15:12:54 UTC (rev 2085)
@@ -7,8 +7,8 @@
real (kind=RKIND), parameter :: omega = 7.29212e-5
real (kind=RKIND), parameter :: gravity = 9.80616
real (kind=RKIND), parameter :: rgas = 287.
- real (kind=RKIND), parameter :: Rv = 461.6
- real (kind=RKIND), parameter :: RvoRd = Rv/rgas
+ real (kind=RKIND), parameter :: rv = 461.6
+ real (kind=RKIND), parameter :: rvord = rv/rgas
real (kind=RKIND), parameter :: cp = 1003.
real (kind=RKIND), parameter :: cv = cp - rgas
real (kind=RKIND), parameter :: cvpm = -cv/cp
</font>
</pre>