遠地の波形

Baba et al. (2016) で提案されているSD+EL+GPの計算はJAGURSで出来ますが，EL (Elastic Loading) を有効にするため，

with_elastic_loading=1

として，パッケージとして配布されている

m_pyfile='PREM_Ggz.nc'

を使うと，GPも含んでいるため，漏れなく gravitational potential change も付いてきます．マニュアルには書かれていませんが，パッケージに入っている PREM_Gz.nc を使うと，ELのみが適用されるようです（確認中）．

海⽔密度効果は，独立に

with_density=1

を設定することで有効に出来ます．

If you specify

with_elastic_loading=1

option to apply elastic loading in your simulation using PREM_Ggz.nc file which is included in JAGURS package, the gravitational potential change is also adopted. Because PREM_Ggz.nc contains not only elastic loading but also gravitational potential change. If you omit GP, PREM_Gz.nc should be set as m_pyfile parameter. 

Error発生：
netcdf err=-51 in file=PREM_Ggz.nc
これは， https://www.unidata.ucar.edu/software/netcdf/docs/nc-error-codes.html によると，"Not a netcdf file"．

どうやら，ライブラリ関係が古かったためだったよう．ただし，詳細は詰めていない．
動いてた環境（MacPorts）：
  hdf5 @1.10.5_0+cxx+hl (active)
  netcdf @4.6.3_0+dap+netcdf4 (active)
  netcdf-fortran @4.4.5_5+gcc7 (active)

駄目だった環境（Source build @linux）：
  hdf5-1.8.19
  netcdf-4.4.1.1
  netcdf-fortran-4.4.4

駄目だった環境（linux）で，hdf5-1.10.5, netcdf-c-4.6.3, netcdf-fortran-4.4.5 を使ったところ，問題なくELが計算できた．

If you have an error of "netcdf err=-51 in file=PREM_Ggz.nc", please update hdf5 and netcdf libraries. I'm not sure the real reason, but I succeeded to execute JAGURS with hdf5-1.10.5, netcdf-c-4.6.3 and netcdf-fortran-4.4.5 instead of hdf5-1.8.19, netcdf-4.4.1.1 and netcdf-fortran-4.4.4.