<p><b>duda</b> 2013-03-07 15:41:19 -0700 (Thu, 07 Mar 2013)</p><p>BRANCH COMMIT<br>
<br>
Two changes to MPAS-A initialization from Joe:<br>
<br>
1) Correct coefficient in terrain smoother to give fourth-order smoothing.<br>
<br>
2) Modify several parameters that influence coordinate surface smoothing.<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/atmos_physics/src/core_init_nhyd_atmos/mpas_init_atm_test_cases.F
===================================================================
--- branches/atmos_physics/src/core_init_nhyd_atmos/mpas_init_atm_test_cases.F        2013-03-07 22:36:53 UTC (rev 2563)
+++ branches/atmos_physics/src/core_init_nhyd_atmos/mpas_init_atm_test_cases.F        2013-03-07 22:41:19 UTC (rev 2564)
@@ -2510,7 +2510,7 @@
                                         * (ter(cellsOnCell(j,iCell))-ter(iCell))
                end do
             endif
-            hs(iCell) = ter(iCell) + 0.25*hs(iCell)
+            hs(iCell) = ter(iCell) + 0.125*hs(iCell)
          end do
 
          do iCell=1,grid %nCells
@@ -2633,7 +2633,7 @@
 
       if (smooth) then
 
-         dzmin = 0.3
+         dzmin = 0.5
 
          do k=2,kz-1
             hx(k,:) = hx(k-1,:)
@@ -2641,9 +2641,9 @@
 
 !            dzmin = max(0.5_RKIND,1.-.5*zw(k)/hm)
 
-            sm = .05*min(0.5_RKIND*zw(k)/hm,1.0_RKIND)
+            sm = .02*min(0.5_RKIND*zw(k)/hm,1.0_RKIND)
           
-            do i=1,50
+            do i=1,30
                do iCell=1,grid % nCells
                   hs1(iCell) = 0.
                   do j = 1,nEdgesOnCell(iCell)

</font>
</pre>