<p><b>sprice@lanl.gov</b> 2011-08-31 13:49:07 -0600 (Wed, 31 Aug 2011)</p><p>Branch Commit (land ice): more corrections to modified basin code (icesheet.F) for extracting ice sheet mesh from periodic hex mesh; thickness and basal topog fileds now included (thck &amp; topg) in output .nc file; mesh partitioning files also written out now. Code still needs some general tidying up.<br>
</p><hr noshade><pre><font color="gray">Modified: branches/land_ice/icesheet/src/icesheet.F
===================================================================
--- branches/land_ice/icesheet/src/icesheet.F        2011-08-31 17:56:45 UTC (rev 969)
+++ branches/land_ice/icesheet/src/icesheet.F        2011-08-31 19:49:07 UTC (rev 970)
@@ -188,9 +188,9 @@
 call write_grid
 
 ! dump graph for partioning
-!write(6,*) ' call write_graph'
-!write(6,*)
-!call write_graph
+write(6,*) ' call write_graph'
+write(6,*)
+call write_graph
 
 ! write OpenDx file
 write(6,*) ' calling write_OpenDX'
@@ -695,7 +695,7 @@
 !                    u_srcNew, &amp;
 !                    uNew, &amp;
 !                    vNew, &amp;
-!                    hNew, &amp;
+                    hNew, &amp;
 !                    rhoNew, &amp;
 !                    temperatureNew, &amp;
 !                    salinityNew, &amp;

Modified: branches/land_ice/icesheet/src/module_write_netcdf.F
===================================================================
--- branches/land_ice/icesheet/src/module_write_netcdf.F        2011-08-31 17:56:45 UTC (rev 969)
+++ branches/land_ice/icesheet/src/module_write_netcdf.F        2011-08-31 19:49:07 UTC (rev 970)
@@ -54,7 +54,7 @@
    integer :: wrVarIDboundaryVertex
 !   integer :: wrVarIDu_src
 !   integer :: wrVarIDv
-!   integer :: wrVarIDh
+   integer :: wrVarIDh
 !   integer :: wrVarIDrho
 !   integer :: wrVarIDtemperature
 !   integer :: wrVarIDsalinity
@@ -241,10 +241,10 @@
 !      dimlist( 1) = wrDimIDnVertLevels
 !      dimlist( 2) = wrDimIDnCells
 !      dimlist( 3) = wrDimIDTime
-!      nferr = nf_def_var(wr_ncid, 'h', NF_DOUBLE,  3, dimlist, wrVarIDh)
-!      dimlist( 1) = wrDimIDnVertLevels
-!      dimlist( 2) = wrDimIDnCells
-!      dimlist( 3) = wrDimIDTime
+      nferr = nf_def_var(wr_ncid, 'h', NF_DOUBLE,  3, dimlist, wrVarIDh)
+      dimlist( 1) = wrDimIDnVertLevels
+      dimlist( 2) = wrDimIDnCells
+      dimlist( 3) = wrDimIDTime
 !      nferr = nf_def_var(wr_ncid, 'rho', NF_DOUBLE,  3, dimlist, wrVarIDrho)
 !      dimlist( 1) = wrDimIDnVertLevels
 !      dimlist( 2) = wrDimIDnCells
@@ -325,7 +325,7 @@
 !                                  u_src, &amp;
 !                                  u, &amp;
 !                                  v, &amp;
-!                                  h, &amp;
+                                  h, &amp;
 !                                  rho, &amp;
 !                                  temperature, &amp;
 !                                  salinity, &amp;
@@ -385,7 +385,7 @@
 !      real (kind=8), dimension(:,:), intent(in) :: u_src
 !      real (kind=8), dimension(:,:,:), intent(in) :: u
 !      real (kind=8), dimension(:,:,:), intent(in) :: v
-!      real (kind=8), dimension(:,:,:), intent(in) :: h
+      real (kind=8), dimension(:,:,:), intent(in) :: h
 !      real (kind=8), dimension(:,:,:), intent(in) :: rho
 !      real (kind=8), dimension(:,:,:), intent(in) :: temperature
 !      real (kind=8), dimension(:,:,:), intent(in) :: salinity
@@ -618,11 +618,11 @@
 !      count3( 3) = 1
 !      nferr = nf_put_vara_double(wr_ncid, wrVarIDv, start3, count3, v)
  
-!      start3(3) = time
-!      count3( 1) = wrLocalnVertLevels
-!      count3( 2) = wrLocalnCells
-!      count3( 3) = 1
-!      nferr = nf_put_vara_double(wr_ncid, wrVarIDh, start3, count3, h)
+      start3(3) = time
+      count3( 1) = wrLocalnVertLevels
+      count3( 2) = wrLocalnCells
+      count3( 3) = 1
+      nferr = nf_put_vara_double(wr_ncid, wrVarIDh, start3, count3, h)
  
 !      start3(3) = time
 !      count3( 1) = wrLocalnVertLevels

</font>
</pre>