<p><b>mperego@fsu.edu</b> 2012-01-10 22:40:48 -0700 (Tue, 10 Jan 2012)</p><p>Branch land_ice: update test_steady_state_greenland.F <br>
</p><hr noshade><pre><font color="gray">Modified: branches/land_ice/icesheet/src/test_steady_state_greenland.F
===================================================================
--- branches/land_ice/icesheet/src/test_steady_state_greenland.F 2012-01-11 00:03:13 UTC (rev 1352)
+++ branches/land_ice/icesheet/src/test_steady_state_greenland.F 2012-01-11 05:40:48 UTC (rev 1353)
@@ -26,7 +26,7 @@
real, dimension(:), pointer :: thckdata, elevdata, betadata, temperaturedata, u, v;
- integer nVerticesSolve, nEdgesSolve, nCellsSolve, nVertices, nEdges, nCells, i, k, iCell, nLevels
+ integer nVerticesSolve, nEdgesSolve, nCellsSolve, nVertices, nEdges, nCells, i, k, iCell, nLevels, nKeptVerteces
logical keep_vertex
@@ -39,7 +39,7 @@
call dmpar_init(dminfo)
! initialize lifev data structure
- call init_lifev(dminfo % comm)
+
call read_namelist(dminfo)
@@ -56,7 +56,7 @@
nCells = domain % blocklist % mesh % nCells
nEdges = domain % blocklist % mesh % nEdges
nVertices = domain % blocklist % mesh % nVertices
- nLevels = 10 !domain % blocklist % mesh % nVertLevels
+ nLevels = domain % blocklist % mesh % nVertLevels
!get elevation and thickness on each cell
@@ -80,16 +80,23 @@
!(discard a vertex if all the cells around it have zero thickness)
allocate(verticesMask(nVertices));
+ nKeptVerteces = 0;
do i = 1,nVertices
keep_vertex = .false.
do k = 1, 3
iCell = domain % blocklist % mesh % cellsOnVertex % array(k,i)
keep_vertex = keep_vertex.or.(thckdata(iCell) > 0.01)
end do
+ if (keep_vertex .and. (i <= nVerticesSolve )) then
+ nKeptVerteces = nKeptVerteces + 1
+ endif
verticesMask(i) = keep_vertex
end do
+
+write(6,*) 'nKeptVerteces', nKeptVerteces
+
+ call init_lifev(dminfo % comm, nKeptVerteces)
-
!build send and receive arrays using exchange_list
call array_from_exchange_list(domain % blocklist % parinfo % verticesToSend, sendVerticesArray)
</font>
</pre>