<p><b>duda</b> 2012-08-01 14:56:02 -0600 (Wed, 01 Aug 2012)</p><p>Update quicksort() to mpas_quicksort() in non-hydro initialization core.<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/omp_blocks/multiple_blocks/src/core_init_nhyd_atmos/mpas_init_atm_test_cases.F
===================================================================
--- branches/omp_blocks/multiple_blocks/src/core_init_nhyd_atmos/mpas_init_atm_test_cases.F        2012-08-01 20:24:12 UTC (rev 2077)
+++ branches/omp_blocks/multiple_blocks/src/core_init_nhyd_atmos/mpas_init_atm_test_cases.F        2012-08-01 20:56:02 UTC (rev 2078)
@@ -4020,7 +4020,7 @@
if (vert_level(k) == 200100.0) sorted_arr(1,k) = 99999.0
sorted_arr(2,k) = fg % t % array(k,iCell)
end do
- call quicksort(config_nfglevels, sorted_arr)
+ call mpas_quicksort(config_nfglevels, sorted_arr)
do k=1,grid%nVertLevels
target_z = 0.5 * (grid % zgrid % array(k,iCell) + grid % zgrid % array(k+1,iCell))
state % theta_m % array(k,iCell) = vertical_interp(target_z, config_nfglevels, sorted_arr, order=1, extrap=1)
@@ -4035,7 +4035,7 @@
if (vert_level(k) == 200100.0) sorted_arr(1,k) = 99999.0
sorted_arr(2,k) = fg % rh % array(k,iCell)
end do
- call quicksort(config_nfglevels, sorted_arr)
+ call mpas_quicksort(config_nfglevels, sorted_arr)
do k=1,grid%nVertLevels
target_z = 0.5 * (grid % zgrid % array(k,iCell) + grid % zgrid % array(k+1,iCell))
state % scalars % array(state % index_qv,k,iCell) = vertical_interp(target_z, config_nfglevels, sorted_arr, order=1, extrap=0)
@@ -4050,7 +4050,7 @@
if (vert_level(k) == 200100.0) sorted_arr(1,k) = 99999.0
sorted_arr(2,k) = fg % z % array(k,iCell)
end do
- call quicksort(config_nfglevels, sorted_arr)
+ call mpas_quicksort(config_nfglevels, sorted_arr)
do k=1,grid%nVertLevels
target_z = 0.5 * (grid % zgrid % array(k,iCell) + grid % zgrid % array(k+1,iCell))
fg % gfs_z % array(k,iCell) = vertical_interp(target_z, config_nfglevels, sorted_arr, order=1, extrap=1)
@@ -4068,7 +4068,7 @@
end if
sorted_arr(2,k) = log(fg % p % array(k,iCell))
end do
- call quicksort(config_nfglevels, sorted_arr)
+ call mpas_quicksort(config_nfglevels, sorted_arr)
do k=1,grid%nVertLevels
target_z = 0.5 * (grid % zgrid % array(k,iCell) + grid % zgrid % array(k+1,iCell))
diag % pressure % array(k,iCell) = exp(vertical_interp(target_z, config_nfglevels, sorted_arr, order=1, extrap=1))
@@ -4086,7 +4086,7 @@
! end if
! sorted_arr(2,k) = log(fg % p % array(k,iCell))
! end do
-! call quicksort(config_nfglevels, sorted_arr)
+! call mpas_quicksort(config_nfglevels, sorted_arr)
! do k=1,grid%nVertLevels+1
! target_z = grid % zgrid % array(k,iCell)
! fg % gfs_p % array(k,iCell) = exp(vertical_interp(target_z, config_nfglevels, sorted_arr, order=1, extrap=1))
@@ -4105,7 +4105,7 @@
if (vert_level(k) == 200100.0) sorted_arr(1,k) = 99999.0
sorted_arr(2,k) = fg % u % array(k,iEdge)
end do
- call quicksort(config_nfglevels, sorted_arr)
+ call mpas_quicksort(config_nfglevels, sorted_arr)
do k=1,grid%nVertLevels
target_z = 0.25 * (grid % zgrid % array(k,cellsOnEdge(1,iEdge)) + grid % zgrid % array(k+1,cellsOnEdge(1,iEdge)) + grid % zgrid % array(k,cellsOnEdge(2,iEdge)) + grid % zgrid % array(k+1,cellsOnEdge(2,iEdge)))
state % u % array(k,iEdge) = vertical_interp(target_z, config_nfglevels, sorted_arr, order=1, extrap=0)
@@ -4146,7 +4146,7 @@
end if
sorted_arr(2,k) = log(fg % p % array(k,iCell))
end do
- call quicksort(config_nfglevels, sorted_arr)
+ call mpas_quicksort(config_nfglevels, sorted_arr)
target_z = grid % zgrid % array(1,iCell)
fg % psfc % array(iCell) = exp(vertical_interp(target_z, config_nfglevels, sorted_arr, order=1, extrap=1))
</font>
</pre>