<p><b>duda</b> 2012-07-27 16:31:48 -0600 (Fri, 27 Jul 2012)</p><p>BRANCH COMMIT<br>
<br>
Minor fix to theta pert in 3-1: need to take z at p-levels as (zgrid(k,iCell)+zgrid(k+1,iCell))/2.<br>
<br>
<br>
M    src/core_init_nhyd_atmos/mpas_init_atm_test_cases.F<br>
</p><hr noshade><pre><font color="gray">Modified: branches/dcmip/src/core_init_nhyd_atmos/mpas_init_atm_test_cases.F
===================================================================
--- branches/dcmip/src/core_init_nhyd_atmos/mpas_init_atm_test_cases.F        2012-07-27 22:25:41 UTC (rev 2072)
+++ branches/dcmip/src/core_init_nhyd_atmos/mpas_init_atm_test_cases.F        2012-07-27 22:31:48 UTC (rev 2073)
@@ -5378,7 +5378,7 @@
                s = widthParm**2.0 / (widthParm**2.0 + sphere_distance(theta_c,                   lambda_c,              &amp;
                                                                       grid%latCell%array(iCell), grid%lonCell%array(iCell), &amp;
                                                                       grid%sphere_radius)**2.0)
-               theta_pert = dTheta * s * sin((2.0_RKIND * pii * zgrid(k,iCell)) / L_z)
+               theta_pert = dTheta * s * sin((2.0_RKIND * pii * 0.5*(zgrid(k,iCell)+zgrid(k+1,iCell))) / L_z)
                diag % theta % array(k,iCell) = diag % theta % array(k,iCell) + theta_pert
             end if
 

</font>
</pre>