RRI on iRIC の計算を実行したところriv(ii, jj) should be 1 ( 6 363 )( 6 364 )とのエラーが発生し、計算実行が中断してしまいました。このエラーの数値がどのファイルの値を示しているのか、またどのような修正を加えれば、エラーが解消されるのか、教えていただけますと大変助かります。
計算結果は以下に示す通りです。
# arguments
%0 = C:\Users\shun-\iRIC\solvers\rri_v1\go.bat
%1 = Case1.cgn
————————————————–
# RRI execute file path is
C:\Users\shun-\iRIC\solvers\rri_v1\
–
————————————————–
# Current directory is
C:\Users\shun-\.iRIC_workspace\67e57c1861259bbf9728225e9894aee4
–
————————————————–
# make “\out”
–
————————————————–
# make “\topo”
–
————————————————–
# Start RRI.exe
2023/10/30_23:58:50.20
————————————————–
# run rri.exe only
OMP: Warning #72: KMP_AFFINITY: affinity only supported for Intel(R) processors.
OMP: Warning #71: KMP_AFFINITY: affinity not supported, using “none”.
rainfile : E:\iRIC_river_analysis\used_RRI_river_analysis\rain_asc-dat\asc3\dat3\rain.dat
demfile : E:\iRIC_river_analysis\used_RRI_river_analysis\ask_file\part2\takasaka_2_elv.asc
accfile : E:\iRIC_river_analysis\used_RRI_river_analysis\ask_file\part2\takasaka_2-2_upg.asc
dirfile : E:\iRIC_river_analysis\used_RRI_river_analysis\ask_file\part2\takasaka_2_dir.asc
utm : 0
eight_dir : 1
lasth : 5
dt : 600
dt_riv : 60
outnum : 5
xllcorner_rain : 139.00000
yllcorner_rain : 35.00000
cellsize_rain_x : 0.01250 cellsize_rain_y : 0.00833
num_of_landuse : 5
dif : 1 1 1 1 1
ns_slope : 1.000 0.400 0.600 0.200 0.100
soildepth : 1.000 1.000 1.000 1.000 1.000
gammaa : 0.471 0.471 0.471 0.471 0.471
ksv : 0.556E-06 0.556E-06 0.000E+00 0.000E+00 0.556E-06
faif : 0.273 0.273 0.000 0.000 0.273
ka : 0.000E+00 0.000E+00 0.100E+00 0.000E+00 0.000E+00
gammam : 0.000 0.000 0.000 0.000 0.000
beta : 4.000 4.000 4.000 4.000 4.000
ksg : 0.000E+00 0.000E+00 0.000E+00 0.000E+00 0.000E+00
gammag : 0.000 0.000 0.000 0.000 0.000
kg0 : 0.000E+00 0.000E+00 0.000E+00 0.000E+00 0.000E+00
fpg : 0.000 0.000 0.000 0.000 0.000
rgl : 0.000E+00 0.000E+00 0.000E+00 0.000E+00 0.000E+00
ns_river : 0.030
riv_thresh : 100
width_param_c : 5.00
width_param_s : 0.35
depth_param_c : 0.95
depth_param_s : 0.20
height_param : 2.00
height_limit_param : 1000
num_of_landuse : 5
dx [m] : 25.0495865360156 dy [m] : 30.8218767174903
num_of_cell : 672968
total area [km2] : 519.581948980443
riv(ii, jj) should be 1 ( 6 363 )( 6 364 )
コメント
多少RRIに詳しいので返信させていただきます。
このエラーはi=6,j=363が河道として設定されているので、その下流のi=6,j=364も河道でないとおかしいですというエラーです。
添付されているascファイルはおそらくaccのファイルかと思いますが、i=6,j=363の次の河道メッシュはi=6,j=364と推察されます。
このaccファイルが正解だとするとおそらくdirファイルのi=6,j=363が128でないといけないのに、1になっているのではないかと思われます。
kon様 ご回答誠にありがとうございます。dirのファイルの方を確認してみます。
追加で大変申し訳ないのですが、iとjが緯度経度のどちらを表しているのか分かりかねているので、お聞きしてもよろしいでしょうか。
accやdir等の地理情報は、左上をi=1,j=1として下にi、右にjとなります。
引き続きありがとうございます。ご指定いただいたdirのセル値を確認したところ、ご提示していただいた通り、(6,363)が128、(6,364)が1と異なる値となっていました。まだ実行はできていませんが、こちらを修正して再度計算を行いたいと思います。お忙しい中、ご回答いただき大変ありがとうございました。
申し訳ありません。追加で確認させていただいてもよろしでしょうか。
追加でお伺いさせていただいてもよろしいでしょうか。