<p><b>mpetersen@lanl.gov</b> 2011-11-10 15:47:15 -0700 (Thu, 10 Nov 2011)</p><p>Updated basin files to include hZLevel variable in ocean.nc initial condition files.<br>
</p><hr noshade><pre><font color="gray">Modified: branches/ocean_projects/basin/src/Makefile
===================================================================
--- branches/ocean_projects/basin/src/Makefile        2011-11-09 21:30:55 UTC (rev 1183)
+++ branches/ocean_projects/basin/src/Makefile        2011-11-10 22:47:15 UTC (rev 1184)
@@ -15,9 +15,9 @@
 # ifort
 FC = ifort
 CC = icc
-FFLAGS = -real-size 64 -g -traceback -check all
-CFLAGS = -g
-LDFLAGS = -g -traceback -check all
+FFLAGS = -real-size 64 #-g -traceback -check all
+CFLAGS = #-g
+LDFLAGS = #-g -traceback -check all
 
 # absoft
 #FC = f90

Modified: branches/ocean_projects/basin/src/basin.F
===================================================================
--- branches/ocean_projects/basin/src/basin.F        2011-11-09 21:30:55 UTC (rev 1183)
+++ branches/ocean_projects/basin/src/basin.F        2011-11-10 22:47:15 UTC (rev 1184)
@@ -116,7 +116,7 @@
 real, allocatable, dimension(:) :: areaCellNew, areaTriangleNew, dcEdgeNew, dvEdgeNew, angleEdgeNew
 real, allocatable, dimension(:,:) :: kiteAreasOnVertexNew, weightsOnEdgeNew, normalsNew
 
-real, allocatable, dimension(:) :: fEdgeNew, fVertexNew, h_sNew
+real, allocatable, dimension(:) :: fEdgeNew, fVertexNew, h_sNew, hZLevel
 real, allocatable, dimension(:,:) :: u_srcNew
 real, allocatable, dimension(:,:,:) :: uNew, vNew, hNew
 real, allocatable, dimension(:,:,:) :: rhoNew, temperatureNew, salinityNew, tracer1New
@@ -324,6 +324,7 @@
        hNew(1,1,:) = 3250.0
        h_sNew(:) = -( hNew(1,1,:) )
    endif
+   hZLevel = hNew(1,:,1)
 
    ! basin-mod
    ! Noise is meant to make the flow unstable at some point 
@@ -544,6 +545,7 @@
 iNoData = 0
 do iCell=1,nCellsNew
   hNew(1,:,iCell) = dz(:)
+  hZLevel = dz
   if(mod(iCell,100).eq.0) write(6,*) 'l_woce t and s',iCell
   rlon = lonCellNew(iCell)/dtr
   rlat = latCellNew(iCell)/dtr
@@ -1052,7 +1054,8 @@
                     salinityNew, &amp;
                     tracer1New, &amp;
                     temperatureRestoreNew, &amp;
-                    salinityRestoreNew &amp;
+                    salinityRestoreNew, &amp;
+                    hZLevel &amp;
                    )
 
 call write_netcdf_finalize
@@ -1329,6 +1332,7 @@
 allocate(uNew(1,nVertLevelsNew,nEdgesNew))
 allocate(vNew(1,nVertLevelsNew,nEdgesNew))
 allocate(hNew(1,nVertLevelsNew,nCellsNew))
+allocate(hZLevel(nVertLevelsNew))
 allocate(rhoNew(1,nVertLevelsNew,nCellsNew))
 allocate(temperatureNew(1,nVertLevelsNew,nCellsNew))
 allocate(salinityNew(1,nVertLevelsNew,nCellsNew))

Modified: branches/ocean_projects/basin/src/module_write_netcdf.F
===================================================================
--- branches/ocean_projects/basin/src/module_write_netcdf.F        2011-11-09 21:30:55 UTC (rev 1183)
+++ branches/ocean_projects/basin/src/module_write_netcdf.F        2011-11-10 22:47:15 UTC (rev 1184)
@@ -62,6 +62,7 @@
    integer :: wrVarIDtracer1
    integer :: wrVarIDtemperatureRestore
    integer :: wrVarIDsalinityRestore
+   integer :: wrVarIDhZLevel
  
    integer :: wrLocalnCells
    integer :: wrLocalnEdges
@@ -221,6 +222,8 @@
       dimlist( 1) = wrDimIDnCells
       nferr = nf_def_var(wr_ncid, 'salinityRestore', NF_DOUBLE,  1, dimlist, wrVarIDsalinityRestore)
       dimlist( 1) = wrDimIDnVertLevels
+      nferr = nf_def_var(wr_ncid, 'hZLevel', NF_DOUBLE,  1, dimlist, wrVarIDhZLevel)
+      dimlist( 1) = wrDimIDnVertLevels
       dimlist( 2) = wrDimIDnEdges
       dimlist( 3) = wrDimIDTime
       nferr = nf_def_var(wr_ncid, 'u', NF_DOUBLE,  3, dimlist, wrVarIDu)
@@ -321,7 +324,8 @@
                                   salinity, &amp;
                                   tracer1, &amp;
                                   temperatureRestore, &amp;
-                                  salinityRestore &amp;
+                                  salinityRestore, &amp;
+                                  hZLevel &amp;
                                  )
  
       implicit none
@@ -381,6 +385,7 @@
       real (kind=8), dimension(:,:,:), intent(in) :: tracer1
       real (kind=8), dimension(:), intent(in) :: temperatureRestore
       real (kind=8), dimension(:), intent(in) :: salinityRestore
+      real (kind=8), dimension(:), intent(in) :: hZLevel
 
  
       integer :: nferr
@@ -596,6 +601,10 @@
       count1( 1) = wrLocalnCells
       nferr = nf_put_vara_double(wr_ncid, wrVarIDsalinityRestore, start1, count1, salinityRestore)
  
+      start1(1) = 1
+      count1( 1) = wrLocalnVertLevels
+      nferr = nf_put_vara_double(wr_ncid, wrVarIDhZLevel, start1, count1, hZLevel)
+ 
       start3(3) = time
       count3( 1) = wrLocalnVertLevels
       count3( 2) = wrLocalnEdges

</font>
</pre>