<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 & 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, &
! uNew, &
! vNew, &
-! hNew, &
+ hNew, &
! rhoNew, &
! temperatureNew, &
! salinityNew, &
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, &
! u, &
! v, &
-! h, &
+ h, &
! rho, &
! temperature, &
! salinity, &
@@ -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>