<p><b>duda</b> 2012-09-04 14:15:16 -0600 (Tue, 04 Sep 2012)</p><p>** Changes for atmosphere cores only **<br>
<br>
With the changes made in r2122, nVertLevelsSolve is no longer set. Rather than<br>
set nVertLevelsSolve to nVertLevels ourselves in MPAS-A cores, just revert to<br>
using nVertLevels for several loop bounds (and leave a comment that these bounds<br>
might be modified if we support vertical decomposition at some future time).<br>
<br>
Note: No other cores (sw, ocean, atmosphere initialization) use<br>
nVertLevelsSolve, so no changes are required there.<br>
<br>
<br>
M    core_hyd_atmos/mpas_atmh_time_integration.F<br>
M    core_nhyd_atmos/mpas_atm_time_integration.F<br>
</p><hr noshade><pre><font color="gray">Modified: trunk/mpas/src/core_hyd_atmos/mpas_atmh_time_integration.F
===================================================================
--- trunk/mpas/src/core_hyd_atmos/mpas_atmh_time_integration.F        2012-09-03 14:00:40 UTC (rev 2139)
+++ trunk/mpas/src/core_hyd_atmos/mpas_atmh_time_integration.F        2012-09-04 20:15:16 UTC (rev 2140)
@@ -276,7 +276,7 @@
                                          block % mesh % areaCell % array (iCell) &amp;
                                        - block % state % time_levs(2) % state % pressure % array (block % mesh % nVertLevels + 1, 1) * &amp;
                                          block % mesh % areaCell % array (iCell)
-             do k=1, block % mesh % nVertLevelsSolve
+             do k=1, block % mesh % nVertLevels   ! Could be nVertLevelsSolve?
                scalar_mass = scalar_mass - block % state % time_levs(2) % state % scalars % array (2,k,iCell) * &amp;
                                            block % state % time_levs(2) % state % h % array (k,iCell) * &amp;
                                            block % mesh % dnw % array (k) * &amp;
@@ -1378,7 +1378,7 @@
         end do
         wdtn(:,nVertLevels+1) = 0.
 
-         do k=1,grid % nVertLevelsSolve
+         do k=1,grid % nVertLevels   ! Could be nVertLevelsSolve?
             do iScalar=1,num_scalars
               scalar_new(iScalar,k,iCell) = (   scalar_old(iScalar,k,iCell)*h_old(k,iCell) &amp;
                     + dt*( scalar_tend(iScalar,k,iCell) -rdnw(k)*(wdtn(iScalar,k+1)-wdtn(iScalar,k)) ) )/h_new(k,iCell)

Modified: trunk/mpas/src/core_nhyd_atmos/mpas_atm_time_integration.F
===================================================================
--- trunk/mpas/src/core_nhyd_atmos/mpas_atm_time_integration.F        2012-09-03 14:00:40 UTC (rev 2139)
+++ trunk/mpas/src/core_nhyd_atmos/mpas_atm_time_integration.F        2012-09-04 20:15:16 UTC (rev 2140)
@@ -1326,7 +1326,7 @@
                   wdtn(iScalar,k) = wwAvg(k,iCell)*(fnm(k)*scalar_new(iScalar,k,iCell)+fnp(k)*scalar_new(iScalar,k-1,iCell))
                 end do
              end do
-             do k=1,grid % nVertLevelsSolve
+             do k=1,grid % nVertLevels   ! Could be nVertLevelsSolve?
                 do iScalar=1,s_old % num_scalars
                   scalar_new(iScalar,k,iCell) = (   scalar_old(iScalar,k,iCell)*h_old(k,iCell) &amp;
                         + dt*( scalar_tend(iScalar,k,iCell) -rdnw(k)*(wdtn(iScalar,k+1)-wdtn(iScalar,k)) ) )/h_new(k,iCell)
@@ -1355,7 +1355,7 @@
                wdtn(iScalar,k) = wwAvg(k,iCell)*(fnm(k)*scalar_new(iScalar,k,iCell)+fnp(k)*scalar_new(iScalar,k-1,iCell))
              enddo
 
-             do k=1,grid % nVertLevelsSolve
+             do k=1,grid % nVertLevels   ! Could be nVertLevelsSolve?
                 do iScalar=1,s_old % num_scalars
                   scalar_new(iScalar,k,iCell) = (   scalar_old(iScalar,k,iCell)*h_old(k,iCell) &amp;
                         + dt*( scalar_tend(iScalar,k,iCell) -rdnw(k)*(wdtn(iScalar,k+1)-wdtn(iScalar,k)) ) )/h_new(k,iCell)
@@ -1383,7 +1383,7 @@
                wdtn(iScalar,k) = wwAvg(k,iCell)*(fnm(k)*scalar_new(iScalar,k,iCell)+fnp(k)*scalar_new(iScalar,k-1,iCell))
              enddo
 
-             do k=1,grid % nVertLevelsSolve
+             do k=1,grid % nVertLevels   ! Could be nVertLevelsSolve?
                 do iScalar=1,s_old % num_scalars
                   scalar_new(iScalar,k,iCell) = (   scalar_old(iScalar,k,iCell)*h_old(k,iCell) &amp;
                         + dt*( scalar_tend(iScalar,k,iCell) -rdnw(k)*(wdtn(iScalar,k+1)-wdtn(iScalar,k)) ) )/h_new(k,iCell)

</font>
</pre>