<p><b>sprice@lanl.gov</b> 2011-08-16 14:57:45 -0600 (Tue, 16 Aug 2011)</p><p>branch commit (land ice): Minor code cleanup; edited / added some readme files; added namelist.input files for use with specific Greenland datasets; corrections to namelist files for grid.nc generation so that no. of hex cells along y direction is correct (i.e., larger than for reg. struct. grid by a factor of 2/sqrt(3)).<br>
</p><hr noshade><pre><font color="gray">Modified: branches/land_ice/icesheet/README
===================================================================
--- branches/land_ice/icesheet/README        2011-08-16 19:40:21 UTC (rev 942)
+++ branches/land_ice/icesheet/README        2011-08-16 20:57:45 UTC (rev 943)
@@ -1,23 +1,28 @@
-Program: greenland
+This dir contains a hacked version of the "basin" code from the "ocean_projects" dir.
+The goal is to do something similar but for ice sheet grids. That is ...
-This is a hacked version of the "basin" code from the "ocean_projects" dir. The goal is to do
-something similar but for a greenland ice sheet grid. That is ...
+1) use periodic_hex code to generate a periodic hex grid w/ appropriate overall dims that contain
+the geometry data from the ice sheet of interest (e.g. that encompass Greenland or Ant.)
-1) use periodic_hex code to generate a periodic hex grid w/ appropriate dims for greenland DEM
-2) use altered version of basin code to extract cells from that grid that have zero thickness
+2) use altered version of basin code (icesheet.F) to extract cells from that grid for which there are
+non-zero geom data (e.g. for which thickness data are non-zero)
-The .nc file containing the original (lat/lon georefed) data for greenland go in the "topo" dir.
+The .nc files containing the original data for ice sheets are in "topo". The "grid.nc" file for the
+generic, MPAS periodic hex mesh are in the "periodic_hex" subdir. See additional notes in the REAME
+there for creating of that "grid.nc" file.
-The info from the original readme file is included below.
+To execute the code, copy or link to the "grid.nc" file in the "periodic_hex" subdir, then execute
+the script "runit". This will build and run the code that does all of the work (reads in the topo
+data and information from the original MPAS hex mesh, compares them, and extracts only cells from the
+MPAS mesh for which data exist).
-sfp 7/13/11
+sfp 8/16/11
+--- original contents of this readme file ---
---- original contents of readme file ---
-
Program: basin
This program reads in a MPAS grid file (grid.nc) and outputs a modified grid file
Added: branches/land_ice/icesheet/periodic_hex/README
===================================================================
--- branches/land_ice/icesheet/periodic_hex/README         (rev 0)
+++ branches/land_ice/icesheet/periodic_hex/README        2011-08-16 20:57:45 UTC (rev 943)
@@ -0,0 +1,24 @@
+
+This dir contains files for use in creating MPAS period hex mesh grids based on ...
+
+1) a generic perodic hex mesh constructed using the "namelist.input" files stored here
+
+2) a given input .nc dataset, stored in "../topo/"
+
+3) the "periodic_grid" executable in "../../grid_gen/periodic_hex/"
+
+To construct the mesh, first link to the necessary executable ...
+
+ln -s ../../grid_gen/periodic_hex/periodic_grid
+
+Then copy one of the "namelist.specific-file-info.input" files to "namelist.input".
+
+Next, execute "periodic_grid", which will produce a generic MPAS periodic hex mesh
+with the dimensions specified in the "namelist.input" file. This will be used along
+with the other code in "../src/" and the .nc files in "../topo/" to build a more
+specific, boundary conforming hex mesh, based on the geometry read from the input
+.nc files in "../topo/".
+
+Additional instructions for these later steps can be found one dir level up.
+
+sfp 8/16/11
Added: branches/land_ice/icesheet/periodic_hex/namelist.input.Greenland.20km
===================================================================
--- branches/land_ice/icesheet/periodic_hex/namelist.input.Greenland.20km         (rev 0)
+++ branches/land_ice/icesheet/periodic_hex/namelist.input.Greenland.20km        2011-08-16 20:57:45 UTC (rev 943)
@@ -0,0 +1,8 @@
+&periodic_grid
+ nx = 76,
+ ny = 163,
+ dc = 20000.,
+ nVertLevels = 1,
+ nTracers = 1,
+ nproc = 2,
+/
Property changes on: branches/land_ice/icesheet/periodic_hex/namelist.input.Greenland.20km
___________________________________________________________________
Added: svn:executable
+ *
Added: branches/land_ice/icesheet/periodic_hex/namelist.input.Greenland.5km
===================================================================
--- branches/land_ice/icesheet/periodic_hex/namelist.input.Greenland.5km         (rev 0)
+++ branches/land_ice/icesheet/periodic_hex/namelist.input.Greenland.5km        2011-08-16 20:57:45 UTC (rev 943)
@@ -0,0 +1,8 @@
+&periodic_grid
+ nx = 301,
+ ny = 648,
+ dc = 5000.,
+ nVertLevels = 3,
+ nTracers = 1,
+ nproc = 2,
+/
Property changes on: branches/land_ice/icesheet/periodic_hex/namelist.input.Greenland.5km
___________________________________________________________________
Added: svn:executable
+ *
Modified: branches/land_ice/icesheet/src/icesheet.F
===================================================================
--- branches/land_ice/icesheet/src/icesheet.F        2011-08-16 19:40:21 UTC (rev 942)
+++ branches/land_ice/icesheet/src/icesheet.F        2011-08-16 20:57:45 UTC (rev 943)
@@ -795,7 +795,7 @@
write(6,*) xCell(1), yCell(1)
write(6,*) ' '
- pause
+! pause
do iCell=1,nCells
@@ -817,13 +817,7 @@
print *, 'y1(ix,iy) = ', y(iy)
print *, 'thck(ix,iy) = ', thkdata
-! if(ix .gt. 1 .and. ix .lt. nx )then
-! if(iy .gt. 1 .and. iy .lt. ny )then
-
-! if(ix .gt. 20 .and. ix .lt. nx-20)then
-! if(iy .gt. 20.and. iy .lt. ny-20)then
-
- if( thkdata .ge. 0 )then
+ if( thkdata .gt. 0 )then
kmt(iCell) = nVertLevelsMod
print *, 'cell kept'
print *, ' '
@@ -833,11 +827,6 @@
print *, ' '
endif
-! endif
-! endif
-
-! pause
-
enddo
! output no. of cells removed
Modified: branches/land_ice/icesheet/src/module_read_topo.F
===================================================================
--- branches/land_ice/icesheet/src/module_read_topo.F        2011-08-16 19:40:21 UTC (rev 942)
+++ branches/land_ice/icesheet/src/module_read_topo.F        2011-08-16 20:57:45 UTC (rev 943)
@@ -30,7 +30,6 @@
! nferr = nf_open('topo/gis_5km.180511.nc', NF_SHARE, rd_ncid)
nferr = nf_open('topo/gis_20km.180511.nc', NF_SHARE, rd_ncid)
write(6,*) ' nferr ncid ', nferr, rd_ncid
-
!
! Get IDs for variable dimensions
</font>
</pre>