<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>