v3.0.1.1
Pre-releaseGitHub TITAN2D/titan2d -b default-grid-resolution-mods repository release v3.0.1.1
Modifications to set the default grid resolutions for titan2d output, to the e-w and n-s resol values specified in the GIS raster map.
Used the IntelliJ IDEA program and verified v3.0.1 release of the GitHub TITAN2D/titan2d repository is identical to the gmfg titan-3.0.0 release except for formatting, minor modifications to configuration.ac, INSTALL updated with additional instructions for when installing titan from the github repository, a new file titan2d.h added, the use of -I$(srcdir) in Makefile.am files, and the removal of the Makefile.in files.
Created a branch of the v3.0.1 release and named this branch default-grid-resolution-mods. Tagged default-grid-resolution-mods modifications to release v3.0.1.1. Performed tests for release v3.0.1.1 on a laptop running the Fedora 20 operating system and on UB's ccr general compute cluster.
Tested release v3.0.1.1 with GIS formats GIS_GRASS and GDAL. As tested, the GDAL format requires the GDAL-GRASS plugin. Please see the GDAL-GRASS plugin section below for notes on how to compile and install the GDAL-GRASS plugin.
Test Environment (07/03/2015):
Fedora 20 operating system
grass-6.4.4
gdal-1.11.0
gdal-grass-1.11.2
proj 4.8.0
hdf5-1.8.11
Tested gmfg titan 3.0.0 with default grid resolution modifications and GitHub TITAN2D/titan2d -b default-grid-resolution-mods release v3.0.1.1.
Downloaded titan2d-3.0.1.1.tar.gz from https://github.com/TITAN2D/titan2d/releases; entered: tar -zxf titan2d-3.0.1.1.tar.gz to extract.
Set export LD_LIBRARY_PATH="/usr/local/lib:/usr/local/src/hdf5-1.8.11/hdf5/lib".
GIS Format: GIS
titan2d build:
./configure --with-mpi-includes=/usr/local/include --with-mpi-libs=/usr/local/lib --with-hdf5=/usr/local/src/hdf5-1.8.11/hdf5 --with-grass=/usr/local/grass-6.4.4 --without-gdal CC=gcc F77=gfortran, make and make install.
Tested with 1 processor; simulation.data configured for GIS and using the following GIS raster map:
/home/Renette/grass.data/grass5 Colima ColimaR ColimaR:
proj: 1
zone: 13
north: 2210030
south: 2109950
east: 700050
west: 599970
cols: 1667
rows: 1667
e-w resol: 60.0359928
n-s resol: 60.0359928
format: -1
compressed: 1
Result from titan_preprocess:
The GIS region is: x = 5.999700e+05 to 7.000500e+05 y = 2.109950e+06 to 2.210030e+06.
The simulation region is: x = 5.999700e+05 to 7.000500e+05 y = 2.109950e+06 to 2.210030e+06.
The simulation region excluding the margin is: x = 6.097778e+05 to 6.902422e+05 y = 2.119758e+06 to 2.200222e+06.
Ran titan for 1000 time steps.
Result:
elevation.grid header:
Nx=1340: X={ 609777.83999976, 690242.15999784}
Ny=1340: Y={ 2119757.8399998, 2200222.1599978}
pileheightrecord header:
Nx=1340: X={ 609777.83999976, 690242.15999784}
Ny=1340: Y={ 2119757.8399998, 2200222.1599978}
GDAL Format: GDAL
titan2d build:
./configure --with-mpi-includes=/usr/local/include --with-mpi-libs=/usr/local/lib --with-hdf5=/usr/local/src/hdf5-1.8.11/hdf5 -with-grass=/usr/local/grass-6.4.4 --with-gdal-config=gdal-config CC=gcc F77=gfortran, make and make install.
Verified gdal-grass installation:
/usr/local/bin/gdalinfo /home/Renette/grass.data/grass5/Colima/ColimaR/cellhd/ColimaR.
Result:
Driver: GRASS/GRASS Database Rasters (5.7+)
Files: /home/Renette/grass.data/grass5/Colima/ColimaR/cellhd/ColimaR
Size is 1667, 1667
Coordinate System is:
PROJCS["UTM Zone 13, Northern Hemisphere",
GEOGCS["clark66",
DATUM["unknown",
SPHEROID["Clarke_1866",6378206.4,294.9786982]],
PRIMEM["Greenwich",0],
UNIT["degree",0.0174532925199433]],
PROJECTION["Transverse_Mercator"],
PARAMETER["latitude_of_origin",0],
PARAMETER["central_meridian",-105],
PARAMETER["scale_factor",0.9996],
PARAMETER["false_easting",500000],
PARAMETER["false_northing",0],
UNIT["Meter",1]]
Origin = (599970.000000000000000,2210030.000000000000000)
Pixel Size = (60.035992801439711,-60.035992801439711)
Corner Coordinates:
Upper Left ( 599970.000, 2210030.000) (104d 2'40.01"W, 19d59' 7.79"N)
Lower Left ( 599970.000, 2109950.000) (104d 2'59.11"W, 19d 4'52.10"N)
Upper Right ( 700050.000, 2210030.000) (103d 5'17.31"W, 19d58'39.95"N)
Lower Right ( 700050.000, 2109950.000) (103d 5'55.50"W, 19d 4'25.62"N)
Center ( 650010.000, 2159990.000) (103d34'12.98"W, 19d31'48.67"N)
Band 1 Block=1667x1 Type=Float32, ColorInterp=Gray
Min=140.000 Max=4219.000
NoData Value=nan
Metadata:
COLOR_TABLE_RULES_COUNT=0
Tested with 1 processor; simulation.data configured for GDAL and the same GIS raster map used for the GIS test:
/home/Renette/grass.data/grass5/Colima/ColimaR/cellhd/ColimaR.
Ran titan for 1000 time steps.
Result:
elevation.grid header:
Nx=1340: X={ 609777.84, 690242.16}
Ny=1340: Y={ 2119757.84, 2200222.16}
pileheightrecord header:
Nx=1340: X={ 609777.84, 690242.16}
Ny=1340: Y={ 2119757.84, 2200222.16}
Test Environment (07/03/2015):
ccr x86_64 GNU/Linux (per uname -a) operating system
grass-6.4.4
gdal-1.11.2
gdal-grass-1.11.2
proj.4-4.9.1
hdf5-1.8.15
Tested TITAN2D/titan2d -b default-grid-resolution-mods release v3.0.1.1.
Downloaded titan2d-3.0.1.1.tar.gz from https://github.com/TITAN2D/titan2d/releases; entered: tar -zxf titan2d-3.0.1.1.tar.gz to extract.
Set export LD_LIBRARY_PATH="/util/academic/intel/impi/5.0.2.044/intel64/lib:/gpfs/user/renettej/TITAN2D/hdf5-1.8.11/hdf5/lib:/projects/gmfg/renettej/grass-6.4.4/local/grass-6.4.4/lib:/util/academic/grass/proj.4-4.9.1/lib:/projects/gmfg/renettej/gdal-1.11.2/local/lib".
GIS Format: GIS
titan2d build:
./configure --with-mpi-includes=/util/academic/intel/impi/5.0.2.044/intel64/include --with-mpi-libs=/util/academic/intel/impi/5.0.2.044/intel64/lib --with-hdf5=/util/academic/grass/hdf5-1.8.15/hdf5 --with-grass=/projects/gmfg/renettej/grass-6.4.4/local/grass-6.4.4 --without-gdal CC=mpicc CXX=mpicxx F77=mpif90, make and make install.
Tested with 1 and 6 processors; simulation.data configured for GIS and using the following GIS raster map:
/gpfs/user/renettej/TITAN2D/grass.data/grass5 Colima ColimaR ColimaR:
proj: 1
zone: 13
north: 2210030
south: 2109950
east: 700050
west: 599970
cols: 1667
rows: 1667
e-w resol: 60.0359928
n-s resol: 60.0359928
format: -1
compressed: 1
Result from titan_preprocess same as laptop result.
Ran titan for 1000 time steps.
Result:
elevation.grid header:
Nx=1340: X={ 609777.83999976, 690242.15999784}
Ny=1340: Y={ 2119757.8399998, 2200222.1599978}
pileheightrecord header:
Nx=1340: X={ 609777.83999976, 690242.15999784}
Ny=1340: Y={ 2119757.8399998, 2200222.1599978}
Display of titan2d output for 6 processor run:
GDAL Format: GDAL
titan2d build:
./configure --with-mpi-includes=/util/academic/intel/impi/5.0.2.044/intel64/include --with-mpi-libs=/util/academic/intel/impi/5.0.2.044/intel64/lib --with-hdf5=/util/academic/grass/hdf5-1.8.15/hdf5 --with-grass=/projects/gmfg/renettej/grass-6.4.4/local/grass-6.4.4 --with-gdal-config=/projects/gmfg/renettej/gdal-1.11.2/local/bin/gdal-config CC=mpicc CXX=mpicxx F77=mpif90, make and make install.
Verified gdal-grass installation:
/projects/gmfg/renettej/gdal-1.11.2/local/bin/gdalinfo /gpfs/user/renettej/TITAN2D/grass.data/grass5/Colima/ColimaR/cellhd/ColimaR.
Result same as laptop result.
Tested with 1 and 6 processors; simulation.data configured for GDAL and the same GIS raster map used for the GIS test:
/gpfs/user/renettej/TITAN2D/grass.data/grass5/Colima/ColimaR/cellhd/ColimaR.
Ran titan for 1000 time steps.
Result:
elevation.grid header:
Nx=1340: X={ 609777.84, 690242.16}
Ny=1340: Y={ 2119757.84, 2200222.16}
pileheightrecord header:
Nx=1340: X={ 609777.84, 690242.16}
Ny=1340: Y={ 2119757.84, 2200222.16}
Display of titan2d output for 6 processor run:
GDAL-GRASS plugin
The GDAL-GRASS plugin allows direct access to GRASS raster map data from titan2d GDAL enabled software.
References:
http://grasswiki.osgeo.org/wiki/Compile_and_install_GDAL-GRASS_plugin
https://trac.osgeo.org/gdal/wiki/GRASS
Test Environment (07/01/2015):
Fedora 20 operating system
grass-6.4.4
gdal-1.11.0
gdal-grass-1.11.2
proj 4.8.0
hdf5-1.8.11
Tested gmfg titan 3.0.0 with default grid resolution modifications and GitHub TITAN2D/titan2d -b default-grid-resolution-mods release v3.0.1.1.
Note: The following steps were performed as su. Packages were installed to standard directories.
- Built grass:
In the grass-6.4.4 subdirectory entered: ./configure, make and make install.
- Built gdal without grass support:
In the gdal-1.11.0 subdirectory entered: ./configure -without-grass, make and make install. Grass support is provided via the plugin.
- Built gdal-grass:
Created /etc/ld.so.conf.d/grass6.conf containing the line /usr/local/grass-6.4.4/lib. Entered: ldconfig.
In the gdal-grass-1.11.2 subdirectory, entered: ./configure --with-gdal=/usr/local/bin/gdal-config --with-grass=/usr/local/grass-6.4.4, make and make install.
Test Environment (07/01/2015):
ccr x86_64 GNU/Linux (per uname -a) operating system
grass-6.4.4
gdal-1.11.2
gdal-grass-1.11.2
proj.4-4.9.1
hdf5-1.8.15
Tested GitHub TITAN2D/titan2d -b default-grid-resolution-mods release v3.0.1.1.
Note: on ccr cannot perform steps as su. PROJ.4 not installed to standard directories. Needed to add required paths to the LD_LIBRARY_PATH environment variable instead of using ldconfig. For examples given below, accessing ccr /util/academic/grass subdirectories and building in ccr /projects/gmfg/renettej subdirectories.
- Built grass:
In the grass-6.4.4 subdirectory, created the local subdirectory, edited configure and set both occurrences of ac_default_prefix=/projects/gmfg/renettej/grass-6.4.4/local. Needed to explicitly configure with the PROJ.4 directories to prevent external PROJ.4 errors, entered: ./configure --with-proj-includes=/util/academic/grass/proj.4-4.9.1/include --with-proj-libs=/util/academic/grass/proj.4-4.9.1/lib --with-proj-share=/util/academic/grass/proj.4-4.9.1/share/proj, make and make install.
- Built gdal without grass support:
In the gdal-1.11.2 subdirectory, created the local subdirectory and entered: ./configure --prefix=/projects/gmfg/renettej/gdal-1.11.2/local --without-grass, make and make install. Grass support is provided via the plugin.
Set gdal-config, entered: /projects/gmfg/renettej/gdal-1.11.2/local/bin/gdal-config --libs gdal --cflags -I/projects/gmfg/renettej/gdal-1.11.2/local/include.
For a precheck entered: /projects/gmfg/renettej/gdal-1.11.2/bin/gdal_translate --formats | grep -i grass.
Result:
GRASSASCIIGrid(rov): GRASS ASCII Grid
- Built gdal-grass:
Needed to make sure /projects/gmfg/renettej/grass-6.4.4/local/grass-6.4.4/lib and /util/academic/grass/proj.4-4.9.1/lib were exported to the LD_LIBRARY_PATH environment variable. Made sure LD_LIBRARY_PATH path names were colon separated.
In the gdal-grass-1.11.2 subdirectory, entered: ./configure --with-gdal=/projects/gmfg/renettej/gdal-1.11.2/local/bin/gdal-config --with-grass=/projects/gmfg/renettej/grass-6.4.4/local/grass-6.4.4.
(Per http://git.net/gdal-development-gis-osgeo/dsc07579.html#07580 from Frank Waterdam - need to hand edit the Makefile and
add "-I/usr/local/include" or whatever is appropriate to
refer to the PROJ.4 include files.)
Modified the Makefile, added -I/util/academic/grass/proj.4-4.9.1/include to CPPFLAGS. Also in the Makefile, set GRASSTABLES_DIR to /projects/gmfg/renettej/gdal-1.11.2/local/share/gdal/grass; the gdal-grass plugin copies files from the /projects/gmfg/renettej/grass-6.4.4/local/grass-6.4.4 subdirectories to this directory. The gdal-grass plugin copies generated gdal_GRASS.so and ogr_GRASS.so files to the
/projects/gmfg/renettej/gdal-1.11.2/local/lib/gdalplugins directory. Entered make and make install.
For check entered: /projects/gmfg/renettej/gdal-1.11.2/bin/gdal_translate --formats | grep -i grass.
Result:
GRASS (r0): GRASS Database Rasters (5.7+)
GRASSASCIIGrid(rov): GRASS ASCII Grid