Takuya Miyashita
This site
Web
Start:
* fgmax メモ [#m531e0ae]
fg は fixed grid monitoring のこと.~
[[Fixed grid monitoring - Clawpack 5>https://www.clawpack...
fixed grid の中で,水位・流速等の最大値の空間分布の情報が...
** v5.7.0 [#c39b725d]
v5.7.0より,ファイルによる定義がなくなり,setrun.py 上で...
以下は定義の例.
#codeprettify{{
from clawpack.geoclaw import fgmax_tools
}}
#codeprettify{{
# Domain 1
fg = fgmax_tools.FGmaxGrid()
fg.point_style = 2 # uniform rectangular x-y grid
fg.x1 = -20790.0
fg.x2 = 28710.0
fg.y1 = -82690.0
fg.y2 = -19690.0
fg.dx = 90.0 # desired resolution of fgmax ...
fg.min_level_check = 1 # which levels to monitor max on
fg.tstart_max = 0.0 # just before wave arrives
fg.tend_max = 1.e10 # when to stop monitoring max ...
fg.dt_check = 60.0 # how often to update max values
fg.interp_method = 0 # 0 ==> pw const in cells, rec...
rundata.fgmax_data.fgmax_grids.append(fg) # written ...
# Domain 2
fg = fgmax_tools.FGmaxGrid()
fg.point_style = 2 # uniform rectangular x-y grid
fg.x1 = -11730.0
fg.x2 = 7770.0
fg.y1 = -49630.0
fg.y2 = -19930.0
fg.dx = 30.0 # desired resolution of fgmax ...
fg.min_level_check = 1 # which levels to monitor max on
fg.tstart_max = 0.0 # just before wave arrives
fg.tend_max = 1.e10 # when to stop monitoring max ...
fg.dt_check = 60.0 # how often to update max values
fg.interp_method = 0 # 0 ==> pw const in cells, rec...
rundata.fgmax_data.fgmax_grids.append(fg) # written ...
# num_fgmax_val
rundata.fgmax_data.num_fgmax_val = 5
}}
~
地形ファイルの範囲をそのまま使用したいときは,clawpack.ge...
#codeprettify{{
# Domain 4
topo_file = topotools.Topography(os.path.join('topofi...
fg = fgmax_tools.FGmaxGrid()
fg.point_style = 2 # uniform rectangular x-y grid
fg.dx = topo_file.delta[0] # desired resolutio...
fg.x1 = topo_file.x[0]
fg.x2 = topo_file.x[-1]
fg.y1 = topo_file.y[0]
fg.y2 = topo_file.y[-1]
fg.min_level_check = 1 # which levels to monitor max on
fg.arrival_tol = 1.0e-1
fg.tstart_max = 0.0 # just before wave arrives
fg.tend_max = 1.e10 # when to stop monitoring max ...
fg.dt_check = 1.0 # how often to update max values
fg.interp_method = 0 # 0 ==> pw const in cells, rec...
rundata.fgmax_data.fgmax_grids.append(fg) # written ...
}}
** v5.6.1 [#y11ad48d]
読み取りのカラムについては,GEOLIB の fgmax_finalize.f90 ...
#codeprettify(lang-fortran){{
do k=1,fg%npts
do mv=1,FG_NUM_VAL
if (abs(fg%valuemax(mv,k)) .lt. 1.d-90) t...
fg%valuemax(mv,k) = 0.d0
endif
enddo
write(FG_UNIT,111) fg%x(k),fg%y(k), fg%levelm...
(fg%valuemax(mv,k), mv=1,FG_NUM_VAL), &
(fg%tmax(mv,k), mv=1,FG_NUM_VAL), fg%ar...
111 format(2e17.8,i4,21e17.8)
enddo
}}
~
fort.FGX.aux1 には,1, 2 列目は x, y, 3列目以降に各 AMR l...
同じく GEOLIB の fgmax_finalize.f90 に,その出力形式が記...
#codeprettify(lang-fortran){{
do ma=1,FG_NUM_AUX
cma = char(ichar('0') + ma)
fname = 'fort.FG' // cfg // '.aux' // cma
print *, 'Writing to file ', fname
open(unit=FG_UNIT,file=trim(fname),status='un...
do k=1,fg%npts
write(FG_UNIT,112) fg%x(k),fg%y(k), &
(fg%aux(level,ma,k), level=1,mxnest)
112 format(2e17.8,20e17.8)
enddo
close(FG_UNIT)
enddo
}}
~
chile2010 の topo ファイルを例に,fgmax を全体に定義し,f...
解像度が細かい Level 3 の黒塗りの部分は,-0.99999000E+99 ...
抽出対象とする最小レベルを指定することができるが,対象時...
#htmlinsert(svg/fgmax_bath.svg)
~
fort.FG1.valuemax に吐き出される h は全水深のため,fort.F...
初期潮位があるときはさらに補正が必要.~
~
これらのファイルを読み取り図にしてみた.最大水位,最大流...
#htmlinsert(svg/FGmax4vars.svg)
End:
* fgmax メモ [#m531e0ae]
fg は fixed grid monitoring のこと.~
[[Fixed grid monitoring - Clawpack 5>https://www.clawpack...
fixed grid の中で,水位・流速等の最大値の空間分布の情報が...
** v5.7.0 [#c39b725d]
v5.7.0より,ファイルによる定義がなくなり,setrun.py 上で...
以下は定義の例.
#codeprettify{{
from clawpack.geoclaw import fgmax_tools
}}
#codeprettify{{
# Domain 1
fg = fgmax_tools.FGmaxGrid()
fg.point_style = 2 # uniform rectangular x-y grid
fg.x1 = -20790.0
fg.x2 = 28710.0
fg.y1 = -82690.0
fg.y2 = -19690.0
fg.dx = 90.0 # desired resolution of fgmax ...
fg.min_level_check = 1 # which levels to monitor max on
fg.tstart_max = 0.0 # just before wave arrives
fg.tend_max = 1.e10 # when to stop monitoring max ...
fg.dt_check = 60.0 # how often to update max values
fg.interp_method = 0 # 0 ==> pw const in cells, rec...
rundata.fgmax_data.fgmax_grids.append(fg) # written ...
# Domain 2
fg = fgmax_tools.FGmaxGrid()
fg.point_style = 2 # uniform rectangular x-y grid
fg.x1 = -11730.0
fg.x2 = 7770.0
fg.y1 = -49630.0
fg.y2 = -19930.0
fg.dx = 30.0 # desired resolution of fgmax ...
fg.min_level_check = 1 # which levels to monitor max on
fg.tstart_max = 0.0 # just before wave arrives
fg.tend_max = 1.e10 # when to stop monitoring max ...
fg.dt_check = 60.0 # how often to update max values
fg.interp_method = 0 # 0 ==> pw const in cells, rec...
rundata.fgmax_data.fgmax_grids.append(fg) # written ...
# num_fgmax_val
rundata.fgmax_data.num_fgmax_val = 5
}}
~
地形ファイルの範囲をそのまま使用したいときは,clawpack.ge...
#codeprettify{{
# Domain 4
topo_file = topotools.Topography(os.path.join('topofi...
fg = fgmax_tools.FGmaxGrid()
fg.point_style = 2 # uniform rectangular x-y grid
fg.dx = topo_file.delta[0] # desired resolutio...
fg.x1 = topo_file.x[0]
fg.x2 = topo_file.x[-1]
fg.y1 = topo_file.y[0]
fg.y2 = topo_file.y[-1]
fg.min_level_check = 1 # which levels to monitor max on
fg.arrival_tol = 1.0e-1
fg.tstart_max = 0.0 # just before wave arrives
fg.tend_max = 1.e10 # when to stop monitoring max ...
fg.dt_check = 1.0 # how often to update max values
fg.interp_method = 0 # 0 ==> pw const in cells, rec...
rundata.fgmax_data.fgmax_grids.append(fg) # written ...
}}
** v5.6.1 [#y11ad48d]
読み取りのカラムについては,GEOLIB の fgmax_finalize.f90 ...
#codeprettify(lang-fortran){{
do k=1,fg%npts
do mv=1,FG_NUM_VAL
if (abs(fg%valuemax(mv,k)) .lt. 1.d-90) t...
fg%valuemax(mv,k) = 0.d0
endif
enddo
write(FG_UNIT,111) fg%x(k),fg%y(k), fg%levelm...
(fg%valuemax(mv,k), mv=1,FG_NUM_VAL), &
(fg%tmax(mv,k), mv=1,FG_NUM_VAL), fg%ar...
111 format(2e17.8,i4,21e17.8)
enddo
}}
~
fort.FGX.aux1 には,1, 2 列目は x, y, 3列目以降に各 AMR l...
同じく GEOLIB の fgmax_finalize.f90 に,その出力形式が記...
#codeprettify(lang-fortran){{
do ma=1,FG_NUM_AUX
cma = char(ichar('0') + ma)
fname = 'fort.FG' // cfg // '.aux' // cma
print *, 'Writing to file ', fname
open(unit=FG_UNIT,file=trim(fname),status='un...
do k=1,fg%npts
write(FG_UNIT,112) fg%x(k),fg%y(k), &
(fg%aux(level,ma,k), level=1,mxnest)
112 format(2e17.8,20e17.8)
enddo
close(FG_UNIT)
enddo
}}
~
chile2010 の topo ファイルを例に,fgmax を全体に定義し,f...
解像度が細かい Level 3 の黒塗りの部分は,-0.99999000E+99 ...
抽出対象とする最小レベルを指定することができるが,対象時...
#htmlinsert(svg/fgmax_bath.svg)
~
fort.FG1.valuemax に吐き出される h は全水深のため,fort.F...
初期潮位があるときはさらに補正が必要.~
~
これらのファイルを読み取り図にしてみた.最大水位,最大流...
#htmlinsert(svg/FGmax4vars.svg)
Page:
Edit with a page name which already exists