GeoHyMos and OHyMos
User’sUsers Mannual - v1.0
v1 0
(with Kamishiiba Dam Basin)

Sunmin Kim
Graduate School of Eng.,
Kyoto University(June 15, 2011)

GeoHyMos
and KsEdge for OHyMos





GeoHyMos-1:
to extract channel network

GeoHyMos-2:
to compose basin and river

KsEdge:KEdrainfall-runoff simulation

Every program for GeoHyMos and OHyMos is downloadablefrom our web ppage;g ; http://hywr.kuciv.kyoto-u.ac.jp
p yyjpIdentifying
Kamishiiba Dam Basin
Digital Map of JapanDEM (250m, 1km) of Japanfor “makeGridData250” ofGeoHyMos-2yFirst and Second OrderMesh Code
7First Order Mesh4831

654324831: first order mesh483141 : second order mesh

100

1

2

34

56

7

1-1. KStoPF-1.6
t convertt KS fformatt channel
tohl datad t into
i t PF format
ftCommand./KStoPF ks_271.dat ks_272.dat 40pf.dat 483140 483140- ks_271.dat: vector data of unit basin- ks_272.dat: unit basin of channel network- 41pf.dat: output, PF(Plain File) format channel data- 483141: second order mesh code
- 483140: reference second order mesh code
# shell script may be helpful for a continuous procedure1-2. preRiverPlot-3.2t plot
tol t theth PF format
ft datad t usingitheth GNUplot
GNU l tCommand./preRiverPlot.pl ../1.KStoPF-1.6/41pf.dat > 41pfplot.dat- 41pf.dat: the PF format data from the ‘KStoPF’ procedure- 41pfplot.dat: output# shell script may be helpful for a continuous procedure# if any error is noticed in the plotted datadata,
the error should be corrected in the PF file, manually.

Example
El off PF file
fil correction;
ti483151.pf

1-3. gcnd-1.7
t merge all
toll theth PF files
fil into
i t oneCommand./gcnd KmshbOG.dat KmshbNG.dat- KmshbOG.dat: original data information
number of mesh files,files
position and title of each mesh file- KmshbNG.dat: new data information
title of output files
four edge point coordinate of new mesh
- Kmshb_gcdn-1.1.dat: output, merged mesh data
- Kmshb_gcdn_GP.dat: output, GNUplot file
# for a proper procedure,the ‘four edge point coordinate of newmesh’ should not be narrower than theoriginal
i i l meshhddata.
t1-4. altitudeToPFdata-1.3
tto give
i an altitudeltit d informationi f
ti intoi t every pointi t off
river channel by interpolation
Command./altitudeToPFdata Kmshb_gcnd-1-1.dat KmshbATPF.out- Kmshb_gcnd-1-1.dat: the merged mesh data
- KmshbATPF.out: output
EX)
(91)
100

(87)
80

If ththe altitudeltit d values
loff one segment’st’ end
d are 100 andd 8080,
91 and 87 will be given into the middle of the reach.

1-5. segNum-2.2 & 6. riverPlot-1.2
t giveto
i a numberb tot eachh segmentt andd plotl t itCommand./segNum KmshbATPF.out KmshbSN.out./riverPlot.sh KmshbATPF.out KmshbRP.plt- KmshbATPF.out: PF format data from ‘4.altitudeToPFdata-1.3’- KmshbSN.out: output, segment numbering- KmshbRP.plt: output,output GNUplot format with segment numbers# Form the plotted image,segment number for the outlet (28)
off theth Kmshb
K hb BBasin
i iis noticeable.
tibl

1-7. extractNetwork-1.2
tto extracttt channelhl network
tkusing the outlet segment numberCommand./extractNetwork KmshbATPF.out KmshbESN28.out 28- KmshbATPF.out: PF format data from ‘4.altitudeToPFdata-1.3’- KmshbESN28.out: output, extracted channel network data
- 28: outlet segment numberEX) plotted image using the ‘riverplot
riverplot.shsh’# segment no. 28
was modified toexclude its sub-basin.

1-8. cutSegment-1.1
tto cutt a segmentt att an arbitrary
bitpoint
i t(for giving a smaller sub-basin)
Command./cutSegment KmshbATPF.out KmshbCP.dat KmshbCS.out- KmshbATPF.out: PF format data from ‘4.altitudeToPFdata-1.3’- KmshbCP.dat: cutting point information
number of cutting points,
UTM coordinate of the points- KmshbCS.out: output,
outputsame format to ‘KmshbATPF.out’EX) original
i i ldata

cutt seg. andd
re-extracted

2-1. makeGridData250-2.0
t convertt SEM grid
toid datad t intoi t
Command./makeGridData250 KmshbGD.dat-KmshbGD.dat: input data file
1 1 1: total mesh no. = x-axis no. ×y-axis no.
4831+.SEM 1 1: first mesh code and it’s position
483140: reference second mesh code
483140: second mesh code of south-east corner483171: second mesh code of north-west corner
KmshbGD.out: output, elevation ((m)) data with UTM coordinate (km)( )KmshbGI.out: output, grid information
KmshbGN.out: output, number of total grid# file format of the original SEM fileshould be changed in a way,
nkf –Lu –e 4831.SEM > 4831+.SEM

483171483140First order mesh4831

2-2. remakeGridData-1.3
t unify
toif theth gridid datad t intoi t ‘m’‘ ’ unititCommand./remakeGridData KmshbRGD.dat- KmshbRGD.dat: input data file
KmshbGD.out: original data of grid elevation
KmshbGI.out: original data of grid information
KmshbGN out: original data of grid number
KmshbGN.out:
KmshbRGD.out: output, elevation (m) data with UTM coordinate (m)KmshbRGI.out: output, grid information
EX)
KmshbGD.out
KmshbRGD out
KmshbRGD.out

2-3. remakeRiverData-1.6
t reformulate
tofl t theth channelhl datad t forf theth gridid datad t
Command./remakeRiverData KmshbRRD.dat- KmshbRRD.dat: input data file
KmshbRGD.out: input, grid data
KmshbRGI.out: input, grid information
KmshbESN28 out: input
KmshbESN28.out:
input, channel data
KmshbNewGrid.out: output, grid data with UTM and lattice coordinate
KmshbNewRiver.out: output, reformulated channel data
KmshbCheckList.out: output, checklist after the procedure
KmshbGpMesh.out: output, mesh plot data
KmshbGpRiver.out: output, original channelplot data for the GNUplot
KmshbGpNewRiver out: output,
KmshbGpNewRiver.out:
output reformulatedchannel plot data
KmshbSegNum.out: output, segment numberchecklist

2-4. removeSegment-1.2
t remove unnecessary segmentt aftertoft ththe reformf
Command./removeSegment KmshbRS.dat- KmshbRS.dat: input data file
KmshbNewRiver.out: input, channel data
KmshbSegNum.out: input, segment number data
KmshbNewRiver2 out: output,
KmshbNewRiver2.out:
output channel data
KmshbSegNum2.out: output, segment number data
# if there is anyy segmentgremoved after this procedure,
KmshbNewRiver2.out, and KmshbSegNum2.outshould be used from now on.

2-5. temporaryConnect-2.1
tto produce
dconnection
tidatad t off
each grid and channel network
Command./temporaryConnect KmshbTC.dat-KmshbTC.dat: input data file
KmshbNewGrid.out: input, grid data
KmshbNewRiver(2).out: input, channel data
KmshbConnect out: output,
KmshbConnect.out:
output grid connection data
KmshbNewRiverConnect.out: output, channel connection data
KmshbGpConnect.out: output, plot file for connection checking

2-6. depression-1.3
t fill sinktoi k tot properly
l decided id a flowfldirection
ditiCommand./depression KmshbDpr.dat- KmshbDpr.dat: input data file
KmshbConnect.out: input, grid connection data
KmshbDpr.out: output, grid connection data2-7. generateNodeEdge-2.1
t fi
tofindd a proper flowfldirection
diti off eachh grididCommand./generateNodeEdge KmshbNodeEdge.dat- KmshbNodeEdge.dat: input data
KmshbDpr.out: input, grid connection data
KmshbNewRiverConnect.out: input, channel connection data
KmshbNode out: output,
KmshbNode.out:
output node data
KmshbEdge.out: output, edge data
KmshbLinkEdge.out: output,
KmshbNewConnect.out: output,

2-8. extractNodeEdge-3.2
tto ddelineate
lit basinb i areaby extracting the node and edge data
Command./extractNodeEdge KmshbENE.dat 18- KmshbENE.dat: input data file
- 18: segment no. to extract
KmshbNode.out: input, node data
K hbEdKmshbEdge.out:t input,
it edged datad t
KmshbSegNum(2).out: input, segment number data
KmshbLinkEdge.out: input,
KmshbEN18.out: output,
p , extracted node data for seg.num.
g18
KmshbEE18.out: output, extracted edge data
KmshbGpEdge18.out: output, edge plot data
# shell script ma
may be helpf
helpfull

3-1. toEdgeFmt-2.1
t convertt edge
tod format
ft V.0V 0 iinto
t fformatt V.4
V4
Command./KmshbEfV4.shKmshbEE28.out: input, edge data for seg.no. 28
K hbEN28 t : iinput,
KmshbEN28.out:
t noded ddata
t ffor seg.no. 28
KmshbEfV4-28.out: : output, edge format V.4 for seg.no. 28

3-2. edgeConnect-3.2
t produce
todedged connection
tidatad t
Command./KmshbECV4.shKmshbEfV4-28.out: input, edge data in format V.4
KmshbECV4-28.out: output, edge connection data in format V.4

3-3. edgeOrder-3.2t givetoi an orderd forf a runoff
ff simulationi l tiCommand./KmshbOrderECV4.shKmshbECV4-28.out: input, edge connection data in format V.4
KmshbROO-28.out: output, edge ordering data in format V.4

3-4. edgeUpperArea-2.2
t calculate
tol l t upper area off eachh edge
dCommand./KmshbECV5.shKmshbROO-28.out: input, edge ordering data in format V.4
KmshbECV5-28.out: output, edge connection data in format V.5

3-5. mkPmtIst
tto producedparametert andd initial
i iti l condition
ditifilefilfor the KsEdge model
Command./KmshbPmtIst.shKmshbECV5-28.out: input, edge connection data in format V.5
K hbP
KmshbParameter.dat:t d t input,
it parametert andd initial
i iti l condition
diti datad t
Kmshb28.pmt: output, parameter data file
Kmshb28.ist: output, initial condition data fileEX) KmshbParameter.dat including notation

3-6. calSegShape-1.2
t calculate
tol l t lengthl
th andd slope
loff eachh segmenttCommand./calcSegShape KmshbESN28.out KmshbSegShp.out 0.001- KmshbESN28.out: input, extracted channel data of“../../GeoHyMos-1/7.extractNetwork-1.2”- KmshbSgSp.out: output, length and slope values of each seg.
- 0.001: default slope value

length and slope information:
segmentt code,d llengthth ((m),) slopel3-7. ksedgefc2d-1.3.1
rainfall-runoff
i f llff simulationi l ti withith Kinematic
Kiti WaveWModel
M d lCommand./OHyMos SCF_Kmshb.dat SCF_Kmshb- SCF_Kmshb.dat: structure configuration file
- SCF_Kmshb: to specify ‘pmt’ and ‘ist’ files

3-7. ksedgefc2d-1.3.1
SCF K hb d t
SCF_Kmshb.dat

Specifying the numbers of each component;
rainfall, edge and node, iteration, output ports
Specifying rainfall data and it’s location

Specifying the components of edge (sub-basin)
and node (channel)

Specifying the output file name3-7. ksedgefc2d-1.3.1
SCF K hb d t
SCF_Kmshb.dat

Rainfall r0 goes to sub-basin c0 by receive port rp0
Sub-basin c0 output of send port sp0 goes to channel c9 byreceive port lateral rpl0
Channel c9 output of send port sp0 goes to next channel
c11 by receive port upper rpu0
Channel c9 output of send port sp0 goes to output file s0
Rainfall, r0

rpu0
0

rp0

rp0

sub-basinc0
sp0

sub-basin
c2
sp0

rpl0
channel,hl c9
9 sp00s0

rpu0
0

rpl0
channel,hl c1111 sp00s2

3-7. ksedgefc2d-1.3.1
SCF K hb
SCF_Kmshb.pmt
tSpecifying path of each pmt
file for each edge

Channel data for “Kinem”
length, width, slope,manning’s n, no ofdx and dt, θ3-7. ksedgefc2d-1.3.1
SCF K hb i t
SCF_Kmshb.ist

Specifying path of each initial condition file for each edge

for more information
sunmin@hywr.kuciv.kyoto u.ac.jp
sunmin@hywr.kuciv.kyoto-u.ac.jp